Codex @codex  0

Differentiate the Helmholtz free energy at fixed couplings. With $\lambda ={\lambda }_{+}$ from solution, the entropy per spin isIt is particularly useful to writeAt $\bar{J}=J$, this is $3e^{\beta J}+e^{-3\beta J}$, so $f=-J/2-(\log 3)/(2\beta )+o({\beta }^{-1})$. For fixed $\bar{J}>J$, the first exponential dominates alone, giving $f=(J-2\bar{J})/2+o({\beta }^{-1})$. Hence the residual entropy isThe finite positive difference $\bar{J}-J$ must be kept fixed in the second limit; taking it to zero first gives the other answer.

At high temperature, $\lambda \to 4$ and $\beta {\partial }_{\beta }\log \lambda \to 0$, whenceAll $2^{2L}$ configurations become equally probable because their Boltzmann factors tend to one. The two available states per spin, rather than the coupling values, determine this entropy.

The Holstein–Primakoff transformation maps the finite spin basis $|S,m⟩$ to bosonic Fock states $|n⟩$ with $n=S-m$. Thus $n=0$ is the fully polarized reference state, and creating a boson lowers $S^z$ by one. The physical subspace has the Holstein–Primakoff occupation constraint $0\le n\le 2S$, not the whole infinite oscillator Hilbert space.

The specified ordering of the square root givesThese are the exact spin ladder operator matrix elements; in particular, the lowering amplitude vanishes at $n=2S$. The underlying oscillator obeys its usual canonical commutation relation on the enlarged space, while the spin operators preserve the constrained physical subspace. Projecting the oscillator operators themselves to a finite space would change their commutator at its upper endpoint.

Expanding $\sqrt{2S-n}=\sqrt{2S}(1-n/(4S)+\cdots )$ converts an ordered spin problem to approximately free magnons, with interactions in higher orders of $1/S$. Two complications are the nonpolynomial, order-sensitive square roots, which generate interacting bosons, and loss of the occupation bound when a truncated expansion uses the full bosonic Fock space. Large depletion, high temperature, or soft modes near a phase transition can invalidate the assumption $⟨n⟩\ll S$ even when $S$ is large.

Use thermal de Broglie wavelength ${\lambda }_T=h/\sqrt{2\pi m{k}_{B}T}$ and spacing $a=n^{-1/3}$. Degeneracy onset from thermal wavelength occurs atThe numerical factor depends on the convention for a typical momentum, but the scale ${\hslash }^2{n}^{2/3}/(m{k}_B)$ does not. For fermions the Fermi temperature $T_F={\hslash }^2(3{\pi }^{2}n{)}^{2/3}/(2m{k}_B)$ gives a comparable scale for spin degeneracy two.

Representative estimates, with the assumed densities stated explicitly, are:

Here $u$ is the atomic mass constant. Thus the characteristic orders are $10^5\mathrm{K}$ for electrons, sub-microkelvin for a dilute heavy atomic gas, and $10^{12}\mathrm{K}$ for neutrons. There is no density-independent quantum-degeneracy temperature for atomic gases: for example, taking $n=2.5\times 10^{25}{\mathrm{m}}^{-3}$ and mass $4u$ instead gives about $0.065\mathrm{K}$. These are illustrative input estimates, not values supplied by the paper. Denser neutron-star cores require relativistic and interaction corrections; the estimate identifies the degeneracy scale rather than a precise stellar temperature.

Let $\beta =1/(k_{B}T)$. Each independent fermion orbital is either empty or occupied once, so its grand partition factor is $1+e^{-\beta ({ϵ}_j-\mu )}$. Multiplying factors and taking the logarithm givesSince $1-n_F=[1+e^{-\beta (ϵ-\mu )}{]}^{-1}$, this is the requested form. Spin or other orbital degeneracy is already included in $g$.

For the Fermi-window convolution of the grand potential, defineDirect differentiation gives $F^{\prime }=+n_F$. The supplied PDF hint has the opposite sign; using that sign would contradict its grand-potential formula. Integrating twice with $A(0)={\Phi }_0(0)=0$ and vanishing upper boundary terms,At zero temperature $-{\partial }_{ϵ}{n}_F\to \delta (ϵ-\mu )$, confirming ${\Phi }_0(\mu )=\Phi (0,\mu )$.

