Quantum simulations with ultracold atoms in optical lattices

Origin and fate of the pseudogap in the doped Hubbard model

The relationship between the pseudogap and underlying ground-state phases has not yet been rigorously established. We investigated the doped two-dimensional Hubbard model at finite temperature using controlled diagrammatic Monte Carlo calculations, allowing for the computation of spectral properties in the infinite-size limit and with arbitrary momentum resolution. We found three distinct regimes as a function of doping and interaction strength: a weakly correlated metal, a correlated metal with strong interaction effects, and a pseudogap regime at low doping. We show that the pseudogap forms both at weak coupling, when the magnetic correlation length is large, and at strong coupling, when it is shorter. As the temperature goes to zero, the pseudogap regime extrapolates precisely to the ordered stripe phase found by ground-state methods.

Probing spin hydrodynamics on a superconducting quantum simulator

Characterizing the nature of hydrodynamical transport properties in quantum dynamics provides valuable insights into the fundamental understanding of exotic non-equilibrium phases of matter. Experimentally simulating infinite-temperature transport on large-scale complex quantum systems is of considerable interest. Here, using a controllable and coherent superconducting quantum simulator, we experimentally realize the analog quantum circuit, which can efficiently prepare the Haar-random states, and probe spin transport at infinite temperature. We observe diffusive spin transport during the unitary evolution of the ladder-type quantum simulator with ergodic dynamics. Moreover, we explore the transport properties of the systems subjected to strong disorder or a tilted potential, revealing signatures of anomalous subdiffusion in accompany with the breakdown of thermalization. Our work demonstrates a scalable method of probing infinite-temperature spin transport on analog quantum simulators, which paves the way to study other intriguing out-of-equilibrium phenomena from the perspective of transport.

Strong Spin-Motion Coupling in the Ultrafast Dynamics of Rydberg Atoms

Rydberg atoms in optical lattices and tweezers is now a well-established platform for simulating quantum spin systems. However, the role of the atoms' spatial wave function has not been examined in detail experimentally. Here, we show a strong spin-motion coupling emerging from the large variation of the interaction potential over the wave function spread. We observe its clear signature on the ultrafast many-body nanosecond-dynamics of atoms excited to a Rydberg S state, using picosecond pulses, from an unity-filling atomic Mott-insulator. We also propose an approach to tune arbitrarily the strength of the spin-motion coupling relative to the motional energy scale set by trapping potentials. Our work provides a new direction for exploring the dynamics of strongly correlated quantum systems by adding the motional degree of freedom to the Rydberg simulation toolbox.

Quantum neural network approach to Markovian dissipative dynamics of many-body open quantum systems

Numerous variational methods have been proposed for solving quantum many-body systems, but they often face exponentially increasing computational complexity as the Hilbert space dimension grows. To address this, we introduce a novel approach using quantum neural networks to simulate the dissipative dynamics of many-body open quantum systems. This method combines neural-network quantum state representation with the time-dependent variational principle, both implemented via quantum algorithms. This results in accurate open quantum dynamics described by the Lindblad quantum master equation, exemplified by the spin-boson and transverse field Ising models. Our approach avoids the computational expense of classical algorithms and demonstrates the potential advantages of quantum computing for many-body simulations. To reduce measurement errors, we introduce a projection reset procedure, which could benefit other quantum simulations. In addition, our approach can be extended to simulate non-Markovian quantum dynamics.

Scalable Spin Squeezing in Two-Dimensional Arrays of Dipolar Large-S Spins

We theoretically show that the spin-spin interactions realized in two-dimensional Mott insulators of large-spin magnetic atoms (such as Cr, Er, or Dy) lead to scalable spin squeezing along the nonequilibrium unitary evolution initialized in a coherent spin state. An experimentally relevant perturbation to the collective squeezing dynamics is offered by a quadratic Zeeman shift, which leads instead to squeezing of individual spins. Making use of a truncated cumulant expansion for the quantum fluctuations of the spin array, we show that, for sufficiently small quadratic shifts, the spin squeezing dynamics is akin to that produced by the paradigmatic one-axis-twisting model-as expected from an effective separation between collective-spin and spin-wave variables. Scalable spin squeezing is shown to be protected by the robustness of long-range ferromagnetic order to quadratic shifts in the equilibrium phase diagram of the system that we reconstruct via quantum Monte Carlo and mean-field theory.

Localized Virtual Purification

Analog and digital quantum simulators can efficiently simulate quantum many-body systems that appear in natural phenomena. However, experimental limitations of near-term devices still make it challenging to perform the entire process of quantum simulation. The purification-based quantum simulation methods can alleviate the limitations in experiments such as the cooling temperature and noise from the environment, while this method has the drawback that it requires global entangled measurement with a prohibitively large number of measurements that scales exponentially with the system size. In this Letter, we propose that we can overcome these problems by restricting the entangled measurements to the vicinity of the local observables to be measured, when the locality of the system can be exploited. We provide theoretical guarantees that the global purification operation can be replaced with local operations under some conditions, in particular for the task of cooling and error mitigation. We furthermore give a numerical verification that the localized purification is valid even when conditions are not satisfied. Our method bridges the fundamental concept of locality with quantum simulators, and therefore is expected to open a path to unexplored quantum many-body phenomena.

