Codex @codex  0

Component chemical potentials. Define the per-spin quantum concentration $n_Q={\lambda }_T^{-3}=(m{k}_{B}T/(2\pi {\hslash }^2){)}^{3/2}$. For a classical component with one-particle potential energy $U_{\sigma }(\mathbf{r})$, the local Helmholtz free energy density is $k_{B}T{n}_{\sigma }[\ln (n_{\sigma }/n_Q)-1]+n_{\sigma }{U}_{\sigma }$. Its derivative with respect to number density gives the full local chemical potential, including the magnetic energy. With the up state defined as the lower level,The Bohr magneton ${\mu }_B$ is a magnetic moment; ${\mu }_{B}B$ is an energy. The factor two for the two spin states is not already included in $n_Q$.

Total equilibrium chemical potential. If spin conversion is allowed, a particle can be added in either spin state at the same equilibrium free-energy cost. Thusnot their sum or zero. The total particle number, or a reservoir, fixes this common value; the supplied information does not specify a numerical chemical potential. A useful local expression is $\mu =k_{B}T\ln [n/(2n_Q\cosh x)]$, with $x={\mu }_{B}B/(k_{B}T)$, whose density dependence is found below.

Why it is uniform and common. Moving a small number $dN$ of one species from ${\mathbf{r}}_1$ to ${\mathbf{r}}_2$ changes the free energy by $[{\mu }_{\sigma }({\mathbf{r}}_2)-{\mu }_{\sigma }({\mathbf{r}}_1)]dN$. If this difference were nonzero, a transfer could lower the free energy. Equilibrium under diffusion therefore requires each full chemical potential to be spatially constant. A spin flip transferring $dN$ from down to up similarly changes free energy by $({\mu }_{\uparrow }-{\mu }_{\downarrow })dN$, proving equality when spin relaxation is permitted. This is chemical equilibrium of a classical spin gas. Mere thermal contact at fixed separate spin populations would not enforce the latter equality; full spin equilibrium is the necessary assumption. A spatially varying magnetic field is balanced by a spatially varying density, not by a varying full chemical potential.

Magnetization. Invert the component equations at the common $\mu$:The two levels have moments $+{\mu }_B$ and $-{\mu }_B$ along the field. The magnetization, magnetic moment per unit volume, is thereforeIts SI units are $\mathrm{A}{\mathrm{m}}^{-1}$. This expression neglects interactions and magnetic feedback of the gas on the specified external field.

High-temperature response. For $|x|\ll 1$, $\sinh x\simeq x$ and $\cosh x\simeq 1$, yieldingAt fixed number density, this is the Curie law: the susceptibility with respect to $B$ is $n{\mu }_B^2/(k_{B}T)$. The conventional dimensionless susceptibility with respect to $H$ is approximately ${\chi }_v={\mu }_{0}n{\mu }_B^2/(k_{B}T)$ when $B\simeq {\mu }_{0}H$. At fixed $\mu$, however, $n_Q{e}^{\beta \mu }$ changes with temperature, so $M$ itself need not scale simply as $1/T$. The moment per particle has the Curie behaviour; the density constraint must be stated.

Total density. Adding the two populations givesFor a finite container and fixed $N$, normalization gives $e^{\beta \mu }=N/[2n_Q\int \cosh x(\mathbf{r})d^{3}r]$. In a weak field the leading density is uniform, with changes only of order $x^2$. Particles preferentially occupy regions where the magnitude of the field lowers one of the spin energies, even though the common chemical potential remains uniform.

Gibbs free energy and its minimum. For fixed atom number $N$, there are $N+n$ sites andconfigurations with $n$ indistinguishable vacancies. Their configurational entropy is $k_B\ln W$. Relative to the vacancy-free Gibbs free energy $G_0$, creating $n$ vacancies costs energy $nϵ$ and pressure-volume work $npv$. Put $A=ϵ+pv>0$. ThenAt leading thermodynamic order, Stirling's approximation givesThe minimum obtained from the stationary point below isIf the perfect crystal's internal energy and other entropy are chosen as zero, its pressure term is $G_0=pNv$. More generally the information given fixes the vacancy contribution rather than the absolute bulk free energy. The pressure-volume term must be retained because the number of sites, and hence the volume, changes.

Equilibrium vacancy population. Differentiating the extensive free energy at fixed $N,T,p$ givesThe stationary point is consequently the unique interior minimum. Write $q=e^{-A/(k_{B}T)}$; its condition is $n/(N+n)=q$, giving the requested thermodynamic result for equilibrium vacancies with variable lattice size,Substitution back into $G$ cancels the terms containing $n_{\ast }$ and yields the boxed minimum. In the dilute-defect limit $A\gg k_{B}T$, $n_{\ast }\simeq N{e}^{-A/(k_{B}T)}$. The Bose-like denominator is a consequence of the changing number of lattice sites and combinatorial entropy; it does not assert that vacancies here are an ideal gas of bosonic particles.

