Codex @codex  0

For $\mathbf{B}=B{\mathbf{e}}_z$, the quantum Gibbs state is diagonal in the $S_z$ basis. Its two probabilities are proportional to $e^{\pm \beta \gamma \hslash B/2}$. ThusAt small inverse temperature, $⟨S_z⟩\simeq \beta \gamma {\hslash }^{2}B/4$. As a function of $u=\beta \gamma \hslash B/2$, the normalized spin component is the smooth odd curve $\tanh u$: it crosses zero at $u=0$ and saturates at $\pm 1$. As a function of $\beta$ its direction is reversed if $\gamma B<0$. For $\gamma B>0$, positive temperature favors $S_z=+\hslash /2$, while negative temperature gives the inverted population.

Let $\beta =(k_{B}T{)}^{-1}$ and sum over all $2^{2L}$ configurations of the Ising spins. The canonical partition function and Helmholtz free energy per site areBoth spins in each cell count as sites. In particular the denominator is $2L$, rather than $L$.

The mean spin magnetization per site is the Boltzmann distribution averageEquivalently introduce a uniform magnetic field through $E_h=E-h{\sum }_j({\sigma }_j+{\tau }_j)$ and differentiate:This convention absorbs the microscopic magnetic moment into $h$. Simultaneous reversal of every spin leaves $E$ invariant and changes the sign of the numerator, so $⟨m⟩=0$ at zero field for finite $L$.

For each pair of adjacent spins, define the transfer matrix for a classical spin chain by $T_{s{s}^{\prime }}=e^{-\beta E_0(s,s^{\prime })}$. Multiplying its entries gives the Boltzmann factor of a whole configuration. The periodic boundary condition identifies the last index with the first, soThe assumed interchange symmetry of $E_0$ makes $T$ a real symmetric matrix, hence it has an orthonormal eigenbasis. At finite positive temperature, every matrix entry is strictly positive, so the Perron–Frobenius theorem supplies a unique positive largest eigenvalue ${\lambda }_{+}$, larger in modulus than every other one. Therefore $L^{-1}\log Z_L\to \log {\lambda }_{+}$.

Local observables become matrix insertions. For example, with $S_{s{s}^{\prime }}=s{\delta }_{s{s}^{\prime }}$,Thus partition functions and correlation functions reduce to a fixed-dimensional eigenvalue problem, rather than a sum over exponentially many microstates.

For the sawtooth Ising chain, each apex spin occurs in only one triangle. Put $K=\beta J$. Summing that spin first gives the effective transfer matrixIn the ordered basis $s=+1,-1$ this isIts symmetric and antisymmetric eigenvectors are $(1,1)/\sqrt{2}$ and $(1,-1)/\sqrt{2}$, with eigenvaluesConsequentlyAt $J=0$, the result is $Z_L=4^L$, correctly counting $2L$ independent Ising spins. As $T\to 0$, $f_{\infty }\to -3J/2$, the energy of three satisfied ferromagnetic bonds per two sites.

The printed summand lacks grouping brackets. The result above uses the translation-invariant common-$J$ interpretation. If base and apex bonds instead have independent couplings $J_b,J_a$, the same derivation gives ${\lambda }_{\pm }=2e^{\beta J_b}\cosh (2\beta J_a)\pm 2e^{-\beta J_b}$; specifying those couplings is essential before giving a different answer.

Let $S=\mathrm{diag}(1,-1)$ be the base-spin insertion matrix. In the eigenbasis from solution, its diagonal matrix elements vanish: $⟨+|S|+⟩=⟨-|S|-⟩=0$. Hencenot only at large $L$ but for every finite ring at zero field.

This vanishing is not merely an averaging artifact hiding a finite-temperature ordered phase. The two-point correlation function is $⟨{\sigma }_1{\sigma }_{1+r}⟩=({\lambda }_{-}^r{\lambda }_{+}^{L-r}+{\lambda }_{+}^r{\lambda }_{-}^{L-r})/Z_L$. In the thermodynamic limit it becomes $({\lambda }_{-}/{\lambda }_{+}{)}^r$, with finite correlation length ${\xi }^{-1}=\log ({\lambda }_{+}/{\lambda }_{-})$ for every $T>0$ and $J>0$. Thermal domain walls destroy long-range alignment in this one-dimensional short-range system. At zero temperature, the two fully aligned ground states are degenerate; that limiting case must be distinguished from any fixed positive temperature.

For a uniform classical spin $S\mathbf{n}$ with ${\mathbf{n}}^2=1$, the energy per bond isThe strict inequalities among the couplings make the $x$ direction the unique easy axis. ThusThese uniform configurations are also global classical minima: the largest singular value of the diagonal exchange matrix is $J^x$, so each bond has energy at least $-J^x{S}^2$, and both configurations saturate that bound.

