[Image goes here]
Published on

Statistical Physics of Particles/Fields

SPP: Statistical Physics of Particles

SPF: Statistical Physics of Fields

Chapter 1

What actually is a field??? Kardar defines an "average deformation field" u(x)u(x) as an alternative to the interacting particle model of phonons, with a corresponding velocity field u/t\partial u/\partial t.

  • Locality gives us short range interactions between particles, such that we can define a potential energy density Φ\Phi for every xx, where V[u]=dxΦ(u(x),u/x,)V[u] = \int dx \Phi(u(x), \partial u/\partial x, \cdots)
  • Translational Symmetry For a 1-d chain, uniform translation does not change the energy, so the energy density satisfies a periodic Φ(u(x)+c)=Φ(u(x))\Phi(u(x) + c) = \Phi(u(x)). This removes the dependence of Φ\Phi on uu itself and only on its derivatives.
  • Stability Looking at equilibrium, we don't want linear terms in uu (linearizing everything should mean all particles are stable). Then, quadratic term of VV should be pos-definite.

What potential does this give us? The most general one is of the form

V(u)=dx[K2(ux)2+L2(2ux2)2+]V(u) = \int dx[\frac{K}{2} (\frac{\partial u}{\partial x})^2 + \frac{L}{2}(\frac{\partial^2 u}{\partial x^2})^2 + \cdots]

Update: A better answer may be in Zee

Problems

2. (SPF) Ising model of magnetism: Hamiltonian for configuration {σi}\{\sigma_i\} of spins is

H=12i,j=1NJijσiσjhiσi\begin{aligned} \mathcal{H} = \frac{1}{2} \sum_{i, j = 1}^N J_{ij} \sigma_i \sigma_j - h\sum_i \sigma_i \end{aligned}
  • Suppose Jij=J/NJ_{ij} = -J/N, so H=J2Ni,j=1Nσiσjhiσi\mathcal{H} = -\frac{J}{2N} \sum_{i, j = 1}^N \sigma_i \sigma_j - h\sum_i \sigma_i. Take the magnetization to be m=i=1Nσi/N=M/Nm = \sum_{i = 1}^N \sigma_i/N = M/N (average spin direction), so
E(M,h)=H=J2Ni,j=1Nσiσjhiσi=N(Jm2/2+hm)\begin{aligned} E(M, h) = \langle \mathcal{H} \rangle = \langle -\frac{J}{2N} \sum_{i, j = 1}^N \sigma_i \sigma_j - h\sum_i \sigma_i\rangle = -N(Jm^2/2 + hm) \end{aligned}
  • The partition function is Z(h,T)=(σi)exp(βH)Z(h, T) = \sum_{(\sigma_i)} \exp(-\beta \mathcal{H}), where our sum is over all spin configurations (σi)(\sigma_i). However, our Hamiltonian is invariant to specific spin values and only depends on the total net spin, so we can alter our sum to be over the number of positive spins NpN_p, which also fixes a magnetization MM such that H=E(M,h)\mathcal{H} = E(M, h):
Z(h,T)=Np=0N( number of configs with Np up spins )×exp(βE(M,h))=Np=0N(NNp)exp(βE(M,h))=Np=0Nexp(β(1/βln(NNp)))exp(βE(M,h))=Np=0Nexp(β(E(M,h)1/βln(NNp)))\begin{aligned} Z(h, T) & = \sum_{N_p = 0}^N (\text{ number of configs with Np up spins }) \times \exp(-\beta E(M, h)) \\ & = \sum_{N_p = 0}^N {N \choose N_p} \exp(-\beta E(M, h)) \\ & = \sum_{N_p = 0}^N \exp(-\beta(-1/\beta \ln{N \choose N_p}))\exp(-\beta E(M, h)) \\ & = \sum_{N_p = 0}^N \exp(-\beta(E(M, h) -1/\beta \ln{N \choose N_p})) \end{aligned}

But having NpN_p up spins forces NNpN - N_p down spins, so we might as well just sum over M=2NpNM = 2N_p - N instead, which gives us

Z(h,T)=Mexp(β(E(M,h)kBTln(N(M+N)/2)))=Mexp(βF(M,h))F(m,h)=E(M,h)kBTln(N(M+N)/2)\begin{aligned} Z(h, T) & = \sum_{M} \exp(-\beta(E(M, h) - k_B T \ln{N \choose (M+N)/2})) = \sum_{M} \exp(-\beta F(M, h)) \\ & F(m, h) = E(M, h) - k_B T \ln{N \choose (M+N)/2}\\ \end{aligned}
  • Free energy F(h,T)=kBTlnZ(h,T)=minm(F(m,h))F(h, T) = -k_B T \ln Z(h, T) = \min_m(F(m, h))