Physlib

PhyslibAlpha.CondensedMatter.IsingModel.GibbsStates

Equilibrium states of the Ising model

The Gibbs states of the classical and of the quantum Ising model in a finite region, and the states of the infinite lattice that arise from them in the thermodynamic limit.

i. Overview

A finite region `Y` of an Ising system in thermal equilibrium at inverse temperature `β` is in the *Gibbs state*, in which the configuration `c` of its spins has probability `exp (- β E_Y(c)) / Z_Y`. The expectation value of an observable depends on the model:

* in the classical Ising model an observable is a function `f` of the configuration and its expectation value is `∑ c, f c exp (- β E_Y(c)) / Z_Y`; * in the quantum Ising model an observable is a self-adjoint operator `M` on the spins, the Hamiltonian is diagonal in the `σᶻ` basis with the eigenvalue `E_Y(c)` on `c`, and the expectation value is `tr (exp (- β H_Y) M) / Z_Y = ∑ c, M (c, c) exp (- β E_Y(c)) / Z_Y`.

For a function `f` on configurations, which is a diagonal operator, the two expectation values agree: the classical model is the `σᶻ` sector of the quantum model.

The finite-volume Gibbs states are generally not compatible under restriction: boundary interactions change when the region grows. They nevertheless define a *thermodynamic limit* along a chosen ultrafilter containing the tails of the finite regions. Compactness gives limits of local expectations, and those limits are compatible, hence define a state of the net. This construction proves existence of a limit point. It does not prove uniqueness, convergence along all finite regions, translation invariance, or the occurrence of a phase transition.

ii. Key results

- `IsingModel.gibbsClassical` : the Gibbs state of the classical Ising model in a region. - `IsingModel.gibbs` : the Gibbs state of the quantum Ising model in a region. - `IsingModel.thermodynamicFamily` : the Gibbs states of the finite regions, exhausting the lattice. - `IsingModel.thermodynamicLimitClassical` : a limit of classical Gibbs states. - `IsingModel.tendsto_thermodynamicLimitClassical` : convergence of local classical expectations along the chosen ultrafilter. - `IsingModel.thermodynamicLimit` : the corresponding construction for quantum spins.

iii. Table of contents

  • A. Gibbs states of the quantum Ising model
  • B. Thermodynamic limits

iv. References

* O. Bratteli & D. Robinson, Operator Algebras and Quantum Statistical Mechanics 2, 2nd ed., Section 6.2. * L. D. Landau & E. M. Lifshitz, Statistical Physics, Part 1, Chapter 14. [ref: landau_statphys1]

A. Gibbs states of the quantum Ising model

B. Thermodynamic limits

5 declarations

definition

Gibbs state of the quantum Ising model in region YY

For a quantum Ising model at inverse temperature β\beta, the Gibbs state in a finite region Y⊂SY \subset S is a unital positive linear map on the algebra of observables AY\mathcal{A}_Y (self-adjoint operators on the spin space). For an observable M∈AYM \in \mathcal{A}_Y, its expectation value is given by the weighted sum over all spin configurations c∈Cfg(Y)c \in \text{Cfg}(Y): ωY,β(M)=∑c∈Cfg(Y)Pβ(c)Re(⟨c∣M∣c⟩) \omega_{Y, \beta}(M) = \sum_{c \in \text{Cfg}(Y)} P_\beta(c) \text{Re}(\langle c | M | c \rangle) where Pβ(c)P_\beta(c) is the Gibbs probability of configuration cc defined by Pβ(c)=e−βEY(c)Zβ(Y)P_\beta(c) = \frac{e^{-\beta E_Y(c)}}{Z_\beta(Y)}, and ⟨c∣M∣c⟩\langle c | M | c \rangle denotes the diagonal element of the operator MM in the σz\sigma^z basis. This corresponds to the quantum mechanical trace Tr(ρM)\text{Tr}(\rho M) where the density matrix ρ=e−βHYZY\rho = \frac{e^{-\beta H_Y}}{Z_Y} is diagonal in the configuration basis.

theorem

The expectation value of the quantum Gibbs state ωY,β(M)\omega_{Y, \beta}(M) equals ∑cPβ(c)Re(⟨c∣M∣c⟩)\sum_{c} P_\beta(c) \text{Re}(\langle c | M | c \rangle)

For a finite region YY and a quantum observable MM (a self-adjoint operator on the spin space of YY), the expectation value of MM in the Gibbs state at inverse temperature β\beta is equal to the sum over all spin configurations c∈Cfg(Y)c \in \text{Cfg}(Y) of the product of the Gibbs probability Pβ(c)P_\beta(c) and the real part of the diagonal matrix element of MM corresponding to configuration cc: ωY,β(M)=∑c∈Cfg(Y)Pβ(c)Re(⟨c∣M∣c⟩) \omega_{Y, \beta}(M) = \sum_{c \in \text{Cfg}(Y)} P_\beta(c) \text{Re}(\langle c | M | c \rangle)

definition

Thermodynamic family of finite-volume Gibbs states ωY,β\omega_{Y, \beta}

For an Ising model II at inverse temperature β\beta on a lattice SS, the thermodynamic family is a collection of Gibbs states (ωY,β)Y∈F(S)(\omega_{Y, \beta})_{Y \in \mathcal{F}(S)} indexed by the set of finite regions F(S)\mathcal{F}(S). For each finite region YY, the associated state is the Gibbs state ωY,β\omega_{Y, \beta} of the quantum Ising model. This family is equipped with an ultrafilter U\mathcal{U} on F(S)\mathcal{F}(S) that refines the filter of finite sets directed by inclusion (the "at top" filter), which is used to define the thermodynamic limit of local expectations as the regions exhaust the lattice.

definition

Thermodynamic limit state of the quantum Ising model

For a quantum Ising model II on a lattice SS at inverse temperature β\beta, the thermodynamic limit ωβ\omega_\beta is a state of the infinite lattice (represented as a `NetState`). This state is constructed by taking the limit of the expectations of the thermodynamic family of finite-volume quantum Gibbs states ωY,β\omega_{Y, \beta} along a chosen ultrafilter U\mathcal{U} on the finite subsets Y⊆SY \subseteq S. For any local observable MM, the expectation value in the thermodynamic limit state is given by ωβ(M)=lim⁡Y→UωY,β(M) \omega_\beta(M) = \lim_{Y \to \mathcal{U}} \omega_{Y, \beta}(M) where ωY,β\omega_{Y, \beta} is the Gibbs state for the finite region YY.

theorem

Expectation Values Converge to the Thermodynamic Limit State: lim⁡Y→UωY,β(M)=ωβ(M)\lim_{Y \to \mathcal{U}} \omega_{Y, \beta}(M) = \omega_\beta(M)

For a quantum Ising model on a lattice SS at inverse temperature β\beta, let X⊆SX \subseteq S be a finite region and MM be a local observable associated with XX. Let ωY,β\omega_{Y, \beta} denote the Gibbs state of the system in a finite region Y⊆SY \subseteq S, and let ωβ\omega_\beta denote the thermodynamic limit state constructed via an ultrafilter U\mathcal{U} on the finite subsets of SS. The expectation value of the local observable MM in the finite-volume Gibbs states converges along U\mathcal{U} to the expectation value of MM in the thermodynamic limit state: lim⁡Y→UωY,β(M)=ωβ(M)\lim_{Y \to \mathcal{U}} \omega_{Y, \beta}(M) = \omega_\beta(M) where ωY,β(M)\omega_{Y, \beta}(M) is the expectation value of MM (embedded into the algebra of region YY) when X⊆YX \subseteq Y.