Differences

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

--- un:bose-einstein-statistics [2025/11/30 07:16] – asad
+++ un:bose-einstein-statistics [2025/11/30 07:37] (current) – asad
@@ Line 1: / Line 1: @@
 ====== Bose–Einstein statistics ======
-Bose–Einstein statistics describes the equilibrium occupation of single-particle energy levels by indistinguishable bosons with integer spin. Because bosons are not constrained by the Pauli exclusion principle, many particles can occupy the same quantum state. Photons, some atomic species in ultracold gases, and various quasiparticles such as phonons all obey Bose–Einstein statistics. In high-occupancy regimes, this leads to strong enhancement of low-energy states compared with classical expectations.
+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.
-In thermal and diffusive equilibrium at temperature \( T \) and chemical potential \( \mu \), the probability of finding a given number of particles in a state is determined by the grand canonical ensemble. Maximizing the entropy under constraints of fixed average particle number and energy leads to a characteristic form for the mean occupation number of each energy level.
+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.
-===== Bose–Einstein occupation number =====
+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.
-Consider an ideal gas of non-interacting bosons with single-particle energies \( \epsilon_i \). The mean occupation number \( \bar{n}_i \) of the level \( \epsilon_i \) in thermal and chemical equilibrium is
+===== Bose–Einstein phase-space distribution =====
+The statistical mechanics of bosons gives the number of particles per phase-space cell of volume \( h^{3} \) as
 $$
-\bar{n}_i = \frac{1}{\exp\!\left(\frac{\epsilon_i - \mu}{kT}\right) - 1},
+f = \frac{2}{h^{3}} \, \frac{1}{e^{h\nu/(kT)} - 1}.
 $$
-where \( k \) is Boltzmann’s constant and \( \mu \) is the chemical potential. The total particle number and energy are obtained by summing over all single-particle states,
+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 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:
 $$
-N = \sum_i \bar{n}_i, \qquad
+f = \frac{g}{h^{3}} \, \frac{1}{e^{(\epsilon-\mu)/(kT)} - 1}.
-E = \sum_i \epsilon_i \bar{n}_i.
 $$
-In the continuum limit, one replaces the sum by an integral over energy with density of states \( g(\epsilon) \),
+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
 $$
-n = \frac{N}{V} = \int_0^\infty \frac{g(\epsilon)}{\exp\!\left(\frac{\epsilon - \mu}{kT}\right) - 1} \, d\epsilon.
+\bar{n}_i = \frac{1}{e^{(\epsilon_i - \mu)/(kT)} - 1}.
 $$
-For nonrelativistic particles of mass \( m \) in three dimensions, \( g(\epsilon) \propto V\,\epsilon^{1/2} \), so low-energy states contribute strongly to the total number and energy.
+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
+E = \sum_i \epsilon_i \bar{n}_i.
+$$
+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.
+$$
-Because the denominator can become small as \( \epsilon_i \rightarrow \mu \), the mean occupation \( \bar{n}_i \) can become very large. For an ideal Bose gas in equilibrium one always has \( \mu \leq \epsilon_0 \), where \( \epsilon_0 \) is the ground-state energy. Approaching this limit leads to qualitatively new behavior.
+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 =====
-In the dilute, high-temperature regime where
+When the exponent is large,
 $$
-\exp\!\left(\frac{\epsilon_i - \mu}{kT}\right) \gg 1,
+e^{(\epsilon-\mu)/(kT)} \gg 1,
 $$
-the “\(-1\)” in the denominator of the Bose–Einstein distribution can be neglected. The mean occupation number then reduces to the classical Maxwell–Boltzmann form,
+the Bose–Einstein distribution reduces to the Maxwell–Boltzmann form
 $$
-\bar{n}_i \approx \exp\!\left(-\frac{\epsilon_i - \mu}{kT}\right).
+\bar{n}_i \approx e^{-(\epsilon_i-\mu)/(kT)},
 $$
-In this limit, the average occupation of each quantum state satisfies \( \bar{n}_i \ll 1 \), and the discrete nature of the energy levels becomes unimportant. Quantum statistics becomes necessary when the phase-space density approaches unity. This condition can be expressed using the thermal de Broglie wavelength,
+and the phase-space distribution becomes
 $$
-\lambda_{\mathrm{th}} = \sqrt{\frac{2\pi \hbar^2}{m k T}},
+f \approx \frac{g}{h^{3}} \, e^{-(\epsilon-\mu)/(kT)}.
 $$
-so that quantum degeneracy sets in when
+Quantum statistical effects appear when the thermal de Broglie wavelength
 $$
-n \lambda_{\mathrm{th}}^3 \gtrsim 1.
+\lambda_{\mathrm{th}}
+= \sqrt{\frac{2\pi\hbar^{2}}{m kT}}
 $$
