Differences

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

--- un:bose-einstein-statistics [2025/11/30 07:13] – created asad
+++ un:bose-einstein-statistics [2025/11/30 07:37] (current) – asad
@@ Line 1: / Line 1: @@
-====== Bose–Einstein Statistics ======
+====== Bose–Einstein statistics ======
-Bose–Einstein (BE) statistics describes the equilibrium occupation of single–particle energy levels by indistinguishable bosons (integer spin), which are not constrained by the Pauli exclusion principle. A given energy level with energy <math>\epsilon_i</math> can therefore host an arbitrarily large number of particles.
+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.
-For an ideal Bose gas in thermal and diffusive equilibrium at temperature <math>T</math> and chemical potential <math>\mu</math>, the mean occupation number of the level <math>\epsilon_i</math> is
+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.
-<math>
-\bar{n}_i \;=\; \frac{1}{\exp\!\left(\dfrac{\epsilon_i-\mu}{k_{\mathrm B}T}\right) - 1},
-</math>
-where <math>k_{\mathrm B}</math> is Boltzmann’s constant. The total particle number and total energy are
-<math>
-N \;=\; \sum_i \bar{n}_i,
-\qquad
-E \;=\; \sum_i \epsilon_i \bar{n}_i.
-</math>
-In the continuum limit with density of states <math>g(\epsilon)</math>, the number density is
+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.
-<math>
-n \;=\; \frac{N}{V}
-\;=\;
-\int_0^\infty
-\frac{g(\epsilon)}{\exp\!\left(\dfrac{\epsilon-\mu}{k_{\mathrm B}T}\right) - 1}\,d\epsilon.
-</math>
-For non–relativistic particles in three dimensions, <math>g(\epsilon) \propto V\,\epsilon^{1/2}</math>, so low–energy states have a strong influence on the thermodynamics.
-Because the denominator can become small as <math>\epsilon_i \to \mu</math>, the mean occupation <math>\bar{n}_i</math> can become very large. In an ideal gas, the chemical potential satisfies <math>\mu \le \epsilon_0</math>, the ground–state energy. As <math>T</math> decreases at fixed density, <math>\mu</math> approaches <math>\epsilon_0</math> from below and a macroscopic population of the ground state appears: Bose–Einstein condensation.
+===== Bose–Einstein phase-space distribution =====
-In the dilute, high–temperature limit,
+The statistical mechanics of bosons gives the number of particles per phase-space cell of volume \( h^{3} \) as
-<math>
-\exp\!\left(\dfrac{\epsilon_i-\mu}{k_{\mathrm B}T}\right) \gg 1
-</math>
-and the “–1” in the denominator can be neglected, giving the classical Maxwell–Boltzmann distribution
-<math>
-\bar{n}_i \;\approx\; \exp\!\left(-\dfrac{\epsilon_i-\mu}{k_{\mathrm B}T}\right).
-</math>
-Quantum degeneracy is important when the phase–space density
-<math>
-n \lambda_{\mathrm{th}}^3 \gtrsim 1,
-\qquad
-\lambda_{\mathrm{th}} = \sqrt{\frac{2\pi \hbar^2}{m k_{\mathrm B}T}}
-</math>
-(thermal de Broglie wavelength) is of order unity, signalling that BE statistics is required.
-In astrophysics and cosmology, BE statistics underlies the distribution of photons (with <math>\mu = 0</math>) and is important for bosonic quasiparticles (e.g. phonons) and certain dark–matter candidates in high–occupation–number regimes.
+$$
+f = \frac{2}{h^{3}} \, \frac{1}{e^{h\nu/(kT)} - 1}.
+$$
-===== Plotly visualization of the Bose–Einstein occupation function =====
+Here the factor of 2 accounts for two photon polarization states. For massive bosons this factor is replaced by the appropriate spin degeneracy \( g \). This function expresses a key difference between classical and quantum statistics: the denominator can become small, allowing large occupation of low-energy states.
-For visualization it is convenient to introduce the dimensionless variables
+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:
-<math>
-x \;=\; \frac{\epsilon}{k_{\mathrm B}T},
+$$
+f = \frac{g}{h^{3}} \, \frac{1}{e^{(\epsilon-\mu)/(kT)} - 1}.
+$$
+The form of \( f \) therefore unifies the physics of blackbody photons and ideal Bose gases.
+===== Mean occupation of quantum states =====
+The phase-space distribution function determines the mean number of bosons occupying each single-particle quantum state. Integrating \( f \) over a phase-space cell associated with state \( i \) leads to the Bose–Einstein occupation number
+$$
+\bar{n}_i = \frac{1}{e^{(\epsilon_i - \mu)/(kT)} - 1}.
+$$
+This quantity plays the same role for massive bosons that the photon distribution function plays for radiation. Summing over all states gives
+$$
+N = \sum_i \bar{n}_i,
 \qquad
