Simulation of Hubbard model physics in WSe2/WS2 moire superlattices

Study of WSe2/WS2 moire superlattices reveals the phase diagram of the triangular-lattice Hubbard model, including a Mott insulating state at half-filling and a possible magnetic quantum phase transition near 0.6 filling. The Hubbard model, formulated by physicist John Hubbard in the 1960s(1), is a simple theoretical model of interacting quantum particles in a lattice. The model is thought to capture the essential physics of high-temperature superconductors, magnetic insulators and other complex quantum many-body ground states(2,3). Although the Hubbard model provides a greatly simplified representation of most real materials, it is nevertheless difficult to solve accurately except in the one-dimensional case(2,3). Therefore, the physical realization of the Hubbard model in two or three dimensions, which can act as an analogue quantum simulator (that is, it can mimic the model and simulate its phase diagram and dynamics(4,5)), has a vital role in solving the strong-correlation puzzle, namely, revealing the physics of a large number of strongly interacting quantum particles. Here we obtain the phase diagram of the two-dimensional triangular-lattice Hubbard model by studying angle-aligned WSe2/WS2 bilayers, which form moire superlattices(6) because of the difference between the lattice constants of the two materials. We probe the charge and magnetic properties of the system by measuring the dependence of its optical response on an out-of-plane magnetic field and on the gate-tuned carrier density. At half-filling of the first hole moire superlattice band, we observe a Mott insulating state with antiferromagnetic Curie-Weiss behaviour, as expected for a Hubbard model in the strong-interaction regime(2,3,7-9). Above half-filling, our experiment suggests a possible quantum phase transition from an antiferromagnetic to a weak ferromagnetic state at filling factors near 0.6. Our results establish a new solid-state platform based on moire superlattices that can be used to simulate problems in strong-correlation physics that are described by triangular-lattice Hubbard models.