-At this point, either Bose–Einstein or Fermi–Dirac statistics must be used depending on the particle spin.
+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 =====
-As a Bose gas is cooled at fixed particle density, the chemical potential \( \mu \) rises toward the ground-state energy \( \epsilon_0 \). For a uniform ideal gas, \( \mu \) cannot exceed \( \epsilon_0 \), and at a critical temperature \( T_{\mathrm{c}} \) the excited states can no longer accommodate all particles. The excess particles accumulate in the ground state, producing a macroscopic occupation \( \bar{n}_0 \).
+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 phenomenon is known as Bose–Einstein condensation. Below \( T_{\mathrm{c}} \), a finite fraction of the total particle number resides in the lowest-energy state even in the thermodynamic limit. The condensed component behaves as a coherent matter wave and exhibits properties such as long-range phase coherence and superfluidity in interacting systems.
+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 astrophysics and cosmology, Bose–Einstein statistics is essential for describing photon fields such as the cosmic microwave background, where the chemical potential \( \mu = 0 \), and for modeling bosonic dark-matter candidates in certain theoretical scenarios.
+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 =====
-It is often convenient to introduce a dimensionless energy variable
+To analyze how the Bose–Einstein distribution behaves, it is useful to define dimensionless variables
 $$
-x = \frac{\epsilon}{kT},
+x = \frac{\epsilon}{kT}, \qquad
+a = -\frac{\mu}{kT},
 $$
-and a dimensionless chemical potential parameter
+so that
 $$
-a = -\frac{\mu}{kT} \geq 0.
+\bar{n}(x; a) = \frac{1}{e^{x+a} - 1}.
 $$
-In these variables the Bose–Einstein occupation number becomes
+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.
-$$
-\bar{n}(x; a) = \frac{1}{\exp(x + a) - 1}.
-$$
-For fixed \( a \), the curve \( \bar{n}(x; a) \) rises sharply at small \( x \), indicating strong occupation of low-energy states, and then falls off approximately as \( \exp(-x) \) at large \( x \). Smaller values of \( a \) (corresponding to larger \( \mu \)) enhance the low-energy occupation even further. In the special case of photons, such as in [[Blackbody radiation]], one has \( \mu = 0 \) and therefore \( a = 0 \).
-The interactive Plotly figure below shows \( \bar{n}(x; a) \) as a function of \( x = \epsilon/(kT) \) for different values of \( a = -\mu/(kT) \). Moving the slider changes the value of \( a \), and the curve deforms accordingly.
 <html>
@@ Line 110: / Line 134: @@
 <div style="margin-top:0.3em; text-align:center; font-family:sans-serif;">
-  Parameter \(a = -\mu/(kT)\):
+  Parameter a = -μ/(kT):
   <input type="range" id="aslider" min="0.0" max="3.0" step="0.1" value="1.0"
          style="width:400px;">
@@ Line 116: / Line 140: @@
 </div>
-<script type="text/javascript">
-window.PlotlyConfig = {MathJaxConfig: 'local'};
-</script>
 <script src="https://cdn.plot.ly/plotly-latest.min.js"></script>
@@ Line 124: / Line 145: @@
 //<![CDATA[
-// Dimensionless Bose–Einstein occupation:
+// Linear x-grid for BE occupation
-// nbar(x; a) = 1 / (exp(x + a) - 1),
-// with x = epsilon / (kT) and a = -mu / (kT) >= 0.
 const Npts = 400;
 const xMin = 0.05;
 const xMax = 8.0;