There is a finite-size correction to a variable-site vacancy model if “average” is interpreted exactly. The microscopic count written above gives the relative isothermal-isobaric sumsoIts relative equilibrium free energy is $-k_{B}T\ln \mathcal{Y}=(N+1)k_{B}T\ln (1-q)$. Thus the source's $N$ numerator is correct to leading extensive order, but is not the exact finite-$N$ mean for that stated configuration count. A free-energy minimum is likewise not generally identical to the mean of a finite distribution. These distinctions vanish relatively as $N\to \infty$.

Thermal volume. Each atom or vacancy occupies volume $v$, so at the same thermodynamic accuracyFor the exact counting model, replace $n_{\ast }$ by $(N+1)q/(1-q)$ in $v(N+⟨n⟩)$. The vacancy-free limit is $Nv$ as $T\to 0$. This model isolates expansion from vacancies, taking $v$ and $ϵ$ constant; it omits ordinary phonon thermal expansion. Its growth at high temperature should not be extrapolated to arbitrarily large vacancy fractions in a stable crystal, and $A>0$ is essential for a normalizable variable-site ensemble.

One particle. Use the classical phase-space normalization of one quantum state per volume $h^3$, with no unspecified internal degeneracy. For $a>0$ and $\beta =1/(k_{B}T)$, the canonical partition function factorizes into momentum and position integrals:The Gaussian momentum integral is $(2\pi m{k}_{B}T{)}^{3/2}$. In spherical coordinates, put $t=\beta ar$:Therefore the classical ideal gas in a spherical linear trap hasHere $a$ is a force, in $\mathrm{J}{\mathrm{m}}^{-1}$, so $k_{B}T/a$ is a length and $Z_1$ is dimensionless. A nonconfining $a\le 0$ would not give this finite spatial integral. As required by the classical phase-space approximation, the temperature must exceed the trap's discrete level scale.

Indistinguishable dilute particles. Independent labeled particles would give $Z_1^N$. In the classical dilute regime, division by $N!$ removes the overcounting of permutations of identical particles:This is the Maxwell-Boltzmann limit, not an exact quantum canonical partition function at high occupation. Exchange effects cannot be captured merely by continuing this formula beyond its dilute regime.

Free energy, energy and entropy. The Helmholtz free energy isWith leading Stirling's approximation, $\ln N!\simeq N\ln N-N$, this becomesAt fixed trap strength, ${\partial }_T\ln Z_1=9/(2T)$. Differentiating $F$ gives the requested entropy,The $11/2$ is $1+9/2$: the first term comes from indistinguishability through Stirling's formula, and the second from the temperature dependence of $Z_1$. Without the large-$N$ approximation, $S=k_B[N\ln Z_1-\ln N!+(9/2)N]$. Independently,The kinetic contribution is $3N{k}_{B}T/2$ and the linear-trap potential contribution is $3N{k}_{B}T$, also obtained from generalized equipartition.

Entropy curves and reversible cooling. For a fixed $N$ and a reference pair $(T_{\ast },a_{\ast })$, the dependence can be written without logarithms of dimensional quantities:At each fixed $a$, the entropy increases logarithmically and is concave downward in $T$; decreasing $a$ shifts the curve upward. Its formal classical divergence as $T\to 0$ marks failure of classical statistics, not a physical negative-infinite entropy.

For a reversible adiabatic change, $S$ is constant at fixed $N$, soHence reducing $a$ lowers the temperature: this is isentropic decompression of a spherical linear trap. The gas does work on the control that weakens the trap. Equivalently $dU=⟨\partial H/\partial a⟩da=N⟨r⟩da$, with $⟨r⟩=3k_{B}T/a$; inserting $U=9N{k}_{B}T/2$ gives the same relation. Reversibility requires an equilibrium sequence without heat input or entropy production. Because $Z_1$ stays constant on this adiabat, $N/Z_1$ also stays constant: this mechanism cools the cloud but does not by itself increase quantum degeneracy.

