Casey Chu
- Statistical Mechanics
- I took these notes in winter 2017 for Physics 212, taught by Professor Stephen Shenker.
- Ising model
- The Ising model is a lattice of spins with Hamiltonian where indexes over neighboring sites and . Thus the partition function is where we’ve nondimensionalized by setting and . In the following, let be the number of bonds, and be the number of sites.
- In the high-/low-temperature regime, the system is ordered, dominated by states in which most spins are pointing in the same direction. In the low-/high-temperature regime, the system is disordered, and many states are equally important.
- Useful diagostics of order are magnetization in the symmetry-breaking limit and the spin-spin correlation function .
- High-temperature expansion ()
- In the following, set . Note that for , the following identity holds: We can use this to rewrite the partition function as The product expands into a power series in . Noticing that we see that terms corresponding to open paths on the lattice vanish: but terms corresponding to closed paths on the lattice survive: More abstractly, where the sum indices over closed (but perhaps disconnected) paths in the lattice.
- For the 2D lattice, for example, we obtain
- The magnetization is clearly , and the spin-spin correlation function can be seen as times the above expression for the partition function that counts closed paths, except that the sites and are “pre-covered.” Thus this counts paths from to (as well as disconnected closed paths, which tend to cancel). The first contribution is the path of length directly from to , so .
- Low-temperature expansion ()
- In the low-temperature regime, we can explicitly list the most important Boltzmann factors. For the 2D lattice, the all-up/all-down states dominate, followed by the states in which only spin differs from the rest, followed by the states in which two adjacent spins differ from the rest, and so on, so that The is the Boltzmann weight for the ordered all up/all down state, and the is the cost of breaking bonds. For example, in the state in which only one spin differs from the rest, incident bonds are broken.
- In fact, the broken bonds surrounding an island of flipped spins can be thought of as, again, closed paths, but this time in the dual lattice. Abstractly, we write This is Kramers-Wannier duality. We can define an involution with the relation that identifies a lattice at with its dual lattice at . More symmetrically, this relation is . Since the 2D lattice is its own dual, if we assume that there is one critical point, we can solve for by setting .
- On the 2D lattice, we can calculate the magnetization directly from the expression for the partition function: where the factor of is for symmetry-breaking. For the spin-spin correlation function, the first two factors that don’t cancel correspond to the totally ordered state and the one where spins are flipped, giving .
- Mean field theory/effective field theory
- Now coarse-grain by coalescing the individual microstates into macrostates determined by a magnetization field , roughly the average spin in a block of spins at a location . An effective Hamiltonian is given by and the partition function can now be written as the functional integral
- The Landau-Ginzburg Hamiltonian is defined by where:
- the integral over space arises because we’d like to express as a sum of local effects,
- the expansion in powers of is justified because non-analyticity is washed out by the coarse-graining process (think central limit theorem),
- the gradient terms arise to capture local interactions,
- odd powers of vanish via the symmetry when ,
- mixed derivative terms vanish via rotational symmetry,
- the highest power of has positive coefficent (in this case, ) to ensure zero probability for large .
- The parameters , , , etc. are phenomenological and depend on the specific model and properties like temperature and pressure.
- We make the mean-field approximation and examine only the most probable configuration, that is, the one that minimizes . When is minimized, the field is constant (), and the most probable configuration can be viewed as the minimum of the function Assume for now. For , the minimum is only at , but for , the minimum goes as . If is an increasing differentiable function of the temperature , we conclude as well that for and for , at least locally around .
- Thus we see that a phase transition occurs at . Further (higher-order) phase transitions occur when, for example, both and , and can be investigated by including more terms in the power series. The mean-field approximation is extremely accurate for the Ising model in dimensions .
- For , we can consider small fluctuations around the minimum of . That is, we write , so that where , and the inner product is the inner product. This is an example of a Gaussian integral. We’re ignoring the divergent factor of . On a lattice, the Laplacian becomes the discrete Laplacian.
- Viewed as a function of , the free energy is the moment-generating function of . In particular, the correlation function is in , where is the solution to the PDE , found using, for example, the Fourier transform. (In momentum space with momenta , is diagonal with eigenvalues , allowing us to easily invert . On a square lattice, the eigenvalues of become instead of , erasing rotational symmetry at high momenta/small scales.) At the critical point, the correlation scales as a power law.
- Renormalization group