Codex @codex  0

Put $\rho ={\phi }^{†}\phi$. The heavy-field equation of motion isThis solution describes the low-energy sector without an independent incoming heavy-particle wave. It is valid when all external momentum invariants and light masses squared are small compared with $M^2$.

After an integration by parts, the terms involving the heavy field are $-\frac{1}{2}\varphi (\square +M^2)\varphi -g\rho \varphi$. Substituting the solution into both terms, rather than the source term alone, givesHence, with the question's sign convention for the local quartic interaction,This is integrating out a field followed by a derivative expansion. The negative quartic means attraction. It also makes the truncated potential unbounded at large fields; minimizing the original potential over $\varphi$ reveals the same negative quartic. Thus neither displayed model, without extra stabilizing interactions, supplies a globally stable nonperturbative vacuum. The matching describes perturbative low-energy scattering about the local vacuum. The simple equation-of-motion matching fixes the leading tree-level coefficient; precision quantum matching requires consistent loop renormalization and the additional effective operators at the requested accuracy.

Use the unitary discrete Fourier transformThe transformed operators obey $[a_k,a_q^{†}]={\delta }_{kq}$. The hopping terms contribute $2A\cos k$ to the number coefficient. In the pairing sum, the odd sine part cancels between $k$ and $-k$. ThuswhereThe factor one half in the full momentum sum avoids counting an anomalous pair twice. For $H_f>0$, ${\alpha }_k>|{\Delta }_k|$ for every allowed momentum. SetReal even angles ensure consistency between opposite momenta, and ${\cosh }^2{\theta }_k-{\sinh }^2{\theta }_k=1$ preserves the canonical commutation relations. The anomalous coefficient becomes ${\Delta }_k\cosh 2{\theta }_k-{\alpha }_k\sinh 2{\theta }_k=0$. Normal ordering then giveswithThe field-polarized planar ferromagnetic spin wave has energy gap $ϵ(0)=S\sqrt{H_f(H_f+2J)}$. The momenta run over the whole Brillouin zone once, including the self-inverse momenta $0$ and, for even $L$, $\pi$; these use single-mode squeezing rather than introducing an extra independent partner.

The zero-field case has a zero-mode obstruction to bosonic Bogoliubov diagonalization. Its nonzero-momentum modes have $ϵ(k)=2\sqrt{2}JS|\sin (k/2)|$, but ${\alpha }_0={\Delta }_0=JS$, so no finite squeezing angle exists for the uniform mode. With $Q_0=(a_0+a_0^{†})/\sqrt{2}$, its Hamiltonian is $JS{Q}_0^2-JS/2$, not an identically zero oscillator. For example, its matrix element between the vacuum and the two-boson number state is $JS/\sqrt{2}$, whereas $ϵ(0)b_0^{†}{b}_0$ would vanish. Consequently the requested oscillator diagonalization is literal for nonzero field; at zero field treat the collective orientation separately or introduce a field regulator before taking a limit.

Choose the current convention corresponding to $\delta \phi =i\alpha \phi$. Noether's theorem givesFor the requested terms with no $\sigma$ and $f=\sqrt{2}v$, the Lagrangian is $\frac{1}{2}(\partial \pi {)}^2$ and the current becomes $J^{\mu }=f{\partial }^{\mu }\pi$. The supplied covariant oscillator normalization cancels the invariant measure when the annihilation operator contracts with $a^{†}(p)$, soThe momentum is on the massless shell, $p^2=0$. Differentiating once more therefore givesThis Goldstone current matrix element in a phase-only scalar theory shows that a conserved current can have a nonzero matrix element proportional to a null momentum. Thus current conservation is fully consistent with the current creating a massless one-particle state from the vacuum; the nonzero coefficient $f$ is its Goldstone-current coupling. It is the divergence that vanishes, not the current matrix element. Reversing the current convention reverses the nonzero matrix element but leaves this zero divergence unchanged.

Varying the vector field gives the Proca equation,Its divergence uses antisymmetry of $F^{\mu \nu }$ to cancel the first term. Since $m_A\ne 0$,For an on-shell momentum this imposes $p_{\nu }{ϵ}^{\nu }=0$ on the four polarization components, leaving three independent solutions. In the rest frame, $p=(m_A,0)$ makes ${ϵ}^0=0$, leaving precisely the three spatial polarizations of a massive spin-one particle.