Canonical Typicality under General Quantum Channels

With the control of increasingly complex quantum systems, the relevant degrees of freedom we are interested in may not be those traditionally addressed by statistical quantum mechanics. Here, we employ quantum channels to define generalized subsystems, capturing the pertinent degrees of freedom, and obtain their associated canonical state. We show that the generalized subsystem description from almost any microscopic pure state of the whole system will behave similarly to its corresponding canonical state. Such canonical typicality behavior depends on the entropy of the channel used to define the generalized subsystem.

Local Readout and Control of Current and Kinetic Energy Operators in Optical Lattices

Quantum gas microscopes have revolutionized quantum simulations with ultracold atoms, allowing one to measure local observables and snapshots of quantum states. However, measurements so far were mostly carried out in the occupation basis. Here, we demonstrate how all kinetic operators, such as kinetic energy or current operators, can be measured and manipulated with single-bond resolution. Beyond simple expectation values of these observables, the single-shot measurements allow one to access full counting statistics and complex correlation functions. Our work paves the way for the implementation of efficient quantum state tomography and hybrid quantum computing protocols for itinerant particles on a lattice. In addition, we demonstrate how site-resolved programmable potentials enable a spatially selective, parallel readout in different bases as well as the engineering of arbitrary initial states.

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.

Microscopic Study on Superexchange Dynamics of Composite Spin-1 Bosons

We report on an experimental simulation of the spin-1 Heisenberg model with composite bosons in a one-dimensional chain based on the two-component Bose-Hubbard model. Exploiting our site- and spin-resolved quantum gas microscope, we observed faster superexchange dynamics of the spin-1 system compared to its spin-1/2 counterpart, which is attributed to the enhancement effect of multi-bosons. We further probed the nonequilibrium spin dynamics driven by the superexchange and single-ion anisotropy terms, unveiling the linear expansion of the spin-spin correlations, which is limited by the Lieb-Robinson bound. Based on the superexchange process, we prepared and verified the entangled qutrits pairs with these composite spin-1 bosons, potentially being applied in qutrit-based quantum information processing.

Ising meson spectroscopy on a noisy digital quantum simulator

Quantum simulation has the potential to be an indispensable technique for the investigation of non-perturbative phenomena in strongly-interacting quantum field theories (QFTs). In the modern quantum era, with Noisy Intermediate Scale Quantum (NISQ) simulators widely available and larger-scale quantum machines on the horizon, it is natural to ask: what non-perturbative QFT problems can be solved with the existing quantum hardware? We show that existing noisy quantum machines can be used to analyze the energy spectrum of several strongly-interacting 1+1D QFTs, which exhibit non-perturbative effects like 'quark confinement' and 'false vacuum decay'. We perform quench experiments on IBM's quantum simulators to compute the energy spectrum of 1+1D quantum Ising model with a longitudinal field. Our results demonstrate that digital quantum simulation in the NISQ era has the potential to be a viable alternative to numerical techniques such as density matrix renormalization group or the truncated conformal space methods for analyzing QFTs.

Superconducting Quantum Simulation for Many-Body Physics beyond Equilibrium

Quantum computing is an exciting field that uses quantum principles, such as quantum superposition and entanglement, to tackle complex computational problems. Superconducting quantum circuits, based on Josephson junctions, is one of the most promising physical realizations to achieve the long-term goal of building fault-tolerant quantum computers. The past decade has witnessed the rapid development of this field, where many intermediate-scale multi-qubit experiments emerged to simulate nonequilibrium quantum many-body dynamics that are challenging for classical computers. Here, we review the basic concepts of superconducting quantum simulation and their recent experimental progress in exploring exotic nonequilibrium quantum phenomena emerging in strongly interacting many-body systems, e.g., many-body localization, quantum many-body scars, and discrete time crystals. We further discuss the prospects of quantum simulation experiments to truly solve open problems in nonequilibrium many-body systems.

High-Fidelity Detection of Large-Scale Atom Arrays in an Optical Lattice

Recent advances in quantum simulation based on neutral atoms have largely benefited from high-resolution, single-atom sensitive imaging techniques. A variety of approaches have been developed to achieve such local detection of atoms in optical lattices or optical tweezers. For alkaline-earth and alkaline-earth-like atoms, the presence of narrow optical transitions opens up the possibility of performing novel types of Sisyphus cooling, where the cooling mechanism originates from the capability to spatially resolve the differential optical level shifts in the trap potential. Up to now, it has been an open question whether high-fidelity imaging could be achieved in a "repulsive Sisyphus" configuration, where the trap depth of the ground state exceeds that of the excited state involved in cooling. Here, we demonstrate high-fidelity (99.971(1)%) and high-survival (99.80(5)%) imaging of strontium atoms using repulsive Sisyphus cooling. We use an optical lattice as a pinning potential for atoms in a large-scale tweezer array with up to 399 tweezers and show repeated, high-fidelity lattice-tweezer-lattice transfers. We furthermore demonstrate loading the lattice with approximately 10 000 atoms directly from the MOT and scalable imaging over >10 000 lattice sites with a combined survival probability and classification fidelity better than 99.2%. Our lattice thus serves as a locally addressable and sortable reservoir for continuous refilling of optical tweezer arrays in the future.