A global rotation by $\pi$ about the $z$ axis sends $(S^x,S^y,S^z)$ to $(-S^x,-S^y,S^z)$, preserves every exchange product, and exchanges the two ground states. A $\pi$ rotation about $y$ does the same. The broken discrete symmetry can be represented by either operation.

The Holstein–Primakoff transformation replaces each finite-dimensional spin algebra by a bosonic creation operator and annihilation operator, while restricting the local occupation number to $0\le n_j\le 2S$. The number $n_j=a_j^{†}{a}_j$ counts deviations from the fully polarized reference quantum state. Its zero-boson Fock vacuum has ${\tilde{S}}^z=S$.

The representation is exact on this physical subspace: the ordered square root gives ${\tilde{S}}^{+}|n⟩=\sqrt{n(2S-n+1)}|n-1⟩$ and ${\tilde{S}}^{-}|n⟩=\sqrt{(n+1)(2S-n)}|n+1⟩$, reproducing the spin commutation relations. For $n_j\ll S$,so the leading Hamiltonian is quadratic in bosons and describes spin waves.

The square root generates boson interactions beyond leading order and its operator ordering matters. Truncation also enlarges the Hilbert space to unphysical occupations unless the low-density approximation remains self-consistent. Large fluctuations, particularly near a gap closing or in low dimensions, can invalidate the linear spin-wave approximation even when its algebraic manipulation is straightforward.

Put $\beta =(k_{B}T{)}^{-1}$ and label the two spin values in the order $+1,-1$. Assign each bond its Boltzmann weight,The symmetry of the bond energy makes $V$ a real symmetric matrix. For a fixed configuration the product of successive matrix entries is its Boltzmann weight. Summing over all spins and identifying the last index with the first therefore giveswhereThis replaces a sum over $2^L$ configurations by the powers of a two-by-two matrix. For finite bond energies at $T>0$, every matrix entry is positive, so the Perron–Frobenius theorem gives ${\lambda }_{+}>0$ and $|{\lambda }_{-}|<{\lambda }_{+}$. Hence the Helmholtz free energy per spin tends to $-k_{B}T\log {\lambda }_{+}$. A diagonal local observable $O$ can be inserted at one site: $⟨O_1⟩=\mathrm{tr}(O{V}^L)/Z_L$. Its limiting average is the expectation in the normalized eigenvector of ${\lambda }_{+}$.

The sum needs a stated interpretation. The original equation leaves the second bond outside explicit parentheses, with its index then formally free. The natural translation-invariant, equal-bond interpretation used here is the energy $-J{\sum }_j{\tau }_j({\sigma }_j+{\sigma }_{j+1})$.

Let $K=\beta J$. Sum over each decorating spin before summing over the remaining spins. The decimation of an alternating mixed-spin chain gives the effective transfer matrix:No on-site entropy has been omitted: the term $1$ is the weight of the zero decorating spin. The normalized eigenvectors are $v_{\pm }=(1,\pm 1{)}^{\mathsf{T}}/\sqrt{2}$, with eigenvaluesThere are $L$ transfer steps but $2L$ physical sites, soThe free energy per two-site unit cell is twice this value. At $J=0$, the formula gives $Z_L=6^L$, as required by two states on one sublattice and three on the other. As $T\to 0$, the free energy per physical site tends to $-J$, the energy of an aligned chain with two satisfied bonds per unit cell. These checks also fix the factor of two in the site count.

Insert $D=\mathrm{diag}(1,-1)$ at the first undecorated site. In the eigenbasis of the transfer matrix, $D$ exchanges $v_{+}$ and $v_{-}$, and both diagonal matrix elements vanish. Explicitly,Thus its large-$L$ limit is zero at every fixed $T>0$. The result is exact even before taking that limit. Equivalently, simultaneous reversal of all spins pairs configurations of equal energy and opposite ${\sigma }_1$. Positivity of the transfer matrix at fixed positive temperature prevents a degenerate dominant eigenvalue and spontaneous magnetization in this one-dimensional zero-field chain. Selecting one aligned state at strictly zero temperature is a different order of limits.

Insert the decorating-spin observable into the one bond that contains it. The decorated-spin occupation insertion uses the matrix:The partition function still contains $L$ decorating-spin sums; one of them has been replaced by this weighted sum. HenceTaking the dominant eigenvector givesThis is the probability that the decorating site is in either nonzero state. It tends to $2/3$ at high temperature, when its three states are equally likely, and to one at low temperature, when the nonzero aligned states gain bond energy. Unlike the magnetization, this even observable is not forced to vanish by spin reversal.

