Principal component analysis of absorbing state phase transitions
Muzzi C., Cortes R.S., Bhakuni D.S., Jelić A., Gambassi A., We perform a principal component analysis (PCA) of two one-dimensional lattice models belonging to distinct nonequilibrium universality classes - directed bond percolation and branching and annihilating random walks with an even number of offspring. We find that the uncentered PCA of data sets storing various system's configurations can be successfully used to determine the critical properties of these nonequilibrium phase transitions. In particular, in both cases, we obtain good estimates of the critical point and the dynamical critical exponent of the models. For directed bond percolation, we are furthermore able to extract critical exponents associated with the correlation length and the order parameter. We discuss the relation of our analysis with low-rank approximations of data sets.
Entanglement - Nonstabilizerness separation in hybrid quantum circuits
Fux G.E., Tirrito E., Nonstabilizerness describes the distance of a quantum state to its closest stabilizer state. It is - like entanglement - a necessary resource for a quantum advantage over classical computing. We study nonstabilizerness, quantified by stabilizer entropy, in a hybrid quantum circuit with projective measurements and a controlled injection of non-Clifford resources. We discover a phase transition between a power law and constant scaling of nonstabilizerness with system size controlled by the rate of measurements. The same circuit also exhibits a phase transition in entanglement that appears, however, at a different critical measurement rate. This mechanism shows how, from the viewpoint of a quantum advantage, hybrid circuits can host multiple distinct transitions where not only entanglement, but also other nonlinear properties of the density matrix come into play.
Diagnosing quantum transport from wave function snapshots
Bhakuni D.S., Verdel R., Muzzi C., Andreoni R., Aidelsburger M., We study nonequilibrium quantum dynamics of spin chains by employing principal component analysis on data sets of wave function snapshots and examine how information propagates within these data sets. The quantities we employ are derived from the spectrum of the sample second moment matrix, built directly from data sets. Our investigations on several interacting spin chains featuring distinct spin or energy transport reveal that the growth of data information spreading follows the same dynamical exponents as that of the underlying quantum transport of spin or energy. Specifically, our approach enables an easy, data-driven, and, importantly, interpretable diagnostic to track energy transport with a limited number of samples, which is usually challenging without any assumption on the Hamiltonian form. These observations are obtained at a modest finite-size and evolution time, which aligns with experimental and numerical constraints. Our framework directly applies to experimental quantum simulator data sets of dynamics in higher-dimensional systems, where classical simulation methods usually face significant limitations and apply equally to both near- and far-from-equilibrium quenches.
Nonstabilizerness versus entanglement in matrix product states
Frau M., Tarabunga P.S., Collura M., In this paper, we investigate the relationship between entanglement and nonstabilizerness (also known as magic) in matrix product states (MPSs). We study the relation between magic and the bond dimension used to approximate the ground state of a many-body system in two different contexts: full state of magic and mutual magic (the nonstabilizer analog of mutual information, thus free of boundary effects) of spin-1 anisotropic Heisenberg chains. Our results indicate that obtaining converged results for nonstabilizerness is typically considerably easier than entanglement. For full state magic at critical points and at sufficiently large volumes, we observe convergence with 1/χ2, with χ being the MPS bond dimension. At small volumes, magic saturation is so quick that, within error bars, we cannot appreciate any finite-χ correction. Mutual magic also shows a fast convergence with bond dimension, whose specific functional form is however hindered by sampling errors. As a byproduct of our study, we show how Pauli-Markov chains (originally formulated to evaluate magic) resets the state of the art in terms of computing mutual information for MPS. We illustrate this last fact by verifying the logarithmic increase of mutual information between connected partitions at critical points. By comparing mutual information and mutual magic, we observe that, for connected partitions, the latter is typically scaling much slower - if at all - with the partition size, while for disconnected partitions, both are constant in size.
Nonstabilizerness via Matrix Product States in the Pauli Basis
Tarabunga P.S., Tirrito E., Bañuls M.C., Nonstabilizerness, also known as "magic,"stands as a crucial resource for achieving a potential advantage in quantum computing. Its connection to many-body physical phenomena is poorly understood at present, mostly due to a lack of practical methods to compute it at large scales. We present a novel approach for the evaluation of nonstabilizerness within the framework of matrix product states (MPSs), based on expressing the MPS directly in the Pauli basis. Our framework provides a powerful tool for efficiently calculating various measures of nonstabilizerness, including stabilizer Rényi entropies, stabilizer nullity, and Bell magic, and enables the learning of the stabilizer group of an MPS. We showcase the efficacy and versatility of our method in the ground states of Ising and XXZ spin chains, as well as in circuits dynamics that has recently been realized in Rydberg atom arrays, where we provide concrete benchmarks for future experiments on logical qubits up to twice the sizes already realized.
Network science: Ising states of matter
Sun H., Panda R.K., Verdel R., Rodriguez A., Network science provides very powerful tools for extracting information from interacting data. Although recently the unsupervised detection of phases of matter using machine learning has raised significant interest, the full prediction power of network science has not yet been systematically explored in this context. Here we fill this gap by providing an in-depth statistical, combinatorial, geometrical, and topological characterization of 2D Ising snapshot networks (IsingNets) extracted from Monte Carlo simulations of the 2D Ising model at different temperatures, going across the phase transition. Our analysis reveals the complex organization properties of IsingNets in both the ferromagnetic and paramagnetic phases and demonstrates the significant deviations of the IsingNets with respect to randomized null models. In particular percolation properties of the IsingNets reflect the existence of the symmetry between configurations with opposite magnetization below the critical temperature and the very compact nature of the two emerging giant clusters revealed by our persistent homology analysis of the IsingNets. Moreover, the IsingNets display a very broad degree distribution and significant degree-degree correlations and weight-degree correlations demonstrating that they encode relevant information present in the configuration space of the 2D Ising model. The geometrical organization of the critical IsingNets is reflected in their spectral properties deviating from the one of the null model. This work reveals the important insights that network science can bring to the characterization of phases of matter. The set of tools described hereby can be applied as well to numerical and experimental data.
Wave-Function Network Description and Kolmogorov Complexity of Quantum Many-Body Systems
Mendes-Santos T., Schmitt M., Angelone A., Rodriguez A., Scholl P., Williams H.J., Barredo D., Lahaye T., Browaeys A., Heyl M., Programmable quantum devices are now able to probe wave functions at unprecedented levels. This is based on the ability to project the many-body state of atom and qubit arrays onto a measurement basis which produces snapshots of the system wave function. Extracting and processing information from such observations remains, however, an open quest. One often resorts to analyzing low-order correlation functions - that is, discarding most of the available information content. Here, we introduce wave-function networks - a mathematical framework to describe wave-function snapshots based on network theory. For many-body systems, these networks can become scale-free - a mathematical structure that has found tremendous success and applications in a broad set of fields, ranging from biology to epidemics to Internet science. We demonstrate the potential of applying these techniques to quantum science by introducing protocols to extract the Kolmogorov complexity corresponding to the output of a quantum simulator and implementing tools for fully scalable cross-platform certification based on similarity tests between networks. We demonstrate the emergence of scale-free networks analyzing experimental data obtained with a Rydberg quantum simulator manipulating up to 100 atoms. Our approach illustrates how, upon crossing a phase transition, the simulator complexity decreases while correlation length increases - a direct signature of buildup of universal behavior in data space. Comparing experiments with numerical simulations, we achieve cross-certification at the wave-function level up to timescales of 4 μs with a confidence level of 90% and determine experimental calibration intervals with unprecedented accuracy. Our framework is generically applicable to the output of quantum computers and simulators with in situ access to the system wave function and requires probing accuracy and repetition rates accessible to most currently available platforms.
Quantifying nonstabilizerness through entanglement spectrum flatness
Tirrito E., Tarabunga P.S., Lami G., Chanda T., Leone L., Oliviero S.F.E., Nonstabilizerness, also colloquially referred to as magic, is a resource for advantage in quantum computing and lies in the access to non-Clifford operations. Developing a comprehensive understanding of how nonstabilizerness can be quantified and how it relates to other quantum resources is crucial for studying and characterizing the origin of quantum complexity. In this work, we establish a direct connection between nonstabilizerness and entanglement spectrum flatness for a pure quantum state. We show that this connection can be exploited to efficiently probe nonstabilizerness even in the presence of noise. Our results reveal a direct connection between nonstabilizerness and entanglement response, and define a clear experimental protocol to probe nonstabilizerness in cold atom and solid-state platforms.
Topological Kolmogorov complexity and the Berezinskii-Kosterlitz-Thouless mechanism
Vitale V., Mendes-Santos T., Rodriguez A., Topology plays a fundamental role in our understanding of many-body physics, from vortices and solitons in classical field theory to phases and excitations in quantum matter. Topological phenomena are intimately connected to the distribution of information content that, differently from ordinary matter, is now governed by nonlocal degrees of freedom. However, a precise characterization of how topological effects govern the complexity of a many-body state, i.e., its partition function, is presently unclear. In this paper, we show how topology and complexity are directly intertwined concepts in the context of classical statistical mechanics. We concretely present a theory that shows how the Kolmogorov complexity of a classical partition function sampling carries unique, distinctive features depending on the presence of topological excitations in the system. We confront two-dimensional Ising, Heisenberg, and XY models on several topologies and study the corresponding samplings as high-dimensional manifolds in configuration space, quantifying their complexity via the intrinsic dimension. While for the Ising and Heisenberg models the intrinsic dimension is independent of the real-space topology, for the XY model it depends crucially on temperature: across the Berezkinskii-Kosterlitz-Thouless (BKT) transition, complexity becomes topology dependent. In the BKT phase, it displays a characteristic dependence on the homology of the real-space manifold, and, for g-torii, it follows a scaling that is solely genus dependent. We argue that this behavior is intimately connected to the emergence of an order parameter in data space, the conditional connectivity, which displays scaling behavior. Our approach paves the way for an understanding of topological phenomena emergent from many-body interactions from the perspective of Kolmogorov complexity.
Data-driven discovery of statistically relevant information in quantum simulators
Verdel R., Vitale V., Panda R.K., Donkor E.D., Rodriguez A., Lannig S., Deller Y., Strobel H., Oberthaler M.K., Quantum simulators offer powerful means to investigate strongly correlated quantum matter. However, interpreting measurement outcomes in such systems poses significant challenges. Here, we present a theoretical framework for information extraction in synthetic quantum matter, illustrated for the case of a quantum quench in a spinor Bose-Einstein condensate experiment. Employing nonparametric unsupervised learning tools that provide different measures of information content, we demonstrate a theory-agnostic approach to identify dominant degrees of freedom. This enables us to rank operators according to their relevance, akin to effective field theory. To characterize the corresponding effective description, we then explore the intrinsic dimension of data sets as a measure of the complexity of the dynamics. This reveals a simplification of the data structure, which correlates with the emergence of time-dependent universal behavior in the studied system. Our assumption-free approach can be immediately applied in a variety of experimental platforms.
Spectral properties of the critical (1+1)-dimensional Abelian-Higgs model
Chanda T., The presence of gauge symmetry in 1+1 dimensions is known to be redundant, since it does not imply the existence of dynamical gauge bosons. As a consequence, in the continuum, the Abelian-Higgs model (i.e., the theory of bosonic matter interacting with photons) just possesses a single phase, as the higher-dimensional Higgs and Coulomb phases are connected via nonperturbative effects. However, recent research published in Phys. Rev. Lett. 128, 090601 (2022)0031-900710.1103/PhysRevLett.128.090601 has revealed an unexpected phase transition when the system is discretized on the lattice. This transition is described by a conformal field theory with a central charge of c=3/2. In this paper, we aim to characterize the two components of this c=3/2 theory - namely the free Majorana fermionic and bosonic parts - through equilibrium and out-of-equilibrium spectral analyses.
Scalable, ab initio protocol for quantum simulating SU(N)×U(1) Lattice Gauge Theories
Surace F.M., Fromholz P., Scazza F., We propose a protocol for the scalable quantum simulation of SU(N)×U(1) lattice gauge theories with alkaline-earth like atoms in optical lattices. The protocol exploits the combination of naturally occurring SU(N) pseudo-spin symmetry and strong inter-orbital interactions that is unique to such atomic species. A detailed ab initio study of the microscopic dynamics shows how gauge invariance emerges in an accessible parameter regime, and allows us to identify the main challenges in the simulation of such theories. We provide rigorous estimates about the requirements in terms of experimental stability in relation to observing gauge invariant dynamics in both one- and two-dimensional systems, a key element for a deeper analysis on the functioning of such class of theories in both quantum simulators and computers.
Measurement induced transitions in non-Markovian free fermion ladders
Tsitsishvili M., Poletti D., Recently there has been an intense effort to understand measurement induced transitions, but we still lack a good understanding of non-Markovian effects on these phenomena. To that end, we consider two coupled chains of free fermions, one acting as the system of interest, and one as a bath. The bath chain is subject to Markovian measurements, resulting in an effective non-Markovian dissipative dynamics acting on the system chain which is still amenable to numerical studies in terms of quantum trajectories. Within this setting, we study the entanglement within the system chain, and use it to characterize the phase diagram depending on the ladder hopping parameters and on the measurement probability. For the case of pure state evolution, the system is in an area law phase when the internal hopping of the bath chain is small, while a non-area law phase appears when the dynamics of the bath is fast. The non-area law exhibits a logarithmic scaling of the entropy compatible with a conformal phase, but also displays linear corrections for the finite system sizes we can study. For the case of mixed state evolution, we instead observe regions with both area, and non-area scaling of the entanglement negativity. We quantify the non-Markovianity of the system chain dynamics and find that for the regimes of parameters we study, a stronger non-Markovianity is associated to a larger entanglement within the system.
Many-Body Magic Via Pauli-Markov Chains - From Criticality to Gauge Theories
Tarabunga P.S., Tirrito E., Chanda T., We introduce a method to measure many-body magic in quantum systems based on a statistical exploration of Pauli strings via Markov chains. We demonstrate that sampling such Pauli-Markov chains gives ample flexibility in terms of partitions where to sample from: in particular, it enables the efficient extraction of the magic contained in the correlations between widely separated subsystems, which characterizes the nonlocality of magic. Our method can be implemented in a variety of situations. We describe an efficient sampling procedure using tree tensor networks, that exploit their hierarchical structure leading to a modest O(log N) computational scaling with system size. To showcase the applicability and efficiency of our method, we demonstrate the importance of magic in many-body systems via the following discoveries: (a) for one-dimensional systems, we show that long-range magic displays strong signatures of conformal quantum criticality (Ising, Potts, and Gaussian), overcoming the limitations of full state magic; (b) in two-dimensional Z2 lattice gauge theories, we provide conclusive evidence that magic is able to identify the confinement-deconfinement transition, and displays critical scaling behavior even at relatively modest volumes. Finally, we discuss an experimental implementation of the method, which relies only on measurements of Pauli observables.
Complexity of spin configuration dynamics due to unitary evolution and periodic projective measurements
Casagrande H.P., Xing B., 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.
First-order photon condensation in magnetic cavities: A two-leg ladder model
Bacciconi Z., Andolina G.M., Chanda T., Chiriacò G., Schirò M., We consider a model of free fermions in a ladder geometry coupled to a non-uniform cavity mode via Peierls substitution. Since the cavity mode generates a magnetic field, no-go theorems on spontaneous photon condensation do not apply, and we indeed observe a phase transition to a photon condensed phase characterized by finite circulating currents, alternatively referred to as the equilibrium superradiant phase. We consider both square and triangular ladder geometries, and characterize the transition by studying the energy structure of the system, light-matter entanglement, the properties of the photon mode, and chiral currents. The transition is of first order and corresponds to a sudden change in the fermionic band structure as well as the number of its Fermi points. Thanks to the quasi-one dimensional geometry we scrutinize the accuracy of (mean field) cavity-matter decoupling against large scale density-matrix renormalization group simulations. We find that light-matter entanglement is essential for capturing corrections to matter properties at finite sizes and for the description of the correct photon state. The latter remains Gaussian in the the thermodynamic limit both in the normal and photon condensed phases.
Diagrammatic method for many-body non-Markovian dynamics: Memory effects and entanglement transitions
Chiriacò G., Tsitsishvili M., Poletti D., Fazio R., We study the quantum dynamics of a many-body system subject to coherent evolution and coupled to a non-Markovian bath. We propose a technique to unravel the non-Markovian dynamics in terms of quantum jumps, a connection that was so far only understood for single-body systems. We develop a systematic method to calculate the probability of a quantum trajectory and formulate it in a diagrammatic structure. We find that non-Markovianity renormalizes the probability of realizing a quantum trajectory and that memory effects can be interpreted as a perturbation on top of the Markovian dynamics. We show that the diagrammatic structure is akin to that of a Dyson equation and that the probability of the trajectories can be calculated analytically. We then apply our results to study the measurement-induced entanglement transition in random unitary circuits. We find that non-Markovianity does not significantly shift the transition but stabilizes the volume law phase of the entanglement by shielding it from transient strong dissipation.
Emergence of non-Abelian SU(2) invariance in Abelian frustrated fermionic ladders
Beradze B., Tsitsishvili M., Tirrito E., We consider a system of interacting spinless fermions on a two-leg triangular ladder with π/2 magnetic flux per triangular plaquette. Microscopically, the system exhibits a U(1) symmetry corresponding to the conservation of total fermionic charge and a discrete Z2 symmetry - a product of parity transformation and chain permutation. Using bosonization, we show that, in the low-energy limit, the system is described by the quantum double-frequency sine-Gordon model. On the basis of this correspondence, a rich phase diagram of the system is obtained. It includes trivial and topological band insulators for weak interactions, separated by a Gaussian critical line, whereas at larger interactions a strongly correlated phase with spontaneously broken Z2 symmetry sets in, exhibiting a net charge imbalance and nonzero total current. At the intersection of the three phases, the system features a critical point with an emergent SU(2) symmetry. This non-Abelian symmetry, absent in the microscopic description, is realized at low energies as a combined effect of the magnetic flux, frustration, and many-body correlations. The criticality belongs to the SU(2)1 Wess-Zumino-Novikov-Witten universality class. The critical point bifurcates into two Ising critical lines that separate the band insulators from the strong-coupling symmetry broken phase. We establish an analytical connection between the low-energy description of our model around the critical bifurcation point on one hand and the Ashkin-Teller model and a weakly dimerized XXZ spin-1/2 chain on the other. We complement our field-theory understanding via tensor network simulations, providing compelling quantitative evidences of all bosonization predictions. Our findings are of interest to up-to-date cold atom experiments utilizing Rydberg dressing that have already demonstrated correlated ladder dynamics.
Classification and emergence of quantum spin liquids in chiral Rydberg models
Tarabunga P.S., Giudici G., Chanda T., We investigate the nature of quantum phases arising in chiral interacting Hamiltonians recently realized in Rydberg atom arrays. We classify all possible fermionic chiral spin liquids with U(1) global symmetry using parton construction on the honeycomb lattice. The resulting classification includes six distinct classes of gapped quantum spin liquids: the corresponding variational wavefunctions obtained from two of these classes accurately describe the Rydberg many-body ground state at 1/2 and 1/4 particle density. Complementing this analysis with tensor network simulations, we conclude that both particle filling sectors host a spin liquid with the same topological order of a ν=1/2 fractional quantum Hall effect. At density 1/2, our results clarify the phase diagram of the model, while at density 1/4, they provide an explicit construction of the ground-state wavefunction with almost unit overlap with the microscopic one. These findings pave the way to the use of parton wavefunctions to guide the discovery of quantum spin liquids in chiral Rydberg models.
Ab Initio Derivation of Lattice-Gauge-Theory Dynamics for Cold Gases in Optical Lattices
Surace F.M., Fromholz P., Oppong N.D., We introduce a method for quantum simulation of U(1) lattice gauge theories coupled to matter, utilizing alkaline-earth(-like) atoms in state-dependent optical lattices. The proposal enables the study of both gauge and fermionic matter fields without integrating out one of them in one and two dimensions. We focus on a realistic and robust implementation that utilizes the long-lived metastable clock state available in alkaline-earth(-like) atomic species. Starting from an ab initio modeling of the experimental setting, we systematically carry out a derivation of the target U(1) gauge theory. This approach allows us to identify and address conceptual and practical challenges for the implementation of lattice gauge theories that - while pivotal for a successful implementation - have never been rigorously addressed in the literature: those include the specific engineering of lattice potentials to achieve the desired structure of Wannier functions and the subtleties involved in realizing the proper separation of energy scales to enable gauge-invariant dynamics. We discuss realistic experiments that can be carried out within such a platform using the fermionic isotope 173Yb, addressing via simulations all key sources of imperfections, and provide concrete parameter estimates for relevant energy scales in both one- and two-dimensional settings.