Codex @codex  0

Bose-Einstein condensation is macroscopic occupation of the lowest one-particle state. It is a change in occupation statistics, not ordinary liquid formation and not a requirement that all atoms occupy one spatial point. For ideal bosons,Set the lowest one-particle energy to zero. Then $\mu \le 0$, and the population of excited states in a uniform three-dimensional box isHere $g_s(z)={\sum }_{j\ge 1}{z}^j/j^s$. As $\mu$ approaches zero, the excited-state population saturates at $V\zeta (3/2)/{\lambda }^3$. Any excess particles occupy the ground state. In the thermodynamic limit,Thus the onset occurs when the thermal de Broglie wavelength is comparable to the mean interparticle spacing. At finite $N$, the chemical potential remains slightly below the ground-state energy; setting it equal to zero is a thermodynamic-limit statement.

The 1995 rubidium experiment used magnetic confinement and evaporative cooling of atoms. Absorption images after trap manipulation revealed a narrow, anisotropically expanding component superposed on the broader thermal distribution; the new component grew as cooling proceeded. The thermal component supplied an internal comparison through the same imaging optics. Trap decompression and release geometry mattered, and finite imaging resolution and interactions limited direct comparison with an ideal-gas ground state.^[ref]

A narrow cloud alone is not sufficient evidence. One should check its emergence across a controlled cooling sequence, atom-number and temperature calibration, the thermal background and the predicted expansion anisotropy. These are tests of an occupation transition rather than a simple focusing artifact. The box formula above is not the exact critical-temperature or condensate-fraction formula for the trapped experiment: the trap changes the density of states, and interactions also matter. Condensation and superfluid response are related but distinct properties; the former definition is ground-state occupation.

Excited-state saturation. Consider first a uniform three-dimensional noninteracting gas of identical spinless bosons in a large box. Set the one-particle ground energy to zero. The Bose-Einstein distribution is $\bar{n}(ϵ)=[z^{-1}{e}^{\beta ϵ}-1{]}^{-1}$, with quantum gas fugacity $z=e^{\beta \mu }\le 1$. For dispersion $ϵ=p^2/(2m)$, the excited density of states isUsing the thermal de Broglie wavelength ${\lambda }_T=\sqrt{2\pi {\hslash }^2/(m{k}_{B}T)}$, integration givesThe polylogarithm is ${\mathrm{Li}}_s(z)={\sum }_{\ell \ge 1}{z}^{\ell }/{\ell }^s$. Its finite value at $z=1$ implies a maximum excited population $V\zeta (3/2)/{\lambda }_T^3$, with $\zeta (3/2)\simeq 2.612$. Once this is less than $N$, the excess particles occupy the ground state macroscopically. This is Bose-Einstein condensation, rather than occupation of a state by merely one or a few atoms.

Transition and condensate fraction. Setting the saturation number equal to $N$ gives the homogeneous ideal Bose gas condensation temperatureBelow this temperature, $\mu$ approaches zero from below in the thermodynamic limit, and particle conservation givesAbove $T_c$ there is no macroscopic condensate and $\mu <0$ adjusts to keep $N$ fixed. At finite $N$, the chemical potential remains slightly below the ground energy and the transition is rounded.

Energy and the role of geometry. The excited-state internal energy isThe zero-energy condensate contributes no term to this energy in the ideal model. Below $T_c$,Condensation depends on low-energy state counting: for a uniform quadratic-dispersion gas, a finite-temperature bulk condensate requires spatial dimension greater than two. A trap changes the density of states and therefore the temperature law; a three-dimensional harmonic trap gives a cubic, rather than $3/2$, condensate-fraction exponent. Finally, an ideal-gas condensate should not automatically be called a superfluid: its free quadratic spectrum has zero critical speed by the Landau criterion. Interactions matter for the robust collective-flow properties of real dilute condensates.

Particles within each species are indistinguishable, while the two species are physically distinguishable and their numbers are fixed. Independence makes the canonical partition function factor by species:Since $Z_1\propto T^{3/2}$ at fixed volume, $U=-{\partial }_{\beta }\ln Z=3N{k}_{B}T/2$. The canonical entropy is thereforeFor large species populations, Stirling's approximation gives the equivalent Sackur-Tetrode equation formFor structureless particles with the same translational partition function, $Z_1=V/{\lambda }^3$, where $\lambda$ is their common thermal de Broglie wavelength. More generally the supplied common constant factor in $Z_1$ should be retained in the entropy.

Writing $x_i=N_i/N$ separates the total-number contribution from the composition term:The last term is the entropy of mixing relative to separated ideal components at the same temperature and pressure. This is the entropy of a two-species classical ideal gas. Artificially labeling otherwise identical chemical atoms would not justify dividing by separate factorials or adding a mixing entropy. The pure-species limit is obtained by $x\ln x\to 0$.