There is no fourth independent oscillator hidden in $A^0$. Its canonical momentum is zero because the Lagrangian has no ${\dot{A}}^0$, and the temporal equation constrains $A^0$ in terms of the spatial data. Equivalently the primary and secondary second-class constraints remove two variables from the eight-dimensional phase space, leaving six phase-space variables, or three physical degrees of freedom.

Take $p^{\mu }=(e,0,0,k)$ with $e^2=k^2+m_A^2$. A real orthonormal basis that satisfies the source's stated bilinear normalization isEach obeys $p\cdot {ϵ}_i=0$. The two transverse norms are $-1$, while ${ϵ}_L^2=(k^2-e^2)/m_A^2=-1$, and all cross products vanish. For momentum in the opposite direction, one consistent longitudinal choice is ${ϵ}_L(-k)=(k,0,0,-e)/m_A$ for $k>0$; an overall sign of any polarization is immaterial.

Here $+$ and $-$ label the two linear transverse basis vectors. If those symbols are intended to denote circular helicities, choose ${ϵ}_{\pm }^{\mathrm{hel}}=({ϵ}_x\pm i{ϵ}_y)/\sqrt{2}$ instead. The corresponding orthonormality relation must contain complex conjugation, ${ϵ}_i^{\ast }\cdot {ϵ}_j=-{\delta }_{ij}$. Without that conjugation, a circular polarization has ${ϵ}_{+}^{\mathrm{hel}}\cdot {ϵ}_{+}^{\mathrm{hel}}=0$, a concrete counterexample to interpreting the printed bilinear condition as a helicity normalization.

A spherical three-momentum ultraviolet cutoff givesFor $m>0$ and $\Lambda \gg m$, the asymptotic expansion isIn particular, the leading ultraviolet divergence is quartic. For $m=0$, $\mathcal{E}={\Lambda }^4/(16{\pi }^2)$ exactly.

Write $\phi =⟨\psi ⟩$. Under the proposed identification, the scalar potential from a field-dependent momentum cutoff is $V_{\mathrm{ind}}(\phi )=\mathcal{E}(|\phi |)$. It is even and nonnegative, and its derivative with respect to $|\phi |$ is $|\phi {|}^2\sqrt{{\phi }^2+m^2}/(4{\pi }^2)>0$ away from zero. Consequently its only global minimum is at $\phi =0$: this contribution alone does not cause spontaneous symmetry breaking. At large $|\phi |$ it grows quartically. At small $|\phi |$ with fixed positive $m$, however,Thus the quartic sketch is an asymptotic shape, not the exact small-field behavior for a massive original field. Identifying a regulator with a new field is an additional modeling prescription; an actual effective potential also depends on the new field's dynamics and any other allowed terms.

Use one common ultraviolet cutoff for all physical degrees of freedom, with positive sign for each bosonic degree of freedom and negative sign for each fermionic degree of freedom. For the mass-moment cancellation of vacuum energy divergences, define the weighted supertrace $\mathrm{Str}f(m)={\sum }_{B}f(m_B)-{\sum }_{F}f(m_F)$. The large-momentum expansion isAfter multiplication by the radial measure $p^{2}dp$, these three terms produce quartic, quadratic and logarithmic ultraviolet divergences. Therefore the necessary and sufficient cancellations for the displayed unsubtracted finite-species vacuum energy areThey respectively mean equal numbers of bosonic and fermionic physical degrees of freedom, equal sums of squared masses, and equal sums of fourth powers of masses. Once they hold, the remaining radial integrand falls as $p^{-3}$ and is integrable at infinity. Nonnegative masses introduce no lower-end divergence.

Pairing every bosonic degree of freedom with a fermionic one of the same mass, as in an unbroken supersymmetric spectrum, is sufficient and cancels the entire expression. It is stronger than these three moment conditions. For a concrete non-paired example, take bosonic squared masses $\{0,3,3\}{m}_{\ast }^2$ and fermionic squared masses $\{1,1,4\}{m}_{\ast }^2$. Both have three degrees of freedom, squared-mass sum $6m_{\ast }^2$ and fourth-power-mass sum $18m_{\ast }^4$, but their spectra differ. If a vacuum energy counterterm is instead permitted, a renormalized finite answer need not satisfy them; that is a different prescription from finiteness of the source's bare sum.

