N spins s1,s2,…,sN, each =±1 on a regular grid; only neighbors interact. 1D = chain, 2D = grid.

Energy

E(s1,s2,…)=−h∑i=1Nsi−J∑i=1N−1sisi+1

Studied in the canonical ensemble. Define h^=h/kT, J^=J/kT. Boltzmann factor:

e−E/kT=eh^∑si+J^∑sisi+1

Partition function:

Z(h^,J^)=∑s1=±1,s2=±1,…eh^∑si+J^∑sisi+1

Single spin (J=0, N=1)

Z0(h^)=∑s1=±1e−h^s1=e−h^+eh^=cosh⁡(h^)

Average spin:

⟨s1⟩=1Z0(h^)∑s1=±1s1e−h^s1=12cosh⁡h^(−e−h^+eh^)=sinh⁡h^cosh⁡h^=tanh⁡(h^)

Response saturates at ⟨s⟩=±1 for large |h^|.


Non-interacting spins (J=0, N spins)

Boltzmann factor factorizes → partition function factorizes:

Z(h^,N)=Z0(h^)N=cosh⁡(h^)N

Magnetization density m=1N∑isi, and by the same argument per spin:

⟨m⟩=tanh⁡h^

Strong coupling limit (J→∞)

J^→∞ locks spins into alignment; anti-aligned configurations have Boltzmann factor →0. System behaves as a single dipole of strength N:

⟨m⟩J→∞=tanh⁡(Nh^)

As N→∞: sharp step — any nonzero h fully magnetizes the system. (J→∞ equivalent to T→0 at fixed J.)


Define li=sisi+1=±1 (aligned vs. anti-aligned). Energy becomes:

E=−J∑i=1N−1sisi+1=−J∑i=1N−1li

The N−1 link variables are independent — this is just a non-interacting Ising model for the links with external field J. This is a duality: the interacting spin model maps exactly onto a non-interacting link model.


General result (N→∞, transfer matrix)

⟨m⟩=sinh⁡h^sinh2⁡h^+e−4J^

Recovers tanh⁡h^ at J^=0 and sharp step as J^→∞. No finite-T phase transition in 1D for any J>0. The 2D Ising model does have a finite critical Tc>0.