A cold-atom Fermi–Hubbard antiferromagnet

Magnetic Phase Diagram of the Three-Dimensional Doped Hubbard Model

We establish the phase diagram of the Hubbard model on a cubic lattice for a wide range of temperatures, dopings, and interaction strengths, considering both commensurate and incommensurate magnetic orders. We use the dynamical mean-field theory together with an efficient method to compute the free energy which enable the determination of the correct ordering vectors. Besides an antiferromagnetic state close to half filling, we identify a number of different magnetic spiral phases with ordering vectors (q,π,π), (q,q,π), and (q,q,q), as well as a region with close competition between them, hinting at spatial phase separation or at the onset of a stripe phase. Additionally, we extensively study several thermodynamic properties with direct relevance to cold-atom experiments: the entropy, energy, and double occupancy.

Homogeneous Fermionic Hubbard Gases in a Flattop Optical Lattice

Fermionic atoms in a large-scale, homogeneous optical lattice provide an ideal quantum simulator for investigating the fermionic Hubbard model, yet achieving this remains challenging. Here, by developing a hybrid potential that integrates a flat-top optical lattice with an optical box trap, we successfully realize the creation of three-dimensional, homogeneous fermionic Hubbard gases across approximately 8×10^{5} lattice sites. This homogeneous system enables us to capture a well-defined energy band occupation that aligns perfectly with the theoretical calculations for a zero-temperature, ideal fermionic Hubbard model. Furthermore, by employing novel radio-frequency spectroscopy, we precisely measure the doublon fraction D as a function of interaction strength U and temperature T, respectively. The crossover from metal to Mott insulator is detected, where D smoothly decreases with increasing U. More importantly, we observe a nonmonotonic temperature dependence in D, revealing the Pomeranchuk effect and the development of extended antiferromagnetic correlations.

Formation of individual stripes in a mixed-dimensional cold-atom Fermi-Hubbard system

The relation between d-wave superconductivity and stripes is fundamental to the understanding of ordered phases in high-temperature cuprate superconductors1-6. These phases can be strongly influenced by anisotropic couplings, leading to higher critical temperatures, as emphasized by the recent discovery of superconductivity in nickelates7-10. Quantum simulators with ultracold atoms provide a versatile platform to engineer such couplings and to observe emergent structures in real space with single-particle resolution. Here we show, to our knowledge, the first signatures of individual stripes in a cold-atom Fermi-Hubbard quantum simulator using mixed-dimensional (mixD) settings. Increasing the energy scale of hole-hole attraction to the spin exchange energy, we access the interesting crossover temperature regime in which stripes begin to form11. We observe extended, attractive correlations between hole dopants and find an increased probability of forming larger structures akin to individual stripes. In the spin sector, we study correlation functions up to the third order and find results consistent with stripe formation. These observations are interpreted as a precursor to the stripe phase, which is characterized by interleaved charge and spin density wave ordering with fluctuating lines of dopants separating domains of opposite antiferromagnetic order12-14.

Symmetry: A Fundamental Resource for Quantum Coherence and Metrology

We introduce a new paradigm for the preparation of deeply entangled states useful for quantum metrology. We show that, when the quantum state is an eigenstate of an operator A, observables G which are completely off diagonal with respect to A have purely quantum fluctuations, as quantified by the quantum Fisher information, namely, F_{Q}(G)=4⟨G^{2}⟩. This property holds regardless of the purity of the quantum state, and it implies that off-diagonal fluctuations represent a metrological resource for phase estimation. In particular, for many-body systems such as quantum spin ensembles or bosonic gases, the presence of off-diagonal long-range order (for a spin observable or for bosonic operators) directly translates into a metrological resource, provided that the system remains in a well-defined symmetry sector. The latter is defined, e.g., by one component of the collective spin or by its parity in spin systems; and by the particle number for bosons. Our results establish the optimal use for metrology of arbitrarily non-Gaussian quantum correlations in a large variety of many-body systems.

Quantum Turnstiles for Robust Measurement of Full Counting Statistics

We present a scalable protocol for measuring full counting statistics (FCS) in experiments or tensor-network simulations. In this method, an ancilla in the middle of the system acts as a turnstile, with its phase keeping track of the time-integrated particle flux. Unlike quantum gas microscopy, the turnstile protocol faithfully captures FCS starting from number-indefinite initial states or in the presence of noisy dynamics. In addition, by mapping the FCS onto a single-body observable, it allows for stable numerical calculations of FCS using approximate tensor-network methods. We demonstrate the wide-ranging utility of this approach by computing the FCS of the transferred magnetization in a Floquet Heisenberg spin chain, as studied in a recent experiment with superconducting qubits, as well as the FCS of charge transfer in random circuits.

Many-body physics of ultracold alkaline-earth atoms with SU(N)-symmetric interactions

Symmetries play a crucial role in understanding phases of matter and the transitions between them. Theoretical investigations of quantum models with SU(N) symmetry have provided important insights into many-body phenomena. However, these models have generally remained a theoretical idealization, since it is very difficult to exactly realize the SU(N) symmetry in conventional quantum materials for largeN. Intriguingly however, in recent years, ultracold alkaline-earth-atom (AEA) quantum simulators have paved the path to realize SU(N)-symmetric many-body models, whereNis tunable and can be as large as 10. This symmetry emerges due to the closed shell structure of AEAs, thereby leading to a perfect decoupling of the electronic degrees of freedom from the nuclear spin. In this work, we provide a systematic review of recent theoretical and experimental work on the many-body physics of these systems. We first discuss the thermodynamic properties and collective modes of trapped Fermi gases, highlighting the enhanced interaction effects that appear asNincreases. We then discuss the properties of the SU(N) Fermi-Hubbard model, focusing on some of the major experimental achievements in this area. We conclude with a compendium highlighting some of the significant theoretical progress on SU(N) lattice models and a discussion of some exciting directions for future research.