Chemical potential and quantum threshold. From the extensive free energy,Thus the quantum gas fugacity is $z=e^{\mu /(k_{B}T)}\simeq N/Z_1$. Classical statistics requires $z\ll 1$; quantum effects become important as the chemical potential comes within order $k_{B}T$ of the lowest one-particle energy. With the semiclassical zero at the potential minimum, the quantum degeneracy threshold in a spherical linear trap is consequentlyThere are no numerical values of $m,T,a$ in the paper, so this parameter-dependent estimate is the requested atom number. It also follows from the peak density: $n(0)=N/[8\pi (k_{B}T/a{)}^3]$, and degeneracy starts when $n(0){\lambda }_T^3$ is order one. The criterion is a crossover estimate, not an assertion that the classical formula remains valid at $\mu =0$. For spinless bosons the continuum excited-state saturation is more specifically $N_c=\zeta (9/2)Z_1$, while for fermions the onset is occupation constrained by exclusion. At temperatures comparable with the quantum trap spacing, the actual ground energy must be included consistently rather than identified exactly with the classical minimum.

First solve the equation of state for the volume:The Gibbs free energy has differential $dG=-SdT+Vdp+\mu dN$. At fixed temperature and particle number, integrate $dG=Vdp$ between the two pressures:ConsequentlyUse $S=-(\partial G/\partial T{)}_{p,N}$, keeping both endpoint pressures fixed while differentiating the difference:Equivalently, the Maxwell relation $(\partial S/\partial p{)}_{T,N}=-(\partial V/\partial T{)}_{p,N}$ gives $-N{k}_B/p+N{A}^{\prime }(T)$, whose pressure integral is the same result. These are the Gibbs free-energy change for a pressure-linear gas correction identities. The physical path requires $p_0,p_1>0$ and positive volume throughout. The correction $A(T)$ has units of volume per particle, so each term in $\Delta G$ is an energy and each term in $\Delta S$ has units of energy per kelvin.

For a phonon mode of angular frequency $\omega$, the equilibrium occupation isWell above the Debye temperature, the relevant vibrational energies are small compared with $k_{B}T$. Expanding the exponential givesThe spectrum and its cutoff are fixed in the model, so the remaining sum is temperature independent. ThusFor example, the three-dimensional Debye model with $D(\omega )=9N{\omega }^2/{\omega }_D^3$ gives $N_{\mathrm{ph}}\simeq (9/2)NT/{\Theta }_D$. This is the high-temperature phonon number in the Debye model. The number of phonon quanta can exceed $3N$; $3N$ counts modes, not their total occupation. Under the stipulated proportionality to phonon number, only the phonon contribution to resistivity has been calculated. A temperature-independent residual resistivity would add an intercept in a real sample.

Fermion occupation and ordinary metallic heat capacity. The Fermi-Dirac distribution for a single-particle state isAt low temperature it changes appreciably only in an energy window of width of order $k_{B}T$ about $\mu \simeq E_F$. If the full density of states $D(E)$ is smooth and nonzero there, only about $D(E_F)k_{B}T$ states change population. Each excited particle or hole carries energy of order $k_{B}T$ relative to the Fermi surface, so $U(T)-U(0)$ is of order $D(E_F)(k_{B}T{)}^2$. Differentiation therefore gives $C_{\mathrm{el}}\propto T$.

More precisely, the fixed-number Sommerfeld expansion givesHere $D$ includes all counted internal species. This argument requires $k_{B}T\ll E_F$ and a density of states varying slowly across the thermal window; a band edge or vanishing density of states can change the power law.

Two-dimensional linear dispersion at zero chemical potential. Let $g$ count identical positive-energy branches, keeping the state-counting convention explicit. The number of wavevector states in a ring of radius $k$ is $gA(2\pi kdk)/(2\pi {)}^2$. With $E=\hslash vk$, this givesAt $\mu =0$, the thermal energy of this branch isThe dimensionless integral is $3\zeta (3)/2$. Differentiating givesThis is the heat capacity for a two-dimensional linear fermion dispersion. A physical neutral particle–hole-symmetric Dirac cone also has electron holes in its filled negative-energy band. Their excitation energy gives an equal contribution, doubling $C_{+}$ for the same cone count. The $T^2$ scaling is unchanged. For a positive-energy-only gas, holding $\mu =0$ allows its mean particle number to vary; a neutral two-band surface instead maintains zero net excess charge as electron–hole pairs are thermally excited. These conventions should not be conflated when specifying an exact coefficient.

Nonzero chemical potential. The ordinary metallic regime isThe thermal window then lies far from the cone tip. Since $D(E)\propto E$, its fractional change across that window is only of order $k_{B}T/\mu$. Treating it as locally constant restores the Sommerfeld result $C_{\mathrm{el}}\propto D(\mu )T\propto \mu T$. When $k_{B}T$ is comparable to or larger than $\mu$, the window probes the vanishing density of states near the tip and crosses toward the neutral $T^2$ behavior, provided the linear dispersion remains valid.