Classical translational statistics. Write the number density as $n=N/V$ and the thermal de Broglie wavelength as $\lambda =\sqrt{2\pi {\hslash }^2/(m{k}_{B}T)}$. Quantum exchange is negligible when wave packets rarely overlap:The equivalent state-counting statement is that the mean occupation of each thermally relevant one-particle state is much smaller than one. The equations in this question count one internal species; with $g$ equally populated spin states the relevant density per species is $n/g$.

Let $z=e^{\beta \mu }$ be the dimensionless quantum gas fugacity, with the lowest one-particle energy chosen as zero. Both quantum distributions can be writtenIf $z\ll 1$, then $n_{\eta }(E)\simeq z{e}^{-\beta E}\ll 1$ for all $E\ge 0$, recovering Maxwell–Boltzmann statistics. Conversely, low occupation even at the spectrum's bottom requires $z/(1-z)\ll 1$ for bosons or $z/(1+z)\ll 1$ for fermions; either condition impliesA single high-energy state can be sparsely occupied even in a degenerate gas, so the condition must hold across the relevant spectrum, not just for one arbitrarily chosen energy. Also, classical translational statistics does not imply that all internal excitation energies or other quantum effects can be ignored.

Grand partition function, density and pressure. Noninteracting particles have independent occupations of each one-particle state. Put $x_{\mathbf{k}}=z{e}^{-\beta {\epsilon }_{\mathbf{k}}}$. Summing the allowed occupations givesThereforeThe grand potential is ${\Omega }_G=-k_{B}T\ln \Xi =-pV$ for this homogeneous gas. Replacing the sum by $V\int d^{3}k/(2\pi {)}^3$ and using ${\epsilon }_{\mathbf{k}}={\hslash }^2{k}^2/(2m)$ givesThe mean number is $N=z{\partial }_z\ln \Xi$, givingFor bosons, the grand sum requires $z<1$ when the ground energy is zero. A macroscopic condensate at the limiting $z\to 1$ must be separated from the continuum integral in the condensed phase. The following small-$z$ calculation is in the normal, nondegenerate regime, so that complication does not enter.

Density expansion and inversion. Expand the occupation as a geometric series:The Gaussian wavevector integral for its $\ell$th term isIt follows thatFor the inverse series, put $y=n{\lambda }^3$ and try $z=y+c{y}^2+O(y^3)$. Substitution gives $0=(c+\eta /2^{3/2})y^2+O(y^3)$, henceThe leading equality $n{\lambda }^3\simeq z$ directly connects the density and occupancy criteria for classical behavior.

Pressure and its first quantum correction. Expanding the logarithm in the pressure expression gives one extra factor $1/\ell$ compared with the density series:Insert the inverse series. To second order in $y=n{\lambda }^3$,ThusThe first quantum correction at fixed density and temperature isIt is negative for bosons and positive for fermions, and vanishes in the classical density limit.

Virial coefficient. Define the number-density virial expansion byThe coefficient $B_j$ has units ${\mathrm{volume}}^{j-1}$ in this convention. Comparing terms givesThis is the exchange contribution to the second virial coefficient. It arises without a pair-interaction potential; the underlying gas is noninteracting. A molar-density convention would multiply this $B_2$ by the Avogadro constant. For $g$ equally populated internal species, exchange occurs within each species and the coefficient in terms of total density is smaller by a factor $g$.

Why the two signs differ. Bosonic exchange symmetry enhances the probability of nearby identical bosons, whereas fermionic exchange symmetry suppresses nearby identical fermions in the same internal state. In the dilute regime, the normalized equal-time pair detection correlation isThe exchange part follows by squaring the Gaussian one-particle coherence factor $e^{-\pi r^2/{\lambda }^2}$. At zero separation this gives $g^{(2)}(0)=2$ for bosons and $0$ for one-component fermions. Integrating the exchange pair weight provides a direct coefficient check:Bunching therefore acts like an effective statistical attraction in the dilute equation of state, lowering the pressure needed at fixed density. Exclusion acts like an effective statistical repulsion, raising it. Neither description asserts an actual interparticle force. Different internal fermion species do not exclude each other from a common spatial location in the same way, explaining the importance of the species convention.

Put $\beta =1/(k_{B}T)$ and define the thermal de Broglie wavelength of species $\alpha$ byIgnoring the internal degrees of freedom as stipulated, each molecule's translational canonical partition function is $V/{\lambda }_{\alpha }^3$. The two species do not interact, so their partition functions multiply; the treatment of indistinguishable particles requires a separate factorial for each species, not a factorial of the total particle number. ThusUsing $F=-k_{B}T\ln Z$ and Stirling's approximation gives the extensive Helmholtz free energyThese are ideal gas expressions, valid for $n_{\alpha }{\lambda }_{\alpha }^3\ll 1$, with $n_{\alpha }=N_{\alpha }/V$. The energy zero has been chosen so that separated CO and oxygen molecules have zero internal reference energy.

