Differences

This shows you the differences between two versions of the page.

--- un:blackbody-radiation [2025/11/29 10:10] – [Blackbody spectra] asad
+++ un:blackbody-radiation [2025/12/06 10:58] (current) – asad
@@ Line 1: / Line 1: @@
 ====== Blackbody radiation ======
-Blackbody radiation is the equilibrium electromagnetic field that arises when photons and matter interact frequently enough to share a common temperature. In such environments, photons scatter repeatedly and exchange energy with particles until the entire system reaches thermal equilibrium. Once equilibrium is established, the radiation field depends only on temperature, not on the material composition of the emitting region. Stellar interiors, deep layers of stellar atmospheres, and the early Universe before photon decoupling all approximate these conditions.
+Blackbody radiation refers to the electromagnetic radiation emitted by a non-reflecting body held at a constant temperature. It represents a state of thermal equilibrium between matter and radiation, often realized in astrophysics within "optically thick" environments such as stellar interiors or the early universe. This radiation field is isotropic and unpolarized, and its properties depend exclusively on the temperature $T$ of the source.
-When matter and radiation are in equilibrium, the average photon energy matches the characteristic kinetic energy of nonrelativistic particles. This condition can be expressed as
+===== Specific intensity =====
-$$
+The fundamental descriptor of this radiation is the specific intensity, denoted as $I(\nu, T)$. This quantity measures the radiant power passing through a unit area, per unit solid angle, per unit frequency interval. For a blackbody, the specific intensity is universal and is described by the Planck function.
-\frac{3}{2} kT = \left( \frac{1}{2} mv^2 \right)_{av} = h \nu_{\mathrm{av}}.
-$$
-A hypothetical container with walls and gas at temperature \( T \) would produce this equilibrium spectrum exactly. Although astrophysical environments are more complex than an ideal cavity, sufficiently high density and optical thickness allow photons to reach the same statistical equilibrium state.
-===== Planck spectrum =====
-The distribution of blackbody radiation follows from the Bose–Einstein statistics of photons. The photon distribution function describes how many photons occupy a phase-space cell of volume \( h^3 \) and is given by
-$$
-f = \frac{2}{h^3} \frac{1}{e^{h\nu/(kT)} - 1}.
-$$
-From this distribution one obtains the specific intensity,
-$$
-I_\nu(T) = \frac{2 h \nu^3}{c^2} \frac{1}{e^{h\nu/(kT)} - 1},
-$$
-which is known as the Planck function. It gives the radiant energy per unit area, time, solid angle, and frequency interval. The spectrum rises as \( \nu^3 \) at low frequencies, reaches a well-defined maximum, and then declines rapidly at high frequencies. Expressed in wavelength units, the Planck function becomes
-$$
-I_\lambda(T) = \frac{2 h c^2}{\lambda^5} \frac{1}{e^{hc/(\lambda kT)} - 1}.
-$$
-While these two forms describe the same radiation field, their peaks occur at different numerical locations because frequency and wavelength depend on each other nonlinearly.
-===== Rayleigh–Jeans and Wien regimes =====
-At low frequencies where \( h\nu \ll kT \), the exponential in the denominator can be expanded using a Taylor series:
-$$
-e^{h\nu/(kT)} \approx 1 + \frac{h\nu}{kT}.
-$$
-Substituting this into the Planck function yields the Rayleigh–Jeans approximation,
-$$
-I_\nu(T) \approx \frac{2 \nu^2 kT}{c^2},
-$$
-which grows quadratically with frequency and linearly with temperature. This form applies to radio and microwave wavelengths where photon energies are much smaller than the thermal energy scale.
-At high frequencies where \( h\nu \gg kT \), the exponential term dominates and Planck’s law reduces to the Wien approximation,
-$$
-I_\nu(T) \approx \frac{2 h \nu^3}{c^2} e^{-h\nu/(kT)}.
-$$
-This approximation describes the exponential falloff at short wavelengths. The location of the maximum of \( I_\nu(T) \) is determined by setting the derivative with respect to frequency equal to zero, which gives
-$$
-\nu_{\mathrm{peak}} = 5.88 \times 10^{10} T.
-$$
-Thus the peak frequency scales directly with temperature.
-===== Displacement law and peak =====
-Because \( I_\nu \) and \( I_\lambda \) have different functional forms, the peak of the wavelength spectrum does not occur at \( \lambda = c/\nu_{\mathrm{peak}} \). Differentiating the wavelength form of the Planck function yields the Wien displacement law,
-$$
-T \lambda_{\mathrm{peak}} = 2.898 \times 10^{-3} \ \mathrm{m\,K}.
-$$
-Increasing the temperature therefore shifts the peak emission to shorter wavelengths. Cold sources of a few kelvin peak in the millimeter regime, while hotter bodies with temperatures of thousands of kelvin peak in the visible or ultraviolet.
-===== Flux and luminosity =====
-Integrating the Planck function over all frequencies gives the total radiation energy density. The number density of photons varies as \( T^3 \), while the radiation energy density varies as \( T^4 \). The total flux emerging from a unit surface area is given by the Stefan–Boltzmann law,
-\begin{equation}
-F = \sigma T^4.
-\end{equation}
-A spherical star of radius \( R \) radiates a luminosity
-$$
-L = 4\pi R^2 \sigma T^4.
-$$
-This connects the effective temperature of a radiating surface to its total output of energy and allows estimates of stellar radii and luminosities when the flux and temperature are known.
-===== Blackbody spectra =====
-When blackbody spectra are plotted as functions of frequency, \(I(\nu)\), several consistent trends appear. At low frequencies the curves rise approximately as \(\nu^2\), following the Rayleigh–Jeans form that emerges from the classical limit of the Planck law. At high frequencies the intensity drops exponentially, producing the Wien tail whose slope steepens rapidly with increasing \(\nu\). As the temperature increases, the entire curve shifts upward and the peak moves systematically to higher frequencies, illustrating the strong temperature dependence of the radiation field. The total emitted power scales as \(T^4\), so hotter blackbodies produce far greater integrated intensity even if their spectral shapes remain similar.
 <html>