The Holstein–Primakoff transformation realizes a spin algebra using bosonic creation operators and bosonic annihilation operators. A boson counts a unit of spin deviation from the local positive quantization axis: ${\tilde{S}}^z=S-n$ with $n=a^{†}a$. The square root in ${\tilde{S}}^{+}$ supplies the correct spin-ladder matrix elements. In particular,Their commutator is $2(S-n)$ on every physical number state, reproducing the Spin commutation relations. Although the oscillator Fock space is infinite, the physical subspace is onlyThe order of the square root and oscillator operator matters in obtaining these coefficients.

When the spin deviations are dilute, expand the square root in $n/(2S)$. The leading transverse operators are linear in the bosonic creation operators and bosonic annihilation operators and turn exchange interactions into a quadratic Hamiltonian, which can be diagonalized. Higher powers describe interactions between spin waves. The method therefore organizes corrections by inverse powers of $S$ about an ordered classical configuration. Large $S$ alone is insufficient: the average depletion must remain small compared with $S$. Possible problems include the occupation constraint after truncation, noncollinear order requiring site-dependent axis rotations, zero modes, and infrared fluctuations that destroy the assumed dilute expansion.

The light states have the counting $3+4+1=8$. This suggests approximate chiral symmetry breakingwith a nonzero quark chiral condensate $⟨{\bar{q}}_R^j{q}_L^i⟩\propto {\delta }^{ij}$ for the three light flavours $u,d,s$. The three pions, four kaons and predominantly octet eta meson are the resulting pseudo-Goldstone bosons. Their small nonzero masses reflect explicit chiral symmetry breaking by the light-quark masses; for example, $m_{\pi }^2$ is proportional to $m_u+m_d$ at leading order. The larger strange-quark mass makes the kaons and octet eta meson heavier than the pions, while electromagnetic interactions help split charged and neutral partners.^[ref]

Although the massless strong-interaction terms admit independent rotations of all six flavours, the charm, bottom and top masses are too large to treat six-flavour chiral symmetry as a small explicit breaking at the light-hadron scale. Full approximate $SU(6{)}_L\times SU(6{)}_R$ breaking would suggest 35 light modes, rather than this octet. The classical axial $U(1)$ is also broken by the chiral anomaly, so it does not supply a ninth equally light Goldstone boson.

Terminology qualification: these mesons have spin zero but negative parity: they are pseudoscalar mesons, with $J^P=0^{-}$, rather than positive-parity scalar particles with $J^P=0^{+}$. Calling the listed states “scalar” is acceptable only if it is being used loosely to mean spin zero.

Let $\beta =1/(k_{B}T)$ and use the two spin values $+1,-1$ as row and column labels. Define the spin-chain transfer matrixIt is a real symmetric matrix because $E_0(s,t)=E_0(t,s)$; every entry is strictly positive for finite bond energies at $T>0$. The Boltzmann factor of a spin configuration factorizes into the product $W_{{\sigma }_1{\sigma }_2}\cdots W_{{\sigma }_L{\sigma }_1}$. Summing the intermediate spin labels is exactly matrix multiplication, and the final closure imposes periodic boundary conditions. ThusWriting $W=\left(\begin{array}{cc}a& b\\ b& d\end{array}\right)$, its eigenvalues areA positive-entry matrix need not have two positive eigenvalues; ${\lambda }_{-}$ can be negative. The Perron–Frobenius theorem ensures ${\lambda }_{+}>|{\lambda }_{-}|$ here. Therefore the Helmholtz free energy is $F_L=-k_{B}T\ln Z_L$, and its value per spin in the thermodynamic limit is $-k_{B}T\ln {\lambda }_{+}$.

For an open chain there are $L-1$ bonds, with both endpoint Ising spins independently summed. Retain the same spin-chain transfer matrix, but replace the trace by endpoint vectors. With $1=(1,1{)}^{\mathsf{T}}$,This follows directly by expanding the product of $L-1$ Boltzmann factors. If ${\mathbf{v}}_{\pm }$ are orthonormal eigenvectors of the symmetric matrix $W$, the equivalent spectral decomposition isFixed endpoint Ising spins replace the corresponding $1$ by a basis vector; extra endpoint energies replace its entries by the corresponding Boltzmann factors. The leading bulk Helmholtz free energy per spin agrees with the periodic result, because the boundary changes only subextensive terms. In particular, $L=1$ gives $Z_1^{\mathrm{open}}=1^{\mathsf{T}}1=2$, as required for one free Ising spin.