Learning Fermionic Correlations by Evolving with Random Translationally Invariant Hamiltonians

Schemes of classical shadows have been developed to facilitate the readout of digital quantum devices, but similar tools for analog quantum simulators are scarce and experimentally impractical. In this Letter, we provide a measurement scheme for fermionic quantum devices that estimates second and fourth order correlation functions by means of free fermionic, translationally invariant evolutions-or quenches-and measurements in the mode occupation number basis. We precisely characterize what correlation functions can be recovered and equip the estimates with rigorous bounds on sample complexities, a particularly important feature in light of the difficulty of getting good statistics in reasonable experimental platforms, with measurements being slow. Finally, we demonstrate how our procedure can be approximately implemented with just nearest-neighbor, translationally invariant hopping quenches, a very plausible procedure under current experimental requirements and requiring only random time evolution with respect to a single native Hamiltonian. On a conceptual level, this Letter brings the idea of classical shadows to the realm of large scale analog quantum simulators.

Quantum algorithms for scientific computing

Quantum computing promises to provide the next step up in computational power for diverse application areas. In this review, we examine the science behind the quantum hype, and the breakthroughs required to achieve true quantum advantage in real world applications. Areas that are likely to have the greatest impact on high performance computing (HPC) include simulation of quantum systems, optimization, and machine learning. We draw our examples from electronic structure calculations and computational fluid dynamics which account for a large fraction of current scientific and engineering use of HPC. Potential challenges include encoding and decoding classical data for quantum devices, and mismatched clock speeds between classical and quantum processors. Even a modest quantum enhancement to current classical techniques would have far-reaching impacts in areas such as weather forecasting, aerospace engineering, and the design of 'green' materials for sustainable development. This requires significant effort from the computational science, engineering and quantum computing communities working together.

Conformational Tuning of Magnetic Interactions in Coupled Nanographenes

Phenalenyl (C13H9) is an open-shell spin-1/2 nanographene. Using scanning tunneling microscopy (STM) inelastic electron tunneling spectroscopy (IETS), covalently bonded phenalenyl dimers have been shown to feature conductance steps associated with singlet-triplet excitations of a spin-1/2 dimer with antiferromagnetic exchange. Here, we address the possibility of tuning the magnitude of the exchange interactions by varying the dihedral angle between the two molecules within a dimer. Theoretical methods ranging from density functional theory calculations to many-body model Hamiltonians solved within different levels of approximation are used to explain STM-IETS measurements of phenalenyl dimers on a hexagonal boron nitride (h-BN)/Rh(111) surface, which exhibit signatures of twisting. By means of first-principles calculations, we also propose strategies to induce sizable twist angles in surface-adsorbed phenalenyl dimers via functional groups, including a photoswitchable scheme. This work paves the way toward tuning magnetic couplings in carbon-based spin chains and two-dimensional lattices.

Altermagnetic Anomalous Hall Effect Emerging from Electronic Correlations

While altermagnetic materials are characterized by a vanishing net magnetic moment, their symmetry in principle allows for the existence of an anomalous Hall effect. Here, we introduce a model with altermagnetism in which the emergence of an anomalous Hall effect is driven by interactions. This model is grounded in a modified Kane-Mele framework with antiferromagnetic spin-spin correlations. Quantum Monte Carlo simulations show that the system undergoes a finite temperature phase transition governed by a primary antiferromagnetic order parameter accompanied by a secondary one of Haldane type. The emergence of both orders turns the metallic state of the system, away from half-filling, to an altermagnet with a finite anomalous Hall conductivity. A mean field ansatz corroborates these results, which pave the way into the study of correlation induced altermagnets with finite Berry curvature.

Antiferromagnetic phase transition in a 3D fermionic Hubbard model

The fermionic Hubbard model (FHM)1 describes a wide range of physical phenomena resulting from strong electron-electron correlations, including conjectured mechanisms for unconventional superconductivity. Resolving its low-temperature physics is, however, challenging theoretically or numerically. Ultracold fermions in optical lattices2,3 provide a clean and well-controlled platform offering a path to simulate the FHM. Doping the antiferromagnetic ground state of a FHM simulator at half-filling is expected to yield various exotic phases, including stripe order4, pseudogap5, and d-wave superfluid6, offering valuable insights into high-temperature superconductivity7-9. Although the observation of antiferromagnetic correlations over short10 and extended distances11 has been obtained, the antiferromagnetic phase has yet to be realized as it requires sufficiently low temperatures in a large and uniform quantum simulator. Here we report the observation of the antiferromagnetic phase transition in a three-dimensional fermionic Hubbard system comprising lithium-6 atoms in a uniform optical lattice with approximately 800,000 sites. When the interaction strength, temperature and doping concentration are finely tuned to approach their respective critical values, a sharp increase in the spin structure factor is observed. These observations can be well described by a power-law divergence, with a critical exponent of 1.396 from the Heisenberg universality class12. At half-filling and with optimal interaction strength, the measured spin structure factor reaches 123(8), signifying the establishment of an antiferromagnetic phase. Our results provide opportunities for exploring the low-temperature phase diagram of the FHM.

Realizing Altermagnetism in Fermi-Hubbard Models with Ultracold Atoms

Altermagnetism represents a type of collinear magnetism, that is in some aspects distinct from ferromagnetism and from conventional antiferromagnetism. In contrast to the latter, sublattices of opposite spin are related by spatial rotations and not only by translations and inversions. As a result, altermagnets have spin-split bands leading to unique experimental signatures. Here, we show theoretically how a d-wave altermagnetic phase can be realized with ultracold fermionic atoms in optical lattices. We propose an altermagnetic Hubbard model with anisotropic next-nearest neighbor hopping and obtain the Hartree-Fock phase diagram. The altermagnetic phase separates in a metallic and an insulating phase and is robust over a large parameter regime. We show that one of the defining characteristics of altermagnetism, the anisotropic spin transport, can be probed with trap-expansion experiments.

