Spectroscopic signatures of localization with interacting photons in superconducting qubits
P. Roushan, C. Neill, Jirawat Tangpanitanon, V. M. Bastidas, A. Megrant, R. Barends, Yu Chen, Z. Chen, B. Chiaro, A. Dunsworth, Austin G. Fowler, Brooks Foxen, Marissa Giustina, Evan R. Jeffrey, J. Kelly, Erik A. Lucero, J. Mutus, M. Neeley, Chris Quintana, D. Sank, A. Vainsencher, J. Wenner, T. White, Hartmut Neven, Dimitris G. Angelakis, John M. Martinis
Quantized eigenenergies and their associated wave functions provide extensive information for predicting the physics of quantum many-body systems. Using a chain of nine superconducting qubits, we implement a technique for resolving the energy levels of interacting photons. We benchmark this method by capturing the main features of the intricate energy spectrum predicted for two-dimensional electrons in a magnetic field-the Hofstadter butterfly. We introduce disorder to study the statistics of the energy levels of the system as it undergoes the transition from a thermalized to a localized phase. Our work introduces a many-body spectroscopy technique to study quantum phases of matter.
Quantized eigenenergies and their associated wave functions provide extensive information for predicting the physics of quantum many-body systems. Using a chain of nine superconducting qubits, we implement a technique for resolving the energy levels of interacting photons. We benchmark this method by capturing the main features of the intricate energy spectrum predicted for two-dimensional electrons in a magnetic field-the Hofstadter butterfly. We introduce disorder to study the statistics of the energy levels of the system as it undergoes the transition from a thermalized to a localized phase. Our work introduces a many-body spectroscopy technique to study quantum phases of matter.
A promising approach to study quantum materials is to simulate them on an engineered quantum platform. However, achieving the accuracy needed to outperform classical methods has been an outstanding challenge. Here, using eighteen superconducting qubits, we provide an experimental blueprint for a programmable and accurate quantum matter simulator and demonstrate how to probe fundamental electronic properties. We illustrate the underlying method by reconstructing the single-particle band-structure of a one-dimensional wire. We demonstrate nearly complete mitigation of decoherence and readout errors and arrive at an accuracy in measuring energy eigenvalues of this wire with an error of ~0.01 radians, whereas typical energy scales are of order 1 radian. Insight into this unprecedented algorithm fidelity is gained by highlighting robust properties of a Fourier transform, including the ability to resolve eigenenergies with a statistical uncertainty of 10−4 radians. Furthermore, we synthesize magnetic flux and disordered local potentials, two key tenets of a condensed-matter system. When sweeping the magnetic flux, we observe avoided level crossings in the spectrum, a detailed fingerprint of the spatial distribution of local disorder. Combining these methods, we reconstruct electronic properties of the eigenstates where we observe persistent currents and a strong suppression of conductance with added disorder. Our work describes an accurate method for quantum simulation and paves the way to study novel quantum materials with superconducting qubits.
Abstract A promising approach to study condensed-matter systems is to simulate them on an engineered quantum platform1–4. However, the accuracy needed to outperform classical methods has not been achieved so far. Here, using 18 superconducting qubits, we provide an experimental blueprint for an accurate condensed-matter simulator and demonstrate how to investigate fundamental electronic properties. We benchmark the underlying method by reconstructing the single-particle band structure of a one-dimensional wire. We demonstrate nearly complete mitigation of decoherence and readout errors, and measure the energy eigenvalues of this wire with an error of approximately 0.01 rad, whereas typical energy scales are of the order of 1 rad. Insight into the fidelity of this algorithm is gained by highlighting the robust properties of a Fourier transform, including the ability to resolve eigenenergies with a statistical uncertainty of 10−4 rad. We also synthesize magnetic flux and disordered local potentials, which are two key tenets of a condensed-matter system. When sweeping the magnetic flux we observe avoided level crossings in the spectrum, providing a detailed fingerprint of the spatial distribution of local disorder. By combining these methods we reconstruct electronic properties of the eigenstates, observing persistent currents and a strong suppression of conductance with added disorder. Our work describes an accurate method for quantum simulation5, 6 and paves the way to study new quantum materials with superconducting qubits.
The question of how the low-energy properties of disordered quantum systems may be connected to exotic localization phenomena at high energy is a key open question in the context of quantum glasses and many-body localization. In the preceding Letter [L. Villa, S. J. Thomson, and L. Sanchez-Palencia, preceding Letter, Phys. Rev. A 104, L021301 (2021)] we have shown that key features of the excitation spectrum of a disordered system can be efficiently probed from out-of-equilibrium dynamics following a quantum quench, providing distinctive signatures of the various phases. Here we extend this work by providing a more-in-depth study of the behavior of the quench spectral functions associated with different observables and investigating an extended parameter regime. We provide a detailed introduction to quench spectroscopy for disordered systems and show how spectral properties can be probed using both local operators and two-point correlation functions. We benchmark the technique using the one-dimensional Bose-Hubbard model in the presence of a random external potential, focusing on the low-lying excitations, and demonstrate that quench spectroscopy can distinguish the Mott insulator, superfluid, and Bose glass phases. We then explicitly reconstruct the zero-temperature phase diagram of the disordered Bose-Hubbard at fixed filling using two independent methods, experimentally accessible via both time-of-flight imaging and quantum gas microscopy, respectively, and demonstrate that quench spectroscopy can give valuable insights into the distribution of rare regions within disordered systems.
The study of interacting quantum many-body systems is constrained by the exponential growth of the Hilbert space with system size, which limits both numerical simulations and theoretical understanding. This thesis explores how the interplay between disorder and interactions shapes localization and ergodicity in quantum many-body systems, using the exact diagonalization method, and how machine learning can serve as an alternative approach for studying these systems. We first investigate the many-body localization problem in a one-dimensional SU(2)-invariant random antiferromagnetic Heisenberg chain using exact and sparse matrix diagonalization methods. By analyzing energy level gap-ratio statistics, participation ratio, entanglement entropy and the entanglement spectral parameter as a function of disorder strength under di!erent random coupling distributions, we identify clear qualitative di!erences from the widely studied random field Heisenberg chain. For system sizes up to L = 24, signatures of localized behaviour are much less pronounced than in the case of the random field model. Motivated by the limitations of conventional numerical methods, we examine machine learning as a scalable paradigm for many-body physics. Using the Bose- Hubbard model in one and two dimensions as a benchmark, we test HUBBARDNET, a deep neural network architecture, by replacing standard stochastic gradient descent with adaptive moment estimation (ADAM), adjusting the learning rate for improved optimization, and introducing physics-informed output activations that can resolve extremely small wave-function amplitudes. This approach yields ground-state energy errors reduced by orders of magnitude and wave function fidelities exceeding 99%, remaining robust across four decades of interaction strength and over superfluid, Mott-insulating, and quantum chaotic regimes. Beyond energy prediction, the trained models accurately reproduce generalized inverse participation ratios and multifractal dimensions for ground and excited states, capturing localization, delocalization, and multifractality across the spectrum, with excited-state training explicitly guided by fractal dimensions that encode the multifractal properties of the wave functions.
Simulating quantum many-body systems on a classical computer generally requires a computational cost that grows exponential with the number of particles. This computational complexity has been the main obstacle to understanding various fundamental emergent phenomena in condensed matters such as high-Tc superconductivity and the fractional quantum-Hall effect. The difficulty arises because even the simplest models that are proposed to capture those phenomena cannot be simulated on a classical computer. Recognizing this problem in 1981, Richard Feynman envisioned a quantum simulator, an entirely new type of machine that exploits quantum superposition and operates by individually manipulating its constituting quantum particles and their interactions. Recent advances in various experimental platforms from cold atoms in optical lattices, trapped ions, to solid-state systems have brought the idea of Feynman to the realm of reality. Among those, interacting photons in superconducting circuits has been one of the promising platforms thanks to their local controllability and long coherence times. Early theoretical proposals have shown possibilities to realize quantum many-body phenomena of light using coupled cavity arrays such as Mott to superfluid transitions and fractional quantum Hall states. Start-of-the-art experiments include realization of interacting chiral edge states and stroboscopic signatures of localization of interacting photons in a three-site and a nine-site superconducting circuit, respectively. Interacting photons also serve as a natural platform to simulate driven-dissipative quantum many-body phenomena. A 72-site superconducting circuit has also recently been fabricated to study a dissipative phase transition of light.
Simulating quantum many-body systems on a classical computer generally requires a computational cost that grows exponential with the number of particles. This computational complexity has been the main obstacle to understanding various fundamental emergent phenomena in condensed matters such as high-Tc superconductivity and the fractional quantum-Hall effect. The difficulty arises because even the simplest models that are proposed to capture those phenomena cannot be simulated on a classical computer. Recognizing this problem in 1981, Richard Feynman envisioned a quantum simulator, an entirely new type of machine that exploits quantum superposition and operates by individually manipulating its constituting quantum particles and their interactions. Recent advances in various experimental platforms from cold atoms in optical lattices, trapped ions, to solid-state systems have brought the idea of Feynman to the realm of reality. Among those, interacting photons in superconducting circuits has been one of the promising platforms thanks to their local controllability and long coherence times. Early theoretical proposals have shown possibilities to realize quantum many-body phenomena of light using coupled cavity arrays such as Mott to superfluid transitions and fractional quantum Hall states. Start-of-the-art experiments include realization of interacting chiral edge states and stroboscopic signatures of localization of interacting photons in a three-site and a nine-site superconducting circuit, respectively. Interacting photons also serve as a natural platform to simulate driven-dissipative quantum many-body phenomena. A 72-site superconducting circuit has also recently been fabricated to study a dissipative phase transition of light.