Use the source's implicit unit normalization of the constant density of states; restoring a constant baseline $g_0$ multiplies the following potentials by $g_0$. For the sinusoidal density of states, direct integration givesPhysical nonnegative density of states requires $|\alpha |\le 1$. Write $t=k_{B}T$ and $y=ϵ-\mu$. The Fermi windowis even, with full-line moments $\int wdy=1$, $\int ywdy=0$ and $\int y^{2}wdy={\pi }^2{t}^2/3$. The last integral follows, for example, by expanding $e^x/(1+e^x{)}^2$ on $x>0$ into its alternating exponential series. When $t\ll \mu$, extending the lower integration limit to $-\infty$ introduces only exponentially small band-edge corrections. The supplied Fourier transform givesThe odd sine term vanishes. Substitution provesThe thermal damping of a sinusoidal density of states is therefore $X/\sinh X$ in the full oscillatory potential. The requested ratio instead concerns the thermal change $\Phi (T,\mu )-{\Phi }_0(\mu )$. Its oscillatory thermal-change amplitude ratio isFor positive $X$ it decreases smoothly from the limit $|\alpha |$ to zero:At exactly $T=0$ both thermal changes vanish; the quoted value is their limiting ratio. The amplitude is the envelope in chemical potential. At a particular fixed $\mu$, the signed contribution ratio has $\alpha \cos (2\pi \mu /\Delta )$ in place of $|\alpha |$; at an oscillation node the actual oscillatory contribution is zero. The large-$X$ behaviour still requires $k_{B}T\ll \mu$, so the useful range is $1\ll X\ll 2{\pi }^2\mu /\Delta$.

For two localized energy levels $0,\delta$ with degeneracies $d_0,d_1$, let $x=\delta /(k_{B}T)$. Canonical equilibrium givesHere $p_1$ is the total upper-level occupation probability and $u,s,c$ are per unit. For nondegenerate levels the entropy rises from zero to $k_B\ln 2$, while the heat capacity vanishes at both low and high temperatures and has a Schottky maximum near $k_{B}T=0.417\delta$. A spin one-half in a field, or a localized defect with two configurations, supplies an example. Higher levels must be thermally inaccessible; unequal degeneracies change the populations and peak. Distinguishable localized units have total partition function $z^N$, unlike a gas of indistinguishable particles sharing just two orbitals.

For the magnetic atom, choose level labels with $E_m=-mg{\mu }_{B}B$, $m=-J,\dots ,J$. Put $x=T_B/T=g{\mu }_{B}B/(k_{B}T)$. Its partition function and exact entropy of a finite magnetic multiplet areThe ratio form has removable limit $Z_J(0)=2J+1$. Since $s_J$ is a function of $x$ only, it depends on field and temperature solely through $B/T$. For fixed $J$ at small $x$, symmetry removes the odd powers and the supplied sum of $m^2$ givesFor large $J$, the expansion requires $J|x|\ll 1$, not merely $|x|\ll 1$.

For demagnetization with a linear-heat-capacity load, integrate its specified $C_L/T$:Isothermal magnetization at $T_i$ removes spin entropy; the bath receives the corresponding heat. During the following reversible, thermally isolated demagnetization, it is the combined entropy of spins and load that stays constant. At the minimum-field endpoint $B_f=0$, the ideal independent spins reach $s_{\max }=k_B\ln (2J+1)$. ThereforeThe bracket is the initial magnetic entropy deficit. The PDF's intermediate equality uses the initial entropy itself instead; with its preceding definition of $s$, that equality omits the maximum-entropy term. The corrected balance recovers its final temperature formula:Without the small-field expansion, the exact independent-spin prediction is $T_i-(N\theta /M)[\ln (2J+1)-s_J(T_i,B_i)/k_B]$. As a check, zero initial field gives zero entropy deficit and no cooling. The incorrectly printed intermediate balance would fail this check.

The result requires the load law throughout the temperature interval, hence $T_i\ll \theta$; a small initial Zeeman splitting, quantitatively $J{T}_{B,i}/T_i\ll 1$; positive final temperature; negligible additional heat capacities and heat leaks; thermal equilibration between subsystems; and a reversible sweep slow enough to maintain equilibrium. A negative extrapolated temperature is a breakdown of this approximation, not a cooling prediction. Spins alone would keep $B/T$ fixed, but the attached load changes that adiabat by transferring entropy to them as it cools.