@@ Line 97: / Line 12: @@
 #bb_nu, #bb_lambda {
     width: 100%;
-    height: 420px;
+    height: 440px;
     position: relative;
     margin-bottom: 1.2em;
@@ Line 117: / Line 32: @@
 <div style="margin-top:0.3em; text-align:center; font-family:sans-serif;">
   Temperature:
-  <input type="range" id="Tslider" min="1000" max="50000" step="100" value="5800"
+  <input type="range" id="Tslider" min="200" max="20000" step="100" value="5800"
          style="width:400px; accent-color:#cc0000;">
   <span id="Tval" style="margin-left:6px; font-weight:bold;">5800</span> K
@@ Line 123: / Line 38: @@
 <script type="text/javascript">
-  window.PlotlyConfig = {MathJaxConfig: 'local'};
+window.PlotlyConfig = {MathJaxConfig: 'local'};
 </script>
 <script src="https://cdn.plot.ly/plotly-latest.min.js"></script>
@@ Line 129: / Line 44: @@
 <script type="text/javascript">
 //<![CDATA[
-const h=6.626e-34,k=1.381e-23,c=3e8,N=400;
+const h=6.626e-34, k=1.381e-23, c=3e8;
-const Bnu=(n,T)=>(2*h*n**3/c**2)/(Math.exp(h*n/(k*T))-1);
-const Blambda=(l,T)=>(2*h*c**2/l**5)/(Math.exp(h*c/(l*k*T))-1);
+const N=400;
 const logspace=(a,b,n)=>Array.from({length:n},(_,i)=>10**(a+(b-a)*i/(n-1)));