Symmetry-Broken Perturbation Theory to Large Orders in Antiferromagnetic Phases

We introduce a spin-symmetry-broken extension of the connected determinant algorithm [Riccardo Rossi, Determinant diagrammatic Monte Carlo algorithm in the thermodynamic limit, Phys. Rev. Lett. 119, 045701 (2017).PRLTAO0031-900710.1103/PhysRevLett.119.045701]. The resulting systematic perturbative expansions around an antiferromagnetic state allow for numerically exact calculations directly inside a magnetically ordered phase. We show new precise results for the magnetic phase diagram and thermodynamics of the three-dimensional cubic Hubbard model at half-filling. With detailed computations of the order parameter in the low to intermediate-coupling regime, we establish the Néel phase boundary. The critical behavior in its vicinity is shown to be compatible with the O(3) Heisenberg universality class. By determining the evolution of the entropy with decreasing temperature through the phase transition we identify the different physical regimes at U/t=4. We provide quantitative results for several thermodynamic quantities deep inside the antiferromagnetic dome up to large interaction strengths and investigate the crossover between the Slater and Heisenberg regimes.

Antiferromagnetic Bosonic t-J Models and Their Quantum Simulation in Tweezer Arrays

The combination of optical tweezer arrays with strong interactions-via dipole exchange of molecules and Van der Waals interactions of Rydberg atoms-has opened the door for the exploration of a wide variety of quantum spin models. A next significant step will be the combination of such settings with mobile dopants. This will enable one to simulate the physics believed to underlie many strongly correlated quantum materials. Here, we propose an experimental scheme to realize bosonic t-J models via encoding the local Hilbert space in a set of three internal atomic or molecular states. By engineering antiferromagnetic (AFM) couplings between spins, competition between charge motion and magnetic order similar to that in high-T_{c} cuprates can be realized. Since the ground states of the 2D bosonic AFM t-J model we propose to realize have not been studied extensively before, we start by analyzing the case of two dopants-the simplest instance in which their bosonic statistics plays a role-and compare our results to the fermionic case. We perform large-scale density matrix renormalization group calculations on six-legged cylinders, and find a strong tendency for bosonic holes to form stripes. This demonstrates that bosonic, AFM t-J models may contain similar physics as the collective phases in strongly correlated electrons.

Constants of motion characterizing continuous symmetry-broken phases

We present a theory characterizing the phases emerging as a consequence of continuous symmetry-breaking in quantum and classical systems. In symmetry-breaking phases, dynamics is restricted due to the existence of a set of conserved charges derived from the order parameter of the phase transition. Their expectation values are determined by the privileged direction appearing in the ordered phase as a consequence of symmetry breaking, and thus they can be used to determine whether this direction is well defined or has quantum fluctuations. Our theory is numerically exemplified via the two-dimensional limit of the vibron model, a fully connected system invariant under a rotation operator which generates the continuous symmetry-breaking.

Engineering of a Low-Entropy Quantum Simulator for Strongly Correlated Electrons Using Cold Atoms with SU(N)-Symmetric Interactions

An advanced cooling scheme, incorporating entropy engineering, is vital for isolated artificial quantum systems designed to emulate the low-temperature physics of strongly correlated electron systems. This study theoretically demonstrates a cooling method employing multicomponent Fermi gases with SU(N)-symmetric interactions, focusing on the case of ^{173}Yb atoms in a two-dimensional optical lattice. Adiabatically introducing a nonuniform state-selective laser gives rise to two distinct subsystems: a central low-entropy region, exclusively composed of two specific spin components, acts as a quantum simulator for strongly correlated electron systems, while the surrounding N-component mixture retains a significant portion of the entropy of the system. The total particle numbers for each component are good quantum numbers, creating a sharp boundary for the two-component region. The cooling efficiency is assessed through extensive finite-temperature Lanczos calculations. The results lay the foundation for quantum simulations of two-dimensional systems of Hubbard or Heisenberg type, offering crucial insights into intriguing low-temperature phenomena in condensed-matter physics.

Tunable Superdiffusion in Integrable Spin Chains Using Correlated Initial States

Although integrable spin chains host only ballistically propagating particles, they can still feature diffusive charge transfer. This diffusive charge transfer originates from quasiparticle charge fluctuations inherited from the initial state's magnetization Gaussian fluctuations. We show that ensembles of initial states with quasi-long-range correlations lead to superdiffusive charge transfer with a tunable dynamical exponent. We substantiate our prediction with numerical simulations and discuss how finite time and finite size effects might cause deviations.

Cold-Atom Elevator: From Edge-State Injection to the Preparation of Fractional Chern Insulators

Optical box traps offer new possibilities for quantum-gas experiments. Building on their exquisite spatial and temporal control, we propose to engineer system-reservoir configurations using box traps, in view of preparing and manipulating topological atomic states in optical lattices. First, we consider the injection of particles from the reservoir to the system: this scenario is shown to be particularly well suited to activating energy-selective chiral edge currents, but also to prepare fractional Chern insulating ground states. Then, we devise a practical evaporative-cooling scheme to effectively cool down atomic gases into topological ground states. Our open-system approach to optical-lattice settings provides a new path for the investigation of ultracold quantum matter, including strongly correlated and topological phases.

Stable quantum-correlated many-body states through engineered dissipation

Engineered dissipative reservoirs have the potential to steer many-body quantum systems toward correlated steady states useful for quantum simulation of high-temperature superconductivity or quantum magnetism. Using up to 49 superconducting qubits, we prepared low-energy states of the transverse-field Ising model through coupling to dissipative auxiliary qubits. In one dimension, we observed long-range quantum correlations and a ground-state fidelity of 0.86 for 18 qubits at the critical point. In two dimensions, we found mutual information that extends beyond nearest neighbors. Lastly, by coupling the system to auxiliaries emulating reservoirs with different chemical potentials, we explored transport in the quantum Heisenberg model. Our results establish engineered dissipation as a scalable alternative to unitary evolution for preparing entangled many-body states on noisy quantum processors.