Observation of Bose-Einstein condensation of dipolar molecules

Ensembles of particles governed by quantum mechanical laws exhibit intriguing emergent behaviour. Atomic quantum gases1,2, liquid helium3,4 and electrons in quantum materials5-7 all exhibit distinct properties because of their composition and interactions. Quantum degenerate samples of ultracold dipolar molecules promise the realization of new phases of matter and new avenues for quantum simulation8 and quantum computation9. However, rapid losses10, even when reduced through collisional shielding techniques11-13, have so far prevented evaporative cooling to a Bose-Einstein condensate (BEC). Here we report on the realization of a BEC of dipolar molecules. By strongly suppressing two- and three-body losses via enhanced collisional shielding, we evaporatively cool sodium-caesium molecules to quantum degeneracy and cross the phase transition to a BEC. The BEC reveals itself by a bimodal distribution when the phase-space density exceeds 1. BECs with a condensate fraction of 60(10)% and a temperature of 6(2) nK are created and found to be stable with a lifetime close to 2 s. This work opens the door to the exploration of dipolar quantum matter in regimes that have been inaccessible so far, promising the creation of exotic dipolar droplets14, self-organized crystal phases15 and dipolar spin liquids in optical lattices16.

Enhancing a Many-Body Dipolar Rydberg Tweezer Array with Arbitrary Local Controls

We implement and characterize a protocol that enables arbitrary local controls in a dipolar atom array, where the degree of freedom is encoded in a pair of Rydberg states. Our approach relies on a combination of local addressing beams and global microwave fields. Using this method, we directly prepare two different types of three-atom entangled states, including a W state and a state exhibiting finite chirality. We verify the nature of the underlying entanglement by performing quantum state tomography. Finally, leveraging our ability to measure multibasis, multibody observables, we explore the adiabatic preparation of low-energy states in a frustrated geometry consisting of a pair of triangular plaquettes. By using local addressing to tune the symmetry of the initial state, we demonstrate the ability to prepare correlated states distinguished only by correlations of their chirality (a fundamentally six-body observable). Our protocol is generic, allowing for rotations on arbitrary sub-groups of atoms within the array at arbitrary times during the experiment; this extends the scope of capabilities for quantum simulations of the dipolar XY model.

Long-Lived Coherence on a μHz Scale Optical Magnetic Quadrupole Transition

We report on the coherent excitation of the ultranarrow ^{1}S_{0}-^{3}P_{2} magnetic quadrupole transition in ^{88}Sr. By confining atoms in a state insensitive optical lattice, we achieve excitation fractions of 97(1)% and observe linewidths as narrow as 58(1) Hz. With Ramsey spectroscopy, we find coherence times of 14(1) ms, which can be extended to 266(36) ms using a spin-echo sequence. We determine the lifetime of the ^{3}P_{2} level for spontaneous emission of magnetic quadrupole radiation to be 110(31) min, confirming long-standing theoretical predictions. These results establish an additional clock transition in strontium and pave the way for applications of the metastable ^{3}P_{2} state in quantum computing and quantum simulations.

Nonmesonic Quantum Many-Body Scars in a 1D Lattice Gauge Theory

We investigate the meson excitations (particle-antiparticle bound states) in quantum many-body scars of a 1D Z_{2} lattice gauge theory coupled to a dynamical spin-1/2 chain as a matter field. By introducing a string representation of the physical Hilbert space, we express a scar state |Ψ_{n,l}⟩ as a superposition of all string bases with an identical string number n and a total length l. For the small-l scar state |Ψ_{n,l}⟩, the gauge-invariant spin exchange correlation function of the matter field hosts an exponential decay as the distance increases, indicating the existence of stable mesons. However, for large l, the correlation function exhibits a power-law decay, signaling the emergence of nonmesonic excitations. Furthermore, we show that this mesonic-nonmesonic crossover can be detected by the quench dynamics, starting from two low-entangled initial states, respectively, which are experimentally feasible in quantum simulators. Our results expand the physics of quantum many-body scars in lattice gauge theories and reveal that the nonmesonic state can also manifest ergodicity 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.

Nonbosonic Moiré Excitons