The given binding energy $ϵ$ belongs to the formation of two carbon-dioxide molecules. Relative to the reference used in (a), the internal energy offset of one carbon-dioxide molecule is therefore $-ϵ/2$, not $-ϵ$. Its mass is $m_c=m_{\mathrm{CO}}+m_{{\mathrm{O}}_2}/2$ and its thermal de Broglie wavelength is ${\lambda }_c=h/\sqrt{2\pi m_c{k}_{B}T}$, where the subscript $c$ denotes $\mathrm{C}{\mathrm{O}}_2$.

The single-molecule canonical partition function is $(V/{\lambda }_c^3)e^{\beta ϵ/2}$. For $N_c$ indistinguishable noninteracting molecules,The corresponding Helmholtz free energy, using Stirling's approximation, isThe negative sign in the last term expresses the stabilizing binding energy and ensures that low temperature favors association.

At fixed temperature, volume and the numbers of all other species, the chemical potential is ${\mu }_{\alpha }=(\partial F/\partial N_{\alpha }{)}_{T,V,N_{{\alpha }^{\prime }\ne \alpha }}$. For the combined noninteracting mixture, add the free energies from (a) and (b). Since the thermal de Broglie wavelengths are independent of particle number at fixed temperature, differentiation givesEach logarithm has a dimensionless argument, and each chemical potential has units of energy per molecule. The binding offset must be kept consistently when comparing species.

Substituting the chemical potentials into the equilibrium relation givesExponentiation and the positive square root give the law of mass action in number-density form:All three thermal de Broglie wavelengths have the same $T^{-1/2}$ dependence. Define the temperature-independent mass factorIf $f_{\alpha }=n_{\alpha }/n$, this exact equilibrium relation can be writtenIt already has the requested dimensional structure: $n^{3/2}{\lambda }_c^{3/2}$ has units ${\mathrm{m}}^{-3}$. Its coefficient, however, depends on the conserved elemental composition through the species fractions.

For the intended case of a gas obtained by dissociating initially pure carbon dioxide, conservation of the elemental ratio gives $n_{\mathrm{CO}}=2n_{{\mathrm{O}}_2}$ throughout. Henceand the exact relation within this ideal-mixture model isAt high temperature with $k_{B}T\gg ϵ$ and in the dilute ideal gas regime, $n{\lambda }_c^3\ll 1$, this equation is self-consistently dissociation dominated:Thus $n-n_c\simeq n$, and to leading orderUsing nominal molecular masses $28,32,44$ in common mass units gives $A\simeq 0.963$, confirming that it is of order unity. The exponential has been retained as the leading binding correction, although it tends to one as $T\to \infty$.

The printed question does not explicitly fix the elemental composition. Without the initially-pure-$\mathrm{C}{\mathrm{O}}_2$ assumption, the universal result is the boxed mass-action relation, and $A=B{f}_{\mathrm{CO}}\sqrt{f_{{\mathrm{O}}_2}}$ approaches a composition-dependent constant when both dissociation products have nonzero high-temperature fractions. It need not be of order unity for an arbitrarily unbalanced mixture. Also, $k_{B}T\gg ϵ$ alone does not justify replacing $n-n_c$ by $n$ at arbitrarily high density; the dilute classical regime is part of the model.

The dimensionless occupancy scale $n{\lambda }_T^3$, with number density $n$ and thermal de Broglie wavelength ${\lambda }_T=h/\sqrt{2\pi m{k}_{B}T}$. Increasing it approaches the regime where quantum statistics matters. This use of degeneracy must not be confused with degeneracy of an individual energy eigenspace.

The translational concentration scale for a nonrelativistic particle in three dimensions, per internal state. It is the inverse cube of the thermal de Broglie wavelength. A classical component with number density $n$ and external one-particle energy $U$ has chemical potential $\mu =k_{B}T\ln (n/n_Q)+U$. An additional internal degeneracy must be included explicitly; it is not hidden in this definition. The dilute classical condition is $n/n_Q\ll 1$ for each occupied internal component.

An ideal gas whose particles move in an area $A$ with two independent translational momentum components. In its dilute classical limit the one-particle translational partition function is $q=A/{\Lambda }^2$, where $\Lambda$ is the thermal de Broglie wavelength, and $Z_N=q^N/N!$. The translational internal energy is $N{k}_{\mathrm{B}}T$, and the surface pressure conjugate to area is $N{k}_{\mathrm{B}}T/A$. A mobile dilute adsorbed gas is an example; a fully occupied array of fixed adsorption sites has different state counting.

