Differences

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

--- un:blackbody-radiation [2025/11/30 07:44] – [Bose–Einstein occupation curves] asad
+++ un:blackbody-radiation [2025/12/06 10:58] (current) – asad
@@ Line 1: / Line 1: @@
-====== Bose–Einstein statistics ======
+====== Blackbody radiation ======
-Bose–Einstein statistics describes the equilibrium distribution of indistinguishable bosons, particles with integer spin that are not subject to the Pauli exclusion principle. Because any number of bosons may occupy a single quantum state, their statistical behavior differs profoundly from that of fermions. Photons, phonons, and ultracold bosonic atoms are the most important physical examples.
+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.
-In equilibrium, the distribution of bosons in phase space depends only on temperature and (where allowed) chemical potential. For photons in thermal radiation fields, chemical equilibrium with matter forces the chemical potential to vanish, \( \mu = 0 \). For massive bosons, \( \mu \) varies with temperature and density but always satisfies \( \mu \le \epsilon_{0} \), the ground-state energy.
+===== Specific intensity =====
-The fundamental description of a bosonic gas begins not with the occupation of discrete energy levels, but with the **phase-space distribution function**, which specifies how many particles occupy each differential volume of phase space.
+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.
-===== Bose–Einstein phase-space distribution =====
+<html>
-The statistical mechanics of bosons gives the number of particles per phase-space cell of volume \( h^{3} \) as
+<style>
+#bb_nu, #bb_lambda {
+    width: 100%;
+    height: 440px;
+    position: relative;
+    margin-bottom: 1.2em;
+}
+#bb_nu .plot-container,
+#bb_nu .svg-container,
+#bb_nu svg.main-svg,
+#bb_lambda .plot-container,
+#bb_lambda .svg-container,
+#bb_lambda svg.main-svg {
+    width: 100% !important;
+    height: 100% !important;
+}
+</style>
-$$
+<div id="bb_nu"></div>
-f = \frac{2}{h^{3}} \, \frac{1}{e^{h\nu/(kT)} - 1}.
+<div id="bb_lambda"></div>
-$$
-Here the factor of 2 accounts for two photon polarization states. For massive bosons, this factor is replaced by the appropriate spin degeneracy. This function expresses a key difference between classical and quantum statistics: the denominator can become small, allowing **large occupation of low-energy states**.
+<div style="margin-top:0.3em; text-align:center; font-family:sans-serif;">
+  Temperature:
+  <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
+</div>
-For photons, the variable \( \nu \) labels frequency, and this same distribution underlies the Planck spectrum of [[Blackbody radiation]]. For massive bosons, the energy is instead \( \epsilon = p^{2}/2m \), but the statistical form is identical:
+<script type="text/javascript">
+window.PlotlyConfig = {MathJaxConfig: 'local'};
+</script>
+<script src="https://cdn.plot.ly/plotly-latest.min.js"></script>
-$$
+<script type="text/javascript">
-f = \frac{g}{h^{3}} \frac{1}{e^{(\epsilon-\mu)/(kT)} - 1},
+//<![CDATA[
-$$
+const h=6.626e-34, k=1.381e-23, c=3e8;
-with \( g \) the spin degeneracy. The form of \( f \) therefore unifies the physics of blackbody photons and ideal Bose gases.
+const N=400;
+const logspace=(a,b,n)=>Array.from({length:n},(_,i)=>10**(a+(b-a)*i/(n-1)));
-===== Mean occupation of quantum states =====
+const nu_Hz  = logspace(11,17,N);
+const lambda_m  = logspace(-8,-3,N);
-The phase-space distribution function determines the **mean number of bosons occupying each single-particle quantum state**. Integrating \(f\) over a cell associated with the state \(i\) leads directly to the Bose–Einstein occupation number:
+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 => ({
-\bar{n}_i = \frac{1}{\exp\!\left(\frac{\epsilon_i - \mu}{kT}\right) - 1}.
+  x: nu_Hz,
-$$
+  y: nu_Hz.map(n => Bnu(n,T)),
+  mode:'lines',
+  line:{color:'#cc0000'}
+});
-This quantity plays the same role for massive bosons that the photon distribution function plays for radiation. Summing over all states,
+const makeLambdaTrace = T => ({
+  x: lambda_m,
+  y: lambda_m.map(l => Blambda(l,T)),
+  mode:'lines',
+  line:{color:'#0044cc'}
+});
-$$
+// BIG equation text
-N = \sum_i \bar{n}_i,
+const eq_nu = {
-\qquad
+  x: 0.97, y: 0.97, xref:'paper', yref:'paper',
-E = \sum_i \epsilon_i \bar{n}_i.
+  text:'$I(\\nu)=\\frac{2h\\nu^{3}}{c^{2}}\\,\\frac{1}{e^{h\\nu/kT}-1}$',
-$$
+  showarrow:false,
+  font:{size:24},
+  align:'right'
+};
-In the continuum limit with density of states \(g(\epsilon)\),
+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 = {
-n = \frac{N}{V}
+  margin:{l:95,r:20,t:10,b:60},
-  = \int_0^\infty \frac{g(\epsilon)}{\exp[(\epsilon-\mu)/(kT)] - 1} \, d\epsilon.
+  xaxis:{title:'Frequency (Hz)', type:'log'},
-$$
+  yaxis:{title:'I(ν) [W m⁻² Hz⁻¹ sr⁻¹]', automargin:true, titlefont:{size:18}},
+  annotations:[eq_nu]
+};
-The crucial feature is the “\(-1\)” in the denominator; its absence in the classical Boltzmann limit suppresses large occupation at low energy.
+const layout_lambda = {
+  margin:{l:95,r:20,t:10,b:60},
+  xaxis:{title:'Wavelength (m)', type:'log'},
+  yaxis:{title:'I(λ) [W m⁻² m⁻¹ sr⁻¹]', automargin:true, titlefont:{size:18}},
+  annotations:[eq_lambda]
+};
-===== Classical limit and quantum degeneracy =====
+Plotly.newPlot('bb_nu', [makeNuTrace(5800)], layout_nu);
+Plotly.newPlot('bb_lambda', [makeLambdaTrace(5800)], layout_lambda);
-When the exponent is large,
+const slider=document.getElementById('Tslider');
+const Tval=document.getElementById('Tval');
-$$
+slider.oninput=function(){
-e^{(\epsilon-\mu)/(kT)} \gg 1,
+  const T = +this.value;
-$$
+  Tval.textContent = T;
+  Plotly.react('bb_nu', [makeNuTrace(T)], layout_nu);
+  Plotly.react('bb_lambda', [makeLambdaTrace(T)], layout_lambda);
+};
+//]]>
+</script>
-the Bose–Einstein distribution reduces to the Maxwell–Boltzmann form,
+</html>
-$$
+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:
-\bar{n}_i \approx \exp\!\left(-\frac{\epsilon_i-\mu}{kT}\right),
+$$f = \frac{2}{h^3} \frac{1}{e^{E/kT} - 1}$$
-$$
+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$).
-and the phase-space distribution becomes
+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:
+$$I(\nu, T) = \frac{2h\nu^3}{c^2} \frac{1}{e^{h\nu/kT} - 1}$$
+This equation forms the basis for all subsequent characteristics of blackbody radiation.
-$$
+===== Rayleigh-Jeans approximation =====
-f \approx \frac{g}{h^{3}} e^{-(\epsilon-\mu)/(kT)}.
-$$
-Quantum statistical effects appear when the thermal de Broglie wavelength,
+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 =====
-\lambda_{\mathrm{th}}
-= \sqrt{\frac{2\pi\hbar^{2}}{m kT}},
-$$
-becomes comparable to the mean particle spacing. The transition from classical to quantum behavior occurs when
+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) =====
-n\lambda_{\mathrm{th}}^{3} \gtrsim 1.
-$$
-Below this threshold, the large-occupation low-energy behavior of bosons becomes unavoidable.
+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.
-===== Bose–Einstein condensation =====
+===== Wien's displacement law (wavelength) =====
-For a uniform ideal Bose gas, the chemical potential satisfies \( \mu \le \epsilon_{0} \). As the gas is cooled at fixed density, \( \mu \rightarrow \epsilon_{0} \), and the excited states are no longer capable of holding all particles. The excess particles accumulate in the ground state, producing macroscopic occupation:
+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 =====
-\bar{n}_{0} \rightarrow \text{large}.
-$$
-This is **Bose–Einstein condensation**, a striking manifestation of the Bose–Einstein distribution function. The condensation phenomenon is therefore the direct consequence of the form of \(f\) and the associated occupation numbers \(\bar{n}_i\): large phase-space density implies large ground-state occupation.
+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.
-In contrast, for photons the chemical potential is fixed to zero, and instead of condensing, photons redistribute themselves to form the blackbody spectrum.
+===== Spectral energy density =====
-===== Bose–Einstein occupation curves =====
+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.
-To analyze how the Bose–Einstein distribution behaves, it is useful to define dimensionless variables,
+===== 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:
-x = \frac{\epsilon}{kT}, \qquad
+$$u(T) = a T^4$$
-a = -\frac{\mu}{kT},
+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.
-$$
-so that
+===== 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$:
-\bar{n}(x; a) = \frac{1}{e^{x+a} - 1}.
+$$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.
-Small \(a\) (i.e., large \( \mu \)) enhances the low-energy peak of the distribution. The special case \(a=0\) corresponds exactly to the photon distribution function above.
+===== Total number density =====
-<html>
+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.
-<style>
+===== Average photon energy =====
-#be_occ {
-    width: 100%;
-    height: 440px;
-    position: relative;
-    margin-bottom: 1.2em;
-}
-#be_occ .plot-container,
-#be_occ .svg-container,
-#be_occ svg.main-svg {
-    width: 100% !important;
-    height: 100% !important;
-}
-</style>
-<div id="be_occ"></div>
+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.
-<div style="margin-top:0.3em; text-align:center; font-family:sans-serif;">
+===== Inquiries =====
-  Parameter \(a = -\mu/(kT)\):
+  - 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$.
-  <input type="range" id="aslider" min="0.0" max="3.0" step="0.1" value="1.0"
+  - 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$.
-         style="width:400px;">
+  - 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.
-  <span id="Aval" style="margin-left:6px; font-weight:bold;">1.0</span>
+  - 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.
-</div>
+  - 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.
-<script type="text/javascript">
-window.PlotlyConfig = {MathJaxConfig: 'local'};
-</script>
-<script src="https://cdn.plot.ly/plotly-latest.min.js"></script>
-<script type="text/javascript">
-//<![CDATA[
-const Npts = 400;
-const xMin = 0.05;
-const xMax = 8.0;
-const x_vals = Array.from({length:Npts},(_,i)=>xMin+(xMax-xMin)*i/(Npts-1));
-function nbar(x,a){
-  const e = Math.exp(x+a);
-  return (e<=1)?Infinity:1/(e-1);
-}
-function makeBETrace(a){
-  return {
-    x:x_vals,
-    y:x_vals.map(x=>nbar(x,a)),
-    mode:'lines',
-    line:{color:'#cc0000'}
-  };
-}
-// CHANGED: plain-text equation (no MathJax escapes)
-const eq_be = {
-  x:0.97, y:0.97, xref:'paper', yref:'paper',
-  text:'n̄(x; a) = 1 / (exp(x + a) - 1)',
-  showarrow:false, font:{size:24}, align:'right'
-};
-const layout_be = {
-  margin:{l:95,r:20,t:10,b:60},
-  // your fixed xlabel – left untouched
-  xaxis:{title:'$x = (\\epsilon/(kT))$', automargin:true},
-  // CHANGED: plain-text ylabel (no MathJax)
-  yaxis:{title:'n̄(x; a)', type:'log', automargin:true},
-  annotations:[eq_be]
-};
-Plotly.newPlot('be_occ',[makeBETrace(1.0)],layout_be);
-const aslider=document.getElementById('aslider');
-const Aval=document.getElementById('Aval');
-aslider.oninput=function(){
-  const a=parseFloat(this.value);
-  Aval.textContent=a.toFixed(1);
-  Plotly.react('be_occ',[makeBETrace(a)],layout_be);
-};
-//]]>
-</script>
-</html>
-===== Insights =====
-  * The phase-space distribution \( f = g h^{-3} [\,e^{(\epsilon-\mu)/(kT)} - 1\,]^{-1} \) is the foundational expression from which all Bose–Einstein behavior follows.
-  * Large low-energy occupation arises because the denominator may become small as \( \mu \to \epsilon_{0} \).
-  * Classical statistics emerges only when \( e^{(\epsilon-\mu)/(kT)} \gg 1 \), suppressing large occupation.
-  * Quantum degeneracy appears when \( n\lambda_{\mathrm{th}}^{3} \gtrsim 1 \), forcing the use of Bose–Einstein statistics.
-  * Bose–Einstein condensation results directly from the inability of excited states to hold all particles as \( \mu \to \epsilon_0 \).
-  * For photons, the condition \( \mu = 0 \) transforms the Bose–Einstein distribution into the Planck spectrum of blackbody radiation.
-===== Inquiries =====
-  * Derive the phase-space distribution \( f = g h^{-3} [e^{(\epsilon-\mu)/(kT)} - 1]^{-1} \) from the grand canonical ensemble.
-  * Show how the occupation number \( \bar{n}_i \) emerges from integrating \(f\) over a state’s phase-space cell.
-  * Demonstrate explicitly how the Maxwell–Boltzmann limit arises from the Bose–Einstein distribution.
-  * Using \( n\lambda_{\mathrm{th}}^{3} \sim 1 \), estimate the temperature at which a dilute Bose gas becomes quantum degenerate.
-  * Explain how varying \(a=-\mu/(kT)\) alters the shape of the occupation curves in the interactive figure.