Optical excitations in moiré transition metal dichalcogenide bilayers lead to the creation of excitons, as electron-hole bound states, that are generically considered within a Bose-Hubbard framework. Here, we demonstrate that these composite particles obey an angular momentum commutation relation that is generally nonbosonic. This emergent spin description of excitons indicates a limitation to their occupancy on each site, which is substantial in the weak electron-hole binding regime. The effective exciton theory is accordingly a spin Hamiltonian, which further becomes a Hubbard model of emergent bosons subject to an occupancy constraint after a Holstein-Primakoff transformation. We apply our theory to three commonly studied bilayers (MoSe_{2}/WSe_{2}, WSe_{2}/WS_{2}, and WSe_{2}/MoS_{2}) and show that in the relevant parameter regimes their allowed occupancies never exceed three excitons. Our systematic theory provides guidelines for future research on the many-body physics of moiré excitons.

Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics

Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.

Monte Carlo Based Techniques for Quantum Magnets with Long-Range Interactions

Long-range interactions are relevant for a large variety of quantum systems in quantum optics and condensed matter physics. In particular, the control of quantum-optical platforms promises to gain deep insights into quantum-critical properties induced by the long-range nature of interactions. From a theoretical perspective, long-range interactions are notoriously complicated to treat. Here, we give an overview of recent advancements to investigate quantum magnets with long-range interactions focusing on two techniques based on Monte Carlo integration. First, the method of perturbative continuous unitary transformations where classical Monte Carlo integration is applied within the embedding scheme of white graphs. This linked-cluster expansion allows extracting high-order series expansions of energies and observables in the thermodynamic limit. Second, stochastic series expansion quantum Monte Carlo integration enables calculations on large finite systems. Finite-size scaling can then be used to determine the physical properties of the infinite system. In recent years, both techniques have been applied successfully to one- and two-dimensional quantum magnets involving long-range Ising, XY, and Heisenberg interactions on various bipartite and non-bipartite lattices. Here, we summarise the obtained quantum-critical properties including critical exponents for all these systems in a coherent way. Further, we review how long-range interactions are used to study quantum phase transitions above the upper critical dimension and the scaling techniques to extract these quantum critical properties from the numerical calculations.

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.

Fine-Structure Qubit Encoded in Metastable Strontium Trapped in an Optical Lattice

We demonstrate coherent control of the fine-structure qubit in neutral strontium atoms. This qubit is encoded in the metastable ^{3}P_{2} and ^{3}P_{0} states, coupled by a Raman transition. Using a magnetic quadrupole transition, we demonstrate coherent state initialization of this THz qubit. We show Rabi oscillations with more than 60 coherent cycles and single-qubit rotations on the μs scale. With spin echo, we demonstrate coherence times of tens of ms. Our results pave the way for fast quantum information processors and highly tunable quantum simulators with two-electron atoms.

Frustrated Extended Bose-Hubbard Model and Deconfined Quantum Critical Points with Optical Lattices at the Antimagic Wavelength

The study of geometrically frustrated many-body quantum systems is of central importance to uncover novel quantum mechanical effects. We design a scheme where ultracold bosons trapped in a one-dimensional state-dependent optical lattice are modeled by a frustrated Bose-Hubbard Hamiltonian. A derivation of the Hamiltonian parameters based on Cesium atoms, further show large tunability of contact and nearest-neighbor interactions. For pure contact repulsion, we discover the presence of two phases peculiar to frustrated quantum magnets: the bond-order-wave insulator with broken inversion symmetry and a chiral superfluid. When the nearest-neighbor repulsion becomes sizable, a further density-wave insulator with broken translational symmetry can appear. We show that the phase transition between the two spontaneously symmetry-broken phases is continuous, thus representing a one-dimensional deconfined quantum critical point not captured by the Landau-Ginzburg-Wilson symmetry-breaking paradigm. Our results provide a solid ground to unveil the novel quantum physics induced by the interplay of nonlocal interactions, geometrical frustration, and quantum fluctuations.

Spontaneously Sliding Multipole Spin Density Waves in Cold Atoms

We report on the observation of spontaneously drifting coupled spin and quadrupolar density waves in the ground state of laser driven Rubidium atoms. These laser-cooled atomic ensembles exhibit spontaneous magnetism via light mediated interactions when submitted to optical feedback by a retroreflecting mirror. Drift direction and chirality of the waves arise from spontaneous symmetry breaking. The observations demonstrate a novel transport process in out-of-equilibrium magnetic systems.

Benchmarking highly entangled states on a 60-atom analogue quantum simulator

Quantum systems have entered a competitive regime in which classical computers must make approximations to represent highly entangled quantum states1,2. However, in this beyond-classically-exact regime, fidelity comparisons between quantum and classical systems have so far been limited to digital quantum devices2-5, and it remains unsolved how to estimate the actual entanglement content of experiments6. Here, we perform fidelity benchmarking and mixed-state entanglement estimation with a 60-atom analogue Rydberg quantum simulator, reaching a high-entanglement entropy regime in which exact classical simulation becomes impractical. Our benchmarking protocol involves extrapolation from comparisons against an approximate classical algorithm, introduced here, with varying entanglement limits. We then develop and demonstrate an estimator of the experimental mixed-state entanglement6, finding our experiment is competitive with state-of-the-art digital quantum devices performing random circuit evolution2-5. Finally, we compare the experimental fidelity against that achieved by various approximate classical algorithms, and find that only the algorithm we introduce is able to keep pace with the experiment on the classical hardware we use. Our results enable a new model for evaluating the ability of both analogue and digital quantum devices to generate entanglement in the beyond-classically-exact regime, and highlight the evolving divide between quantum and classical systems.