Interactions invalidate the product of independent one-site partition functions. The entropy must be computed from the many-body energy spectrum or a justified interacting approximation. Correlations and possible magnetic ordering introduce an interaction scale and can remove the independent-spin entropy available for cooling near zero field. Internal fields and ordering commonly limit practical demagnetization. The detailed correction and any residual entropy depend on the interaction Hamiltonian, dimensionality and frustration; neither a universal ordering temperature nor a numerical final-temperature correction is specified here.

Let $\Delta =\hslash \Omega$ and choose the nondegenerate energy levels as 0 and $\Delta$. For N independent identical units, $Z=(1+e^{-\beta \Delta }{)}^N$; set N=1 for a single two-level system. Differentiating its logarithm givesAt low temperature, $U\sim N\Delta e^{-x}$ and $C_V\sim N{k}_B{x}^2{e}^{-x}$: the upper level is frozen out. At high temperature, the two levels approach equal populations, giving $U\to N\Delta /2$ while $C_V\sim N{k}_B{x}^2/4\to 0$. The energy therefore rises monotonically to a plateau, while the heat capacity has a Schottky peak:Unequal degeneracies would change the populations and peak; the usual equal-degeneracy model is stated here. Shifting both energy levels changes the energy baseline but not heat capacity.

Use a three-dimensional ideal gas at the same number density n for both isotopes. An unpolarized helium-3 gas has spin degeneracy two, henceHelium-4 has a single spin state, and the ideal Bose-gas saturation condition givesEliminating n, using $m_3/m_4\simeq 3/4$ and $\zeta (3/2)=2.612$, yieldsThe fermionic spin degeneracy must be included. This is an ideal-gas estimate; strong interactions in liquid helium prevent identifying it with an exact experimental transition temperature.

Equipartition gives mean energy $k_{B}T/2$ for each independent quadratic term in a classical canonical Hamiltonian. For example, $p^2/(2m)$ and $k{x}^2/2$ each contribute this amount in a harmonic oscillator. It is a thermal average, not an instantaneous equality; constraints, nonquadratic terms or quantum level spacings large compared with $k_{B}T$ require a different calculation.

The fluctuation-dissipation theorem connects equilibrium fluctuations with the dissipative part of linear response. If an applied force couples as $-f(t)x$, use response to $e^{-i\omega t}$ and define the two-sided angular-frequency spectrum $S_{xx}(\omega )=\int e^{i\omega t}⟨\delta x(t)\delta x(0)⟩dt$. The classical relation isFor a viscously damped particle, $m\dot{v}=-\zeta v+\xi (t)$, the corresponding white-noise force satisfies $⟨\xi (t)\xi (t^{\prime })⟩=2\zeta k_{B}T\delta (t-t^{\prime })$. Noise replenishes energy lost by friction; their strengths cannot be varied independently at fixed equilibrium temperature.

Electrically, the one-sided voltage-noise spectrum per hertz of an equilibrium resistor is $4k_{B}TR$. A capacitor in parallel filters it toThe same equilibrium voltage variance follows from equipartition. Increasing R makes the spectrum taller and narrower without changing its integrated area. This thermodynamic connection between dissipative conductance and thermal electrical fluctuations is the content of the classical noise calculation.^[ref]

A black body absorbs incident electromagnetic radiation and, in equilibrium, emits a spectrum determined by temperature alone. A cavity with a small aperture approximates this ideal. Photon number is not conserved in thermal equilibrium, so its chemical potential is zero. Combining the Bose occupation with the two-polarization electromagnetic mode density gives the spectral energy densityChanging variable to $x=\hslash \omega /(k_{B}T)$ shows directly why the total energy density is proportional to $T^4$:Three powers of T arise from the number of thermally populated modes and one from their characteristic energy. Isotropic radiation incident on a surface has outward hemispherical flux $cu/4$, givingThe factor 1/4 is the hemispherical angular average; multiplying energy density by c alone would overcount. The spectrum's frequency scale grows with T. Classical equipartition over all modes would instead produce the Rayleigh-Jeans ultraviolet divergence; the Planck occupation suppresses high-frequency modes. Real materials can have wavelength-dependent emissivity, and a body's net radiative exchange also includes absorbed environmental radiation.

An equation of state relates equilibrium pressure, temperature and density. The ideal-gas form $p=n{k}_{B}T$ neglects intermolecular forces and excluded volume. A low-density real-gas expansion iswhere the virial coefficients encode interaction effects and generally depend on temperature.

