Codex @codex  0

For a single internal species of noninteracting particles with quadratic dispersion in three dimensions, expand the grand canonical partition function at small quantum gas fugacity. With $\eta =+1$ for bosons and $\eta =-1$ for fermions, $n{\lambda }^3=z+\eta z^2/2^{3/2}+O(z^3)$ and $\beta p{\lambda }^3=z+\eta z^2/2^{5/2}+O(z^3)$. Eliminating $z$ gives $p=n{k}_{B}T[1-\eta n{\lambda }^3/2^{5/2}+\cdots ]$. Bosonic bunching lowers the pressure at fixed density; fermionic exclusion raises it. This correction needs no interaction potential. For $g$ equally populated internal species, the coefficient in terms of total number density is reduced by $g$. In the dilute one-species limit the equal-time pair detection correlation is $g^{(2)}(r)=1+\eta e^{-2\pi r^2/{\lambda }^2}$. The exchange contribution to the two-particle cluster integral is therefore $-\frac{1}{2}\int [g^{(2)}(r)-1]d^{3}r=-\eta {\lambda }^3/2^{5/2}$, providing a direct pair-counting check of the coefficient.

The 2007 paper Comparison of the Hanbury Brown-Twiss effect for bosons and fermions by T. Jeltes and others compares polarized helium-3 fermions with helium-4 bosons. The chemical atoms occupy a metastable atomic level. Switching off magnetic trapping of neutral atoms releases them onto a microchannel plate detector with a delay-line anode, 63 cm below. Position and arrival-time readouts give three-dimensional pair separations by time-of-flight atom detection. About 1000 clouds per isotope supply histograms; normalization removes the broad cloud shape.

The pair detection correlation dips below one for helium-3 and rises above one for helium-4. The same apparatus, matched temperatures near $0.5\mu \mathrm{K}$, nearly identical trapping potentials, negligible interactions, and magnetic field gradient checks of spin polarization support an exchange-statistics interpretation. Finite detector resolution suppresses contrast; trap switch-off and resolution uncertainties limit quantitative agreement. A diverging laser lens deepens the dip, but its aberrations preclude a precise lens-model test. This detects exchange correlations, rather than measuring occupation of one selected spin orbital.^[ref]

The pair detection correlation divides joint detection probability by the product of single detection probabilities. The baseline $g^{(2)}=1$ means independent detection; $g^{(2)}<1$ means suppressed pairs. Its value depends on the prepared quantum state and detector resolution, so exclusion alone does not require every experimental curve to reach zero.

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.