-const nu=logspace(11,15,N),lambda=logspace(-7,-4,N);
-const makeNuTrace=T=>({x:nu,y:nu.map(n=>Bnu(n,T)),mode:'lines',line:{color:'#cc0000'}});
+const nu_Hz  = logspace(11,17,N);
-const makeLambdaTrace=T=>({x:lambda,y:lambda.map(l=>Blambda(l,T)),mode:'lines',line:{color:'#0044cc'}});
+const lambda_m  = logspace(-8,-3,N);
+const Bnu = (nu,T)=>(2*h*nu**3/c**2)/(Math.exp(h*nu/(k*T))-1);
+const Blambda = (lam,T)=>(2*h*c**2/lam**5)/(Math.exp(h*c/(lam*k*T))-1);
+const makeNuTrace = T => ({
+  x: nu_Hz,
+  y: nu_Hz.map(n => Bnu(n,T)),
+  mode:'lines',
+  line:{color:'#cc0000'}
+});
+const makeLambdaTrace = T => ({
+  x: lambda_m,
+  y: lambda_m.map(l => Blambda(l,T)),
+  mode:'lines',
+  line:{color:'#0044cc'}
+});
+// BIG equation text
+const eq_nu = {
+  x: 0.97, y: 0.97, xref:'paper', yref:'paper',
+  text:'$I(\\nu)=\\frac{2h\\nu^{3}}{c^{2}}\\,\\frac{1}{e^{h\\nu/kT}-1}$',
+  showarrow:false,
+  font:{size:24},
+  align:'right'
+};
+const eq_lambda = {
+  x: 0.97, y: 0.97, xref:'paper', yref:'paper',
+  text:'$I(\\lambda)=\\frac{2hc^{2}}{\\lambda^{5}}\\,\\frac{1}{e^{hc/(\\lambda kT)}-1}$',
+  showarrow:false,
+  font:{size:24},
+  align:'right'
+};
-const layout_nu={
+const layout_nu = {
   margin:{l:95,r:20,t:10,b:60},
-  xaxis:{type:'log',title:'Frequency (Hz)',tickformat:'.1e'},
+  xaxis:{title:'Frequency (Hz)', type:'log'},
-  yaxis:{
+  yaxis:{title:'I(ν) [W m⁻² Hz⁻¹ sr⁻¹]', automargin:true, titlefont:{size:18}},
-    type:'log',
+  annotations:[eq_nu]
-    tickformat:'.1e',
-    title:'I(ν)   [W m⁻² Hz⁻¹ sr⁻¹]',
-    titlefont:{size:18},
-    automargin:true
-  }
 };
-const layout_lambda={
+const layout_lambda = {
   margin:{l:95,r:20,t:10,b:60},
-  xaxis:{type:'log',title:'Wavelength (m)',tickformat:'.1e'},
+  xaxis:{title:'Wavelength (m)', type:'log'},
-  yaxis:{
+  yaxis:{title:'I(λ) [W m⁻² m⁻¹ sr⁻¹]', automargin:true, titlefont:{size:18}},
-    type:'log',
+  annotations:[eq_lambda]
-    tickformat:'.1e',
-    title:'I(λ)   [W m⁻² m⁻¹ sr⁻¹]',
-    titlefont:{size:18},
-    automargin:true
-  }
 };
-Plotly.newPlot('bb_nu',[makeNuTrace(5800)],layout_nu);
+Plotly.newPlot('bb_nu', [makeNuTrace(5800)], layout_nu);
-Plotly.newPlot('bb_lambda',[makeLambdaTrace(5800)],layout_lambda);
+Plotly.newPlot('bb_lambda', [makeLambdaTrace(5800)], layout_lambda);
+const slider=document.getElementById('Tslider');
+const Tval=document.getElementById('Tval');
-const slider=document.getElementById('Tslider'),Tval=document.getElementById('Tval');
 slider.oninput=function(){
-  const T=+this.value;
+  const T = +this.value;
-  Tval.textContent=T;
+  Tval.textContent = T;
-  Plotly.react('bb_nu',[makeNuTrace(T)],layout_nu);
+  Plotly.react('bb_nu', [makeNuTrace(T)], layout_nu);
-  Plotly.react('bb_lambda',[makeLambdaTrace(T)],layout_lambda);
+  Plotly.react('bb_lambda', [makeLambdaTrace(T)], layout_lambda);
 };
 //]]>
@@ Line 177: / Line 117: @@
 </html>
-In the wavelength representation, \(I(\lambda)\), the same underlying physics appears but with an important difference: the peak occurs at a different location from the frequency plot. Because frequency and wavelength are related by \(\nu=c/\lambda\), equal intervals in frequency do not correspond to equal intervals in wavelength. As a result, the wavelength-domain curve rises steeply at long wavelengths and falls sharply at short wavelengths, but the maximum occurs at a longer wavelength than one might infer from the frequency-domain peak. This difference directly reflects the Jacobian factor between the two representations, since \(I(\nu)\,d\nu \neq I(\lambda)\,d\lambda\). Wien’s displacement law, \(\lambda_{\max}T \approx 2.9\times10^{-3}\,\mathrm{m\,K}\), becomes particularly clear in the wavelength plot: increasing the temperature shifts the peak systematically toward shorter wavelengths across the visible and ultraviolet regions. Together, the two plots provide complementary perspectives on the same blackbody radiation field and illustrate how spectral peaks depend on the choice of independent variable.
+The derivation of the blackbody spectrum relies on the Bose-Einstein distribution function, which describes the statistical behavior of integer-spin particles (bosons) such as photons. In a phase space defined by position coordinates $\mathbf{x}$ and momentum vectors $\mathbf{p}$, the distribution function $f(\mathbf{x}, \mathbf{p})$ represents the occupation number of quantum states. For photons in thermal equilibrium, this function is given by:
-===== Insights =====
+$$f = \frac{2}{h^3} \frac{1}{e^{E/kT} - 1}$$
-  * Blackbody radiation arises when photons and matter reach full thermal equilibrium, producing a spectrum that depends only on temperature.
+Here, $h$ is Planck's constant, $E = h\nu$ is the energy of the photon, and the factor of 2 arises from the two independent polarization states of the photon. The factor $h^3$ represents the volume of a single cell in six-dimensional phase space ($d^3x d^3p$).
-  * The Planck function determines the intensity at each frequency and follows from the photon distribution in phase space.
-  * The Rayleigh–Jeans and Wien limits describe the low- and high-frequency behavior of the spectrum and explain its characteristic shape.
+The specific intensity $I(\nu, T)$ is directly related to this phase-space density. Since photons travel at speed $c$, the intensity is the energy carried by photons passing through a surface, which equates to the product of the energy per photon $h\nu$, the phase space density $f$, the factor $c/4\pi$ related to isotropic solid angle integration, and the density of states factor $p^2/h^3$. Specifically, the Planck function relates to $f$ via the density of states in frequency space, yielding:
-  * The Wien displacement law shows that higher temperatures shift the peak emission to shorter wavelengths, linking temperature to observed colors.
+$$I(\nu, T) = \frac{2h\nu^3}{c^2} \frac{1}{e^{h\nu/kT} - 1}$$
-  * The Stefan–Boltzmann law connects temperature to total radiative flux, enabling luminosities and radii of stars to be inferred.
+This equation forms the basis for all subsequent characteristics of blackbody radiation.
-  * Dense, optically thick astrophysical regions such as stellar atmospheres and the early Universe closely approximate blackbody conditions.
-  * Blackbody sources serve as calibration standards because their emission depends solely on temperature and follows universal laws.
+===== Rayleigh-Jeans approximation =====
+In the low-frequency limit, where the photon energy is much smaller than the thermal energy ($h\nu \ll kT$), the exponential term in the denominator of the Planck function can be approximated using the Taylor series expansion $e^x \approx 1 + x$. Substituting this into the specific intensity equation cancels the $h$ terms, resulting in the Rayleigh-Jeans law:
+$$I(\nu, T) \approx \frac{2\nu^2 kT}{c^2}$$
+This approximation was historically significant as it corresponds to the classical prediction that failed at high frequencies (the "ultraviolet catastrophe"). In this regime, the intensity is directly proportional to the temperature $T$, a property extensively used in radio astronomy. Radio telescopes measure the power received and convert it into a "brightness temperature" or "antenna temperature" using this linear relation, providing a convenient metric for source intensity even for non-thermal emitters.
+===== Wien approximation =====
+At the opposite end of the spectrum, where frequencies are high and photon energy exceeds the thermal kinetic energy ($h\nu \gg kT$), the exponential term $e^{h\nu/kT}$ becomes very large. The $-1$ in the denominator becomes negligible, allowing the Planck function to be approximated as:
+$$I(\nu, T) \approx \frac{2h\nu^3}{c^2} e^{-h\nu/kT}$$
+This is the Wien approximation. It describes the steep cutoff of the spectrum in the ultraviolet and X-ray regimes for typical stellar temperatures. The exponential factor dominates the $\nu^3$ term, ensuring that the total energy integrated over all frequencies remains finite, thus resolving the divergence predicted by classical physics.
+===== Wien's displacement law (frequency) =====
+The spectral distribution of blackbody radiation peaks at a specific frequency $\nu_{peak}$ that shifts with temperature. Differentiating the Planck function with respect to $\nu$ and solving for the maximum yields a transcendental equation, the solution of which shows that the peak frequency is linearly proportional to temperature:
+$$h\nu_{peak} = 2.82 kT$$
+This relationship implies that hotter objects emit the bulk of their energy at higher frequencies. For example, the Cosmic Microwave Background (CMB) at 2.73 K peaks at approximately 160 GHz, whereas the Sun, at roughly 5800 K, peaks in the visible spectrum.
+===== Wien's displacement law (wavelength) =====
+When the specific intensity is expressed per unit wavelength interval, $I_\lambda(\lambda, T)$, the condition for the spectral peak changes due to the nonlinear relationship between frequency and wavelength differentials ($d\nu = -c/\lambda^2 d\lambda$). The peak wavelength $\lambda_{peak}$ follows the familiar form of Wien's displacement law:
+$$\lambda_{peak} T = 2.898 \times 10^{-3} \text{ K m}$$
+It is important to recognize that the peak of the spectrum in wavelength space does not correspond to the same photon energy as the peak in frequency space. The wavelength peak occurs at a frequency roughly 1.76 times higher than the frequency peak derived from $I(\nu, T)$.
+===== Stefan-Boltzmann law =====
+The total radiant flux $\mathcal{F}$ emitted from the surface of a blackbody is obtained by integrating the specific intensity over all frequencies ($0 \to \infty$) and over the outward-facing hemisphere of solid angles. The integration over solid angle contributes a factor of $\pi$, and the frequency integral yields a dependence on the fourth power of temperature:
+$$\mathcal{F} = \sigma T^4$$
+where $\sigma \approx 5.67 \times 10^{-8} \text{ W m}^{-2} \text{ K}^{-4}$ is the Stefan-Boltzmann constant. For a spherical star of radius $R$, this flux definition leads to the luminosity equation $L = 4\pi R^2 \sigma T^4$. This law defines the effective temperature of stars, relating their total intrinsic power output to their physical size.
+===== Spectral energy density =====
+The spectral energy density $u_\nu(\nu, T)$ measures the amount of radiant energy contained within a unit volume of space per unit frequency interval. For isotropic radiation, the energy density is related to the specific intensity by a factor of $4\pi/c$. This accounts for radiation traveling in all directions within the volume:
+$$u_\nu(\nu, T) = \frac{4\pi}{c} I(\nu, T) = \frac{8\pi h \nu^3}{c^3} \frac{1}{e^{h\nu/kT} - 1}$$
+This function describes the "photon gas" within a cavity. In standard cosmology, this represents the energy density of the CMB photons that permeate the universe.
+===== Total energy density =====
+Integrating the spectral energy density over the entire frequency range yields the total energy density $u(T)$, representing the total joules of radiation energy per cubic meter. Like the flux, this quantity scales with the fourth power of the temperature:
+$$u(T) = a T^4$$
+The radiation constant $a$ is related to the Stefan-Boltzmann constant by $a = 4\sigma/c$. This relationship is critical in cosmology and stellar interior physics, as it governs the energy content of the radiation field which can dominate over matter energy density at sufficiently high temperatures.
+===== Spectral number density =====
+The spectral number density $n_\nu(\nu, T)$ defines the number of photons per unit volume per unit frequency interval. It is derived by dividing the spectral energy density $u_\nu$ by the energy of a single photon, $h\nu$:
+$$n_\nu(\nu, T) = \frac{u_\nu(\nu, T)}{h\nu} = \frac{8\pi \nu^2}{c^3} \frac{1}{e^{h\nu/kT} - 1}$$
+Physically, this can be interpreted as the product of the density of quantum states in phase space (proportional to $\nu^2$) and the Bose-Einstein occupation number (the exponential term). This distribution determines the number of photons available to interact with matter at a specific energy.
+===== Total number density =====
+The total number density of photons $n(T)$ is found by integrating $n_\nu$ over all frequencies. The result depends on the third power of the temperature, rather than the fourth:
+$$n(T) \propto T^3$$
+Specifically, $n(T) \approx 2.03 \times 10^7 T^3 \text{ photons m}^{-3}$. This scaling has important cosmological implications; as the universe expands and temperature drops, the number density of photons decreases effectively as $1/R^3$ (where $R$ is the scale factor), conserving the total number of photons in a comoving volume, whereas the energy density drops as $1/R^4$ due to the additional redshift of photon energy.
+===== Average photon energy =====
+The average energy of a photon in a blackbody radiation field is the ratio of the total energy density to the total number density ($u/n$). Since $u \propto T^4$ and $n \propto T^3$, the average energy is directly proportional to $T$:
+$$h\nu_{avg} \approx 2.70 kT$$
+This value provides a useful rule of thumb for characteristic photon energies. It is slightly lower than the energy at the frequency peak ($2.82 kT$) but provides a representative energy for interactions such as ionization or scattering within a thermal plasma.
 ===== Inquiries =====
-  * How does the equilibrium distribution of photons lead to the mathematical form of the Planck function?
+  - Demonstrate how the Planck function reduces to the Rayleigh-Jeans approximation in the low-frequency limit ($h\nu \ll kT$) using the Taylor series expansion $e^x \approx 1+x$.
-  * Why do the Rayleigh–Jeans and Wien approximations follow from the limiting behavior of the exponential term in Planck’s law?
+  - Explain why the spectral peak in frequency space ($\nu_{peak}$) corresponds to a different photon energy than the spectral peak in wavelength space ($\lambda_{peak}$), specifically referencing the relationship between the differentials $d\nu$ and $d\lambda$.
-  * How does expressing the spectrum in frequency versus wavelength naturally lead to different peak locations?
+  - Using the scaling relations provided, determine the factor by which the total luminosity of a spherical star increases if its temperature doubles while its radius remains constant.
+  - Contrast the temperature dependence of the total photon energy density ($u \propto T^4$) with that of the total photon number density ($n \propto T^3$) and explain what this implies for the average energy per photon as temperature changes.
+  - State the equation of state relating radiation pressure $P$ to total energy density $u$ for an isotropic blackbody field and identify the factor distinguishing it from the kinetic pressure of a non-relativistic gas.