The phenomenon of particle production for quantum field theories in curved spacetimes is crucial to understand the large-scale structure of a universe from an inflationary epoch. In contrast to the free and fixed-background case, the production of particles with strong interactions and back reaction is not completely understood, especially in situations that require going beyond perturbation theory. In this work, we present advances in this direction by focusing on a self-interacting field theory of Dirac fermions in an expanding Friedmann-Robertson-Walker universe. By using a Hamiltonian lattice regularization with continuous conformal time and rescaled fields, this model becomes amenable to either a cold-atom analogue-gravity quantum simulation, or a dynamical variational approach. Leveraging a family of variational fermionic Gaussian states, we investigate how dynamical mass generation and the formation of fermion condensates associated to certain broken symmetries modify some well-known results of the free field theory. In particular, we study how the non-perturbative condensates arise and, more importantly, how their real-time evolution has an impact on particle production. Depending on the Hubble expansion rate, we find an interesting interplay of interactions and particle production, including a non-trivial back reaction on the condensates and a parity-breaking spectrum of produced particles.

Lattice gauge theories (LGTs) describe a broad range of phenomena in condensed matter and particle physics. A prominent example is confinement, responsible for bounding quarks inside hadrons such as protons or neutrons. When quark-antiquark pairs are separated, the energy stored in the string of gluon fields connecting them grows linearly with their distance, until there is enough energy to create new pairs from the vacuum and break the string. While such phenomena are ubiquitous in LGTs, simulating the resulting dynamics is a challenging task. Here, we report the observation of string breaking in synthetic quantum matter using a programmable quantum simulator based on neutral atom arrays. We show that a (2+1)D LGT with dynamical matter can be efficiently implemented when the atoms are placed on a Kagome geometry, with a local U(1) symmetry emerging from the Rydberg blockade, while long-range Rydberg interactions naturally give rise to a linear confining potential for a pair of charges, allowing us to tune both their masses as well as the string tension. We experimentally map out the corresponding phase diagram by adiabatically preparing the ground state of the atom array in the presence of defects, and observe substructure of the confined phase, distinguishing regions dominated by fluctuating strings or by broken string configurations. Finally, by harnessing local control over the atomic detuning, we quench string states and observe string breaking dynamics exhibiting a many-body resonance phenomenon. Our work paves a way to explore phenomena in high-energy physics using programmable quantum simulators.

Many-body fermionic systems can be simulated in a hardware-efficient manner using a fermionic quantum processor. Neutral atoms trapped in optical potentials can realize such processors, where nonlocal fermionic statistics are guaranteed at the hardware level. Implementing quantum error correction in this setup is, however, challenging, due to the atom-number superselection present in atomic systems, that is, the impossibility of creating coherent superpositions of different particle numbers. In this Letter, we overcome this constraint and present a blueprint for an error-corrected fermionic quantum processor that can be implemented using current experimental capabilities. To achieve this, we first consider an ancillary set of fermionic modes and design a fermionic reference, which we then use to construct superpositions of different numbers of referenced fermions. This allows us to build logical fermionic modes that can be error corrected using standard atomic operations. Here, we focus on phase errors, which we expect to be a dominant source of errors in neutral-atom quantum processors. We then construct logical fermionic gates, and show their implementation for the logical particle-number conserving processes relevant for quantum simulation. Finally, our protocol is illustrated with a minimal fermionic circuit, where it leads to a quadratic suppression of the logical error rate.

Fermionic atoms in optical lattices provide a native implementation of Fermi-Hubbard (FH) models that can be used as analog quantum simulators of many-body fermionic systems. Recent experimental advances include the time-dependent local control of chemical potentials and tunnelings, and thus enable us to operate this platform digitally as a programmable quantum simulator. Here, we explore these opportunities and develop ground-state preparation algorithms for different fermionic models, based on the ability to implement both single-particle and many-body high-fidelity fermionic gates, as provided by the native FH Hamiltonian. In particular, we first design variational precompiled quantum circuits to prepare the ground state of the natively implemented FH model, with significant speedups relative to competing adiabatic protocols. Besides, the versatility of this variational approach enables us to target extended FH models, i.e., including terms that are not natively realized on the platform. As an illustration, we include next-nearest-neighbor tunnelings at finite dopings, relevant in the context of d-wave superconductivity. Furthermore, we discuss how to approximate the imaginary-time evolution using variational fermionic circuits, both as an alternative state-preparation strategy and as a subroutine for the quantum Lanczos algorithm to further improve the energy estimation. We benchmark our protocols for ladder geometries, though they can be readily applied to two-dimensional experimental setups to address regimes beyond the capabilities of current classical methods. These results pave the way for more efficient and comprehensive explorations of relevant many-body phases with existing programmable fermionic quantum simulators.

