Print version ISSN 1806-1117
Rev. Bras. Ensino Fís. vol.34 no.2 São Paulo Apr./June 2012
Radiação de corpo negro como função da frequência e comprimento de onda: uma abordagem orientada experimentalmente
Ademir L. Xavier Jr.1; Sergio Celaschi
Centro de Pesquisas Avançadas Wernher von Braun, Campinas, SP, Brasil
It is well known that the peak of black body radiation depends on the adopted description, whether as a function of frequency or wavelength. This paper is a pedagogical exposition of how to write and estimate Planck's distribution spectrum on two distinct spectrograph scales. We show that the emitted power of the Black Body maximum is unique and does not coincide with the peak measured whether in frequency or wavelength. We then make a generalization with other distribution functions in order to further clarify the apparent issue.
Keywords: Planck's law, Wien's law, black body, electromagnetic radiation.
Sabe-se que o pico de radiação de corpo negro depende da descrição adotada, a saber, se em frequência ou em comprimento de onda. Este artigo traz uma exposição pedagógica sobre como escrever e estimar o espectro da distribuição de Planck em duas escalas espectrográficas distintas. Mostramos que a potência emitida do máximo de corpo negro é única e não coincide com o pico medido, seja em frequência ou comprimento de onda. Então, fazemos uma generalização com outras funções de distribuição com o objetivo de clarificar ainda mais a questão.
Palavras-chave: lei de Planck, lei de Wien, corpo negro, radiação eletromagnética.
Every student learning the principles of Quantum Mechanics and following a historical approach is naturally impressed by Max Planck pioneering work . He faced with heroism unsurmonable problems at his time to explain the Black Body (BB) radiation properties [2-5]. He successfully managed to derive a law and solve one of the greatest anomalies that paved pre-quantum physics in the beginning of the 20th century.
A dedicated student might be impressed to know that BB distribution law gives rise to two well distinct maxima if the BB spectrum is written in terms of frequency or wavelength intervals. Many books in Statistical Physics [6-13] prefer the frequency form. This is not surprising since Planck's distribution law is obtained from Bose-Eistein statistics of electromagnetic oscillator modes. Almost none of these books present any values at all for the resulting temperature displacement of the maximum or Wien's displacement law. An exception is Landsberg . On p. 214 (section 13.3), after deducing the frequency law, the wavelength form is obtained. However, historical accounts of the matter show that the spectrum was originally measured on a wavelength scale using a spectrograph [15,16]. Most of the time, derivations come from exclusively mathematical reasoning (Section 2), with scarce or no reference to the experimental setup required to register the BB spectrum.
It is well known that grating laws give rise to a multitude of BB maxima . The subject has also served to motivate speculations on the causes of human vision sensitivity function . The subject is further motivated by the dependence of the BB temperature with the maximum of the spectrum. In other words, the position of this maximum is a direct measurement of the BB temperature which has several applications from material science to astrophysics or even cosmology . If the peak changes position with the representation, which one should be chosen as a reference of the BB temperature? However, spectral measurements dependent strongly on experimental details that are many times absent in reference works. This paper intends to fill this gap following an experimental reasoning (that is, we numerically estimate the BB spectrum as seen by different experimental arragements, while no real experimental measurement was attempt). The approach describes the energy distribution of BB photons on spectrographic scales where wavelength are registered. In fact, there seems to exist no experimental >linear scale for wavelength or frequency on which the BB spectrum maxima would coincide with the theoretical distribution. Indeed, we will show here that at the spectrograph there is no 'ambiguity' in the maximum peak position. In Section 2 we briefly present the problem as it is often stated without any reference to the experimental problem. In Section 3 we discuss some experimental arrangements that are employed to obtain the BB spectrum in wavelength. In Section 4 we generalize our results and make a simple but elucidating application that we hope can further clarify the problem. Finally, Section 5 presents some conclusions.
2. BB spectrum
Planck's distribution function  gives the intensity of the emitted radiation (emitted power per unit area), per unit of a physical quantity interval (frequency, wavelength or wavenumber), from a BB at absolute temperature T. It takes two forms in frequency (ν) and wavelength (λ), respectively
where h = 6.62× 10-34 Js, c = 3× 108 m Hz and k = 1.38 ×10-23 JK-1. Equations (1) and (2) are such that the total energy flux in both descriptions are the same
for all T. These integrals give rise to the well known Stefan's law for the total energy emitted by a BB . The relation between frequency and wavelength is given by λν = c. Imposing conditions for the maxima of Eqs. (1) and (2), two distinct characteristic equations are found for frequency and wavelength respectively
with x = hν/kT and y = hc/λkT. Equations (4) and (5) are similar in form, but result in distinct roots or distinct displacement laws as a function of frequency and wavelength, respectively
As a numerical example, for T = 5500 K (i.e. an approximate value for the Sun's photosphere temperature), the maxima in frequency is νmax = 325.2 THz and in wavelength, λmax = 532.5 nm. Applying the frequency wavelength relationship to the maximum in frequency, the equivalent maximum in wavelength (λ) is located at 922.5 nm. While the maximum in wavelength is in the visible part of the spectrum, the maximum in frequency is located in the near infrared region. As we will see, such maxima occur at different wavelengths depending on the used 'dispersion law'. Regardless of the particular dispersion rule chosen, all possible BB distributions still represent the same physical BB emitted spectrum.
3. Experimental setup for measuring BB spectrum
The relationship between emitted and measured BB spectra is governed by a transfer function that depends on a variety of physical and mathematical properties. These include, among others, transmission scattering, material dispersion and absorption rules as well as quantum responsivity of the detector as a function of frequency or wavelength. To experimentally register the BB spectrum, we need a way to separate frequency or wavelength intervals. In general, light intensity can be measured at distinct wavelengths within a small bandwidth by employing a spectrograph and a bolometer . In the following, we present numerical estimates for the BB spectrum through different spectrograph arrangements.
3.1. Diffraction grating spectrograph
Fig. 1 shows a schematic arrangement of this type using a diffraction grating with normal incidence. A light source S is positioned in front of a slit with width Δw. The light beam is collimated by a lens L1 with focal distance f and diffracted by a grating G. At the detector, an image of the slit is formed. In what follows, we derive the BB spectrum distribution for a diffraction grating. The spectrum of S is scanned by rotating the detector arm as shown in Fig. 1 by an angle θ in relation to the undeflected direction. For reflection gratings with normal incidence, there is a straightforward and simple relationship between the wavelength and the deflection angle at first order 
where є is the groove spacing (lines per inch or mm). We only analyse the spectrum of a single grating order. According to Eq. (8), the grating theoretical range goes from λmin = 0 to λmax = є. If f is the collimator focal distance, then the wavelength resolution will be given by
The spectrograph measures light intensity as a function of the body emissivity and angular resolution δθ = Δw/f at the position θ, or a function proportional to Iθ(θ, T)δθ for which
In reality the sensor photocurrent will be proportional to a function F(IGδθ) that depends on many other details such as: detector response function, grating efficiency curve , sensor operational area and its distance from the grating (given by the focal distance of lens L2). For instance, the detector area can be taken to be lf2δθ where l is an orthogonal dimension (sensor height) and f2 is the focal distance of lens L2. For simplicity, we will assume that the detector response is linear in intensity and that all other details may be represented by a proportionality parameter α < 1. Starting with Eqs. (2), (8) and (9), the detector response (photocurrent) will be proportional to the energy flux within dθ given by
The characteristic equation for the maximum of Eq. (11) is given by
and ζ = hc/єkT. To find a solution zmax, first we note that Eq. (12) can be written as
In the limit of high temperatures, ζ/z << 1, the root zmax approximates the root of Eq. (7)
In this limit, T is such that θmax is between 0 and π/2 (grazing reflectance). In general, the grating spectrograph scale will not give a simple relationship for Wien's law as in Eq. (16). Wien's peak for the spectrum given by Eq. (11) does not coincide with the wavelength peak as would be expected since the new scale introduces a nonlinear relationship for the wavelength. In other words, to get the right peak position, it is not sufficient to simply rescale the spectrograph according to Eq. (8). As a numerical example, we take a grating with 1000 grooves per millimeter so that λmax = 1 µm. The resulting BB spectrum is seen in Fig. (2). For T1 = 9000 K and T2 = 5500 K, the solutions of Eq. (11) are z1 = 5.0782 and z2 = 5.296 corresponding to θmax(T1) = 18.35º and θmax(T2) = 29.6º. Such values correspond to λmax(T1) = 314.7 nm and λmax(T2) = 495.6 nm. Again, for the sun's photosphere temperature, the maximum is in the green part of the spectrum for the grating spectrograph.
We can also scale the spectrograph to work in frequency. In this case, Eq. (10) will be equivalent to
The frequency is given by
If we multiply and divide Eq. (1) - the frequency representation - by ν2 we obtain Eq. (11) again.
Since the condition for the maxium depends on ζ, it is not possible to have a simple relationship as Eqs. (7) and (6) for Wien's law. As a function of T, solutions θmax will be given by the approximate formula
Eq. (20) is of order 10-2 and is valid for a 1000 grooves/ mm grating under the assumed approximations. Other gratings will have different laws with the same funcional dependence on temperature T.
3.2. Prism spectrograph
The experimental setup for a the prism spectrograph  is represented in Fig. 3. However, in this arrangement, the deviation angle is governed by refraction of light in a prism made of an optical material. Again, we assume a prism P with an internal angle A between faces. The light beam hits the entrance face with angle φ1 (in relation to the face normal). The total deviation angle will be given by
where φ2 is the refracted angle in relation to the exit face normal. Function φ2 depends in fact on the wavelength according to
with n(λ) being the prism dispersion relation. To get the material optical behavior, we use an approximate formula for Sellmeier's dispersion relation 
with a = 0.72, γ = 1.91 and λ0 = 0.085 µm for Quartz. The absolute error of this aproximation in the interval 200 nm to 1800 nm is less than 0.5% as compared to the 'exact' Sellmeier's relation. For λ in the interval 400 nm to 600 nm (corresponding to the visible region) n(λ) goes from 1.469 to 1.457. The relation between λ and the total deviation angle θ is more complicated than the equivalent Eq. (8) for the grating spectrograph, but it can be summarized by the following relations
We here write the BB spectrum distribution in the prism reference as determined by the dispersion relation Eq. (24). As usual, we start with the wavelength form Eq. (5) by noting that dθ = dΘ, we have
Fig. 4 is the numerical estimate of the BB spectrum distribution as given by Eq. (26) for a prism with A = 60º and entrance angle φ1 = 45º. To obtain the new peak position at a certain temperature T, one must again impose dIP /dΘ = 0 to Eq. (26). To avoid unnecessary complications, we provide numerical values only. As shown in Fig. 3, a BB peak for T = 5500 K is at θmax = 21.27º. Using the prism dispersion relation, Eq. (24), the corresponding wavelength is λmax = 1671.5 nm which is again in the infrared part of the spectrum.
The results obtained so far can be generalized. Suppose two functions f1(x) and f2(y) are given such that
together with the 'dispersion relation'
where g'(x) = dy/dx. Imposing the condition for the maximum at
that is clearly not zero, unless f2(g()) = 0 or g''() = 0 (that is, this is not the case for the relationship between frequency and wavelength). Suppose that we further make a transformation for a new variable z so that
We call z the 'experimental scale' (as in the case of frequency or wavelength scales). The following relations are valid
Then, in terms of f1(x)
In terms of f2(y)
Given Eq. (30), both ways of calculating the maximum of f3(z), Eqs. (36) and (37) lead to the same value. Eqs. (36) and (37) also show that the condition df3()/dz = 0 does not led to df2()/dy = 0 and f'1() = 0. The derivatives at the extremum = g() can be found by imposing df3/dz = 0 in Eqs. (36) and (37) and extracting both f'1() and df2()/dy
A somewhat more prosaic example is the search for the maximum position through a new 'gauge', say, y = x2, knowing that
To be more specific, suppose that a distribution of bar lengths was found to be in accordance to the normal rule of Eq. (38) with = µ. At this point, the maximum value of Eq. (38) is reached, that is df1()/dx = 0. In the new 'squared' scale, the distribution will be
which has as characteristic equation df2(y)/dy = 0
The distribution given by Eq. (39) diverges as y → 0. If, for example, µ = 10 and σ = 2, then + = 91.82 corresponds to the peak at the 'intuitive' maximum position = 100. Such maximum in y depends on the variance σ2, and only when σ → 0, → 2. If we disregard the region y << 1, for µ < 2σ, there is no real extrema for the distribuition of Eq. (39). At µ = 2σ,
there is only one extremum at = µ2/4. Such simple example shows that when the scale changes acording to some nonlinear function, the extrema chage position and are not described by the simple nonlinear relationship between scales.
Now, if we use another scale (z ) with x = 10z, the resulting distribution f3(z) will be
The characteristic equation, df3(z)/dz = 0, in this case is
Again, a real solution can be written as
For the numerical example given above (µ = 10 and σ = 2), = 1.016 which neither coincides with ( = 10) = 1, as expected by the x = 10z law, or ( = 91.82) = 0.981 as expected by the z = log( p z) relationship. Fig. 5 represents plots of functions f1(x), f2(y) and f3(z) for two distinct values of σ. For f2(y), the maximum position shifts to the left as σ increases, while the opposite happens to f3(z).
In general, the maxima will hardly coincide, although the transformation to a third 'experimental' function shows only one maximum (if there is only one extremum). In the same sense, the spectrum of Eqs. (2) and (1) can only be measured if an experimental device with linear scales for λ or ν becomes available.
This paper addresses and answers the question: 'where is the Wien's peak?' . The emitted power density maxima are derived from Planck's distribution law described either in frequency or wavelength intervals. We showed that such emitted power of a BB body at absolute temperature T is unique and does not coincide with the measured spectra (in fact we present theoretical estimates of measurable values). Estimated spectral shifts are caused by nonlinear frequency or wavelength 'gauge' relations to the experimentally accessible parameter (an intensity within an interval of such parameter). In general, there is no simple relation between the temperature and the maximum position (Wien's law) in the spectrograph scale given by this parameter. For prism spectrographs, the highly non linear dispersion relation complicates the task of obtaining a simple relationship between the maximum position and the temperature. We have found that the maximum position (in fact, the spectrum distribution) depends on the way that light interact with the dispersing element (i.e., prism or grating). The more linear the dispersion relation is, the closer we get to the Planck's distribution, Eq. (5).
All these conclusions are valid under ideal assumptions, which means to assume achromatic response for the lenses and sensing elements. If such assumptions are relaxed, then Eq. (11) will not be valid. The spectrum will be modulated by the transfer function of the combined spectrograph elements and the maximum will be found againg by imposing the condition dF(θ)/dθ = 0 where F(θ) is the output spectrum.
In spite of the nonlinear dependency of the BB radiation spectrum function on ν or λ, the shift observed for the peak can be easily demonstrated with other distribution fuctions under nonlinear scales.
 M. Planck, Annalen der Physik. 4, 553 (1901). [ Links ]
 O. Lummer and E. Prigsheim, Verh. d. D. Phys. Gesell. 1, 215 (1899). [ Links ]
 O. Lummer and E. Prigsheim, Verh. d. D. Phys. Gesell. 2, 163 (1900). [ Links ]
 W.H. Cropper, The Quantum Physicists (Oxford, New York, 1970). [ Links ]
 A. Sommerfeld, Thermodynamics and Statistical Mechanics - Lectures on Theoretical Physics, v. 5 (Ney York Academid Press, New York, 1956), Eq. (20a) on p. 144 [ Links ]
 R.K. Pathria, Statistical Mechanics (Butterworth- Heinemann, Oxford, 1996), 2nd ed., Eq. (8) on p. 170. [ Links ]
 C. Kittel, Elementary Statisticlal Physics (John Wiley and Sons, New York, 1958), Eq. (22.7) on p. 104. [ Links ]
 S.R.A. Salinas, Introduction to Statistical Physics (Springer, New York, 1997), Eq. (10.91) on p. 206. [ Links ]
 B.M. Smirnov, Principles of Statistical Physics - Distributions, Structures, Phenomena, Kinetics of Atomic Systems (Wiley VCH, Weinheim, 2006), Eq. (3.3) on p. 28. [ Links ]
 L.D. Landau and E.M. Lifshitz, Course of Theoretical Physics - Statistical Physics (Pergamon Press, Oxford, 1980), v. 5, Part I, 3rd ed., Eq. (64.9) on p. 184. [ Links ]
 F. Reif, Fundamentals of Statistical and Thermal Physics (McGraw-Hill Book Company, New York, 1963), Eq. (9.13.10) and (9.13.13) on p. 375 and p. 376, respectively. [ Links ]
 J. Jeans, The Dynamical Theory of Gases (Dover Publication, New York, 1925), 4th ed., Eq. (933) on p. 372 and Eq. (930) on p. 371. [ Links ]
 P.T. Landsberg, Thermodynamics and Statistical Mechanics (Dover Publications, New York, 1978), 2nd ed., p. 214, Section 13.3. [ Links ]
 C.E. Mendenhall and F.A. Saunders. Astrophys. Journal, 13, 25 (1901). [ Links ]
 I. Kaplan, Nuclear Physics (Addison-Wesley Pub., New York, 1964), 2nd ed. [ Links ]
 M.A. Heald, Am. J. Phys. 71, 1322 (2003). [ Links ]
 J.M. Overduin,Am. J. Phys.71, 216 (2003). [ Links ]
 D.J. Fixsen, E.S. Cheng, J.M. Gales, J.C. Mather, R.A. Shafer and E.L. Wright, Astrophysical Journal 473, 576 (1996). [ Links ]
 R.C. Jones, J. Opt. Soc. Am. 43, 1 (1953). [ Links ]
 W. Sellmeir, Ann. d. Phy. u. Ch. 219, 271 (1871). [ Links ]
 J. James, Spectrograph Design Fundamentals (Cambridge Univ. Press, Cambridge, 2007). [ Links ]
 A.H. Wilson, Thermodynamics and Statistical Mechanics (Cambridge University Press, Cambridge, 1957), Eq. (6.52.3) on p. 184. [ Links ]
 C. Palmer. Diffraction grating Handbook (Thermo RGL, Rochester, 2002), 5th ed. [ Links ]
 J. Ma, J. Yang and J. Nie. Lat. Am. J. Phys. Educ. 3, 566 (2009). [ Links ]
Recebido em 3/8/2011; Aceito em 8/8/2011; Publicado em 20/4/2012