Abstract:The Hubbard model is one of the primary models for understanding the essential many-body physics in condensed matter systems such as Mott insulators and cuprate high-Tc superconductors. Recent advances in atomically precise fabrication in silicon using scanning tunneling microscopy (STM) have made possible atom-by-atom fabrication of single and few-dopant quantum dots and atomic-scale control of tunneling in dopant-based devices. However, the complex fabrication requirements of multi-component devices have meant that emulating two-dimensional (2D) Fermi-Hubbard physics using these systems has not been demonstrated. Here, we overcome these challenges by integrating the latest developments in atomic fabrication and demonstrate the analog quantum simulation of a 2D extended Fermi-Hubbard Hamiltonian using STM-fabricated 3x3 arrays of single/few-dopant quantum dots. We demonstrate low-temperature quantum transport and tuning of the electron ensemble using in-plane gates as efficient probes to characterize the many-body properties, such as charge addition, tunnel coupling, and the impact of disorder within the array. By controlling the array lattice constants with sub-nm precision, we demonstrate tuning of the hopping amplitude and long-range interactions and observe the finite-size analogue of a transition from Mott insulating to metallic behavior in the array. By increasing the measurement temperature, we simulate the effect of thermally activated hopping and Hubbard band formation in transport spectroscopy. We compare the analog quantum simulations with numerically simulated results to help understand the energy spectrum and resonant tunneling within the array. The results demonstrated in this study serve as a launching point for a new class of engineered artificial lattices to simulate the extended Fermi-Hubbard model of strongly correlated materials.

Quantum Monte Carlo study of mass-imbalanced Hubbard models

Benchmark Study of the Two-Dimensional Hubbard Model with Auxiliary-Field Quantum Monte Carlo Method

A Monte Carlo Study of the Half-Filled Two-Dimensional Hubbard Model

Determinant Quantum Monte Carlo for the Half-Filled Hubbard Model with Nonlocal Density-Density Interactions

Efficient continuous-time quantum Monte Carlo algorithm for fermionic lattice models

Magnetic, thermodynamic and dynamical properties of the three-dimensional fermionic Hubbard model: a comprehensive Monte Carlo study

Quantum phase diagram of fermion mixtures with population imbalance in one-dimensional optical lattices

Quantum magnetic properties of the S U (2 N) Hubbard model in the square lattice: A quantum Monte Carlo study

Quantum Monte Carlo simulations of thermodynamic properties of<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>SU</mml:mi><mml:mo>(</mml:mo><mml:mn>2</mml:mn><mml:mi>N</mml:mi><mml:mo>)</mml:mo></mml:math>ultracold fermions in optical lattices

Boosting Determinant Quantum Monte Carlo with Submatrix Updates: Unveiling the Phase Diagram of the 3D Hubbard Model

Quantum Monte Carlo Simulations Of Thermodynamic Properties Of Su(2n) Ultracold Fermions In Optical Lattices

Simulating the Fermi-Hubbard model with long-range hopping on a quantum computer

Quantum Monte Carlo Simulation of the two-dimensional ionic Hubbard model

Quantum Magnetic Properties Of The Su(2n) Hubbard Model In The Square Lattice: A Quantum Monte Carlo Study

Solving the Fermion Sign Problem in Quantum Monte Carlo Simulations by Majorana Representation

The Sign Problem in Quantum Monte Carlo Simulations of the Hubbard Model in a Weak Interaction Region

Numerical Results on the Short-Range Spin Correlation Functions in the Ground State of the Two-Dimensional Hubbard Model

Quantum Simulation of an Extended Fermi-Hubbard Model Using a 2D Lattice of Dopant-based Quantum Dots

Perturbative solution of fermionic sign problem in lattice Quantum Monte Carlo

Quench Dynamics in the One-Dimensional Mass-Imbalanced Ionic Hubbard Model

Solutions of the two-dimensional Hubbard model: benchmarks and results from a wide range of numerical algorithms