Dynamics of magnetization at infinite temperature in a Heisenberg spin chain

Understanding universal aspects of quantum dynamics is an unresolved problem in statistical mechanics. In particular, the spin dynamics of the one-dimensional Heisenberg model were conjectured as to belong to the Kardar-Parisi-Zhang (KPZ) universality class based on the scaling of the infinite-temperature spin-spin correlation function. In a chain of 46 superconducting qubits, we studied the probability distribution of the magnetization transferred across the chain's center, [Formula: see text]. The first two moments of [Formula: see text] show superdiffusive behavior, a hallmark of KPZ universality. However, the third and fourth moments ruled out the KPZ conjecture and allow for evaluating other theories. Our results highlight the importance of studying higher moments in determining dynamic universality classes and provide insights into universal behavior in quantum systems.

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.

Measurement-induced entanglement and teleportation on a noisy quantum processor

Measurement has a special role in quantum theory1: by collapsing the wavefunction, it can enable phenomena such as teleportation2 and thereby alter the 'arrow of time' that constrains unitary evolution. When integrated in many-body dynamics, measurements can lead to emergent patterns of quantum information in space-time3-10 that go beyond the established paradigms for characterizing phases, either in or out of equilibrium11-13. For present-day noisy intermediate-scale quantum (NISQ) processors14, the experimental realization of such physics can be problematic because of hardware limitations and the stochastic nature of quantum measurement. Here we address these experimental challenges and study measurement-induced quantum information phases on up to 70 superconducting qubits. By leveraging the interchangeability of space and time, we use a duality mapping9,15-17 to avoid mid-circuit measurement and access different manifestations of the underlying phases, from entanglement scaling3,4 to measurement-induced teleportation18. We obtain finite-sized signatures of a phase transition with a decoding protocol that correlates the experimental measurement with classical simulation data. The phases display remarkably different sensitivity to noise, and we use this disparity to turn an inherent hardware limitation into a useful diagnostic. Our work demonstrates an approach to realizing measurement-induced physics at scales that are at the limits of current NISQ processors.

Non-Abelian braiding of graph vertices in a superconducting processor

Indistinguishability of particles is a fundamental principle of quantum mechanics¹. For all elementary and quasiparticles observed to date-including fermions, bosons and Abelian anyons-this principle guarantees that the braiding of identical particles leaves the system unchanged^2,3. However, in two spatial dimensions, an intriguing possibility exists: braiding of non-Abelian anyons causes rotations in a space of topologically degenerate wavefunctions^4-8. Hence, it can change the observables of the system without violating the principle of indistinguishability. Despite the well-developed mathematical description of non-Abelian anyons and numerous theoretical proposals^9-22, the experimental observation of their exchange statistics has remained elusive for decades. Controllable many-body quantum states generated on quantum processors offer another path for exploring these fundamental phenomena. Whereas efforts on conventional solid-state platforms typically involve Hamiltonian dynamics of quasiparticles, superconducting quantum processors allow for directly manipulating the many-body wavefunction by means of unitary gates. Building on predictions that stabilizer codes can host projective non-Abelian Ising anyons^9,10, we implement a generalized stabilizer code and unitary protocol²³ to create and braid them. This allows us to experimentally verify the fusion rules of the anyons and braid them to realize their statistics. We then study the prospect of using the anyons for quantum computation and use braiding to create an entangled state of anyons encoding three logical qubits. Our work provides new insights about non-Abelian braiding and, through the future inclusion of error correction to achieve topological protection, could open a path towards fault-tolerant quantum computing.

Formation of robust bound states of interacting microwave photons

Systems of correlated particles appear in many fields of modern science and represent some of the most intractable computational problems in nature. The computational challenge in these systems arises when interactions become comparable to other energy scales, which makes the state of each particle depend on all other particles1. The lack of general solutions for the three-body problem and acceptable theory for strongly correlated electrons shows that our understanding of correlated systems fades when the particle number or the interaction strength increases. One of the hallmarks of interacting systems is the formation of multiparticle bound states2-9. Here we develop a high-fidelity parameterizable fSim gate and implement the periodic quantum circuit of the spin-½ XXZ model in a ring of 24 superconducting qubits. We study the propagation of these excitations and observe their bound nature for up to five photons. We devise a phase-sensitive method for constructing the few-body spectrum of the bound states and extract their pseudo-charge by introducing a synthetic flux. By introducing interactions between the ring and additional qubits, we observe an unexpected resilience of the bound states to integrability breaking. This finding goes against the idea that bound states in non-integrable systems are unstable when their energies overlap with the continuum spectrum. Our work provides experimental evidence for bound states of interacting photons and discovers their stability beyond the integrability limit.

Noise-resilient edge modes on a chain of superconducting qubits

Inherent symmetry of a quantum system may protect its otherwise fragile states. Leveraging such protection requires testing its robustness against uncontrolled environmental interactions. Using 47 superconducting qubits, we implement the one-dimensional kicked Ising model, which exhibits nonlocal Majorana edge modes (MEMs) with ℤ 2 parity symmetry. We find that any multiqubit Pauli operator overlapping with the MEMs exhibits a uniform late-time decay rate comparable to single-qubit relaxation rates, irrespective of its size or composition. This characteristic allows us to accurately reconstruct the exponentially localized spatial profiles of the MEMs. Furthermore, the MEMs are found to be resilient against certain symmetry-breaking noise owing to a prethermalization mechanism. Our work elucidates the complex interplay between noise and symmetry-protected edge modes in a solid-state environment.