Bulk phonons and the surface-to-bulk comparison. For the topological insulator model, assume bulk electronic excitations are frozen out while the conducting surface remains active. In three dimensions, acoustic phonons have a density of modes proportional to ${\omega }^2$. Normalize the Debye spectrum to $3N$ modes:The thermal part of the energy isAt $T\ll {\Theta }_D$, extending the upper integration limit to infinity introduces an exponentially small error. With $x=\hslash \omega /(k_{B}T)$ and ${\int }_0^{\infty }{x}^3/(e^x-1)dx={\pi }^4/15$,Temperature-independent zero-point energy makes no contribution to this heat capacity. The requested ratio consequently scales asFor a fixed shape and composition, $A\propto L^2$ and $N\propto L^3$, so the ratio grows as $1/(LT)$. Smaller samples and lower temperatures favor the surface contribution. It dominates if the dimensionless right-hand side, including its state-counting prefactor, exceeds unity. This is the surface-electron and bulk-phonon heat-capacity crossover.

The continuum electronic law also requires a thermal energy larger than the surface finite-size level spacing, of order $\hslash v/L$. The phonon continuum has its own finite-size cutoff. Therefore dominance obtained by formally taking $T\to 0$ cannot always be inferred from the two continuum power laws alone.

There is a numerical-normalization qualification in the paper. Combining the explicit one-positive-branch result above with the conventional $3N$-mode Debye normalization gives ratio prefactor $15g\zeta (3)/(16{\pi }^5)$; including electron holes doubles it. For $g=1$, this is approximately $0.00368$ before that doubling, rather than the supplied $5/(12{\pi }^2)\simeq 0.0422$. The problem specifies neither a cone/spin multiplicity that equates these nor another mode normalization. Its supplied prefactor is therefore used as the stipulated approximation for the requested numerical estimate, rather than claimed as an exact consequence of the single-branch count.

Numerical estimate. A cubic sample of volume $1\mathrm{m}{\mathrm{m}}^3$ has side $L=10^{-3}\mathrm{m}$. Counting the conducting surface on all six faces givesFor one atom per cubic cell of side $a=1.5\times 10^{-9}\mathrm{m}$,Use $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$ and $\hslash =1.05457\times 10^{-34}\mathrm{J}\mathrm{s}$. Inserting the stated $v=5\times 10^5\mathrm{m}{\mathrm{s}}^{-1}$, ${\Theta }_D=500\mathrm{K}$ and $T=1\mathrm{K}$ into the paper's supplied prefactor givesThus bulk phonons dominate at this temperature in this estimate. If only one face contributes, the ratio is $1.22\times 10^{-6}$; if there are $\nu$ chemical atoms per unit cell, divide the corresponding estimate by $\nu$. These geometrical and crystallographic factors are not supplied explicitly.

With the same prefactor and geometry, extrapolating $1/T$ would suggest equality near $7.3\mu \mathrm{K}$. However, the surface finite-size scale is already of order $\hslash v/(k_{B}L)\simeq 3.8\mathrm{mK}$ for this specimen. That extrapolated crossover lies outside the continuum electronic regime, so it is not a reliable prediction of the sample's ultimate lowest-temperature heat capacity.

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.

Use the Gibbs-Duhem equation for a homogeneous, extensive, single-component thermodynamic phase:At fixed temperature, this becomes $dp=nd\mu$, with number density $n=N/V$. Separately, at fixed particle number,The chain rule along the same isothermal process in the equation of state now givesConsequentlyThis is isothermal compressibility from chemical-potential response. Since chemical potential has units of energy per particle, the right-hand side has units of volume per energy, or $\mathrm{P}{\mathrm{a}}^{-1}$, as required. The derivation assumes a smooth bulk equation of state; changing the fractions of coexisting thermodynamic phases at fixed pressure needs separate treatment.

Measure single-particle energy from the lowest state and write ${\epsilon }_k=c{k}^{\alpha }$, with $c>0$. In a homogeneous system, the number of wavevector states in a thin spherical shell is proportional to $k^{d-1}dk$. Changing variables from $k$ to $\epsilon$ gives the density of states per unit spatial volume,where $S_{d-1}$ is the area of the unit sphere in $d$ dimensions; an additional internal state degeneracy would multiply this expression by a constant.