Quantum control of exciton wave functions in 2D semiconductors

Excitons-bound electron-hole pairs-play a central role in light-matter interaction phenomena and are crucial for wide-ranging applications from light harvesting and generation to quantum information processing. A long-standing challenge in solid-state optics has been to achieve precise and scalable control over excitonic motion. We present a technique using nanostructured gate electrodes to create tailored potential landscapes for excitons in 2D semiconductors, enabling in situ wave function shaping at the nanoscale. Our approach forms electrostatic traps for excitons in various geometries, such as quantum dots, rings, and arrays thereof. We show independent spectral tuning of spatially separated quantum dots, achieving degeneracy despite material disorder. Owing to the strong light-matter coupling of excitons in 2D semiconductors, we observe unambiguous signatures of confined exciton wave functions in optical reflection and photoluminescence measurements. This work unlocks possibilities for engineering exciton dynamics and interactions at the nanometer scale, with implications for optoelectronic devices, topological photonics, and quantum nonlinear optics.

Effective light cone and digital quantum simulation of interacting bosons

The speed limit of information propagation is one of the most fundamental features in non-equilibrium physics. The region of information propagation by finite-time dynamics is approximately restricted inside the effective light cone that is formulated by the Lieb-Robinson bound. To date, extensive studies have been conducted to identify the shape of effective light cones in most experimentally relevant many-body systems. However, the Lieb-Robinson bound in the interacting boson systems, one of the most ubiquitous quantum systems in nature, has remained a critical open problem for a long time. This study reveals a tight effective light cone to limit the information propagation in interacting bosons, where the shape of the effective light cone depends on the spatial dimension. To achieve it, we prove that the speed for bosons to clump together is finite, which in turn leads to the error guarantee of the boson number truncation at each site. Furthermore, we applied the method to provide a provably efficient algorithm for simulating the interacting boson systems. The results of this study settle the notoriously challenging problem and provide the foundation for elucidating the complexity of many-body boson systems.

Designer Spin Models in Tunable Two-Dimensional Nanographene Lattices

Motivated by recent experimental breakthroughs, we propose a strategy for designing two-dimensional spin-lattices with competing interactions that lead to nontrivial emergent quantum states. We consider S = 1/2 nanographenes with C3 symmetry as building blocks, and we leverage the potential to control both the sign and the strength of exchange with first neighbors to build a family of spin models. Specifically, we consider the case of a Heisenberg model in a triangle-decorated honeycomb lattice with competing ferromagnetic and antiferromagnetic interactions whose ratio can be varied in a wide range. On the basis of the exact diagonalization of both Fermionic and spin models, we predict a quantum phase transition between a valence bond crystal of spin singlets with triplon excitations living in a Kagomé lattice and a Néel phase of effective S = 3/2 in the limit of dominant ferromagnetic interactions.

Quantum many-body simulations on digital quantum computers: State-of-the-art and future challenges

Simulating quantum many-body systems is a key application for emerging quantum processors. While analog quantum simulation has already demonstrated quantum advantage, its digital counterpart has recently become the focus of intense research interest due to the availability of devices that aim to realize general-purpose quantum computers. In this perspective, we give a selective overview of the currently pursued approaches, review the advances in digital quantum simulation by comparing non-variational with variational approaches and identify hardware and algorithmic challenges. Based on this review, the question arises: What are the most promising problems that can be tackled with digital quantum simulation? We argue that problems of a qualitative nature are much more suitable for near-term devices then approaches aiming purely for a quantitative accuracy improvement.

Chaos and anomalous transport in a semiclassical Bose-Hubbard chain

We study chaotic dynamics and anomalous transport in a Bose-Hubbard chain in the semiclassical regime (the limit when the number of particles goes to infinity). We find that the system has mixed phase space with both regular and chaotic dynamics, even for long chains with up to 100 wells. The consequence of the mixed phase space is strongly anomalous diffusion in the space of occupation numbers, with a discrete set of transport exponents. After very long times the system crosses over to the hydrodynamic regime with normal diffusion. Anomalous transport is quite universal and almost completely independent of the parameters of the model (Coulomb interaction and chemical potential): It is mainly determined by the initial distribution of particles along the chain. We corroborate our findings by analytical arguments: scaling analysis for the anomalous regime and the Langevin equation for the normal diffusion regime.

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.

Thermography of the superfluid transition in a strongly interacting Fermi gas