Constraints on Minute-Scale Transient Astrophysical Neutrino Sources

High-energy neutrino emission has been predicted for several short-lived astrophysical transients including gamma-ray bursts (GRBs), core-collapse supernovae with choked jets, and neutron star mergers. IceCube's optical and x-ray follow-up program searches for such transient sources by looking for two or more muon neutrino candidates in directional coincidence and arriving within 100 s. The measured rate of neutrino alerts is consistent with the expected rate of chance coincidences of atmospheric background events and no likely electromagnetic counterparts have been identified in Swift follow-up observations. Here, we calculate generic bounds on the neutrino flux of short-lived transient sources. Assuming an E^{-2.5} neutrino spectrum, we find that the neutrino flux of rare sources, like long gamma-ray bursts, is constrained to <5% of the detected astrophysical flux and the energy released in neutrinos (100 GeV to 10 PeV) by a median bright GRB-like source is <10^{52.5}  erg. For a harder E^{-2.13} neutrino spectrum up to 30% of the flux could be produced by GRBs and the allowed median source energy is <10^{52}  erg. A hypothetical population of transient sources has to be more common than 10^{-5}  Mpc^{-3}  yr^{-1} (5×10^{-8}  Mpc^{-3}  yr^{-1} for the E^{-2.13} spectrum) to account for the complete astrophysical neutrino flux.

Search for neutrinos from decaying dark matter with IceCube: IceCube Collaboration

With the observation of high-energy astrophysical neutrinos by the IceCube Neutrino Observatory, interest has risen in models of PeV-mass decaying dark matter particles to explain the observed flux. We present two dedicated experimental analyses to test this hypothesis. One analysis uses 6 years of IceCube data focusing on muon neutrino 'track' events from the Northern Hemisphere, while the second analysis uses 2 years of 'cascade' events from the full sky. Known background components and the hypothetical flux from unstable dark matter are fitted to the experimental data. Since no significant excess is observed in either analysis, lower limits on the lifetime of dark matter particles are derived: we obtain the strongest constraint to date, excluding lifetimes shorter than10 28 s at 90% CL for dark matter masses above 10 TeV .

Measurement of the ν μ energy spectrum with IceCube-79: IceCube Collaboration

IceCube is a neutrino observatory deployed in the glacial ice at the geographic South Pole. The ν μ energy unfolding described in this paper is based on data taken with IceCube in its 79-string configuration. A sample of muon neutrino charged-current interactions with a purity of 99.5% was selected by means of a multivariate classification process based on machine learning. The subsequent unfolding was performed using the software Truee. The resulting spectrum covers an E ν -range of more than four orders of magnitude from 125 GeV to 3.2 PeV. Compared to the Honda atmospheric neutrino flux model, the energy spectrum shows an excess of more than 1.9 σ in four adjacent bins for neutrino energies E ν ≥ 177.8 TeV . The obtained spectrum is fully compatible with previous measurements of the atmospheric neutrino flux and recent IceCube measurements of a flux of high-energy astrophysical neutrinos.

Constraints on Ultrahigh-Energy Cosmic-Ray Sources from a Search for Neutrinos above 10 PeV with IceCube

We report constraints on the sources of ultrahigh-energy cosmic rays (UHECRs) above 10^{9}  GeV, based on an analysis of seven years of IceCube data. This analysis efficiently selects very high- energy neutrino-induced events which have deposited energies from 5×10^{5}  GeV to above 10^{11}  GeV. Two neutrino-induced events with an estimated deposited energy of (2.6±0.3)×10^{6}  GeV, the highest neutrino energy observed so far, and (7.7±2.0)×10^{5}  GeV were detected. The atmospheric background-only hypothesis of detecting these events is rejected at 3.6σ. The hypothesis that the observed events are of cosmogenic origin is also rejected at >99% CL because of the limited deposited energy and the nonobservation of events at higher energy, while their observation is consistent with an astrophysical origin. Our limits on cosmogenic neutrino fluxes disfavor the UHECR sources having a cosmological evolution stronger than the star formation rate, e.g., active galactic nuclei and γ-ray bursts, assuming proton-dominated UHECRs. Constraints on UHECR sources including mixed and heavy UHECR compositions are obtained for models of neutrino production within UHECR sources. Our limit disfavors a significant part of parameter space for active galactic nuclei and new-born pulsar models. These limits on the ultrahigh-energy neutrino flux models are the most stringent to date.

Searches for Sterile Neutrinos with the IceCube Detector

The IceCube neutrino telescope at the South Pole has measured the atmospheric muon neutrino spectrum as a function of zenith angle and energy in the approximate 320 GeV to 20 TeV range, to search for the oscillation signatures of light sterile neutrinos. No evidence for anomalous ν_{μ} or ν[over ¯]_{μ} disappearance is observed in either of two independently developed analyses, each using one year of atmospheric neutrino data. New exclusion limits are placed on the parameter space of the 3+1 model, in which muon antineutrinos experience a strong Mikheyev-Smirnov-Wolfenstein-resonant oscillation. The exclusion limits extend to sin^{2}2θ_{24}≤0.02 at Δm^{2}∼0.3  eV^{2} at the 90% confidence level. The allowed region from global analysis of appearance experiments, including LSND and MiniBooNE, is excluded at approximately the 99% confidence level for the global best-fit value of |U_{e4}|^{2}.