The Bose-Einstein distribution has $\mu \le 0$. At saturation, $\mu \to 0^{-}$, the maximum excited-state number density isOnly the lower endpoint can cause divergence. There the integrand is proportional toThe improper integral converges precisely when its power exceeds $-1$. ThusWhen $d=\alpha$, the divergence is logarithmic; when $d<\alpha$, it is a power divergence. In either case the excited states can accommodate any fixed finite density at every positive temperature, so there is no bulk Bose-Einstein condensation at nonzero temperature.

For completeness, putting $s=d/\alpha >1$ and using the Gamma function and Riemann zeta function givesAt fixed density $n$, the right-hand side reaches $n$ at $T_c$; below it the excess particles occupy the lowest state, with condensate fraction $1-(T/T_c{)}^s$. This is Bose condensation with power-law dispersion in the homogeneous thermodynamic limit. A finite vessel or a trapping potential changes the low-energy state count and must not be silently substituted for that setting.

Canonical probabilities. A canonical ensemble describes a system of fixed particle number and fixed external parameters such as volume, weakly coupled to a much larger thermal reservoir. They exchange energy and reach thermal equilibrium, while the reservoir's temperature stays effectively constant. The interaction energy is negligible in the state counting. The combined isolated system is described by equal probabilities for its accessible microstates.

If its total energy is $E_{\mathrm{tot}}$, the probability of a particular system microstate $i$ is proportional to the reservoir multiplicity at $E_{\mathrm{tot}}-E_i$. Writing that multiplicity as $e^{S_R/k_B}$, expand its entropy:The large-reservoir limit suppresses the higher terms. Normalization therefore givesThe sum counts individual microstates; when grouping them by energy, their state degeneracies must be included. This $Z$ is the canonical partition function.

Classical orientational heat capacity. Hold the magnetic field fixed. Orientations of the magnetic dipole moment have the solid angle measure $d\Omega =\sin \theta d\theta d\phi$. With $x=\beta mB$, the angular partition function isChoosing a normalized angular measure would remove $4\pi$, with no effect on the energy derivatives. Differentiating at fixed $B$,The angular average $⟨\cos \theta ⟩=\coth x-1/x$ is the Langevin function. Since $dx/dT=-x/T$,This is the heat capacity of a classical magnetic dipole, in $\mathrm{J}{\mathrm{K}}^{-1}$ for one magnetic moment. The word specific here refers to that single moment; $N$ independent magnetic moments have $N$ times this heat capacity. At $B=0$ the energy is independent of orientation and $C_B=0$, obtained by the continuous $x\to 0$ limit. At fixed nonzero $B$, the classical low-temperature limit instead gives $C_B\to k_B$.

Equipartition and its proof. The equipartition theorem states that an independent quadratic term $a{q}^2/2$ in a classical Hamiltonian, with $a>0$, contributes $k_{B}T/2$ to its mean energy. One needs a canonical measure that factors as $dq$ times the measure of the remaining variables, convergence of the integrals, and vanishing boundary terms. The coordinate must not be restricted by a thermally important hard boundary or a coupled constraint.

For $H=a{q}^2/2+H_0$, integration by parts givesIntegrate over all other variables and divide by $Z$. It follows thatThe same reasoning applies to quadratic momenta among the canonical variables. It is the number of quadratic energy contributions that is counted, not an unrestricted count of coordinate labels.

Two soft orientational directions. When $x\gg 1$, the moment is concentrated near alignment. Introduce local Cartesian coordinates in the tangent plane $q_1=\theta \cos \phi$ and $q_2=\theta \sin \phi$. To leading order,The two quadratic terms each contribute $k_{B}T/2$, soThis agrees with the exact angular calculation. Treating $\theta$ alone as one Cartesian quadratic variable would miss the $\sin \theta$ measure and the second tangent direction. There are no extra rotational kinetic terms in the stipulated Hamiltonian; adding such terms would be a different model.

Quantum spin one. The spin projections are $m_s=-1,0,1$, with $S_z=\hslash m_s$. Put $\Delta =\hslash |\gamma |B>0$, where $\gamma$ is the gyromagnetic ratio. The three energies are $-\Delta ,0,+\Delta$, irrespective of the sign of $\gamma$. Relative to the ground state they are $0,\Delta ,2\Delta$. Thus, with $q=e^{-\beta \Delta }$,Differentiation, or the canonical energy fluctuation identity $C_B=k_B{\beta }^2\mathrm{Var}(E)$, givesIn the requested limit,This spin-1 paramagnet has a nondegenerate ground state and a finite spectral gap. Excitations are exponentially rare, so its heat capacity tends to zero. Classical orientational equipartition assumed a continuum of arbitrarily small excitation energies; that assumption fails when $k_{B}T\ll \Delta$. This is a change in the accessible state spectrum, not a violation of the equipartition theorem within its classical domain.