With no pairing, the Hamiltonian is diagonal in the fermionic occupation numbers. The Pauli exclusion principle allows each momentum mode to be either empty or occupied once. Without a fixed fermion number constraint, occupying a mode lowers the energy precisely when its single-particle energy is negative. Therefore, with $|0⟩$ the Fock vacuum, a ground state isChoose a fixed ordering of the fermionic creation operators; changing it alters only the overall phase. Any zero-energy modes may be occupied or empty without changing the ground-state energy.

To identify the occupied modes, complete the square:Here $\beta$ is the dispersion parameter, not inverse temperature. If ${\beta }^2<4\mu$, every mode has positive energy and the unique ground state is the Fock vacuum, with $E_{\mathrm{GS}}=0$. If ${\beta }^2=4\mu$, there are still no negative-energy modes, but $p=\pm \sqrt{\beta /2}$ have zero energy if they lie on the allowed momentum grid. Each such mode doubles the ground-state degeneracy.

For ${\beta }^2>4\mu$, defineBoth are positive, and the occupied Fermi sea consists of the two intervals $p_{-}<|p|<p_{+}$. With $p=2\pi n/L$, the requested finite-size ground-state energy is explicitlyThe factor two counts opposite momenta. Endpoints are excluded from this sum because they have zero energy; when an endpoint is allowed, its optional occupation adds ground-state degeneracy. A narrow negative-energy interval can contain no allowed finite-size momenta, in which case the sum is empty. The zero-momentum mode is always empty since $ϵ(0)=\alpha \mu >0$.

For real angles write $u_p=\cos {\theta }_p$ and $v_p=\sin {\theta }_p$. For one independent pair $p>0$, the transformation and its conjugate readUsing the canonical anticommutation relations for the $\alpha$ modes givesThe remaining anticommutators vanish for independent momentum pairs. The displayed two-by-two matrix is unitary and its inverse is obtained by changing $i$ to $-i$.

When the same formula is indexed over both signs of momentum, its first row at $-p$ must agree with the conjugate of its second row at $p$. Thus a consistent global convention requiresThis is the odd-momentum fermionic Bogoliubov angle convention. If one checks only the first-row formulas, their cross anticommutator is $i\sin ({\theta }_p+{\theta }_{-p})$; its vanishing gives the condition modulo $\pi$, with possible signs absorbable in the $\alpha$ modes. Requiring both rows to describe the same operators fixes the convention above. Merely checking a separate unitary matrix at each momentum misses this redundancy: choosing ${\theta }_p={\theta }_{-p}=\pi /4$ in both first rows would give $\{c_p,c_{-p}\}=i$.

At $p=0$, the annihilation operator and its conjugate belong to one mode. The full displayed formula is consistent only with $v_0=0$; choose ${\theta }_0=0$ and ${\alpha }_0=c_0$. This mode is not a two-mode pair.

Since every quasiparticle energy is positive, the ground state is the vacuum of all the ${\alpha }_p$ annihilation operators, including ${\alpha }_0=c_0$. Its zero-momentum mode is therefore empty. For each positive momentum, the inverse Bogoliubov transformation givesBoth annihilate the normalized pair state $(u_p+i{v}_p{c}_p^{†}{c}_{-p}^{†})|0⟩$. For example, applying ${\alpha }_p$ produces the two cancelling terms $i{u}_p{v}_p{c}_{-p}^{†}|0⟩$. ThusThe factors have even fermion parity and commute for disjoint momentum pairs. Each normalized factor is a superposition of zero and two fermions, whose pair-occupation probability is $v_p^2$. ThereforeThere is no zero-mode contribution. For fixed $L$, the sum converges since $v_p^2=O(p^{-4})$.

A normalized simple excited state is $|k⟩={\alpha }_k^{†}|\Omega ⟩$ for $k>0$. Its excitation energy is $E_k$, and its total momentum is $k$: both terms in ${\alpha }_k^{†}=u_k{c}_k^{†}+i{v}_k{c}_{-k}$ add momentum $k$. Acting on the $k,-k$ pair givesThe affected pair now contains exactly one fermion instead of mean number $2v_k^2$; all other pairs are unchanged. ConsequentlyThis is the fermion number of a Bogoliubov excitation. Its change in mean fermion number lies between zero and one for $\mu >0$, with the value one at $\Delta =0$. The Bogoliubov quasiparticle coherently mixes adding a fermion and removing the opposite-momentum fermion; it is not an eigen-excitation of the original fermion number operator. Indeed,so fermion number is not conserved when $\Delta \ne 0$, while fermion parity is conserved. The generally noninteger change is an expectation value in a coherent paired background, rather than a fractional eigenvalue of the fermion number operator. In the weak-pairing limit at fixed $k$, it tends to one, as expected for an ordinary added fermion.