Quantum computers offer an intriguing path for a paradigmatic change of computing in the natural sciences and beyond, with the potential for achieving a so-called quantum advantage-namely, a significant (in some cases exponential) speedup of numerical simulations. The rapid development of hardware devices with various realizations of qubits enables the execution of small-scale but representative applications on quantum computers. In particular, the high-energy physics community plays a pivotal role in accessing the power of quantum computing, since the field is a driving source for challenging computational problems. This concerns, on the theoretical side, the exploration of models that are very hard or even impossible to address with classical techniques and, on the experimental side, the enormous data challenge of newly emerging experiments, such as the upgrade of the Large Hadron Collider. In this Roadmap paper, led by CERN, DESY, and IBM, we provide the status of high-energy physics quantum computations and give examples of theoretical and experimental target benchmark applications, which can be addressed in the near future. Having in mind hardware with about 100 qubits capable of executing several thousand two-qubit gates, where possible, we also provide resource estimates for the examples given using error-mitigated quantum computing. The ultimate declared goal of this task force is therefore to trigger further research in the high-energy physics community to develop interesting use cases for demonstrations on near-term quantum computers.

Ultracold Fermi gases of spin-3/2 atoms provide a clean platform to realise SO(5) models of 4-Fermi interactions in the laboratory. By confining the atoms in a two-dimensional Raman lattice, we show how this system can be used as a flexible quantum simulator of Dirac quantum field theories (QFTs) that combine Gross-Neveu and Thirring interactions with a higher-order topological twist. We show that the lattice model corresponds to a regularization of this QFT with an anisotropic twisted Wilson mass. This allows us to access higher-order topological states protected by a hidden SO(5) symmetry, a remnant of the original rotational symmetry of the 4-Fermi interactions that is not explicitly broken by the lattice discretization. Using large-$N$ methods, we show that the 4-Fermi interactions lead to a rich phase diagram with various competing fermion condensates. Our work opens a route for the implementation of correlated higher-order topological states with tunable interactions that has interesting connections to non-trivial relativistic QFTs of Dirac fermions in $D = 2 + 1$ dimensions.

Simulating the real-time dynamics of lattice gauge theories, underlying the Standard Model of particle physics, is a notoriously difficult problem where quantum simulators can provide a practical advantage over classical approaches. In this work, we present a complete Rydberg-based architecture, co-designed to digitally simulate the dynamics of general gauge theories coupled to matter fields in a hardware-efficient manner. Ref. \cite{Gonzalez_2022} showed how a qudit processor, where non-abelian gauge fields are locally encoded and time-evolved, considerably reduces the required simulation resources compared to standard qubit-based quantum computers. Here we integrate the latter with a recently introduced fermionic quantum processor \cite{Gonzalez_2023}, where fermionic statistics are accounted for at the hardware level, allowing us to construct quantum circuits that preserve the locality of the gauge-matter interactions. We exemplify the flexibility of such a fermion-qudit processor by focusing on two paradigmatic high-energy phenomena. First, we present a resource-efficient protocol to simulate the Abelian-Higgs model, where the dynamics of confinement and string breaking can be investigated. Then, we show how to prepare hadrons made up of fermionic matter constituents bound by non-abelian gauge fields, and show how to extract the corresponding hadronic tensor. In both cases, we estimate the required resources, showing how quantum devices can be used to calculate experimentally-relevant quantities in particle physics.

Quantum computers offer an intriguing path for a paradigmatic change of computing in the natural sciences and beyond, with the potential for achieving a so-called quantum advantage, namely a significant (in some cases exponential) speed-up of numerical simulations. The rapid development of hardware devices with various realizations of qubits enables the execution of small scale but representative applications on quantum computers. In particular, the high-energy physics community plays a pivotal role in accessing the power of quantum computing, since the field is a driving source for challenging computational problems. This concerns, on the theoretical side, the exploration of models which are very hard or even impossible to address with classical techniques and, on the experimental side, the enormous data challenge of newly emerging experiments, such as the upgrade of the Large Hadron Collider. In this roadmap paper, led by CERN, DESY and IBM, we provide the status of high-energy physics quantum computations and give examples for theoretical and experimental target benchmark applications, which can be addressed in the near future. Having the IBM 100 x 100 challenge in mind, where possible, we also provide resource estimates for the examples given using error mitigated quantum computing.

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.

