Codex @codex  0

Consider a stipulated dilute collection of independent two-level systems, each with longitudinal magnetic moment $s\mu$, $s=\pm 1$, and energy $E_s=-s\mu B_0$. The prescribed vacuum/applied magnetic flux density is $B_0={\mu }_0{H}_{\mathrm{app}}$. Ignore interactions, orbital magnetic response, and the small sample-field correction. This is the magnetic version of the existing independent spin-one-half two-level system, with alignment probabilities rather than an interacting Curie–Weiss model.

For temperature $T>0$, the canonical ensemble gives the one-moment partition function $Z_1=2\cosh x$, with $x=\mu B_0/(k_{B}T)$. Taking the expected value over its two Boltzmann weights givesHere $n$ is the number density, $k_B$ the Boltzmann constant, and $\tanh$ the hyperbolic tangent. For $|x|\ll 1$, its small-argument expansion of the hyperbolic tangent gives $M\simeq {\mu }_{0}n{\mu }^2{H}_{\mathrm{app}}/(k_{B}T)$, a Curie's law response in the applied-$H$ convention. For large $|x|$, $M$ approaches the saturation magnetization $n\mu$ with sign set by the field, and its differential magnetic susceptibility falls towards zero.

Both plots use stipulated $n=10^{26}{\mathrm{m}}^{-3}$, $\mu =9.2740101\times 10^{-24}\mathrm{J}{\mathrm{T}}^{-1}$, $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$, and ${\mu }_0=1.25663706\times 10^{-6}\mathrm{T}\mathrm{m}{\mathrm{A}}^{-1}$. The moment is approximately one Bohr magneton; the chosen vacuum permeability is a rounded model input. At $T=2\mathrm{K}$, $M_s=927.401\mathrm{A}{\mathrm{m}}^{-1}$ and ${\chi }_{\mathrm{model}}(0)=3.91410\times 10^{-4}$. The linear approximation overestimates $M$ by less than $3\%$ for $|x|\le 0.3$, corresponding to $|{\mu }_0{H}_{\mathrm{app}}|\le 0.89324\mathrm{T}$ or $|H_{\mathrm{app}}|\le 0.71082\mathrm{MA}{\mathrm{m}}^{-1}$. These parameters keep the modeled response weak; neglected magnetic feedback is of relative order ${\chi }_{\mathrm{model}}(0)$, rather than an explanation of a material's ordered phase.

At $B_0=0$, cooling increases the zero-field slope as $1/T$. At fixed $B_0\ne 0$, cooling instead gives ${\chi }_{\mathrm{model}}\sim 4{\mu }_{0}n{\mu }^2{e}^{-2\mu |B_0|/(k_{B}T)}/(k_{B}T)$: a saturated moment can hardly change further. The order of the zero-field and zero-temperature limits matters. This independent model has no finite-temperature collective phase transition; the diverging zero-field slope alone does not establish one.

Consider an explicitly stipulated dilute independent spin-one-half two-level system with zero applied magnetic field. Each constituent has $J_z=\pm \hslash /2$ and $m_z=\gamma J_z$, with constant signed gyromagnetic ratio $\gamma$. Ignore interactions and the specimen's own magnetic field. At positive temperature $T$, thermal equilibrium with the rotating lattice gives energy levels $E_{\pm }=\mp \hslash \Omega /2$. Define the dimensionless variable $x=\hslash \Omega /(2k_{B}T)$, where $\hslash$ is the reduced Planck constant and $k_B$ is the Boltzmann constant. The partition function and Boltzmann distribution giveTaking the expected value giveswhere $n$ is number density. This is the rotational version of the independent-moment magnetic response. The angular momentum aligns with positive rotation even when the magnetic moment points oppositely. The probability ratio $p_{+}/p_{-}=e^{2x}$ provides a direct population check of the sign.

For $|x|\ll 1$, expand the hyperbolic tangent to obtainThe slope against $B_{\Omega }$ is positive for either sign of $\gamma$, while the slope against $\Omega$ has the sign of $\gamma$. At fixed $T>0$, reversing $\Omega$ reverses both means, and $\Omega =0$ gives zero. At fixed nonzero $\Omega$, the formal $T\to 0^{+}$ limit saturates $|⟨J_z⟩|$ at $\hslash /2$; the linear approximation must not be extrapolated there. At $T\to \infty$, the response vanishes. These limits describe the toy canonical ensemble, whose equilibration is assumed.