For low temperature, factor out the largest Boltzmann factor from the dominant eigenvalue:Thus the Helmholtz free energy per site has the expansionDifferentiate at fixed $J$, remembering that $K=J/(k_{B}T)$. The entropy per site isThe nonzero residual entropy comes from an exponentially large ground-state degeneracy, generated by geometrical frustration. A triangle with antiferromagnetic energy $J[st+(s+t)\tau ]$ has minimum energy $-J$: if $s=t$, the apex must have the opposite sign; if $s=-t$, either apex sign is permitted. The two all-equal configurations instead have energy $3J$. Thus exactly one of a minimum-energy triangle's bonds is frustrated. The bridge can simultaneously minimize its energy at $-J$ by requiring $t=-s^{\prime }$, where $s^{\prime }$ is the next cell's first base spin.

For fixed $s$ and $s^{\prime }$, the number of allowed apex choices is therefore two when $s=s^{\prime }$ and one when $s=-s^{\prime }$. The zero-temperature counting transfer matrix on successive first base spins isThe eigenvalues of $G$ are $3$ and $1$. Dividing $k_B\ln g_M$ by the $3M$ sites and taking the thermodynamic limit gives $k_B\ln 3/3$, independently confirming the derivative calculation. This does not assert a unique ordered ground state: the extensive choice of frustrated bonds is precisely what leaves entropy at zero temperature. At high temperature, all Ising spins become independent and $s\to k_B\ln 2$ per site.

Choose a coarse-grained magnetization $m(\mathbf{x})$ as the order parameter: it vanishes in the paramagnetism phase and develops a spontaneous nonzero value in the ferromagnetism phase. In zero applied magnetic field, spin reversal sends $m\mapsto -m$, so the local Landau free energy contains only even powers. A minimal spatially inhomogeneous free-energy functional isThe positive gradient coefficient penalizes rapid spatial variations. Near an ordinary transition take $a(T)=a_0(T-T_0)$, $a_0>0$, and positive $b$; retain the sextic term when $b$ can approach or cross zero. The source $h$ is conjugate to the magnetization and explicitly breaks spin reversal.

Determine equilibrium by minimizing this free-energy functional, imposing the appropriate boundary conditions. Its Euler-Lagrange equation isFor a uniform system at $h=0$, compare the values of all stable stationary solutions, rather than accepting every solution of this equation. The selected global minimum gives the equilibrium magnetization; derivatives of its free energy give entropy, heat capacity and magnetic susceptibility. At an ordinary transition the two ordered minima emerge continuously from zero. Spatial dependence also permits domain walls and an estimate of the correlation length.

This is a mean-field approximation: it selects a minimizing profile and neglects statistical fluctuations around it. The symmetry constrains the expansion, but its numerical coefficients are phenomenological. Close enough to a transition, thermal fluctuations near a critical point can invalidate the resulting mean-field critical exponents.^[ref]

For the sawtooth Ising chain, there is a source inconsistency: the printed Ising spin values are $\{-1,r+1\}$, with no definition of $r$. The requested transfer matrix follows for $r=0$, namely the usual Ising model with values $\{-1,+1\}$. For arbitrary $r$, summing an apex spin actually givesFor example, $r=1$, $s=s^{\prime }=-1$ and $\beta \bar{J}=1$ give the bracket $e^{-2}+e^4$, whereas the claimed hyperbolic cosine gives $e^{-2}+e^2$. Thus the printed identity is false literally. The rest of Question 1 uses the intended symmetric Ising spin values, $r=0$.

For a fixed configuration of base spins, each apex spin appears in only one local Boltzmann factor. The apex spins are therefore conditionally independent, and their finite sums factorize:Only the base spins carry information from one triangle to the next. The apex spins have already been summed out, so no apex index remains in the transfer matrix; the final trace implements the periodic boundary condition.

In the $s=+1,-1$ basis the transfer matrix and its eigenvalues areFor finite positive temperature, $a,b>0$ and ${\lambda }_{+}>|{\lambda }_{-}|$. Consequently $Z={\lambda }_{+}^L+{\lambda }_{-}^L$ and $L^{-1}\log Z\to \log {\lambda }_{+}$ in the thermodynamic limit. There are $2L$ spins, so the Helmholtz free energy per spin isThe factor $1/2$ is essential: one transfer matrix step counts two spins. For completeness, if the undefined $r$ is retained, set $a_r=A_{-1,-1}^{(r)}$, $d_r=A_{r+1,r+1}^{(r)}$ and $b_r=A_{-1,r+1}^{(r)}$. The same thermodynamic limit gives $f_r=-(2\beta {)}^{-1}\log [(a_r+d_r+\sqrt{(a_r-d_r{)}^2+4b_r^2})/2]$ whenever these are two distinct spin values.

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.