Equation of State and Thermometry of the 2D SU(N) Fermi-Hubbard Model

We characterize the equation of state (EoS) of the SU(N>2) Fermi-Hubbard Model (FHM) in a two-dimensional single-layer square optical lattice. We probe the density and the site occupation probabilities as functions of interaction strength and temperature for N=3, 4, and 6. Our measurements are used as a benchmark for state-of-the-art numerical methods including determinantal quantum Monte Carlo and numerical linked cluster expansion. By probing the density fluctuations, we compare temperatures determined in a model-independent way by fitting measurements to numerically calculated EoS results, making this a particularly interesting new step in the exploration and characterization of the SU(N) FHM.

Observation and quantification of the pseudogap in unitary Fermi gases

The microscopic origin of high-temperature superconductivity in cuprates remains unknown. It is widely believed that substantial progress could be achieved by better understanding of the pseudogap phase, a normal non-superconducting state of cuprates1,2. In particular, a central issue is whether the pseudogap could originate from strong pairing fluctuations3. Unitary Fermi gases4,5, in which the pseudogap-if it exists-necessarily arises from many-body pairing, offer ideal quantum simulators to address this question. Here we report the observation of a pair-fluctuation-driven pseudogap in homogeneous unitary Fermi gases of lithium-6 atoms, by precisely measuring the fermion spectral function through momentum-resolved microwave spectroscopy and without spurious effects from final-state interactions. The temperature dependence of the pairing gap, inverse pair lifetime and single-particle scattering rate are quantitatively determined by analysing the spectra. We find a large pseudogap above the superfluid transition temperature. The inverse pair lifetime exhibits a thermally activated exponential behaviour, uncovering the microscopic virtual pair breaking and recombination mechanism. The obtained large, temperature-independent single-particle scattering rate is comparable with that set by the Planckian limit6. Our findings quantitatively characterize the pseudogap in strongly interacting Fermi gases and they lend support for the role of preformed pairing as a precursor to superfluidity.

Interference of holon strings in 2D Hubbard model

The 2D Hubbard model with large repulsion is an important problem in condensed matter physics. At half filling, its ground state is an antiferromagnet (AMF). The dope AMF below half filling is believed to capture the physics of highTcsuperconductors. And the fermion excitation of this dope AMF is theorized as splitting up into holons and spinons that carry charge and spin separately. It is believed that these exotic holons and spinons are the origins of the unusual properties of highTcsuperconductors. Despite the interests in holons and spinons, the direct observations of these excitations remain difficult in solid state experiments. Here, we show that with the rapid advances in the experimental techniques in cold atoms, the direct observation of holons is possible in quantum quench dynamic processes in cold atom settings. We show that the well-known holon-strings generated by the motion of a holon as well as their interferences can be detected by the measurements spin-spin correlations and demonstrate the presence of the Marshall phase associated with a holon string reflecting an underlying AMF background. Moreover, we show that the interferences of the holon strings make a holon propagate anisotropically, with a diffusion pattern clearly distinct from that of spinless fermions. At the same time, we show that these interferences lead to a large suppression in magnetic order in the region swept through by the strings (even to about 95% for some bond). We further demonstrate the Marshall phase of the holon-strings by comparing the dynamics of holon in thetJmodel with that of the so-calledσtJ-model, which is thetJmodel with the Marshall phase removed. The holons in these models propagate entirely differently.

Finite-Temperature Hole-Magnon Dynamics in an Antiferromagnet

Employing the numerically accurate multiple Davydov Ansatz in combination with the thermo-field dynamics approach, we delve into the interplay of the finite-temperature dynamics of holes and magnons in an antiferromagnet, which allows for scrutinizing previous predictions from the self-consistent Born approximation while offering, for the first time, accurate finite-temperature computation of detailed magnon dynamics as a response and a facilitator to the hole motion. The study also uncovers a pronounced temperature dependence of the magnon and hole populations, pointing to the feasibility of potential thermal manipulation and control of hole dynamics. Our methodology can be applied not only to the calculation of steady-state angular-resolved photoemission spectra but also to the simulation of femtosecond terahertz pump-probe and other nonlinear signals for the characterization of antiferromagnetic materials.

Commensurate and incommensurate 1D interacting quantum systems

Single-atom imaging resolution of many-body quantum systems in optical lattices is routinely achieved with quantum-gas microscopes. Key to their great versatility as quantum simulators is the ability to use engineered light potentials at the microscopic level. Here, we employ dynamically varying microscopic light potentials in a quantum-gas microscope to study commensurate and incommensurate 1D systems of interacting bosonic Rb atoms. Such incommensurate systems are analogous to doped insulating states that exhibit atom transport and compressibility. Initially, a commensurate system with unit filling and fixed atom number is prepared between two potential barriers. We deterministically create an incommensurate system by dynamically changing the position of the barriers such that the number of available lattice sites is reduced while retaining the atom number. Our systems are characterised by measuring the distribution of particles and holes as a function of the lattice filling, and interaction strength, and we probe the particle mobility by applying a bias potential. Our work provides the foundation for preparation of low-entropy states with controlled filling in optical-lattice experiments.

Dipolar quantum solids emerging in a Hubbard quantum simulator