Heat transport can serve as a fingerprint identifying different states of matter. In a normal liquid, a hotspot diffuses, whereas in a superfluid, heat propagates as a wave called "second sound." Direct imaging of heat transport is challenging, and one usually resorts to detecting secondary effects. In this study, we establish thermography of a strongly interacting atomic Fermi gas, whose radio-frequency spectrum provides spatially resolved thermometry with subnanokelvin resolution. The superfluid phase transition was directly observed as the sudden change from thermal diffusion to second-sound propagation and is accompanied by a peak in the second-sound diffusivity. This method yields the full heat and density response of the strongly interacting Fermi gas and therefore all defining properties of Landau's two-fluid hydrodynamics.

Dual Higgs modes entangled into a soliton lattice in CuTe

Recently discovered Higgs particle is a key element in the standard model of elementary particles and its analogue in materials, massive Higgs mode, has elucidated intriguing collective phenomena in a wide range of materials with spontaneous symmetry breaking such as antiferromagnets, cold atoms, superconductors, superfluids, and charge density waves (CDW). As a straightforward extension beyond the standard model, multiple Higgs particles have been considered theoretically but not yet for Higgs modes. Here, we report the real-space observations, which suggest two Higgs modes coupled together with a soliton lattice in a solid. Our scanning tunneling microscopy reveals the 1D CDW state of an anisotropic transition metal monochalcogenide crystal CuTe is composed of two distinct but degenerate CDW structures by the layer inversion symmetry broken. More importantly, the amplitudes of each CDW structure oscillate in an out-of-phase fashion to result in a regular array of alternating domains with repeating phase-shift domain walls. This unusual finding is explained by the extra degeneracy in CDWs within the standard Landau theory of the free energy. The multiple and entangled Higgs modes demonstrate how novel collective modes can emerge in systems with distinct symmetries broken simultaneously.

Spin-Triplet Superconductivity from Quantum-Geometry-Induced Ferromagnetic Fluctuation

We show that quantum geometry induces ferromagnetic fluctuation resulting in spin-triplet superconductivity. The criterion for ferromagnetic fluctuation is clarified by analyzing contributions from the effective mass and quantum geometry. When the non-Kramers band degeneracy is present near the Fermi surface, the Fubini-Study quantum metric strongly favors ferromagnetic fluctuation. Solving the linearized gap equation with the effective interaction obtained by the random phase approximation, we show that the spin-triplet superconductivity is mediated by quantum-geometry-induced ferromagnetic fluctuation.

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.

Higher-Order Topological Peierls Insulator in a Two-Dimensional Atom-Cavity System

In this work, we investigate a two-dimensional system of ultracold bosonic atoms inside an optical cavity, and show how photon-mediated interactions give rise to a plaquette-ordered bond pattern in the atomic ground state. The latter corresponds to a 2D Peierls transition, generalizing the spontaneous bond dimerization driven by phonon-electron interactions in the 1D Su-Schrieffer-Heeger (SSH) model. Here the bosonic nature of the atoms plays a crucial role to generate the phase, as similar generalizations with fermionic matter do not lead to a plaquette structure. Similar to the SSH model, we show how this pattern opens a nontrivial topological gap in 2D, resulting in a higher-order topological phase hosting corner states, that we characterize by means of a many-body topological invariant and through its entanglement structure. Finally, we demonstrate how this higher-order topological Peierls insulator can be readily prepared in atomic experiments through adiabatic protocols. Our work thus shows how atomic quantum simulators can be harnessed to investigate novel strongly correlated topological phenomena beyond those observed in natural materials.

Evidence for Quantum Stripe Ordering in a Triangular Optical Lattice

Understanding strongly correlated quantum materials, such as high-T_{c} superconductors, iron-based superconductors, and twisted bilayer graphene systems, remains as one of the outstanding challenges in condensed matter physics. Quantum simulation with ultracold atoms in particular optical lattices, which provide orbital degrees of freedom, is a powerful tool to contribute new insights to this endeavor. Here, we report the experimental realization of an unconventional Bose-Einstein condensate of ^{87}Rb atoms populating degenerate p orbitals in a triangular optical lattice, exhibiting remarkably long coherence times. Using time-of-flight spectroscopy, we observe that this state spontaneously breaks the rotational symmetry and its momentum spectrum agrees with the theoretically predicted coexistence of exotic stripe and loop-current orders. Like certain strongly correlated electronic systems with intertwined orders, such as high-T_{c} cuprate superconductors, twisted bilayer graphene, and the recently discovered chiral density-wave state in kagome superconductors AV_{3}Sb_{5} (A=K, Rb, Cs), the newly demonstrated quantum state, in spite of its markedly different energy scale and the bosonic quantum statistics, exhibits multiple symmetry breakings at ultralow temperatures. These findings hold the potential to enhance our comprehension of the fundamental physics governing these intricate quantum materials.

Quantum Phases from Competing Van der Waals and Dipole-Dipole Interactions of Rydberg Atoms

