We demonstrate the existence, in mean field theory, of an antiferromagnetic transition in the transverse-field Hubbard Model and present the phase diagram.
|Fig. 1: Our gauge choice for the problem. Unit cell shown shaded in blue.|
The Hubbard model is a stylized description of fermions interacting via short ranged interactions. The lattice spacing (set equal to 1 in this report) provides a UV regulator, while the IR effects are regulated by replacing the Coulomb interaction by a delta function in space. Mathematically, the Hamiltonian is given by
We use the t elements described in Fig. 1. Lines with arrow represent n.n. hopping elements of strength t and phase &pi/4. Single lines represent diagonal hops, with strength t', and phase 0 (dashed) or &pi (solid). This seemingly arbitrary choice corresponds physically to introducing a unit of magnetic flux through each plaquette and exploiting gauge freedom to reduce the unit cell to two sites.
In the mean field approximation, we replace the quartic interaction by a quadratic one via
|Fig. 2: The energy dispersion for the transverse-field, U=0 Hubbard model, as a function of x and y momenta. The Brillouin zone can be chosen to be within the black diamond.|
Anticipating antiferromagnetism, we guess the mean field values
where &delta is the antiferromagnetic order parameter. With this approximation in place, computing the unit cell Hamiltonian is straightforward:
|Fig. 3:The phase diagram for our model, in mean field approximation. The code used to produce the coordinates of the boundary is here.|
This Hamiltonian admits the energy eigenvalues
The energy spectrum for U=0 is plotted over the Brillouin zone in Fig.1.
With the mean field hamiltonian in hand, we are in a position to determine the antiferromagnetic order parameter self consistently. To do so, we impose the constraint given by the usual "gap equation" which takes the form
This constraint gives &Delta(t,t',U), and thus &delta(t,t',U). Since the integral is convergent for any finite value of &Delta, there are values of U for which there is no self-consistent solution. Physically, the mean field version of our system exhibits a phase transition from nonmagnetic (no solution) to antiferromagnetic (solution exists) at a finite value of U. By evaluating the gap equation integral numerically, we can reconstruct the phase diagram for our model, as shown in Fig. 3. In particular, we find that antiferromagnetism does not emerge for values of U less than
See the next report for a discussion of the critical properties and emergence of a Goldstone mode near the phase boundary.
© 2007 D. Stanford. The author grants permission to copy, distribute and display this work in unaltered form, with attribution to the author, for noncommercial purposes only. All other rights, including commercial rights, are reserved to the author.