In quantum mechanical many-body systems, long-range and anisotropic interactions promote rich spatial structure and can lead to quantum frustration, giving rise to a wealth of complex, strongly correlated quantum phases1. Long-range interactions play an important role in nature; however, quantum simulations of lattice systems have largely not been able to realize such interactions. A wide range of efforts are underway to explore long-range interacting lattice systems using polar molecules2-5, Rydberg atoms2,6-8, optical cavities9-11 or magnetic atoms12-15. Here we realize novel quantum phases in a strongly correlated lattice system with long-range dipolar interactions using ultracold magnetic erbium atoms. As we tune the dipolar interaction to be the dominant energy scale in our system, we observe quantum phase transitions from a superfluid into dipolar quantum solids, which we directly detect using quantum gas microscopy with accordion lattices. Controlling the interaction anisotropy by orienting the dipoles enables us to realize a variety of stripe-ordered states. Furthermore, by transitioning non-adiabatically through the strongly correlated regime, we observe the emergence of a range of metastable stripe-ordered states. This work demonstrates that novel strongly correlated quantum phases can be realized using long-range dipolar interactions in optical lattices, opening the door to quantum simulations of a wide range of lattice models with long-range and anisotropic interactions.

Proposal for Observing Yang-Lee Criticality in Rydberg Atomic Arrays

Yang-Lee edge singularities (YLES) are the edges of the partition function zeros of an interacting spin model in the space of complex control parameters. They play an important role in understanding non-Hermitian phase transitions in many-body physics, as well as characterizing the corresponding nonunitary criticality. Even though such partition function zeroes have been measured in dynamical experiments where time acts as the imaginary control field, experimentally demonstrating such YLES criticality with a physical imaginary field has remained elusive due to the difficulty of physically realizing non-Hermitian many-body models. We provide a protocol for observing the YLES by detecting kinked dynamical magnetization responses due to broken PT symmetry, thus enabling the physical probing of nonunitary phase transitions in nonequilibrium settings. In particular, scaling analyses based on our nonunitary time evolution circuit with matrix product states accurately recover the exponents uniquely associated with the corresponding nonunitary CFT. We provide an explicit proposal for observing YLES criticality in Floquet quenched Rydberg atomic arrays with laser-induced loss, which paves the way towards a universal platform for simulating non-Hermitian many-body dynamical phenomena.

Fermionic Correlation Functions from Randomized Measurements in Programmable Atomic Quantum Devices

We provide an efficient randomized measurement protocol to estimate two- and four-point fermionic correlations in ultracold atom experiments. Our approach is based on combining random atomic beam splitter operations, which can be realized with programmable optical landscapes, with high-resolution imaging systems such as quantum gas microscopes. We illustrate our results in the context of the variational quantum eigensolver algorithm for solving quantum chemistry problems.

Frustration- and doping-induced magnetism in a Fermi-Hubbard simulator

Geometrical frustration in strongly correlated systems can give rise to a plethora of novel ordered states and intriguing magnetic phases, such as quantum spin liquids1-3. Promising candidate materials for such phases4-6 can be described by the Hubbard model on an anisotropic triangular lattice, a paradigmatic model capturing the interplay between strong correlations and magnetic frustration7-11. However, the fate of frustrated magnetism in the presence of itinerant dopants remains unclear, as well as its connection to the doped phases of the square Hubbard model12. Here we investigate the local spin order of a Hubbard model with controllable frustration and doping, using ultracold fermions in anisotropic optical lattices continuously tunable from a square to a triangular geometry. At half-filling and strong interactions U/t ≈ 9, we observe at the single-site level how frustration reduces the range of magnetic correlations and drives a transition from a collinear Néel antiferromagnet to a short-range correlated 120° spiral phase. Away from half-filling, the triangular limit shows enhanced antiferromagnetic correlations on the hole-doped side and a reversal to ferromagnetic correlations at particle dopings above 20%, hinting at the role of kinetic magnetism in frustrated systems. This work paves the way towards exploring possible chiral ordered or superconducting phases in triangular lattices8,13 and realizing t-t' square lattice Hubbard models that may be essential to describe superconductivity in cuprate materials14.

Tangent Space Approach for Thermal Tensor Network Simulations of the 2D Hubbard Model

Accurate simulations of the two-dimensional (2D) Hubbard model constitute one of the most challenging problems in condensed matter and quantum physics. Here we develop a tangent space tensor renormalization group (tanTRG) approach for the calculations of the 2D Hubbard model at finite temperature. An optimal evolution of the density operator is achieved in tanTRG with a mild O(D^{3}) complexity, where the bond dimension D controls the accuracy. With the tanTRG approach we boost the low-temperature calculations of large-scale 2D Hubbard systems on up to a width-8 cylinder and 10×10 square lattice. For the half-filled Hubbard model, the obtained results are in excellent agreement with those of determinant quantum Monte Carlo (DQMC). Moreover, tanTRG can be used to explore the low-temperature, finite-doping regime inaccessible for DQMC. The calculated charge compressibility and Matsubara Green's function are found to reflect the strange metal and pseudogap behaviors, respectively. The superconductive pairing susceptibility is computed down to a low temperature of approximately 1/24 of the hopping energy, where we find d-wave pairing responses are most significant near the optimal doping. Equipped with the tangent-space technique, tanTRG constitutes a well-controlled, highly efficient and accurate tensor network method for strongly correlated 2D lattice models at finite temperature.

Non-orientable order and non-commutative response in frustrated metamaterials

From atomic crystals to animal flocks, the emergence of order in nature is captured by the concept of spontaneous symmetry breaking^1-4. However, this cornerstone of physics is challenged when broken symmetry phases are frustrated by geometrical constraints. Such frustration dictates the behaviour of systems as diverse as spin ices^5-8, confined colloidal suspensions⁹ and crumpled paper sheets¹⁰. These systems typically exhibit strongly degenerated and heterogeneous ground states and hence escape the Ginzburg-Landau paradigm of phase ordering. Here, combining experiments, simulations and theory we uncover an unanticipated form of topological order in globally frustrated matter: non-orientable order. We demonstrate this concept by designing globally frustrated metamaterials that spontaneously break a discrete [Formula: see text] symmetry. We observe that their equilibria are necessarily heteregeneous and extensively degenerated. We explain our observations by generalizing the theory of elasticity to non-orientable order-parameter bundles. We show that non-orientable equilibria are extensively degenerated due to the arbitrary location of topologically protected nodes and lines where the order parameter must vanish. We further show that non-orientable order applies more broadly to objects that are non-orientable themselves, such as buckled Möbius strips and Klein bottles. Finally, by applying time-dependent local perturbations to metamaterials with non-orientable order, we engineer topologically protected mechanical memories^11-19, achieve non-commutative responses and show that they carry an imprint of the braiding of the loads' trajectories. Beyond mechanics, we envision non-orientability as a robust design principle for metamaterials that can effectively store information across scales, in fields as diverse as colloidal science⁸, photonics²⁰, magnetism⁷ and atomic physics²¹.