Simulating the properties of many-body fermionic systems is an outstanding computational challenge relevant to material science, quantum chemistry, and particle physics.-5.4pc]Please note that the spelling of the following author names in the manuscript differs from the spelling provided in the article metadata: D. González-Cuadra, D. Bluvstein, M. Kalinowski, R. Kaubruegger, N. Maskara, P. Naldesi, T. V. Zache, A. M. Kaufman, M. D. Lukin, H. Pichler, B. Vermersch, Jun Ye, and P. Zoller. The spelling provided in the manuscript has been retained; please confirm. Although qubit-based quantum computers can potentially tackle this problem more efficiently than classical devices, encoding nonlocal fermionic statistics introduces an overhead in the required resources, limiting their applicability on near-term architectures. In this work, we present a fermionic quantum processor, where fermionic models are locally encoded in a fermionic register and simulated in a hardware-efficient manner using fermionic gates. We consider in particular fermionic atoms in programmable tweezer arrays and develop different protocols to implement nonlocal gates, guaranteeing Fermi statistics at the hardware level. We use this gate set, together with Rydberg-mediated interaction gates, to find efficient circuit decompositions for digital and variational quantum simulation algorithms, illustrated here for molecular energy estimation. Finally, we consider a combined fermion-qubit architecture, where both the motional and internal degrees of freedom of the atoms are harnessed to efficiently implement quantum phase estimation as well as to simulate lattice gauge theory dynamics.

Treating the infinite-dimensional Hilbert space of non-Abelian gauge theories is an outstanding challenge for classical and quantum simulations. Here, we employ q-deformed Kogut-Susskind lattice gauge theories, obtained by deforming the defining symmetry algebra to a quantum group. In contrast to other formulations, this approach simultaneously provides a controlled regularization of the infinite-dimensional local Hilbert space while preserving essential symmetry-related properties. This enables the development of both quantum as well as quantum-inspired classical spin-network algorithms for q-deformed gauge theories. To be explicit, we focus on SU(2)_{k} gauge theories with k∈N that are controlled by the deformation parameter q=e^{2πi/(k+2)}, a root of unity, and converge to the standard SU(2) Kogut-Susskind model as k→∞. In particular, we demonstrate that this formulation is well suited for efficient tensor network representations by variational ground-state simulations in 2D, providing first evidence that the continuum limit can be reached with k=O(10). Finally, we develop a scalable quantum algorithm for Trotterized real-time evolution by analytically diagonalizing the SU(2)_{k} plaquette interactions. Our work gives a new perspective for the application of tensor network methods to high-energy physics and paves the way for quantum simulations of non-Abelian gauge theories far from equilibrium where no other methods are currently available.

Quantum information platforms made great progress in the control of many-body entanglement and the implementation of quantum error correction, but it remains a challenge to realize both in the same setup. Here, we propose a mixture of two ultracold atomic species as a platform for universal quantum computation with long-range entangling gates, while providing a natural candidate for quantum error-correction. In this proposed setup, one atomic species realizes localized collective spins of tunable length, which form the fundamental unit of information. The second atomic species yields phononic excitations, which are used to entangle collective spins. Finally, we discuss a finite-dimensional version of the Gottesman-Kitaev-Preskill code to protect quantum information encoded in the collective spins, opening up the possibility to universal fault-tolerant quantum computation in ultracold atom systems.

We propose a protocol to detect topological phase transitions of one-dimensional p-wave superconductors from their harmonic emission spectra in strong fields. Specifically, we identify spectral features due to radiating edge modes, which characterize the spectrum and the density of states in the topological phase and are absent in the trivial phase. These features allow us to define a measurable signature, obtained from emission measurements, that unambiguously differentiates between the two phases. Local probing provides insight into the localized and topologically protected nature of the modes. The presented results establish that high-harmonic spectroscopy can be used as an all-optical tool for the detection of Majorana zero modes.

Engineering long-range interactions in cold-atom quantum simulators can lead to exotic quantum many-body behavior. Fermionic atoms in ultracold atomic mixtures can act as mediators, giving rise to long-range Ruderman-Kittel-Kasuya-Yosida-type interactions characterized by the dimensionality and density of the fermionic gas. Here, we propose several tuning knobs, accessible in current experimental platforms, that allow one to further control the range and shape of the mediated interactions, extending the existing quantum simulation toolbox. In particular, we include an additional optical lattice for the fermionic mediator, as well as anisotropic traps to change its dimensionality in a continuous manner. This allows us to interpolate between power-law and exponential decays, introducing an effective cutoff for the interaction range, as well as to tune the relative interaction strengths at different distances. Finally, we show how our approach allows one to investigate frustrated regimes that were not previously accessible, where symmetry-protected topological phases as well as chiral spin liquids emerge.