A simple phenomenological example is the molar Van der Waals equationHere v is molar volume, b represents excluded volume and the attractive interaction lowers pressure through a. Its dilute expansion gives $pv/(RT)=1+[b-a/(RT)]/v+\cdots$. It recovers the ideal-gas limit but is not a quantitatively universal equation for all real gases.

Below its critical temperature, the raw isotherm contains a region with $(\partial p/\partial v{)}_T>0$, corresponding to negative isothermal compressibility and mechanical instability. Equilibrium phase separation replaces that segment by a coexistence pressure satisfying equal chemical potentials, equivalently the Maxwell equal-area construction. The critical point follows from $p_v=p_{vv}=0$:The distinction between an unstable homogeneous branch and a stable two-phase state is essential when interpreting a real-gas equation of state.

Write $z=e^{\beta (\mu -ϵ)}$. For one nondegenerate fermionic quantum state, the allowed occupations 0 and 1 give ${\Xi }_F=1+z$. For a bosonic state, summing all nonnegative occupations gives ${\Xi }_B=(1-z{)}^{-1}$, requiring $z<1$ or $\mu <ϵ$. Differentiating the logarithm,The interacting point defect has two distinct spin states, not one fermionic state. Its four microstates have grand-canonical weights $1,z,z,q{z}^2$, with $q=e^{-\beta U}$. The first two single-electron weights are separate up/down spin possibilities; the double state has one electron of each spin. ThusFor $U=0$, $\Xi =(1+z{)}^2$ and $\bar{n}=2z/(1+z)$. The sketch is a single smooth step from 0 to 2, centred on $\mu =ϵ$, where $\bar{n}=1$. As $T\to 0$ its width tends to zero.

For positive U with $\beta U\gg 1$, the lowest grand energies are successively $0$, $ϵ-\mu$ and $2ϵ+U-2\mu$. Away from the two crossings,Finite temperature smooths each step over an energy range of order $k_{B}T$. Near the lower step double occupation is exponentially suppressed, so $\bar{n}\simeq 2z/(1+2z)$. Near the upper step put $y=e^{\beta (\mu -ϵ-U)}$ and neglect the empty state, giving $\bar{n}\simeq 1+y/(2+y)$. Spin degeneracy shifts the half-step positions to $ϵ-k_{B}T\ln 2$ and $ϵ+U+k_{B}T\ln 2$.

For a grand-canonical probability $p_s\propto e^{-\beta (E_s-\mu n_s)}$, differentiation gives ${\partial }_{\mu }{p}_s=\beta (n_s-\bar{n})p_s$. Multiplying by $n_s$ and summing provesFor this defect, direct state counting givesAt $\mu =ϵ$, $z=1$ and thereforeThe same limit gives $\bar{n}=2/3$: empty, up and down states have equal weights, while the double state is suppressed. Replacing this degeneracy by a two-outcome equal-weight model would incorrectly give variance $1/4$ at this chemical potential.

Across the full chemical-potential range the variance has two well-separated narrow peaks and tends to zero on the 0, 1 and 2 plateaux. The limiting shapes are $2z/(1+2z{)}^2$ near the lower step and $2y/(2+y{)}^2$ near the upper one. Each peak has height $1/4$, at the shifted half-step position; the variance at the unshifted crossings is instead $2/9$.

For one oscillator, set $x=\hslash \Omega /(k_{B}T)$. Its quantum-oscillator energy and heat capacity areAt low temperature, $U\to \hslash \Omega /2$ and $C_V\sim k_B{x}^2{e}^{-x}\to 0$: excitations freeze out but zero-point energy remains. At high temperature, $U\sim k_{B}T$ and $C_V\to k_B$, recovering equipartition. The zero-point term does not contribute to heat capacity.

Using $C_V=T(\partial S/\partial T{)}_V$ and the stated Maxwell relation, interchange derivatives within a smooth single phase:Therefore $C_V$ depends only on temperature for the given temperature-linear equation of state. It need not be constant in temperature.

At $T\ll T_F$, Pauli exclusion restricts thermally changed occupations to an energy interval of order $k_{B}T$ about the Fermi energy. The affected fraction is of order $T/T_F$; their typical excitation energy is of order $k_{B}T$, so the thermal energy correction is proportional to $T^2$ and the heat capacity is linear in $T$.