Continuous symmetry breaking in a two-dimensional Rydberg array

Spontaneous symmetry breaking underlies much of our classification of phases of matter and their associated transitions1-3. The nature of the underlying symmetry being broken determines many of the qualitative properties of the phase; this is illustrated by the case of discrete versus continuous symmetry breaking. Indeed, in contrast to the discrete case, the breaking of a continuous symmetry leads to the emergence of gapless Goldstone modes controlling, for instance, the thermodynamic stability of the ordered phase4,5. Here, we realize a two-dimensional dipolar XY model that shows a continuous spin-rotational symmetry using a programmable Rydberg quantum simulator. We demonstrate the adiabatic preparation of correlated low-temperature states of both the XY ferromagnet and the XY antiferromagnet. In the ferromagnetic case, we characterize the presence of a long-range XY order, a feature prohibited in the absence of long-range dipolar interaction. Our exploration of the many-body physics of XY interactions complements recent works using the Rydberg-blockade mechanism to realize Ising-type interactions showing discrete spin rotation symmetry6-9.

Accurate holographic light potentials using pixel crosstalk modelling

Arbitrary light potentials have proven to be a valuable and versatile tool in many quantum information and quantum simulation experiments with ultracold atoms. Using a phase-modulating spatial light modulator (SLM), we generate arbitrary light potentials holographically with measured efficiencies between 15 and 40% and an accuracy of [Formula: see text] root-mean-squared error. Key to the high accuracy is the modelling of pixel crosstalk of the SLM on a sub-pixel scale which is relevant especially for large light potentials. We employ conjugate gradient minimisation to calculate the SLM phase pattern for a given target light potential after measuring the intensity and wavefront at the SLM. Further, we use camera feedback to reduce experimental errors, we remove optical vortices and investigate the difference between the angular spectrum method and the Fourier transform to simulate the propagation of light. Using a combination of all these techniques, we achieved more accurate and efficient light potentials compared to previous studies, and generated a series of potentials relevant for cold atom experiments.

Nanoscale covariance magnetometry with diamond quantum sensors

Nitrogen vacancy (NV) centers in diamond are atom-scale defects that can be used to sense magnetic fields with high sensitivity and spatial resolution. Typically, the magnetic field is measured by averaging sequential measurements of single NV centers, or by spatial averaging over ensembles of many NV centers, which provides mean values that contain no nonlocal information about the relationship between two points separated in space or time. Here, we propose and implement a sensing modality whereby two or more NV centers are measured simultaneously, and we extract temporal and spatial correlations in their signals that would otherwise be inaccessible. We demonstrate measurements of correlated applied noise using spin-to-charge readout of two NV centers and implement a spectral reconstruction protocol for disentangling local and nonlocal noise sources.

Nonequilibrium Hole Dynamics in Antiferromagnets: Damped Strings and Polarons

We develop a nonperturbative theory for hole dynamics in antiferromagnetic spin lattices, as described by the t-J model. This is achieved by generalizing the self-consistent Born approximation to nonequilibrium systems, making it possible to calculate the full time-dependent many-body wave function. Our approach reveals three distinct dynamical regimes, ultimately leading to the formation of magnetic polarons. Following the initial ballistic stage of the hole dynamics, coherent formation of string excitations gives rise to characteristic oscillations in the hole density. Their damping eventually leaves behind magnetic polarons that undergo ballistic motion with a greatly reduced velocity. The developed theory provides a rigorous framework for understanding nonequilibrium physics of defects in quantum magnets and quantitatively explains recent observations from cold-atom quantum simulations in the strong coupling regime.

Disorder-assisted assembly of strongly correlated fluids of light

Guiding many-body systems to desired states is a central challenge of modern quantum science, with applications from quantum computation1,2 to many-body physics3 and quantum-enhanced metrology4. Approaches to solving this problem include step-by-step assembly5,6, reservoir engineering to irreversibly pump towards a target state7,8 and adiabatic evolution from a known initial state9,10. Here we construct low-entropy quantum fluids of light in a Bose-Hubbard circuit by combining particle-by-particle assembly and adiabatic preparation. We inject individual photons into a disordered lattice for which the eigenstates are known and localized, then adiabatically remove this disorder, enabling quantum fluctuations to melt the photons into a fluid. Using our platform11, we first benchmark this lattice melting technique by building and characterizing arbitrary single-particle-in-a-box states, then assemble multiparticle strongly correlated fluids. Intersite entanglement measurements performed through single-site tomography indicate that the particles in the fluid delocalize, whereas two-body density correlation measurements demonstrate that they also avoid one another, revealing Friedel oscillations characteristic of a Tonks-Girardeau gas12,13. This work opens new possibilities for the preparation of topological and otherwise exotic phases of synthetic matter3,14,15.

Experimental realization of an extended Fermi-Hubbard model using a 2D lattice of dopant-based quantum dots

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.

Atomically precise artificial lattices of dopant-based quantum dots offer a tunable platform for simulations of interacting fermionic models. By leveraging advances in fabrication and atomic-state control, Wang et al. report quantum simulations of the 2D Fermi-Hubbard model on a 3 × 3 few-dopant quantum dot array.

Two-Dimensional Programmable Tweezer Arrays of Fermions