Higher-order topological crystalline phases in low-dimensional interacting quantum systems represent a challenging and largely unexplored research topic. Here, we derive a Hamiltonian describing fermions interacting through correlated hopping processes that break chiral invariance, but preserve both inversion and time-reversal symmetries. In this way, we show that our one-dimensional model gives rise to an interacting second-order topological insulating phase that supports gapped edge states. The topological nature of such interacting phase turns out to be revealed by both long-range order of a non-local string correlation function and by even degeneracy of the entanglement spectrum. For strong interactions we instead find that the topological crystalline phase is destroyed and replaced by a singlet superconducting phase. The latter, characterized by local fermionic pairing, turns out to appear both in a homogeneous and in a phase separated form. Relevantly, the derived one-dimensional model and the second-order topological insulator can be explored and investigated in atomic quantum simulators.

We consider a mixture of ultracold bosonic atoms on a one-dimensional lattice described by the XXZ-Bose-Hubbard model, where the tunneling of one species depends on the spin state of a second deeply trapped species. We show how the inclusion of antiferromagnetic interactions among the spin degrees of freedom generates a Devil's staircase of symmetry-protected topological phases for a wide parameter regime via a bosonic analog of the Peierls mechanism in electron-phonon systems. These topological Peierls insulators are examples of symmetry-breaking topological phases, where long-range order due to spontaneous symmetry breaking coexists with topological properties such as fractionalized edge states. Moreover, we identify a region of supersolid phases that do not require long-range interactions. They appear instead due to a Peierls incommensurability mechanism, where competing orders modify the underlying crystalline structure of Peierls insulators, becoming superfluid. Our work show the possibilities that ultracold atomic systems offer to investigate strongly-correlated topological phenomena beyond those found in natural materials.

We investigate the topological properties of the bond order wave phase arising in the extended Fermi-Hubbard model. In particular, we uncover a topological sector, which remained elusive in previous finite-size numerical studies due to boundary effects. We first show that, for an infinite system, the bond order wave regime is characterized by two degenerate bulk states corresponding to the trivial and topological sectors. The latter turns out to be indeed characterized by an even degeneracy of the entanglement spectrum and longe-range order of a string correlation function. For finite size systems, we show that the topological sector can be stabilized by imposing a suitable border potential. This therefore provides a concrete protocol for the observation of topologically protected degenerate edge modes in finite-size systems. Furthermore, we show that the bulk of the system is characterized by exotic solitonic solutions interpolating between the trivial and topological sectors. Finally, we propose an implementation and detection scheme of this strongly-correlated topological phase in a quantum simulator based on dipolar Fermi gases in optical lattices.

The combination of topology and quantum criticality can give rise to an exotic mix of counterintuitive effects. Here, we show that unexpected topological properties take place in a paradigmatic strongly correlated Hamiltonian: the 1D extended Bose-Hubbard model. In particular, we reveal the presence of two distinct topological quantum critical points with localized edge states and gapless bulk excitations. Our results show that the topological critical points separate two phases, one topologically protected and the other topologically trivial, both characterized by a long-range ordered string correlation function. The longrange order persists also at the topological critical points and explains the presence of localized edge states protected by a finite charge gap. Finally, we introduce a superresolution quantum gas microscopy scheme for dipolar dysprosium atoms, which provides a reliable route towards the experimental study of topological quantum critical points.

The weak ergodicity breaking induced by quantum many-body scars (QMBSs) represents an intriguing concept that has received great attention in recent years due to its relation to unusual nonequilibrium behavior. Here, we reveal that this phenomenon can occur in a particular regime of a lattice gauge theory, where QMBSs emerge due to the presence of an extensive number of local constraints. In particular, by analyzing the gauged Kitaev model, we provide an example where QMBSs appear in a regime where charges are deconfined. By means of both numerical and analytical approaches, we find a variety of scarred states far away from the regime where the model is integrable. The presence of these states is revealed both by tracing them directly from the analytically reachable limit, as well as by quantum quenches showing persistent oscillations for specific initial states.

Quantum paramagnets are strongly-correlated phases of matter where competing interactions frustrate magnetic order down to zero temperature. In certain cases, quantum fluctuations induce instead topological order, supporting, in particular, fractionalized quasi-particle excitations. In this work, we investigate paradigmatic spin models and show how magnetic frustration can also give rise to higher-order topological properties. We first study the frustrated Heisenberg model in a square lattice, where a plaquette valence bond solid appears through the spontaneous breaking of translational invariance. Despite the amount of effort that has been devoted to study this phase, its topological nature has so far been overlooked. By means of tensor network simulations, we establish how such state belongs to a higher-order symmetry-protected topological phase, where long-range plaquette order and non-trivial topology coexist. This interplay allows the system to support excitations that would be absent otherwise, such as corner-like states in the bulk attached to dynamical topological defects. Finally, we demonstrate how this higher-order topological quantum paramagnet can also be induced by dipolar interactions, indicating the possibility to directly observe this phase using atomic quantum simulators.