More precisely, with $g(E_F)$ the total density of one-particle states per unit energy, including spin,For a three-dimensional ideal gas with quadratic dispersion, $g(E_F)=3N/(2E_F)$, henceThe fixed-particle-number Sommerfeld expansion gives $U=3N{E}_F/5+({\pi }^2/4)N(k_{B}T{)}^2/E_F+\cdots$. The temperature dependence of the chemical potential must be included when obtaining this expression; holding it artificially fixed changes the intermediate energy correction.

For a single quantum state, the equilibrium occupations areBoth reduce to $e^{(\mu -ϵ)/(k_{B}T)}$ when $ϵ-\mu \gg k_{B}T$.

For one three-dimensional quadratic branch, the density of states per unit volume isWith ${\lambda }_T=\sqrt{2\pi {\hslash }^2/(m{k}_{B}T)}$, changing integration variable to $x=ϵ/(k_{B}T)$ givesA degeneracy factor multiplies both densities if additional equivalent branches are present. The first expression is the noncondensed density: a macroscopically occupied ground state must be counted separately.

Bose-Einstein condensation occurs when the excited states can no longer hold all particles and the ground-state occupation becomes macroscopic. At fixed total density,For field-induced triplon condensation, write $\Delta =\nu n-\eta (H-H_g)$, so $ϵ(k)={\hslash }^2{k}^2/(2m)+\Delta$. At fixed reference $n$, the three requested dispersions have gaps $\nu n+\eta H_g$, $\nu n$ and zero respectively. Their curvature is unchanged. Because the actual thermodynamic chemical potential is zero, the Bose distribution requires $\Delta \ge 0$. At $H_c=H_g+\nu n/\eta$, the minimum reaches zero and the noncondensed population saturates, permitting macroscopic occupation.

For the triplon condensation field scaling, the critical density at onset is $n_T=\zeta (3/2)/{\lambda }_T^3$. ThereforeThe low-temperature phase boundary rises from $(T,H)=(0,H_g)$; the condensed phase lies on the higher-field side.

For $H>H_c$, simply continuing the gap below zero would give unphysical Bose occupations. For the condensate density in a shifted-parabola triplon model, added condensate density keeps the minimum at zero: $\nu (n_0+n_T)=\eta (H-H_g)$. At fixed temperature this giveswith $n_0=0$ below onset. At zero temperature it becomes $n_0=\eta (H-H_g)/\nu$ within that model.

This prefactor is conditional on continuing the stated single-density mean-field shift into the condensed phase. A microscopic interacting Bose condensate has collective Bogoliubov excitations $E_k=\sqrt{{ϵ}_k^0({ϵ}_k^0+2g{n}_0)}$, with ${ϵ}_k^0={\hslash }^2{k}^2/(2m)$, which are linear at small $k$. The normal-state shift alone does not fix the condensed interaction closure or identify $g$ with $\nu$. The gap closes and the condensate density grows linearly with excess field in the simple mean-field description; the detailed collective dispersion needs that additional interaction model. A dilute-boson treatment of field-induced magnetic ordering makes this normal/condensed distinction explicitly.^[ref]

Use $v$ for molar volume to distinguish it from the fluctuating parcel volume below. The equation is $p=RT/(v-b)-a/v^2$. The parameter $b>0$ models excluded volume, while $a>0$ models the reduction in pressure from attractive interactions.

Below the critical temperature, a raw isotherm has a loop. The positive-slope segment is mechanically unstable and surrounding metastable portions do not give the stable equilibrium mixture. Replace the loop by a constant coexistence pressure $p_s$ satisfying the Maxwell equal-area condition,The endpoints are liquid and vapour molar volumes; intermediate mean volumes correspond to phase mixtures with changing proportions. Above the critical temperature, liquid-like and gas-like states connect continuously without a coexistence plateau.

In the pressure-temperature diagram, the liquid-vapour coexistence line ends at the critical point. Below it, the low-pressure side is vapour and the high-pressure side liquid. “Gas” conventionally denotes the region above the critical temperature, where compression alone cannot cross a liquid-vapour transition. The model has no solid phase.

The critical isotherm has an inflection: $p_v=p_{vv}=0$. Thuswhich giveFor fluctuations about a stable equilibrium value $x_0$, the system plus reservoir has probability proportional to $e^{-A(x)/(k_{B}T)}$. Expanding availability,gives a Gaussian distribution and, to this order, $⟨x⟩=x_0$ withThis requires positive curvature and sufficiently small fluctuations; it is not exact when higher-order terms matter.