Equating the site chemical potential with that of an internally structureless classical vapour gives $\theta =b_{s}p/(1+b_{s}p)$. The coefficient $b_s$ includes the bound vibrational partition factor, the translational thermal de Broglie wavelength, and the positive binding Boltzmann factor. It has inverse-pressure units. Half coverage occurs at $p=1/b_s$; an unspecified $z_s(T)$ prevents a numerical pressure scale, but not the functional isotherm.

The virial expansion describes the first departures from the ideal gas law as a power series in number density. For a homogeneous, single-species gas in three dimensions at thermal equilibrium, write $n=N/V$, with $V$ the volume, and define the dimensionless compressibility factor $Z$ byHere $p$ is absolute pressure, $T$ is absolute thermodynamic temperature, and $k_B$ is the Boltzmann constant. The second virial coefficient determines the initial slope away from $Z=1$; higher coefficients describe further corrections. The symbol $Z$ here is an equation of state ratio, not a canonical partition function. ^[ref]

The left panel sketches the measurement principle; the right panel uses deliberately chosen polynomials, $Z=1+bx+6x^2$ with $x=n{\ell }^3$ and $b=\pm 2$, where $\ell$ is a reference length. These curves are not measured data. Their initial slopes have opposite signs, but both bend because of the quadratic term. Using a finite secant from $x=0$ to $x=h$ gives $(Z(h)-1)/h=b+6h$, rather than $b$: reducing the fitting range reduces this particular truncation bias. This algebra explains why measuring one nonideal pressure is insufficient to extract a reliable second virial coefficient.

Waxman and Hastings's 1971 Burnett experiment used sample and expansion chambers in a thermostatic bath, with a valve between them and a pressure transducer linked to a piston gauge. They measured successive isothermal expansions: evacuate the expansion chamber, open the valve, record the equilibrated pressure, isolate and evacuate that chamber again, then repeat. This avoids separately measuring the absolute gas amount and chamber volumes. For carbon dioxide at $50{}^{\circ }\mathrm{C}$, their measurements below $35$ bar gave a molar second virial coefficient of approximately $-102.2\mathrm{c}{\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$. ^[ref]

An idealized reduction makes the signal explicit. Let $q=(V_A+V_B)/V_A$ be the calibrated chamber volume ratio, $p_i$ the pre-expansion pressure, and $p_f$ the post-expansion pressure. Particle conservation during one valve opening gives $n_f=n_i/q$. At the same temperature, division of the two equations of state yieldsThus the intercept determines $q$, while the dilute-limit slope determines the second virial coefficient. In the actual experiment, helium runs calibrated the cell constant and chamber deformation corrections were included. Agreement of the carbon dioxide low-pressure intercept with the helium calibration, and the absence of an anomalous low-pressure upturn, checked for gas adsorption on chamber walls. ^[ref]

For a limited dilute-gas fit, hold temperature fixed, allow thermal equilibrium after every expansion, calibrate absolute pressure and the chamber ratio, and check leaks and sample composition. Repeat the analysis with a smaller upper pressure cutoff and an added curvature term; compare the resulting second virial coefficients with their uncertainties. Too wide a range mixes in higher coefficients; too narrow a range makes the nonideal signal hard to distinguish from instrument error. The displayed chamber equation neglects deformation, dead spaces and adsorption, so it is a teaching reduction rather than the complete experimental analysis.

Units matter. With number density in ${\mathrm{m}}^{-3}$, $B_2$ has units ${\mathrm{m}}^3$ and $B_3$ has units ${\mathrm{m}}^6$. For molar density $c=n/N_{\mathrm{A}}$ in $\mathrm{mol}{\mathrm{m}}^{-3}$, the same equation of state isHere $N_{\mathrm{A}}$ is the Avogadro constant and $R=N_{\mathrm{A}}{k}_B$ is the molar gas constant; the molar coefficients have units ${\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$ and ${\mathrm{m}}^6\mathrm{mo}{\mathrm{l}}^{-2}$. A coefficient tabulated per mole must therefore multiply molar density, not number density.

For the classical formulas below, use a nondegenerate gas with thermal de Broglie wavelength ${\lambda }_{\mathrm{th}}$ satisfying $n{\lambda }_{\mathrm{th}}^3\ll 1$, negligible quantum corrections to relative motion, and a stable, sufficiently short-ranged pair potential. Diluteness must also control the omitted interaction terms: $|B_{2}n|\ll 1$ alone is not enough when positive and negative contributions cancel. A finite truncation is a local low-number density approximation, not an equation of state for a dense liquid or across a phase transition. ^{[ref] [ref]}