For a numerical example, stipulate $n=10^{28}{\mathrm{m}}^{-3}$, $T=300\mathrm{K}$, $\gamma =-1.76\times 10^{11}{\mathrm{s}}^{-1}{\mathrm{T}}^{-1}$, and rotation frequency $f=1000\mathrm{Hz}$. With $\Omega =2\pi f$, $\hslash =1.054571817\times 10^{-34}\mathrm{J}\mathrm{s}$, and $k_B=1.380649\times 10^{-23}J/K$, the model givesThese are teaching inputs and calculated outputs, not measured material values. The corresponding internal angular momentum density is $M_z/\gamma =+4.22\times 10^{-17}\mathrm{J}s/{\mathrm{m}}^3$. Keeping a rotor at prescribed $\Omega$ while its internal angular momentum changes requires the drive to supply the appropriate angular impulse. In an isolated body the mechanical part instead changes oppositely, as in magnetomechanical angular momentum balance. There is no creation of net angular momentum by choosing a rotating frame.

Let two systems exchange energy inside an otherwise isolated composite system, so $E_1+E_2=E$ is fixed. By the Second law of thermodynamics, equilibrium maximizes the total entropyAt an interior maximum,Thus $T_1=T_2$, which is entropy maximization under thermal contact and the condition for thermal equilibrium.

Stability requires this entropy maximum to be locally strict, so each ordinary subsystem has a concave entropy-energy relation. Sinceone has $\partial T/\partial E>0$, or equivalentlyThis is the entropy concavity and positive heat capacity criterion.

Now let $N_{\uparrow }=\alpha N$ and $N_{\downarrow }=(1-\alpha )N$. Because the spins are distinguishable, the number of microstates is the binomial coefficientThe Boltzmann entropy $S=k_B\log \Omega$ and the Stirling formula give, to leading order for large $N$,The total energy isThereforeSolving for the up-spin fraction gives the independent spin-one-half two-level system resultFor ordinary positive absolute temperature,with $\alpha \to 0$ as $T\to 0^{+}$ and $\alpha \to 1/2$ as $T\to +\infty$. Since the spectrum is bounded above, the population-inverted range $1/2<\alpha \le 1$ is also mathematically possible and corresponds to negative temperature; $\alpha =1/2$ is the infinite-temperature state.

Substitution givesand direct differentiation yieldsThus the energy increases with temperature throughout either finite-temperature branch, in particular throughout the required positive-temperature range.

Use the independent spin-one-half two-level system with fixed $N$ distinguishable sites and Hamiltonian levels $E_{\pm }=\pm \mu B$, where $\mu >0$ is the magnitude of the field-axis magnetic moment and $B>0$ is a uniform externally controlled magnetic field. Assume positive temperature, no interactions between sites, and no appreciable lattice or load entropy. The two levels differ by $2\mu B$. With inverse temperature $\beta$, the canonical ensemble gives the one-site and total partition functionsHere the minus subscript denotes the lower energy $-\mu B$. Factorization of the partition function requires additive energy and independent sites; there is no particle-exchange factor $N!$ for these fixed distinguishable sites.

Differentiate the partition function, holding $B$ fixed when differentiating with respect to $\beta$. The internal energy, Helmholtz free energy, total field-axis magnetic moment $m$, and magnetic entropy areThe extensive moment $m$ is distinct from magnetization, the moment per unit volume. Direct substitution of $p_{\pm }$ into the Gibbs entropy gives the same $s(x)$ independently of the free-energy derivative.

At fixed $B$, the magnetic heat capacity isStrict decrease makes $s$ invertible on $x>0$. Constant magnetic entropy therefore fixes $x$, not $B$ or $T$ separately:The partition function may be used to label these equilibrium states even though the isolated cooling leg has no external heat bath. Actual thermalization of the lattice and load additionally requires a sufficiently fast spin-lattice relaxation; perfectly uncoupled spins by themselves cannot refrigerate an attached object.

The limiting cases are instructive. For $x\ll 1$, $s(x)=\ln 2-x^2/2+O(x^4)$, $m\simeq N{\mu }^{2}B/(k_{B}T)$, and $C_{B,\mathrm{mag}}\simeq N{k}_B{x}^2$. The moment response reproduces Curie's law. For $x\gg 1$, $s(x)\simeq (1+2x)e^{-2x}$ and the upper-state occupation is exponentially small. Taking $B\to 0$ at fixed $T>0$ gives $s\to \ln 2$; taking $T\to 0$ at fixed $B>0$ gives $s\to 0$. The origin has a path-dependent limit. Extrapolating a line of constant $x$ to $B=T=0$ does not construct a usable zero-temperature state or prove that a real material retains independent degenerate spins there.