We prepare high-filling two-component arrays of tens of fermionic ^{6}Li atoms in optical tweezers, with the atoms in the ground motional state of each tweezer. Using a stroboscopic technique, we configure the arrays in various two-dimensional geometries with negligible Floquet heating. A full spin- and density-resolved readout of individual sites allows us to postselect near-zero entropy initial states for fermionic quantum simulation. We prepare a correlated state in a two-by-two tunnel-coupled Hubbard plaquette, demonstrating all the building blocks for realizing a programmable fermionic quantum simulator.

Scalable Spin Squeezing from Spontaneous Breaking of a Continuous Symmetry

Spontaneous symmetry breaking is a property of Hamiltonian equilibrium states which, in the thermodynamic limit, retain a finite average value of an order parameter even after a field coupled to it is adiabatically turned off. In the case of quantum spin models with continuous symmetry, we show that this adiabatic process is also accompanied by the suppression of the fluctuations of the symmetry generator-namely, the collective spin component along an axis of symmetry. In systems of S=1/2 spins or qubits, the combination of the suppression of fluctuations along one direction and of the persistence of transverse magnetization leads to spin squeezing-a much sought-after property of quantum states, both for the purpose of entanglement detection as well as for metrological uses. Focusing on the case of XXZ models spontaneously breaking a U(1) [or even SU(2)] symmetry, we show that the adiabatically prepared states have nearly minimal spin uncertainty; that the minimum phase uncertainty that one can achieve with these states scales as N^{-3/4} with the number of spins N; and that this scaling is attained after an adiabatic preparation time scaling linearly with N. Our findings open the door to the adiabatic preparation of strongly spin-squeezed states in a large variety of quantum many-body devices including, e.g., optical-lattice clocks.

Evaluating Second-Order Phase Transitions with Diagrammatic Monte Carlo: Néel Transition in the Doped Three-Dimensional Hubbard Model

Diagrammatic Monte Carlo-the technique for the numerically exact summation of all Feynman diagrams to high orders-offers a unique unbiased probe of continuous phase transitions. Being formulated directly in the thermodynamic limit, the diagrammatic series is bound to diverge and is not resummable at the transition due to the nonanalyticity of physical observables. This enables the detection of the transition with controlled error bars from an analysis of the series coefficients alone, avoiding the challenge of evaluating physical observables near the transition. We demonstrate this technique by the example of the Néel transition in the 3D Hubbard model. At half filling and higher temperatures, the method matches the accuracy of state-of-the-art finite-size techniques, but surpasses it at low temperatures and allows us to map the phase diagram in the doped regime, where finite-size techniques struggle from the fermion sign problem. At low temperatures and sufficient doping, the transition to an incommensurate spin density wave state is observed.

Realization of a Fermi-Hubbard Optical Tweezer Array

We use lithium-6 atoms in an optical tweezer array to realize an eight-site Fermi-Hubbard chain near half filling. We achieve single site detection by combining the tweezer array with a quantum gas microscope. By reducing disorder in the energy offsets to less than the tunneling energy, we observe Mott insulators with strong antiferromagnetic correlations. The measured spin correlations allow us to put an upper bound on the entropy of 0.26(4)k_{B} per atom, comparable to the lowest entropies achieved with optical lattices. Additionally, we establish the flexibility of the tweezer platform by initializing atoms on one tweezer and observing tunneling dynamics across the array for uniform and staggered 1D geometries.

Realizing the symmetry-protected Haldane phase in Fermi-Hubbard ladders

Topology in quantum many-body systems has profoundly changed our understanding of quantum phases of matter. The model that has played an instrumental role in elucidating these effects is the antiferromagnetic spin-1 Haldane chain^1,2. Its ground state is a disordered state, with symmetry-protected fourfold-degenerate edge states due to fractional spin excitations. In the bulk, it is characterized by vanishing two-point spin correlations, gapped excitations and a characteristic non-local order parameter^3,4. More recently it has been understood that the Haldane chain forms a specific example of a more general classification scheme of symmetry-protected topological phases of matter, which is based on ideas connected to quantum information and entanglement^5–7. Here, we realize a finite-temperature version of such a topological Haldane phase with Fermi–Hubbard ladders in an ultracold-atom quantum simulator. We directly reveal both edge and bulk properties of the system through the use of single-site and particle-resolved measurements, as well as non-local correlation functions. Continuously changing the Hubbard interaction strength of the system enables us to investigate the robustness of the phase to charge (density) fluctuations far from the regime of the Heisenberg model, using a novel correlator.

A ladder-like arrangement of an ultracold gas of lithium atoms trapped in an optical lattice enables the observation of a symmetry-protected topological phase.

Phase driven topological states in correlated Haldane model on a honeycomb lattice

Using mean field method and random phase approximation, we studied the phase driven topological exotic states in correlated Haldane model on a honeycomb lattice. It is found that topological spin density waves emerge with the phase change of next-nearest-neighbor hopping. We also investigated the topological properties of these spin density waves, including Chern number, edge state and Hall conductivity. Our work provides a new insight for topological phase transitions in correlated quantum anomalous Hall insulators.

Geometrical picture of the electron-electron correlation at the large-D limit

In electronic structure calculations, the correlation energy is defined as the difference between the mean field and the exact solution of the non relativistic Schrödinger equation. Such an error in the different calculations is not directly observable as there is no simple quantum mechanical operator, apart from correlation functions, that correspond to such quantity. Here, we use the dimensional scaling approach, in which the electrons are localized at the large-dimensional scaled space, to describe a geometric picture of the electronic correlation. Both, the mean field, and the exact solutions at the large-D limit have distinct geometries. Thus, the difference might be used to describe the correlation effect. Moreover, correlations can be also described and quantified by the entanglement between the electrons, which is a strong correlation without a classical analog. Entanglement is directly observable and it is one of the most striking properties of quantum mechanics and bounded by the area law for local gapped Hamiltonians of interacting many-body systems. This study opens up the possibility of presenting a geometrical picture of the electron-electron correlations and might give a bound on the correlation energy. The results at the large-D limit and at D = 3 indicate the feasibility of using the geometrical picture to get a bound on the electron-electron correlations.