For a free real scalar field in four spacetime dimensions, $\mathcal{L}=(\partial \varphi {)}^2/2-m^2{\varphi }^2/2$ and $\pi =\dot{\varphi }$. Canonical quantization promotes these canonical fields to Hermitian operator-valued distributions with equal-time canonical commutation relationsExpand the solution of the Klein-Gordon equation into positive- and negative-frequency modes:Hermiticity pairs the two coefficients as adjoints. The field commutators are equivalent toFor example, inserting the expansion in $[\varphi ,\dot{\varphi }]$ yields two equal Fourier integrals, each with coefficient $i/2$, and hence $i{\delta }^{(3)}$. The vacuum obeys $a(\mathbf{k})|0⟩=0$, and repeated creation operators generate the symmetric bosonic Fock space. Normal ordering of the quadratic Hamiltonian removes its vacuum constant and gives $H=\int d{\Pi }_k{\omega }_k{a}^{†}a$; the same computation for the spatial momentum gives the stipulated four-momentum $P^{\mu }=\int d{\Pi }_k{k}^{\mu }{a}^{†}a$.

The oscillator algebra now impliesConsequentlyIndeed, ${\partial }^{\mu }{e}^{\mp ikx}=\mp i{k}^{\mu }{e}^{\mp ikx}$. For $\mu =0$, the result is $[H,\varphi ]=-i\dot{\varphi }$, which is the Heisenberg equation of motion with the correct sign. As usual, the field identities are interpreted after smearing the operator-valued distributions, or on a suitable common finite-particle domain.

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.

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.

Occupation and the two-dimensional Fermi energy. For an independent single-particle state, the Fermi-Dirac distribution isThe chemical potential enforces the chosen mean particle number through $N={\sum }_{j}n({\epsilon }_j)$, including spin states in the sum. It is common to all states in particle-exchange equilibrium and must generally change with temperature or field when the total mean number is held fixed. At zero temperature, the distribution becomes a Heaviside step function at the Fermi energy.

Let the sample area be $\mathcal{S}$ and let $g_s$ be its internal spin spin degeneracy. In the continuum limit, filling the disk $k<k_F$ givesThereforeFor spin-half fermions with both spin states available, $g_s=2$ and ${\epsilon }_F=\pi {\hslash }^{2}N/(m\mathcal{S})$; for a spin-polarized gas $g_s=1$. The question does not specify the area or spin degeneracy, so these conventions must be stated rather than extracting an energy from $N$ alone. Although the original text calls $\mathbf{k}$ momentum, the supplied ${\hslash }^2{k}^2/(2m)$ dispersion makes it a wave vector; physical momentum is $\mathbf{p}=\hslash \mathbf{k}$. The result agrees with the constant two-dimensional free-electron density of states. Discrete level spacings in a finite sample would require a specified boundary spectrum.

One level occupied on average. Write $\delta ={\mu }_{B}B$ and $z=e^{\beta (\epsilon -\mu )}>0$. The two mean spin occupations areImposing $N_{+}+N_{-}=1$ gives $2+z(e^{\beta \delta }+e^{-\beta \delta })=1+z(e^{\beta \delta }+e^{-\beta \delta })+z^2$, hence $z^2=1$. Thus, for every positive temperature,The total magnetic moment is $\mathcal{M}=-{\mu }_B(N_{+}-N_{-})$. ConsequentlyThis is the grand-canonical spin susceptibility of a half-filled single level. Empty, singly occupied and doubly occupied orbital states are allowed in the grand canonical ensemble. At zero field, the singly occupied configurations together have probability $1/2$, explaining the factor of two relative to one permanent spin-half moment. If exactly one electron were required at all times, a canonical two-state partition function would instead give $\mathcal{M}={\mu }_B\tanh (\beta {\mu }_{B}B)$ and $\chi ={\mu }_B^2/(k_{B}T)$. That different answer does not implement the stated average-number constraint.