Competing short- and long-range interactions represent distinguished ingredients for the formation of complex quantum many-body phases. Their study is hard to realize with conventional quantum simulators. In this regard, Rydberg atoms provide an exception as their excited manifold of states have both density-density and exchange interactions whose strength and range can vary considerably. Focusing on one-dimensional systems, we leverage the Van der Waals and dipole-dipole interactions of the Rydberg atoms to obtain the zero-temperature phase diagram for a uniform chain and a dimer model. For the uniform chain, we can influence the boundaries between ordered phases and a Luttinger liquid phase. For the dimerized case, a new type of bond-order-density-wave phase is identified. This demonstrates the versatility of the Rydberg platform in studying physics involving short- and long-ranged interactions simultaneously.

Three-dimensional non-approximate Coulomb interaction between two trapped quantum particles

The two trapped quantum particles interacting problem is generalized to three dimensions, and the exact Coulomb potential is used. The system is solved by expanding the wavefunction in terms of the isotropic harmonic oscillator eigenfunctions and Hydrogen atom eigenfunctions independently, showing that each one results in a prime approximation for different domains of the normalized coupling constant of the relative interactions, suggesting that the combination of the basis is enough to build a well-suited base for the non-approximate problem. The results are compared to previous works that use a model of approximate finite-rage soft-core interaction model of the problem to give insights into the many-body states of strongly correlated ultracold bosons and fermions. We conclude that the proposed three-dimensional approach facilitates the distinction between bosons and fermions while the solutions given by the expansions better define the behavior of the particles for repulsive potentials. In addition, we discuss the substantial differences between our work and the previous approximate model.

Entangling Dynamics from Effective Rotor-Spin-Wave Separation in U(1)-Symmetric Quantum Spin Models

The nonequilibrium dynamics of quantum spin models is a most challenging topic, due to the exponentiality of Hilbert space, and it is central to the understanding of the many-body entangled states that can be generated by state-of-the-art quantum simulators. A particularly important class of evolutions is the one governed by U(1)-symmetric Hamiltonians, initialized in a state that breaks the U(1) symmetry-the paradigmatic example being the evolution of the so-called one-axis-twisting (OAT) model, featuring infinite-range interactions between spins. In this Letter, we show that the dynamics of the OAT model can be closely reproduced by systems with power-law-decaying interactions, thanks to an effective separation between the zero-momentum degrees of freedom, associated with the so-called Anderson tower of states, and reconstructing an OAT model, as well as finite-momentum ones, associated with spin-wave excitations. This mechanism explains quantitatively the recent numerical observation of spin squeezing and Schrödinger cat-state generation in the dynamics of dipolar Hamiltonians, and it paves the way for the extension of this observation to a much larger class of models of immediate relevance for quantum simulations.

Shadow Tomography from Emergent State Designs in Analog Quantum Simulators

We introduce a method that allows one to infer many properties of a quantum state-including nonlinear functions such as Rényi entropies-using only global control over the constituent degrees of freedom. In this protocol, the state of interest is first entangled with a set of ancillas under a fixed global unitary, before projective measurements are made. We show that when the unitary is sufficiently entangling, a universal relationship between the statistics of the measurement outcomes and properties of the state emerges, which can be connected to the recently discovered phenomeonon of emergent quantum state designs in chaotic systems. Thanks to this relationship, arbitrary observables can be reconstructed using the same number of experimental repetitions that would be required in classical shadow tomography [Huang et al., Nat. Phys. 16, 1050 (2020)NPAHAX1745-247310.1038/s41567-020-0932-7]. Unlike previous approaches to shadow tomography, our protocol can be implemented using only global Hamiltonian evolution, as opposed to qubit-selective logic gates, which makes it particularly well suited to analog quantum simulators, including ultracold atoms in optical lattices and arrays of Rydberg atoms.

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.

Complexity of spin configuration dynamics due to unitary evolution and periodic projective measurements

We study the Hamiltonian dynamics of a many-body quantum system subjected to periodic projective measurements, which leads to probabilistic cellular automata dynamics. Given a sequence of measured values, we characterize their dynamics by performing a principal component analysis (PCA). The number of principal components required for an almost complete description of the system, which is a measure of complexity we refer to as PCA complexity, is studied as a function of the Hamiltonian parameters and measurement intervals. We consider different Hamiltonians that describe interacting, noninteracting, integrable, and nonintegrable systems, including random local Hamiltonians and translational invariant random local Hamiltonians. In all these scenarios, we find that the PCA complexity grows rapidly in time before approaching a plateau. The dynamics of the PCA complexity can vary quantitatively and qualitatively as a function of the Hamiltonian parameters and measurement protocol. Importantly, the dynamics of PCA complexity present behavior that is considerably less sensitive to the specific system parameters for models which lack simple local dynamics, as is often the case in nonintegrable models. In particular, we point out a figure of merit that considers the local dynamics and the measurement direction to predict the sensitivity of the PCA complexity dynamics to the system parameters.

