We describe here a quantum simulator of extended bipartite Hubbard model with broken sublattice symmetry. The simulator consists of a structured lateral gate confining two-dimensional electrons in a quantum well into artificial minima arranged in a hexagonal lattice. The sublattice symmetry breaking is generated by forming an artificial triangular graphene quantum dot (ATGQD) with zigzag edges. The resulting extended Hubbard model generates tunable ratio of tunneling strength to electron-electron interactions and of sublattice symmetry with control over shape. The validity of the simulator is confirmed for small systems using mean-field and exact diagonalization many-body approaches which show that the ground state changes from a metallic to an antiferromagnetic (AF) phase by varying the distance between sites or depth of the confining potential. The one-electron spectrum of these triangular dots contains a macroscopically degenerate shell at the Fermi level. The shell persists at the mean-field level for weak interactions (metallic phase) but disappears for strong interactions in the AF phase. We determine the effects of electron-electron interactions on the ground state, the total spin, and the excitation spectrum as a function of filling of the ATGQD. We find that the half-filled charge neutral shell leads to a partially spin-polarized state in both metallic and AF regimes in accordance with Lieb's theorem. In both regimes a relatively large gap separates the spin-polarized ground state to the first excited many-body state at half-filling of the degenerate shell. By adding or removing an electron, this gap drops dramatically, and alternate total spin states emerge with energies nearly degenerate to a spin-polarized ground state.
We describe here a quantum simulator of extended bipartite Hubbard model with broken sublattice symmetry. The simulator consists of a structured lateral gate confining two-dimensional electrons in a quantum well into artificial minima arranged in a hexagonal lattice. The sublattice symmetry breaking is generated by forming an artificial triangular graphene quantum dot (ATGQD) with zigzag edges. The resulting extended Hubbard model generates tunable ratio of tunneling strength to electron-electron interactions and of sublattice symmetry with control over shape. The validity of the simulator is confirmed for small systems using mean-field and exact diagonalization many-body approaches which show that the ground state changes from a metallic to an antiferromagnetic (AF) phase by varying the distance between sites or depth of the confining potential. The one-electron spectrum of these triangular dots contains a macroscopically degenerate shell at the Fermi level. The shell persists at the mean-field level for weak interactions (metallic phase) but disappears for strong interactions in the AF phase. We determine the effects of electron-electron interactions on the ground state, the total spin, and the excitation spectrum as a function of filling of the ATGQD. We find that the half-filled charge neutral shell leads to a partially spin-polarized state in both metallic and AF regimes in accordance with Lieb's theorem. In both regimes a relatively large gap separates the spin-polarized ground state to the first excited many-body state at half-filling of the degenerate shell. By adding or removing an electron, this gap drops dramatically, and alternate total spin states emerge with energies nearly degenerate to a spin-polarized ground state.
Semiconductor artificial graphene nanostructures where the Hubbard model parameter U/t can be of the order of 100, provide a highly controllable platform to study strongly correlated quantum many-particle phases. We use accurate variational and diffusion Monte Carlo methods to demonstrate a transition from antiferromagnetic to metallic phases for an experimentally accessible lattice constant a=50 nm in terms of lattice site radius $\ensuremath{\rho}$, for finite-sized artificial honeycomb structures nanopatterned on GaAs quantum wells containing up to 114 electrons. By analyzing spin-spin correlation functions for hexagonal flakes with armchair edges and triangular flakes with zigzag edges, we show that edge type, geometry, and charge nonuniformity affect the steepness and the crossover $\ensuremath{\rho}$ value of the phase transition. For triangular structures, the metal-insulator transition is accompanied with a smoother edge polarization transition.
Quantum Hall systems have a one-body energy spectrum consisting of dispersion-less Landau levels. Electron-electron interactions thus dominate in partially filled Landau levels, which exhibit a myriad of strongly correlated phases such as quantum hall ferromagnets and fractional quantum Hall phases. We study two examples of these phenomena. In the first project, we explore the ground state of a system with an interface between two semi-infinite regions with fillings ν= 4 and ν= 3 respectively. The width of the interface can be controlled by varying the background potential, which provides an additional tuning parameter. For a certain range of interaction strengths, the ν= 4 bulk is unpolarized whereas the ν= 3 bulk is fully polarized. In the parameter space spanned by the interaction strength and width of the interface, we find two phases at the interface. Phase A has spin as a good quantum number, and the long-wavelength spin edge excitations are gapped. In phase B, spin rotation symmetry is spontaneously broken at the mean-field level. Using symmetry arguments we find the effective theory near the interface of phase B. This effective theory is known to have gapless long-wavelength spin excitations. In the second project, we study the ground state of a tunnel-decoupled double-layer graphene system when both layers are undoped. We find a simple Hamiltonian in the continuum limit from symmetries of the system. Using the Hartree-Fock approximation we find a state with inter-layer coherence with broken layer U(1) symmetry. This phase becomes magnetized in presence of a non-zero Zeeman field. A first-order phase transition can be driven from the ferromagnetic phase to the magnetized inter-layer coherent phase by increasing the Zeeman field. We predict the number of gapless modes in the bulk.
We investigated the honeycomb lattice using a two-orbital (α and β) Hubbard model and many-body exact calculations. The spin gap (Δ s ), charge gap (Δ c ), intrasite ( γ i i α β = ⟨ S → i α · S → i β ⟩ ), intersite ( γ i j α β = ⟨ S → i α · S → j β ⟩ ), spin–spin correlations, and orbital magnetic moments μ iα are analyzed for different band-filling ( n ), which denotes the number of electrons per site. Special attention was set in one-quarter-filled ( n = 1) and half-filled ( n = 2) bands. The results are presented by varying Coulomb repulsion ( U), Hund’s coupling ( J H ), external magnetic field ( B ), and orbital energy splitting (Δϵ). For the one-quarter-filled band, an intrasite magnetic disorder ( γ i i α β ≈ 0 ) begins to form for large values of U / t, which is interpreted as a possible signature of a resonating valence bond state. Full orbital polarization ( γ i i α β = 0.25 ) with intersite antiferromagnetic (AFM) order ( γ i j α β = − 0.25 ) was observed in the half-filled band for U / t ≫ 1, indicating that the two-orbital Hubbard model can be mapped to an AFM Heisenberg model with S = 1. The increase of J H is conductive to insulating and metallic phases for half- and one-quarter-filled bands, respectively. For a fractional n, the honeycomb lattice can exhibit spin-glass-like behavior, predominating as a disordered magnetic system characterized by random values of γ i j α β . The magnetic field induced a reduction in Δ c, which became zero at the critical external magnetic field. Twisted bilayer graphene (TBG) exhibits approximately flat electronic bands at magic angles and interesting correlated phenomena, such as unconventional superconductivity and correlated insulating states around n = 1 and n = 2, respectively. The dominance of electron–electron interactions in TBG makes the Hubbard model, in various versions (extended or multi-orbital), offer a practical, structure-scale account for investigating many-body effects in large moiré unit cells, where first-principles density-functional theory calculations are computationally challenging.
The noninteracting electronic structures of tight-binding models on bipartite lattices with unequal numbers of sites in the two sublattices have a number of unique features, including the presence of spatially localized eigenstates and flat bands. When a uniform on-site Hubbard interaction U is turned on, Lieb proved rigorously that at half-filling ($\ensuremath{\rho}=1$) the ground state has a nonzero spin. In this paper we consider a ``CuO2 lattice'' (also known as ``Lieb lattice,'' or as a decorated square lattice), in which ``d orbitals'' occupy the vertices of the squares, while ``p orbitals'' lie halfway between two d orbitals; both d and p orbitals can accommodate only up to two electrons. We use exact determinant quantum Monte Carlo (DQMC) simulations to quantify the nature of magnetic order through the behavior of correlation functions and sublattice magnetizations in the different orbitals as a function of U and temperature; we have also calculated the projected density of states, and the compressibility. We study both the homogeneous (H) case, Ud=Up, originally considered by Lieb, and the inhomogeneous (IH) case, ${U}_{d}\ensuremath{\ne}{U}_{p}$. For the H case at half-filling, we found that the global magnetization rises sharply at weak coupling, and then stabilizes towards the strong-coupling (Heisenberg) value, as a result of the interplay between the ferromagnetism of like sites and the antiferromagnetism between unlike sites; we verified that the system is an insulator for all U. For the IH system at half-filling, we argue that the case ${U}_{p}\ensuremath{\ne}{U}_{d}$ falls under Lieb's theorem, provided they are positive definite, so we used DQMC to probe the cases Up=0,Ud=U and Up=U,Ud=0. We found that the different environments of d and p sites lead to a ferromagnetic insulator when Ud=0; by contrast, Up=0 leads to to a metal without any magnetic ordering. In addition, we have also established that at density $\ensuremath{\rho}=1/3$, strong antiferromagnetic correlations set in, caused by the presence of one fermion on each d site.
We introduce a local formalism to deal with the Hubbard model on an N × N square lattice (for even N ) in terms of eigenstates of number operators, having well defined point symmetry. For U →0, the low-lying shells of the kinetic energy are filled in the ground state. At half-filling, using the 2 N -2 one-body states of the partially occupied shell hf , we build a set of degenerate unperturbed ground states with S z = 0 which are then resolved by the Hubbard interaction Ŵ = U ∑ r r ↑ r ↓ . We study the many-body eigenstates in hf of the kinetic energy with vanishing eigenvalue of the Hubbard repulsion ( W = 0 states). In the S z = 0 sector, this is an N -times-degenerate multiplet. From the singlet component one obtains the ground state of the Hubbard model for U = 0 + , which is unique, in agreement with a theorem of Lieb. The wave function demonstrates an antiferromagnetic order, a lattice step translation being equivalent to a spin flip. We show that the total momentum vanishes, while the point symmetry is s or d for even or odd N /2, respectively.
Mott physics in topological bands has recently attracted considerable attention, particularly in the context of twisted bilayer graphene (TBG). However, the essential ingredients for stabilizing this physics remain unclear. Here, we demonstrate a quantum Mott semimetal phase as the ground state within a one-dimensional spinful Hubbard model featuring only one orbital per unit cell, protected by inversion and particle-hole symmetries. We start from a two-orbital model where a localized f orbital on the A sublattice hybridizes with a delocalized c orbital on the B sublattice. Projecting the f-orbital Hubbard U onto the active flat band yields a lattice model with Wannier orbitals centered on the B sublattice. Similar to TBG, a momentum-space scale k∗ emerges, setting the interaction range in the projected model to 1/k∗. While the ground state is ferromagnetic with only the Hubbard U, introducing an inter-site antiferromagnetic spin coupling J stabilizes a Mott semimetal∗ phase with a central charge c=3. Using exact diagonalization (ED) and density matrix renormalization group (DMRG) methods, we show that this phase hosts a spinful Dirac fermion coexisting with a neutral spin mode -- analogous to the fractionalized Fermi liquid (FL∗) phase in higher dimensions. Furthermore, breaking particle-hole (PH) symmetry via dispersion transforms the Mott semimetal into a Mott insulator, which is separated from a distinct Mott insulating phase by a continuous transition with a polarization jump of 1/2. Our work provides the first unbiased evidence of a Mott semimetal ground state and demonstrates that this 1D model captures some essential aspects of TBG physics, despite lacking a Wannier obstruction.