In quantum simulation, many-body phenomena are probed in controllable quantum systems. Recently, simulation of Bose-Hubbard Hamiltonians using cold atoms revealed previously hidden local correlations. However, fermionic many-body Hubbard phenomena such as unconventional superconductivity and spin liquids are more difficult to simulate using cold atoms. To date the required single-site measurements and cooling remain problematic, while only ensemble measurements have been achieved. Here we simulate a two-site Hubbard Hamiltonian at low effective temperatures with single-site resolution using subsurface dopants in silicon. We measure quasi-particle tunnelling maps of spin-resolved states with atomic resolution, finding interference processes from which the entanglement entropy and Hubbard interactions are quantified. Entanglement, determined by spin and orbital degrees of freedom, increases with increasing valence bond length. We find separation-tunable Hubbard interaction strengths that are suitable for simulating strongly correlated phenomena in larger arrays of dopants, establishing dopants as a platform for quantum simulation of the Hubbard model.
In quantum simulation, many-body phenomena are probed in controllable quantum systems. Recently, simulation of Bose-Hubbard Hamiltonians using cold atoms revealed previously hidden local correlations. However, fermionic many-body Hubbard phenomena such as unconventional superconductivity and spin liquids are more difficult to simulate using cold atoms. To date the required single-site measurements and cooling remain problematic, while only ensemble measurements have been achieved. Here we simulate a two-site Hubbard Hamiltonian at low effective temperatures with single-site resolution using subsurface dopants in silicon. We measure quasi-particle tunnelling maps of spin-resolved states with atomic resolution, finding interference processes from which the entanglement entropy and Hubbard interactions are quantified. Entanglement, determined by spin and orbital degrees of freedom, increases with increasing valence bond length. We find separation-tunable Hubbard interaction strengths that are suitable for simulating strongly correlated phenomena in larger arrays of dopants, establishing dopants as a platform for quantum simulation of the Hubbard model.
In quantum simulation, many-body phenomena are probed in controllable quantum systems. Recently, single-site measurements in cold-atom quantum simulations of bosons revealed previously hidden local correlations. However, interacting Fermi-Hubbard systems, connected to canonical open many-body physics problems, are more difficult to simulate using cold atoms. To date, only ensemble measurements have been achieved, while the required single-site measurements and cooling remain formidable challenges. Here we demonstrate quantum simulation of a two-impurity Fermi system using subsurface Coulomb-coupled dopant atoms in silicon. Quasiparticle tunneling maps of individually identifiable coupled-spin states were obtained by scanning tunneling spectroscopy, at low effective temperatures, and resolving individual dopants. From the maps we quantified the two-body probability amplitudes and the entanglement of the overlapping (indistinguishable) fermions, enabling quantum simulation of a thought-experiment on a stretched covalent bond. Entanglement, determined by spin and orbital components, was found to increase with increasing dopant displacement, as interactions overcome tunneling, localizing spin. Looking ahead, we extracted the effective Hubbard interaction strength, which follows a hydrogenic, separation-tunable trend appropriate to simulate emergent Fermi-Hubbard phenomena in larger sets of dopants. These results pave the way for quantum simulation of larger scale correlated Fermi-Hubbard systems, exploiting atomic-scale fabrication to engineer interactions atom-by-atom.
The Hubbard model is an essential tool for understanding many-body physics in condensed matter systems. Artificial lattices of dopants in silicon are a promising method for the analog quantum simulation of extended Fermi-Hubbard Hamiltonians in the strong interaction regime. However, complex atom-based device fabrication requirements have meant emulating a tunable two-dimensional Fermi-Hubbard Hamiltonian in silicon has not been achieved. Here, we fabricate 3 × 3 arrays of single/few-dopant quantum dots with finite disorder and demonstrate tuning of the electron ensemble using gates and probe the many-body states using quantum transport measurements. By controlling the lattice constants, we tune the hopping amplitude and long-range interactions and observe the finite-size analogue of a transition from metallic to Mott insulating behavior. We simulate thermally activated hopping and Hubbard band formation using increased temperatures. As atomically precise fabrication continues to improve, these results enable a new class of engineered artificial lattices to simulate interactive fermionic models.
Abstract The challenge of simulating many-body models with analogue physical systems requires both experimental precision and very low operational temperatures. Atomically precise placement of dopants in Si permits the construction of nanowires by design. We investigate the suitability of these interacting electron systems as simulators of a fermionic extended Hubbard model on demand. We describe the single-particle wavefunctions as a linear combination of dopant orbitals (LCDO). The electronic states are calculated within configuration interaction (CI). Due to the peculiar oscillatory behavior of each basis orbital, properties of these chains are strongly affected by the interdonor distance R0, in a non-monotonic way. Ground state (T = 0 K) properties such as charge and spin correlations are shown to remain robust under temperatures up to 4 K for specific values of R0. The robustness of the model against disorder is also tested, allowing some fluctuation of the placement site around the target position. We suggest that finite donor chains in Si may serve as an analog simulator for strongly correlated model Hamiltonians. This simulator is, in many ways, complementary to those based on cold atoms in optical lattices—the trade-off between the tunability achievable in the latter and the survival of correlation at higher operation temperatures for the former suggests that both technologies are applicable for different regimes.
Over the past few decades, quantum simulation with ultracold atoms in optical lattices has proven to be a useful tool for understanding many-body Hamiltonians.One of its primary directions is the experimental study of the Fermi Hubbard Model (FHM), which aims to capture the fundamental properties of cuprate superconductors and has produced a wide range of intriguing phenomena, including the metal-to-Mott insulator transition, antiferromagnetism, strange metallicity, pseudogap physics, and is frequently used in the study of superconductivity. In recent years, new kinds of quantum simulators have emerged, and in particular we focus on two of them: 1) quantum simulators that study high symmetry fermions in optical lattices, and 2) optical tweezer arrays which provide a bottom-up approach to engineering the FHM. The first class of emerging quantum simulators utilizes alkaline-earth-like atoms (AEA’s) – which have large nuclear spin I (such as 173Yb and 87Sr) – to engineer the SU(N = 2I +1) symmetric FHM, where N = 2,..., 10. The second class corresponds to quantum simulators based on optical tweezer arrays which combine deterministic single tweezer ground state loading and tunneling coupled sites. Both classes of quantum simulators are interesting on their own. For example, the SU(N) FHM is generating a lot of interest due to its connections to multiorbital solid state systems and theoretical predictions of rich phase diagrams which exhibit exotic phases of matter, including unusual ones like spin liquids. On the other hand, optical tweezer arrays have gained a lot of attention since they can be used to study arbitrary geometries and access very low-entropy samples. This thesis reports on the numerical study of both quantum simulation architectures. Results of the SU(N) FHM at finite temperature, for different fillings and interaction strengths that span the non-interacting to the strongly interacting limit are obtained using Determinant Quantum Monte Carlo, Exact Diagonalization, and Numerical Linked Cluster Expansion. In particular, we focus on how thermodynamic and magnetic observables, such as the number of on-site pairs, energy, entropy, spin correlations, and structure factors depend on N, the interaction strength U/t, temperature T/t, and the anisotropy of tunneling rates. Our results demonstrate that in a homogeneous square lattice with one particle per site on average, thermodynamic observables as a function of temperature obey a universal scaling with N, and that short-range antiferromagnetic correlations are stronger for larger N and in lower dimensions. In addition, we present theory-experiment comparisons where possible, where we perform thermometry and provide a precise characterization of the equation of state of the SU(N) FHM. Our results emerge as a tool to perform thermometry in experiments from theory-experiment comparison and also provide guidance for future experiments with AEAs in optical lattices. Additionally, we present calculations using numerically precise discrete variable representation methods for two-dimensional stroboscopic tweezer arrays and compare the outcomes with experimental data. We quantify the effects of stroboscopic potentials on Hubbard parameters like the interaction strength U and the tunneling t in optical tweezer arrays and illustrate how heating from the stroboscopic potential relies on strobe frequency. Our calculations enable evaluation and optimization of two-dimensional tweezer array experiments. An important milestone of the results of this thesis corresponds to the record temperatures and entropies achieved with these quantum simulation architectures: 1) the coldest fermions ever created in nature in absolute temperature and in cold atoms, and 2) the lowest entropy per particle fermions ever created, albeit in a so-far small system.
An outstanding problem in physics is how to understand strongly interacting quantum many-body systems such as the quark-gluon plasma, neutron stars, superfluid 4He, and the high-temperature superconducting cuprates. The physics approach to this problem is to reduce these complex systems to minimal models that are believed to retain relevant phenomenology. For example, the Hubbard model — the focus of this thesis — describes quantum particles tunneling between sites of a lattice with on-site interactions. The Hubbard model is conjectured to describe the low-energy charge and spin properties of high-temperature superconducting cuprates. Thus far, there are no analytic solutions to the Hubbard model, and numerical calculations are difficult and even impossible in some regimes (e.g., the Fermi-Hubbard model away from half-filling). Therefore, whether the Hubbard model is a minimal model for the cuprates remains unresolved. In the face of these difficulties, a new approach has emerged — quantum simulation. The premise of quan-tum simulation is to perform experiments on a quantum system that is well-described by the model we are trying to study, has tunable parameters, and is easily probed. Ultracold atoms trapped in optical lattices are an ideal candidate for quantum simulation of the Hubbard models. This thesis describes work on two such systems: a 87Rb (boson) optical lattice experiment in the group of Brian DeMarco at the University
Experiments with ultracold atoms in optical lattices provide outstanding opportunities to realize exotic quantum states due to a high degree of tunability and control. In this thesis, I present experiments that extend this control from global parameters to the level of individual particles. \n\nUsing a quantum gas microscope for 87Rb, we have developed a single-site addressing scheme based on digital amplitude holograms. The system self-corrects for aberrations in the imaging setup and creates arbitrary beam profiles. We are thus able to shape optical potentials on the scale of single lattice sites and control the dynamics of individual atoms. \n\nWe study the role of quantum statistics and interactions in the Bose-Hubbard model on the fundamental level of two particles. Bosonic quantum statistics are apparent in the Hong-Ou-Mandel interference of massive particles, which we observe in tailored double-well potentials. These underlying statistics, in combination with tunable repulsive interactions, dominate the dynamics in single- and two-particle quantum walks. We observe highly coherent position-space Bloch oscillations, bosonic bunching in Hanbury Brown-Twiss interference and the fermionization of strongly interacting bosons.\n\nMany-body states of indistinguishable quantum particles are characterized by large-scale spatial entanglement, which is difficult to detect in itinerant systems. Here, we extend the concept of Hong-Ou-Mandel interference from individual particles to many-body states to directly quantify entanglement entropy. We perform collective measurements on two copies of a quantum state and detect entanglement entropy through many-body interference. We measure the second order Rényi entropy in small Bose-Hubbard systems and detect the buildup of spatial entanglement across the superfluid-insulator transition.\n\nOur experiments open new opportunities for the single-particle-resolved preparation and characterization of many-body quantum states.