Density-of-states formula for the susceptibility. Here $g$ is the zero-field single-spin density of states. Regard its argument $E$ as the full spin-shifted energy. The shifted functions are $g_{\pm }(E)=g(E\mp {\mu }_{B}B)$, soThe total number is even in $B$; thus $\mu (B)=\mu (0)+O(B^2)$ at fixed total number on a smooth branch. More directly, terms from ${\mu }^{\prime }(0)$ cancel in the spin difference because the two densities coincide at zero field. To first order,It follows thatFor a physical energy band, take $g$ to vanish below its lower edge. Provided the endpoint term vanishes, integration by parts also givesThis form displays the positive response and the narrow thermal window around $\mu$. A sharp band edge is included by the corresponding distributional derivative in the first expression. The susceptibility defined here is the derivative of total magnetic moment with respect to $B$, with units $\mathrm{J}{\mathrm{T}}^{-2}$; it is not automatically a dimensionless SI volume susceptibility. A density-of-states convention per unit volume likewise gives a response per unit volume. These are the field conventions for magnetic susceptibility.

Chemical potential at fixed number. Define the integrated single-spin state countThe probability kernel $K_T$ has unit integral, mean $\mu$, and variance $a(k_{B}T{)}^2$. Using integration by parts,Expand $G$ about $\mu$. The linear term integrates to zero, while $G^{\prime \prime }=g^{\prime }$ givesThis is the number form of the Sommerfeld expansion. Since $N/2=G({\epsilon }_F)$ at zero temperature, set $\mu ={\epsilon }_F+\delta \mu$ with $\delta \mu =O(T^2)$:ThereforeThe coefficient can be obtained directly. With $y=(E-\mu )/(k_{B}T)$,For example, the sum follows by expanding $e^{-y}/(1+e^{-y}{)}^2$ on $y>0$, using the evenness of the kernel and integrating each exponential. Thus $a/2={\pi }^2/6$. A density of states increasing with energy gives a downward quadratic chemical-potential shift. For the constant two-dimensional density of states that shift vanishes to this order; a finite lower band edge still allows exponentially small corrections.

Low-temperature susceptibility at fixed number. Apply the same kernel expansion to $g$ itself:The shift of $\mu$ must be retained in the first term. Writing $g_F=g({\epsilon }_F)$ and similarly for its derivatives,ConsequentlyThis is the fixed-number low-temperature Pauli susceptibility. The $g_F^{\prime \prime }$ contribution is thermal broadening at fixed chemical potential; the $-g_F^{\prime 2}/g_F$ contribution enforces fixed total number. Both are needed. If the density of states were instead defined with both spin directions included, the prefactor $2$ would be absorbed into that definition. Only the spin response has been calculated; no orbital magnetic response is included.

The expansion assumes $g_F>0$, sufficient smoothness near the Fermi energy, and a temperature small compared with the energy scale over which $g$ changes or encounters a band edge. It does not apply unmodified at a singular density of states. For the ideal two-dimensional parabolic band, $g_F^{\prime }=g_F^{\prime \prime }=0$, so the leading Pauli response is temperature independent, consistently with the earlier state count.

Let the gas have a small bulk velocity $u_x(y)$, an approximately isotropic local thermal velocity distribution, and a mean free path $\lambda$ short compared with the flow-gradient scale. A molecule transports the mean streamwise momentum acquired near its previous collision. To first order, its nonequilibrium momentum relative to the destination layer is approximately $-m{\tau }_c{v}_{y}d{u}_x/dy$. The resulting transverse momentum flux isThus the dynamic viscosity is $\eta \simeq \rho {\tau }_c⟨v_y^2⟩$. With isotropy, a representative collision time ${\tau }_c\sim \lambda /⟨v⟩$, and a chosen thermal speed average, this has the formThe dimensionless $K$ depends on how speeds and collision times are averaged; a simple one-speed estimate gives $K=1/3$. Assumptions include dilute binary collisions, a locally equilibrated distribution, weak shear, and negligible boundary slip. The kinetic theory identity $p=\rho ⟨v^2⟩/3$ then givesFor the Maxwell speed distribution, $⟨v^2⟩/⟨v{⟩}^2=3\pi /8$, soThis explains both the units of inverse time and the relation to molecular collision frequency, up to the averaging convention.