-// Create a linear grid in x
 const x_vals = Array.from({length: Npts}, (_, i) =>
   xMin + (xMax - xMin) * i / (Npts - 1)
@@ Line 139: / Line 155: @@
 // Bose–Einstein occupation function
 function nbar(x, a) {
-  const arg = x + a;
+  const e = Math.exp(x + a);
-  const expArg = Math.exp(arg);
+  return (e <= 1) ? Infinity : 1 / (e - 1);
-  // Avoid negative or zero denominator numerically
-  const denom = expArg - 1.0;
-  if (denom <= 0.0) {
-    return Number.POSITIVE_INFINITY;
-  }
-  return 1.0 / denom;
 }
+// Create a trace for a given "a"
 function makeBETrace(a) {
-  const y_vals = x_vals.map(x => nbar(x, a));
   return {
     x: x_vals,
-    y: y_vals,
+    y: x_vals.map(x => nbar(x, a)),
     mode: 'lines',
     line: {color: '#cc0000'},
@@ Line 160: / Line 170: @@
 }
-// Equation annotation for the plot
+// Annotation with plain-text equation
 const eq_be = {
   x: 0.97, y: 0.97, xref: 'paper', yref: 'paper',
-  text: '$\\bar{n}(x; a) = \\dfrac{1}{\\exp(x + a) - 1}$',
+  text: 'n̄(x; a) = 1 / (exp(x + a) - 1)',
   showarrow: false,
-  font: {size: 24},
+  font: {size: 20},
   align: 'right'
 };
+// Axis labels with no MathJax
 const layout_be = {
   margin: {l: 95, r: 20, t: 10, b: 60},
   xaxis: {
-    title: 'x = \\(\\epsilon / (kT)\\)',
+    title: 'x = ε / (kT)',
     automargin: true
   },
   yaxis: {
-    title: '\\(\\bar{n}(x; a)\\)',
+    title: 'n̄(x; a)',
     type: 'log',
-    automargin: true,
+    automargin: true
-    titlefont: {size: 18}
   },
   annotations: [eq_be]
 };
-// Initial value of a
+// Initial plot
-let a0 = 1.0;
+Plotly.newPlot('be_occ', [makeBETrace(1.0)], layout_be);
-Plotly.newPlot('be_occ', [makeBETrace(a0)], layout_be);
-// Slider behavior
+// Slider control
 const aslider = document.getElementById('aslider');
 const Aval = document.getElementById('Aval');
@@ Line 204: / Line 213: @@
 ===== Insights =====
-  * The Bose–Einstein occupation number \( \bar{n}_i = [\exp((\epsilon_i - \mu)/(kT)) - 1]^{-1} \) allows arbitrarily large occupation of a single quantum state, reflecting bosonic bunching.
+  * 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.
-  * In the dilute limit where \( \exp((\epsilon_i - \mu)/(kT)) \gg 1 \), the Bose–Einstein distribution reduces to the Maxwell–Boltzmann form \( \bar{n}_i \approx \exp[-(\epsilon_i - \mu)/(kT)] \).
+  * Large low-energy occupation arises because the denominator may become small as \( \mu \to \epsilon_{0} \).
-  * Quantum degeneracy becomes important when the phase-space density satisfies \( n \lambda_{\mathrm{th}}^3 \gtrsim 1 \) with \( \lambda_{\mathrm{th}} = \sqrt{2\pi \hbar^2/(m k T)} \), requiring Bose–Einstein or Fermi–Dirac statistics.
+  * Classical statistics emerges only when \( e^{(\epsilon-\mu)/(kT)} \gg 1 \), suppressing large occupation.
-  * Bose–Einstein condensation occurs when \( \mu \rightarrow \epsilon_0 \) as the gas is cooled, leading to a macroscopic ground-state occupation below the critical temperature \( T_{\mathrm{c}} \).
+  * Quantum degeneracy appears when \( n \lambda_{\mathrm{th}}^{3} \gtrsim 1 \), forcing the use of Bose–Einstein statistics.
-  * For photons in blackbody radiation, the chemical potential vanishes (\( \mu = 0 \)), so the distribution reduces to \( \bar{n}(x) = [\exp(x) - 1]^{-1} \) and underlies the Planck spectrum.
+  * 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 =====
-  * Starting from the grand canonical ensemble, derive the Bose–Einstein occupation number \( \bar{n}_i = [\exp((\epsilon_i - \mu)/(kT)) - 1]^{-1} \) by maximizing the entropy under constraints on \( N \) and \( E \).
+  * Derive the phase-space distribution \( f = g h^{-3} [e^{(\epsilon-\mu)/(kT)} - 1]^{-1} \) from the grand canonical ensemble.
-  * Show explicitly how the Maxwell–Boltzmann distribution \( \bar{n}_i \approx \exp[-(\epsilon_i - \mu)/(kT)] \) arises as the low-occupation limit of the Bose–Einstein distribution.
+  * Show how the occupation number \( \bar{n}_i \) emerges from integrating \( f \) over a state’s phase-space cell.
-  * Using the condition \( n \lambda_{\mathrm{th}}^3 \sim 1 \), estimate the temperature at which a given Bose gas becomes quantum degenerate for specified \( n \) and \( m \).
+  * Demonstrate explicitly how the Maxwell–Boltzmann limit arises from the Bose–Einstein distribution.
-  * Explain qualitatively how varying \( a = -\mu/(kT) \) changes the shape of \( \bar{n}(x; a) \) in the Plotly figure, especially at small \( x \).
+  * Using \( n \lambda_{\mathrm{th}}^{3} \sim 1 \), estimate the temperature at which a dilute Bose gas becomes quantum degenerate.
-  * Discuss why photons in thermal equilibrium must have \( \mu = 0 \) and how this leads to the blackbody spectrum described in the [[Blackbody radiation]] article.
+  * Explain how varying \( a = -\mu/(kT) \) alters the shape of the occupation curves in the interactive figure.