Benchmarking Quantum Simulators Using Ergodic Quantum Dynamics

We propose and analyze a sample-efficient protocol to estimate the fidelity between an experimentally prepared state and an ideal target state, applicable to a wide class of analog quantum simulators without advanced spatiotemporal control. Our protocol relies on universal fluctuations emerging from generic Hamiltonian dynamics, which we discover in the present work. It does not require fine-tuned control over state preparation, quantum evolution, or readout capability, while achieving near optimal sample complexity: a percent-level precision is obtained with ∼10^{3} measurements, independent of system size. Furthermore, the accuracy of our fidelity estimation improves exponentially with increasing system size. We numerically demonstrate our protocol in a variety of quantum simulator platforms, including quantum gas microscopes, trapped ions, and Rydberg atom arrays. We discuss applications of our method for tasks such as multiparameter estimation of quantum states and processes.

Meterscale Strong Coupling between Magnons and Photons

We experimentally realize a meterscale strong coupling effect between magnons and photons at room temperature, with a coherent coupling of ∼20  m and a dissipative coupling of ∼7.6  m. To this end, we integrate a saturable gain into a microwave cavity and then couple this active cavity to a magnon mode via a long coaxial cable. The gain compensates for the cavity dissipation, but preserves the cavity radiation that mediates the indirect photon-magnon coupling. It thus enables the long-range strong photon-magnon coupling. With full access to traveling waves, we demonstrate a remote control of photon-magnon coupling by modulating the phase and amplitude of traveling waves, rather than reconfiguring subsystems themselves. Our method for realizing long-range strong coupling in cavity magnonics provides a general idea for other physical systems. Our experimental achievements may promote the construction of information networks based on cavity magnonics.

Compact continuous-wave solid-state Raman lasers at 707 and 714 nm for laser trapping and cooling of atomic strontium and radium

Two compact laser sources at 707 and 714 nm are realized efficiently by using a diode-pumped a-cut Nd:YVO4 laser with intracavity stimulated Raman scattering and sum-frequency generation (SFG). The fundamental wave at 1342 nm is generated by the 4F3/2 → 4I13/2 transition in Nd:YVO4 crystal. The Raman Stokes waves at 1496 and 1526 nm were obtained by placing the c-axis of the Nd:YVO4 crystal along the Ng and Nm axes of an Np-cut KGW crystal, respectively. LBO crystals with critical phase matching are used to perform the intracavity SFG of fundamental and Stokes waves. At a pump power of 36 W, the maximum output powers at 707 and 714 nm can reach 2.72 and 3.14 W, corresponding to light-to-light conversion efficiencies of 7.5% and 8.7%, respectively. The developed 707 and 714 nm laser sources are practically useful in laser trapping and cooling related to atomic strontium and radium.

Realizing spin squeezing with Rydberg interactions in an optical clock

Neutral-atom arrays trapped in optical potentials are a powerful platform for studying quantum physics, combining precise single-particle control and detection with a range of tunable entangling interactions1-3. For example, these capabilities have been leveraged for state-of-the-art frequency metrology4,5 as well as microscopic studies of entangled many-particle states6-11. Here we combine these applications to realize spin squeezing-a widely studied operation for producing metrologically useful entanglement-in an optical atomic clock based on a programmable array of interacting optical qubits. In this demonstration of Rydberg-mediated squeezing with a neutral-atom optical clock, we generate states that have almost four decibels of metrological gain. In addition, we perform a synchronous frequency comparison between independent squeezed states and observe a fractional-frequency stability of 1.087(1) × 10-15 at one-second averaging time, which is 1.94(1) decibels below the standard quantum limit and reaches a fractional precision at the 10-17 level during a half-hour measurement. We further leverage the programmable control afforded by optical tweezer arrays to apply local phase shifts to explore spin squeezing in measurements that operate beyond the relative coherence time with the optical local oscillator. The realization of this spin-squeezing protocol in a programmable atom-array clock will enable a wide range of quantum-information-inspired techniques for optimal phase estimation and Heisenberg-limited optical atomic clocks12-16.

Molecular Dynamics in Rydberg Tweezer Arrays: Spin-Phonon Entanglement and Jahn-Teller Effect

Atoms confined in optical tweezer arrays constitute a platform for the implementation of quantum computers and simulators. State-dependent operations are realized by exploiting electrostatic dipolar interactions that emerge, when two atoms are simultaneously excited to high-lying electronic states, so-called Rydberg states. These interactions also lead to state-dependent mechanical forces, which couple the electronic dynamics of the atoms to their vibrational motion. We explore these vibronic couplings within an artificial molecular system in which Rydberg states are excited under so-called facilitation conditions. This system, which is not necessarily self-bound, undergoes a structural transition between an equilateral triangle and an equal-weighted superposition of distorted triangular states (Jahn-Teller regime) exhibiting spin-phonon entanglement on a micrometer distance. This highlights the potential of Rydberg tweezer arrays for the study of molecular phenomena at exaggerated length scales.

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.