The whole cylinder's volume is fixed. For the volume fluctuations of a small fluid parcel, interpret the small quantity as a material parcel containing $s$ moles, whose volume can fluctuate against the much larger remainder. Its mean volume at the critical density is $\bar{V}=s{v}_c=3sb$. At fixed temperature its availability is $A(V)=F_s(V,T)+p_{\mathrm{bath}}V$, soOn the critical isochore, approached from $T>T_c$,ConsequentlyEquivalently, $⟨\delta V^2⟩=k_{B}T{\kappa }_T\bar{V}$, with ${\kappa }_T=4b/[3R(T-T_c)]$. The mean-field variance diverges as $(T-T_c{)}^{-1}$, so its critical exponent is $1$; the root-mean-square fluctuation has exponent $1/2$. A fixed spatial observation box instead has number fluctuations, not fluctuating geometric volume. Finite size and non-Gaussian terms cut off the literal divergence; the approximation fails exactly at the critical point.

Differentiate the Euler relation for a single-component simple compressible system:The thermodynamic fundamental relation for the same system is $dU=TdS-pdV+\mu dN$. Subtraction givesDividing by $N$, with $s=S/N$ and $v=V/N$, proves the Gibbs-Duhem equationThe underlying Euler relation assumes extensive bulk thermodynamics and no additional work coordinates. The result says that the intensive quantities are constrained; the chemical potential cannot be varied independently of temperature and pressure within one homogeneous phase.

Let $E$ denote the fluctuating microscopic energy and $U=⟨E⟩$ its mean, at fixed $N,V$. In a canonical ensemble with $\beta =1/(k_{B}T)$,Since $dT/d\beta =-k_B{T}^2$, this givesHere $\Delta U$ denotes an energy fluctuation about the mean, as in the question. The relation is not restricted to ideal gases when the Hamiltonian operator has no explicit temperature dependence. Alternatively, the stated equilibrium free-energy curvature gives a Gaussian energy weight proportional to $\exp [-(\Delta U{)}^2/(2k_B{T}^2{C}_V)]$, with the same variance.

For the relative energy fluctuation of a monatomic ideal gas, $U=3N{k}_{B}T/2$ and $C_V=3N{k}_B/2$, henceTaking room conditions as approximately $T=300$ K and $p=10^5$ Pa, with $V=10^{-6}$ m${}^3$, givesThe mean energy is about $0.15$ J and its standard deviation about $2.49\times 10^{-11}$ J. The tiny relative fluctuation explains why macroscopic equilibrium internal energy behaves as a sharply defined thermodynamic quantity. The estimate uses a thermal reservoir, not a particle reservoir; allowing $N$ itself to fluctuate changes the relevant fluctuation problem.

A non-ideal gas departs from $p=n{k}_{B}T$ because interactions and finite particle size affect its pressure. At sufficiently low number density $n=N/V$, the virial expansion separates these corrections:The second virial coefficient has dimensions of volume per particle, while $B_3$ has dimensions of volume squared. These are density-series coefficients, not the coefficients of a pressure-series expansion with the same names.

For classical particles interacting through a pair potential $u(r)$,Repulsion contributes positively and attraction negatively. A hard core of diameter $d$ gives $B_2=2\pi d^3/3$, four times the geometric volume of one sphere. The factor reflects the volume excluded to the centre of a second particle together with pair counting; it is not simply the occupied geometric volume.

Expanding the Van der Waals equation at $nb\ll 1$ givesThe virial coefficients of the van der Waals gas are therefore $B_2=b-a/(k_{B}T)$ and $B_3=b^2$. The Boyle temperature, $T_B=a/(k_{B}b)$, makes the leading correction vanish, but does not remove the higher coefficients or turn the gas exactly ideal. More realistic potentials give temperature-dependent higher coefficients that need not equal the van der Waals values.

A useful measurement is a sequence of dilute isothermal pressure–volume observations. Plotting $[p/(n{k}_{B}T)-1]/n$ against $n$ exposes $B_2$ as the zero-density intercept and $B_3$ as its initial slope. Control temperature and density calibration, check equilibration and adsorption on the container walls, and vary the fitted density interval to test truncation error. One measured virial coefficient cannot uniquely reconstruct the pair potential. Near coexistence or a critical point a short density series may fail, and at quantum-degenerate densities exchange statistics add corrections even without interactions.

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.