-a \;=\; -\frac{\mu}{k_{\mathrm B}T} \;\ge\; 0,
+E = \sum_i \epsilon_i \bar{n}_i.
-</math>
+$$
+In the continuum limit with density of states \( g(\epsilon) \),
+$$
+n = \frac{N}{V}
+  = \int_0^\infty \frac{g(\epsilon)}{e^{(\epsilon-\mu)/(kT)} - 1} \, d\epsilon.
+$$
+The crucial feature is the “\(-1\)” in the denominator; its absence in the classical Boltzmann limit suppresses large occupation at low energy.
+===== Classical limit and quantum degeneracy =====
+When the exponent is large,
+$$
+e^{(\epsilon-\mu)/(kT)} \gg 1,
+$$
+the Bose–Einstein distribution reduces to the Maxwell–Boltzmann form
+$$
+\bar{n}_i \approx e^{-(\epsilon_i-\mu)/(kT)},
+$$
+and the phase-space distribution becomes
+$$
+f \approx \frac{g}{h^{3}} \, e^{-(\epsilon-\mu)/(kT)}.
+$$
+Quantum statistical effects appear when the thermal de Broglie wavelength
+$$
+\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
+$$
+n \lambda_{\mathrm{th}}^{3} \gtrsim 1.
+$$
+Below this threshold, the large-occupation low-energy behavior of bosons becomes unavoidable.
+===== Bose–Einstein condensation =====
+For a uniform ideal Bose gas, the chemical potential satisfies \( \mu \le \epsilon_{0} \). As the gas is cooled at fixed density, \( \mu \to \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,
+$$
+\bar{n}_{0} \gg 1.
+$$
+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.
+In contrast, for photons the chemical potential is fixed to zero, and instead of condensing, photons redistribute themselves to form the blackbody spectrum.
+===== Bose–Einstein occupation curves =====
+To analyze how the Bose–Einstein distribution behaves, it is useful to define dimensionless variables
+$$
+x = \frac{\epsilon}{kT}, \qquad
+a = -\frac{\mu}{kT},
+$$
 so that
-<math>
-\bar{n}(x) \;=\; \frac{1}{\exp(x + a) - 1}.
-</math>
-The following embeddable Plotly block (HTML + JavaScript) can be pasted directly into DokuWiki (with HTML embedding enabled) to plot <math>\bar{n}(x)</math> for several values of <math>a</math>:
+$$
+\bar{n}(x; a) = \frac{1}{e^{x+a} - 1}.
+$$
+Small \( a \) (that is, large \( \mu \)) enhances the low-energy peak of the distribution. The special case \( a = 0 \) corresponds exactly to the photon distribution function in blackbody radiation.
 <html>
-<script src="https://cdn.plot.ly/plotly-latest.min.js"></script>
-<div id="be_plot" style="max-width: 800px; height: 500px;"></div>
+<style>
+#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>
+<div style="margin-top:0.3em; text-align:center; font-family:sans-serif;">
+  Parameter a = -μ/(kT):
+  <input type="range" id="aslider" min="0.0" max="3.0" step="0.1" value="1.0"
+         style="width:400px;">
+  <span id="Aval" style="margin-left:6px; font-weight:bold;">1.0</span>
+</div>
+<script src="https://cdn.plot.ly/plotly-latest.min.js"></script>
 <script type="text/javascript">
-// Dimensionless Bose–Einstein occupation:
+//<![CDATA[
-// nbar(x; a) = 1 / (exp(x + a) - 1)
-// with x = epsilon / (k_B T) and a = -mu / (k_B T).
+// Linear x-grid for BE occupation
+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 boseEinstein(x, a) {
+// Bose–Einstein occupation function
-    return 1.0 / (Math.exp(x + a) - 1.0);
+function nbar(x, a) {
+  const e = Math.exp(x + a);
+  return (e <= 1) ? Infinity : 1 / (e - 1);
 }
-// Build x-array from 0.05 to 8.0 in uniform steps
+// Create a trace for a given "a"
-var xVals = [];
+function makeBETrace(a) {
-var xMin = 0.05;
+  return {
-var xMax = 8.0;
+    x: x_vals,
-var nPoints = 400;
+    y: x_vals.map(x => nbar(x, a)),
-var dx = (xMax - xMin) / (nPoints - 1);
+    mode: 'lines',
+    line: {color: '#cc0000'},
-for (var i = 0; i < nPoints; i++) {
+    name: 'a = ' + a.toFixed(1)
-    xVals.push(xMin + i * dx);
+  };
 }
-// Parameter sets: different a = -mu / (k_B T)
+// Annotation with plain-text equation
-var params = [
+const eq_be = {
-    {a: 0.0, label: "a = 0 (\u03bc = 0, photons)"},
+  x: 0.97, y: 0.97, xref: 'paper', yref: 'paper',
-    {a: 1.0, label: "a = 1 (\u03bc = -k_B T)"},
+  text: 'n̄(x; a) = 1 / (exp(x + a) - 1)',
-    {a: 2.0, label: "a = 2 (\u03bc = -2 k_B T)"}
+  showarrow: false,
-];
+  font: {size: 20},
+  align: 'right'
+};
-var data = [];
+// Axis labels with no MathJax
+const layout_be = {
+  margin: {l: 95, r: 20, t: 10, b: 60},
+  xaxis: {
+    title: 'x = ε / (kT)',
+    automargin: true
+  },
+  yaxis: {
+    title: 'n̄(x; a)',
+    type: 'log',
+    automargin: true
+  },
+  annotations: [eq_be]
+};
-for (var p = 0; p < params.length; p++) {
+// Initial plot
-    var aVal = params[p].a;
+Plotly.newPlot('be_occ', [makeBETrace(1.0)], layout_be);
-    var yVals = [];
-    for (var j = 0; j < xVals.length; j++) {
+// Slider control
-        var x = xVals[j];
+const aslider = document.getElementById('aslider');
-        var y = boseEinstein(x, aVal);
+const Aval = document.getElementById('Aval');
-        yVals.push(y);
-    }
-    data.push({
+aslider.oninput = function() {
-        x: xVals,
+  const a = parseFloat(this.value);
-        y: yVals,
+  Aval.textContent = a.toFixed(1);
-        mode: "lines",
+  Plotly.react('be_occ', [makeBETrace(a)], layout_be);
-        name: params[p].label
-    });
-}
-var layout = {
-    title: "Bose–Einstein Occupation Number  n\u0304(x) = 1 / [exp(x + a) - 1]",
-    xaxis: {
-        title: "x = \u03b5 / (k_B T)",
-        rangemode: "tozero"
-    },
-    yaxis: {
-        title: "n\u0304(x)",
-        type: "log"
-    },
-    margin: {l: 70, r: 20, t: 50, b: 60},
-    legend: {x: 0.02, y: 0.98}
 };
-Plotly.newPlot("be_plot", data, layout);
+//]]>
 </script>
-</html>
-This reproduces the characteristic strong enhancement of low–energy occupation for bosons, especially when <math>a</math> is small (i.e. when <math>\mu</math> is close to the ground–state energy).
+</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.
-* The Bose–Einstein occupation number <math>\bar{n}_i = [\exp((\epsilon_i-\mu)/k_{\mathrm B}T)-1]^{-1}</math> allows arbitrarily large <math>\bar{n}_i</math>, reflecting bosonic bunching at low energy.
+  * Large low-energy occupation arises because the denominator may become small as \( \mu \to \epsilon_{0} \).
-* In the limit <math>\exp((\epsilon_i-\mu)/k_{\mathrm B}T) \gg 1</math>, BE statistics reduces to the Maxwell–Boltzmann form <math>\bar{n}_i \approx \exp[-(\epsilon_i-\mu)/k_{\mathrm B}T]</math>.
+  * Classical statistics emerges only when \( e^{(\epsilon-\mu)/(kT)} \gg 1 \), suppressing large occupation.
-* Quantum degeneracy appears when <math>n \lambda_{\mathrm{th}}^3 \gtrsim 1</math>, with <math>\lambda_{\mathrm{th}} = \sqrt{2\pi\hbar^2/(m k_{\mathrm B}T)}</math>, signalling the need for BE (or FD) statistics.
+  * Quantum degeneracy appears when \( n \lambda_{\mathrm{th}}^{3} \gtrsim 1 \), forcing the use of Bose–Einstein statistics.
-* As <math>T</math> is lowered at fixed <math>n</math>, the approach of <math>\mu</math> to <math>\epsilon_0</math> leads to Bose–Einstein condensation with macroscopic ground–state occupation.
+  * Bose–Einstein condensation results directly from the inability of excited states to hold all particles as \( \mu \to \epsilon_0 \).
-* For photons with <math>\mu = 0</math>, BE statistics combined with the electromagnetic density of states yields the Planck blackbody spectrum and the shape of the cosmic microwave background.
+  * 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.
-* Derive <math>\bar{n}_i = [\exp((\epsilon_i-\mu)/k_{\mathrm B}T)-1]^{-1}</math> from the grand–canonical ensemble by maximizing entropy with constraints on <math>N</math> and <math>E</math>.
+  * Show how the occupation number \( \bar{n}_i \) emerges from integrating \( f \) over a state’s phase-space cell.
-* Show analytically how the BE distribution reduces to <math>\bar{n}_i \approx \exp[-(\epsilon_i-\mu)/k_{\mathrm B}T]</math> when <math>\bar{n}_i \ll 1</math>, and interpret this limit physically.
+  * Demonstrate explicitly how the Maxwell–Boltzmann limit arises from the Bose–Einstein distribution.
-* For a 3D ideal Bose gas with <math>g(\epsilon) \propto \epsilon^{1/2}</math>, derive the critical temperature <math>T_{\mathrm c}(n,m)</math> of Bose–Einstein condensation.
+  * Using \( n \lambda_{\mathrm{th}}^{3} \sim 1 \), estimate the temperature at which a dilute Bose gas becomes quantum degenerate.
-* Using <math>\bar{n}(x) = [\exp(x + a)-1]^{-1}</math>, explain how increasing <math>a = -\mu/(k_{\mathrm B}T)</math> modifies the low–energy part of the curve in the Plotly plot.
+  * Explain how varying \( a = -\mu/(kT) \) alters the shape of the occupation curves in the interactive figure.
-* Discuss why a photon gas in thermal equilibrium must have <math>\mu = 0</math>, and how this condition, together with BE statistics, leads to the observed Planck spectrum.