It has been demonstrated that small plaquettes of quantum dot spin qubits are capable of simulating condensed matter phenomena that arise from the Hubbard model, such as the collective Coulomb blockade and Nagaoka ferromagnetism. Motivated by recent materials developments, we investigate a bilayer arrangement of quantum dots with four dots in each layer which exhibits a complex ground state behavior. We find using a generalized Hubbard model with long-range Coulomb interactions, several distinct magnetic configurations occur as the Coulomb interaction strength is varied, with possible ground states that are ferromagnetic, antiferromagnetic, or having both one antiferromagnetic and one ferromagnetic layer. We map out the full phase diagram of the system as it depends on the inter- and intralayer Coulomb interaction strengths, and find that for a single layer, a similar but simpler effect occurs. We also predict interesting contrasts among electron, hole, and electron-hole bilayer systems arising from complex correlation physics. Observing the predicted magnetic configuration in already-existing few-dot semiconductor bilayer structures could prove to be an important assessment of current experimental quantum dot devices, particularly in the context of spin-qubit-based analog quantum simulations.
It has been demonstrated that small plaquettes of quantum dot spin qubits are capable of simulating condensed matter phenomena that arise from the Hubbard model, such as the collective Coulomb blockade and Nagaoka ferromagnetism. Motivated by recent materials developments, we investigate a bilayer arrangement of quantum dots with four dots in each layer which exhibits a complex ground state behavior. We find using a generalized Hubbard model with long-range Coulomb interactions, several distinct magnetic configurations occur as the Coulomb interaction strength is varied, with possible ground states that are ferromagnetic, antiferromagnetic, or having both one antiferromagnetic and one ferromagnetic layer. We map out the full phase diagram of the system as it depends on the inter- and intralayer Coulomb interaction strengths, and find that for a single layer, a similar but simpler effect occurs. We also predict interesting contrasts among electron, hole, and electron-hole bilayer systems arising from complex correlation physics. Observing the predicted magnetic configuration in already-existing few-dot semiconductor bilayer structures could prove to be an important assessment of current experimental quantum dot devices, particularly in the context of spin-qubit-based analog quantum simulations.
Abstract : The authors use electrostatically defined quantum dot arrays to study quantum simulation, examining both large arrays to explore the ground state and excitations of the Fermi-Hubbard model, and to examine very small arrays of 3-4 quantum dots as the basis for simulating superexchange, frustration, magnetism, etc. The large arrays were created by defining a grid-shaped gate on a semiconductor, forcing the electrons to see a periodic potential and be confined in a lattice of sites, where the low temperature reachable by the system in a dilution refrigerator allows examining the phase diagram of the model in a regime where spin correlations have an important influence. Initial plans to perform transport measurement were unsuccessful as the clean room realization of contacts to the two layers of the bilayered heterostructure had very low yield, so the authors performed capacitance spectroscopy measurement instead. This provides direct measurement of density of states, and the lithography process involved allowed grids with 150-200nm periods to be successfully written. The small array examination analyzed how a system of four quantum dots on a square and filled with three electrons allowed the authors to study magnetism of few-body systems, examining whether the Nagaoka theorem could be realizable with quantum dots (conclusion: the ferromagnetic-paramagnetic state is in reach experimentally) and examining the experimental techniques to test the spin-structure of the ground state and induce the transition (finding: application of magnetic field can break the ferromagnetic ground state, in contrast to the usual). The authors also studied a linear array of three dots, resulting in publications in Nature Nanotechnology, Applied Physics Letters, and Annalen Der Physik.
Nagaoka showed that ferromagnetism can arise from a basic Hubbard model of certain lattices under special conditions. Motivated by recent experimental advances in simulating such systems using small plaquettes of a few quantum dots, this work investigates what geometrical arrangement of these dots and what number of electrons are necessary to lead to nontrivial ferromagnetic ground states. This includes many situations not covered by Nagaoka's theorem, which is only valid under specific constraints.
Recent progress of quantum simulators provides insight into the fundamental problems of strongly correlated systems. To adequately assess the accuracy of these simulators, the precise modeling of the many-body physics, with accurate model parameters, is crucially important. In this paper, we employed an ab initio exact diagonalization framework to compute the correlated physics of a few electrons in artificial potentials. We apply this approach to a quantum-dot system and study the magnetism of the correlated electrons, obtaining good agreement with recent experimental measurements in a plaquette. Through control of dot potentials and separation, including geometric manipulation of tunneling, we examine the Nagaoka transition and determine the robustness of the ferromagnetic state. While the Nagaoka theorem considers only a single-band Hubbard model, in this work we perform extensive ab initio calculations that include realistic multiorbital conditions in which the level splitting is smaller than the interactions. This simulation complements the experiments and provides insight into the formation of ferromagnetism in correlated systems. More generally, our calculation sets the stage for further theoretical analysis of analog quantum simulators at a quantitative level.
A wide variety of experimental platforms, ranging from semiconductor quantum-dot arrays to moiré materials, have recently emerged as powerful quantum simulators for studying the Hubbard model and its variants. Motivated by these developments, here, we investigate a generalization of the Hubbard model which includes the effects of long-range Coulomb interactions. Working on finite-sized two-dimensional square and triangular lattices, we use exact diagonalization and density-matrix renormalization group calculations to probe the magnetic structure of the ground state in the strong-coupling regime, where U (the onsite repulsion) ≫t (the nearest-neighbor hopping). For small electron dopings above the half-filled antiferromagnet, we numerically uncover a rich variety of magnetically ordered states, and in conjunction with theoretical arguments, infer the phase diagram of the system as a function of doping and interaction strengths. In particular, we find that the inclusion of long-range Coulomb interactions induces an instability of high-spin states$\unicode{x2014}$such as the saturated Nagaoka ferromagnet$\unicode{x2014}$towards phase separation and stripe ordering. We also present proposals for the observation of some of our key findings in experiments that would shed further light on this paradigmatic strongly correlated system.
A phase diagram of two Mott-Hubbard planes interacting with a short-range Coulomb repulsion is presented. Considering the case of equal amount of doping by holes in one layer as electrons in the other, a holon-doublon inter-layer exciton formation is shown to be a natural consequence of Coulomb attraction. Quasiparticle spectrum is gapped and incoherent below a critical doping δc due to the formation of excitons. A spin liquid insulator (SLI) phase is thus realized without the lattice frustration. The critical value δc sensitively depends on the inter-layer interaction strength. In the tJ model description of each layer with the d-wave pairing, δc marks the crossover between SLI and d-wave superconductor. The SLI phase, despite being non-superconducting and charge-gapped, still shows electromagnetic response similar to that of a superfluid due to the exciton transport. Including antiferromagnetic order in the tJ model introduces magnetically ordered phases at low doping and pushes the spin liquid phase to a larger inter-layer interaction strength and higher doping concentrations.
We are entering the era of the second quantum revolution, where we aim to harness the power of quantum mechanics to create new technologies. Quantum technologies have the potential to revolutionize the fields of simulation, computation, communication, sensing, metrology, and many others. Here we discuss analog quantum simulation, which has attracted a lot of attention in the last few years from several platforms. Although arrays of gate-defined quantum dots exhibit significant potential for analog simulation, example experiments have been few and far between. This thesis focuses on simulating the Fermi-Hubbard model using two dimensional (2d) arrays of quantum dots. The first experiment describes the creation and measurement of a 2x2 quantum dot array. Historically, most experiments with quantum dots have been performed with linear arrays due to the relative ease of fabrication. We introduce a bi-layer gate structure, facilitated by the lift-off of sputtered silicon nitride, to create the 2x2 dot array. This gate design enables us to achieve unprecedented tunability of the tunnel coupling between all nearest-neighbor pairs of dots in 2d arrays. We also demonstrate individual control over the chemical potential and the electron occupation of each dot along with accurate measurement of the on-site and inter-site interaction terms. The use of virtual gates significantly aids in the tuning of tunnel coupling and chemical potential. The demonstrated high degree of control of the system along with fast single-shot spin-readout achieved through Pauli spin blockade establish this dot array as a promising simulator of the Fermi-Hubbard model. The 2x2 dot array is used to simulate Nagaoka ferromagnetism in the next experiment. This form of itinerant ferromagnetism arises from the Fermi-Hubbard model, and was first shown analytically in the limit of infinite interaction strengths and infinite lattices by Nagaoka in 1966. Nagaoka ferromagnetism has been a topic of rigorous theoretical studies ever since, but its experimental signature has eluded us for more than five decades. In this experiment, we load the four dot plaquette with three electrons and demonstrate the emergence of spontaneous ferromagnetism by measuring the spin correlation of two out of the three electrons. Changing the topology of the array to an open chain is shown to destroy the ferromagnetic signature, consistent with the Lieb-Mattis theorem. We also show indications that this ferromagnetic ground state can be destroyed by applying a perpendicular magnetic field, unlike most other forms of ferromagnetism. However, this ground state shows striking robustness to the offset in the local potential of any dot. This is the first experimental verification of Nagaoka’s prediction as well as the first simulation of magnetism using quantumdot arrays. The final experiment takes a different approach to simulate the Fermi-Hubbardmodel with a large 2d array of quantum dots. The dot array is created using only three gates in a top-down approach. This allows for only global control over the electron filling and tunnel coupling of the dots, contrary to the previous experiments. The readout is performed with capacitance spectroscopy, which allows us to directly probe the density of states of the two-dimensional electron systems. We measure the disorder levels and optimize both substrates and gating strategies to induce periodic potential, sufficiently stronger than the disorder level, at the 2d electron gas. Although we demonstrate a novel platformfor the realization of artificial lattices of interacting particles, this effort is currently limited by the substrate inhomogeneity.