Collective Resonance of D States in Rubidium Atoms Probed by Optical Two-Dimensional Coherent Spectroscopy

Collective resonance of interacting particles has important implications in many-body quantum systems and their applications. Strong interactions can lead to a blockade that prohibits the excitation of a collective resonance of two or more nearby atoms. However, a collective resonance can be excited with the presence of weak interaction and has been observed for atoms in the first excited state (P state). Here, we report the observation of collective resonance of rubidium atoms in a higher excited state (D state) in addition to the first excited state. The collective resonance is excited by a double-quantum four-pulse excitation sequence. The resulting double-quantum two-dimensional (2D) spectrum displays well-isolated peaks that can be attributed to collective resonances of atoms in P and D states. The experimental one-quantum and double-quantum 2D spectra can be reproduced by a simulation based on the perturbative solutions to the optical Bloch equations, confirming collective resonances as the origin of the measured spectra. The experimental technique provides a new approach for preparing and probing collective resonances of atoms in highly excited states.

Preparation of the Spin-Mott State: A Spinful Mott Insulator of Repulsively Bound Pairs

We observe and study a special ground state of bosons with two spin states in an optical lattice: the spin-Mott insulator, a state that consists of repulsively bound pairs that is insulating for both spin and charge transport. Because of the pairing gap created by the interaction anisotropy, it can be prepared with low entropy and can serve as a starting point for adiabatic state preparation. We find that the stability of the spin-Mott state depends on the pairing energy, and observe two qualitatively different decay regimes, one of which exhibits protection by the gap.

Quantum simulation of quantum many-body systems with ultracold two-electron atoms in an optical lattice

Ultracold atoms in an optical lattice provide a unique approach to study quantum many-body systems, previously only possible by using condensed-matter experimental systems. This new approach, often called quantum simulation, becomes possible because of the high controllability of the system parameters and the inherent cleanness without lattice defects and impurities. In this article, we review recent developments in this rapidly growing field of ultracold atoms in an optical lattice, with special focus on quantum simulations using our newly created quantum many-body system of two-electron atoms of ytterbium. In addition, we also mention other interesting possibilities offered by this novel experimental platform, such as applications to precision measurements for studying fundamental physics and a Rydberg atom quantum computation.

Phonon-Induced Pairing in Quantum Dot Quantum Simulator

Quantum simulations can provide new insights into the physics of strongly correlated electronic systems. A well-studied system, but still open in many regards, is the Hubbard-Holstein Hamiltonian, where electronic repulsion is in competition with attraction generated by the electron-phonon coupling. In this context, we study the behavior of four quantum dots in a suspended carbon nanotube and coupled to its flexural degrees of freedom. The system is described by a Hamiltonian of the Hubbard-Holstein class, where electrons on different sites interact with the same phonon. We find that the system presents a transition from the Mott insulating state to a polaronic state, with the appearance of pairing correlations and the breaking of the translational symmetry. These findings will motivate further theoretical and experimental efforts to employ nanoelectromechanical systems to simulate strongly correlated systems with electron-phonon interactions.

Rotational Resonances and Regge-like Trajectories in Lightly Doped Antiferromagnets

Understanding the nature of charge carriers in doped Mott insulators holds the key to unravelling puzzling properties of strongly correlated electron systems, including cuprate superconductors. Several theoretical models suggested that dopants can be understood as bound states of partons, the analogues of quarks in high-energy physics. However, direct signatures of spinon-chargon bound states are lacking, both in experiment and theory. Here we propose a rotational variant of angle-resolved photo-emission spectroscopy (ARPES) and calculate rotational spectra numerically using the density-matrix renormalization group. We identify long-lived rotational resonances for an individual dopant, which we interpret as a direct indicator of the microscopic structure of spinon-chargon bound states. Similar to Regge trajectories reflecting the quark structure of mesons, we establish a linear dependence of the rotational energy on the superexchange coupling. The rotational peaks we find are strongly suppressed in standard ARPES spectra, but we suggest a multiphoton extension of ARPES which allows us to access rotational spectra. Our findings suggest that multiphoton spectroscopy experiments should provide new insights into emergent universal features of strongly correlated electron systems.

Critical dynamics and phase transition of a strongly interacting warm spin gas

Phase transitions are emergent phenomena where microscopic interactions drive a disordered system into a collectively ordered phase. Near the boundary between two phases, the system can exhibit critical, scale-invariant behavior. Here, we report on a second-order phase transition accompanied by critical behavior in a system of warm cesium spins driven by linearly polarized light. The ordered phase exhibits macroscopic magnetization when the interactions between the spins become dominant. We measure the phase diagram of the system and observe the collective behavior near the phase boundaries, including power-law dependence of the magnetization and divergence of the susceptibility. Out of equilibrium, we observe a critical slowdown of the spin response time by two orders of magnitude, exceeding 5 s near the phase boundary. This work establishes a controlled platform for investigating equilibrium and nonequilibrium properties of magnetic phases.

Quantum Variational Learning of the Entanglement Hamiltonian

Learning the structure of the entanglement Hamiltonian (EH) is central to characterizing quantum many-body states in analog quantum simulation. We describe a protocol where spatial deformations of the many-body Hamiltonian, physically realized on the quantum device, serve as an efficient variational ansatz for a local EH. Optimal variational parameters are determined in a feedback loop, involving quench dynamics with the deformed Hamiltonian as a quantum processing step, and classical optimization. We simulate the protocol for the ground state of Fermi-Hubbard models in quasi-1D geometries, finding excellent agreement of the EH with Bisognano-Wichmann predictions. Subsequent on-device spectroscopy enables a direct measurement of the entanglement spectrum, which we illustrate for a Fermi Hubbard model in a topological phase.