In this work we study a system of interacting fermions on a triangular lattice in the presence of an external magnetic field. We neglect spin and fix a density of one third, with one unit of magnetic flux per particle. The infinite density matrix renormalization group algorithm is used to compute the ground state of this generalized Fermi-Hubbard model. Increasing the strength of the nearest-neighbor repulsion, we find a first order transition between an Integer Quantum Hall phase and a crystalline, generalized Wigner crystal state. The first-order nature of the phase transition is consistent with a Ginzburg-Landau argument. We expect our results to be relevant for moiré heterostructures of two-dimensional materials.
We consider the crystallization of a two-dimensional electron system in a perpendicular magnetic field using composite boson theory. There are three possible states to consider: the Hall liquid, the Wigner crystal, and the Hall crystal (a state with both broken translation symmetry and a quantized Hall response). Within composite boson theory, these states map onto a superconductor, a Mott insulator, and a supersolid of composite bosons respectively. We show that when a
ν=1 Hall liquid has a sufficiently soft roton, there is a first order transition to a triangular lattice Hall crystal. If we continue to decrease the roton mass, there is a continuous transition from the Hall crystal to a Wigner crystal. {When the Hall crystal exhibits the integer quantum Hall effect,} this transition {is} described by a free Dirac fermion and, at the critical point, the coupling to the phonons of the crystal is irrelevant, {in the {renormalization group} sense}. We extend this analysis to fractional
ν=1/m Hall liquids. There, due to kinetic frustration arising from flux attachment, honeycomb lattice Hall crystals are preferred over triangular ones at intermediate interaction strength.
Two-dimensional density-matrix renormalization group method is employed to examine the ground-state phase diagram of the Hubbard model on the triangular lattice at half-filling. The calculation reveals two discontinuities in the double occupancy with increasing the repulsive Hubbard interaction
U at ${U}_{\mathrm{c}1}\ensuremath{\sim}7.8t$ and ${U}_{\mathrm{c}2}\ensuremath{\sim}9.9t$ (
t being the hopping integral), indicating that there are three phases separated by first-order transitions. The absence of any singularity in physical quantities for $0\ensuremath{\le}U<{U}_{\mathrm{c}1}$ implies a metallic phase in this regime. For
U>Uc2, the local spin density induced by an applied pinning magnetic field exhibits a three sublattice feature, which is compatible with the ${120}^{\ensuremath{\circ}}$ N\'eel-ordered state realized in the limit of $U\ensuremath{\rightarrow}\ensuremath{\infty}$. For
Uc1<U<Uc2, a response to the applied pinning magnetic field is comparable to that in the metallic phase with a relatively large spin correlation length, but showing neither valence bond nor chiral magnetic order, which therefore resembles gapless spin liquid. However, the spin structure factor for the intermediate phase exhibits the maximum at the
K and ${\mathrm{K}}^{\ensuremath{'}}$ points in the momentum space, which is not compatible to spin liquid with a large spinon Fermi surface. The calculation also finds that the pairing correlation function monotonically decreases with increasing
U and thus the superconductivity is unlikely in the intermediate phase.
Motivated by experimental studies that have found signatures of a quantum spin liquid phase in organic crystals whose structure is well described by the two-dimensional triangular lattice, we study the Hubbard model on this lattice at half filling using the infinite-system density matrix renormalization group (iDMRG) method. On infinite cylinders with finite circumference, we identify an intermediate phase between observed metallic behavior at low interaction strength and Mott insulating spin-ordered behavior at strong interactions. Chiral ordering from spontaneous breaking of time-reversal symmetry, a fractionally quantized spin Hall response, and characteristic level statistics in the entanglement spectrum in the intermediate phase provide strong evidence for the existence of a chiral spin liquid in the full two-dimensional limit of the model.
The generalized Wigner crystal (GWC) is a novel quantum phase of matter driven by further-range interaction at fractional fillings of a lattice. The role of further-range interaction as the driver for the incompressible state is akin to the Wigner crystal. On the other hand, the significant role of commensurate filling is akin to the Mott insulator. Recent progress in simulator platforms presents unprecedented opportunities to investigate quantum melting in the strongly interacting regime through synergy between theory and experiments. However, the earlier theory literature presents diverging predictions. We study the quantum freezing of GWC through large-scale density matrix renormalization group simulations of a triangular lattice extended Hubbard model. We find a single first-order phase transition between the Fermi liquid and the sqrt[3]×sqrt[3] GWC state. The GWC state shows long-range antiferromagnetic 120° Néel order. Our results present the simplest answers to the question of the quantum phase transition into the GWC phase and the properties of the GWC phase.
The magnetic phase diagram of the quarter-filled generalized Wigner lattice with nearest-neighbor and next-nearest-neighbor hoppings,
t1 and
t2, is explored. We find a region at negative
t2 with fully saturated ferromagnetic ground states that we attribute to kinetic exchange. Such interaction disfavors antiferromagnetism at
t2<0 and stems from virtual excitations across the charge gap of the Wigner lattice, which is much smaller than the Mott-Hubbard gap $\ensuremath{\propto}U$. Remarkably, we find a strong dependence of the charge structure factor on magnetism even in the limit $U\ensuremath{\rightarrow}\ensuremath{\infty}$, in contrast to the expectation that charge ordering in the Wigner lattice regime should be well described by spinless fermions. Our results, obtained using the density-matrix renormalization group and exact diagonalization, can be transparently explained by means of an effective low-energy Hamiltonian.