All Publications

In the context of ground states of quantum many-body systems, the locality of entanglement between connected regions of space is directly tied to the locality of the corresponding entanglement Hamiltonian: the latter is dominated by local, few-body terms. In this work, we introduce the negativity Hamiltonian as the (non-Hermitian) effective Hamiltonian operator describing the logarithm of the partial transpose of a many-body system. This allows us to address the connection between entanglement and operator locality beyond the paradigm of bipartite pure systems. As a first step in this direction, we study the structure of the negativity Hamiltonian for fermionic conformal field theories and a free-fermion chain: in both cases, we show that the negativity Hamiltonian assumes a quasilocal functional form, that is captured by simple functional relations.

We explore the dynamical phases of unitary Clifford circuits with variable-range interactions, coupled to a monitoring environment. We investigate two classes of models, distinguished by the action of the unitary gates, which either are organized in clusters of finite-range two-body gates, or are pair-wise interactions randomly distributed throughout the system with a power-law distribution. We find the range of the interactions plays a key role in characterizing both phases and their measurement-induced transitions. For the cluster unitary gates we find a transition between a phase with volume-law scaling of the entanglement entropy and a phase with area-law entanglement entropy. Our results indicate that the universality class of the phase transition is compatible to that of short range hybrid Clifford circuits. Oppositely, in the case of power-law distributed gates, we find the universality class of the phase transition changes continuously with the parameter controlling the range of interactions. In particular, for intermediate values of the control parameter, we find a non-conformal critical line which separates a phase with volume-law scaling of the entanglement entropy from one with sub-extensive scaling. Within this region, we find the entanglement entropy and the logarithmic negativity present a cross-over from a phase with algebraic growth of entanglement with system size, and an area-law phase.

When a quantum system initialized in a product state is subjected to either coherent or incoherent dynamics, the entropy of any of its connected partitions generically increases as a function of time, signalling the inevitable spreading of (quantum) information throughout the system. Here, we show that, in the presence of continuous symmetries and under ubiquitous experimental conditions, symmetry-resolved information spreading is inhibited due to the competition of coherent and incoherent dynamics: in given quantum number sectors, entropy decreases as a function of time, signalling dynamical purification. Such dynamical purification bridges between two distinct short and intermediate time regimes, characterized by a log-volume and log-area entropy law, respectively. It is generic to symmetric quantum evolution, and as such occurs for different partition geometry and topology, and classes of (local) Liouville dynamics. We then develop a protocol to measure symmetry-resolved entropies and negativities in synthetic quantum systems based on the random unitary toolbox, and demonstrate the generality of dynamical purification using experimental data from trapped ion experiments [Brydges et al., Science 364, 260 (2019)]. Our work shows that symmetry plays a key role as a magnifying glass to characterize many-body dynamics in open quantum systems, and, in particular, in noisy-intermediate scale quantum devices.

In quantum statistical mechanics, finite-temperature phase transitions are typically governed by classical field theories. In this context, the role of quantum correlations is unclear: recent contributions have shown how entanglement is typically very short-ranged, and thus uninformative about long-ranged critical correlations. In this work, we show the existence of finite-temperature phase transitions where a broader form of quantum correlation than entanglement, the entropic quantum discord, can display genuine signatures of critical behavior. We consider integrable bosonic field theories in both two- and three-dimensional lattices, and show how the two-mode Gaussian discord decays algebraically with the distance even in cases where the entanglement negativity vanishes beyond nearest-neighbor separations. Systematically approaching the zero-temperature limit allows us to connect discord to entanglement, drawing a generic picture of quantum correlations and critical behavior that naturally describes the transition between entangled and discordant quantum matter.

Quantum systems evolving unitarily and subject to quantum measurements exhibit various types of non-equilibrium phase transitions, arising from the competition between unitary evolution and measurements. Dissipative phase transitions in steady states of time-independent Liouvillians and measurement induced phase transitions at the level of quantum trajectories are two primary examples of such transitions. Investigating a many-body spin system subject to periodic resetting measurements, we argue that many-body dissipative Floquet dynamics provides a natural framework to analyze both types of transitions. We show that a dissipative phase transition between a ferromagnetic ordered phase and a paramagnetic disordered phase emerges for long-range systems as a function of measurement probabilities. A measurement induced transition of the entanglement entropy between volume law scaling and sub-volume law scaling is also present, and is distinct from the ordering transition. The two phases correspond to an error-correcting and a quantum-Zeno regimes, respectively. The ferromagnetic phase is lost for short range interactions, while the volume law phase of the entanglement is enhanced. An analysis of multifractal properties of wave function in Hilbert space provides a common perspective on both types of transitions in the system. Our findings are immediately relevant to trapped ion experiments, for which we detail a blueprint proposal based on currently available platforms.

We propose an ordered set of experimentally accessible conditions for detecting entanglement in mixed states. The k-th condition involves comparing moments of the partially transposed density operator up to order k. Remarkably, the union of all moment inequalities reproduces the Peres-Horodecki criterion for detecting entanglement. Our empirical studies highlight that the first four conditions already detect mixed state entanglement reliably in a variety of quantum architectures. Exploiting symmetries can help to further improve their detection capabilities. We also show how to estimate moment inequalities based on local random measurements of single state copies (classical shadows) and derive statistically sound confidence intervals as a function of the number of performed measurements. Our analysis includes the experimentally relevant situation of drifting sources, i.e. non-identical, but independent, state copies.

Learning the structure of the entanglement Hamiltonian (EH) is central to characterizing quantum many-body states in analog quantum simulation. We describe a protocol where spatial deformations of the many-body Hamiltonian, physically realized on the quantum device, serve as an efficient variational ansatz for a local EH. Optimal variational parameters are determined in a feedback loop, involving quench dynamics with the deformed Hamiltonian as a quantum processing step, and classical optimization. We simulate the protocol for the ground state of Fermi-Hubbard models in quasi-1D geometries, finding excellent agreement of the EH with Bisognano-Wichmann predictions. Subsequent on-device spectroscopy enables a direct measurement of the entanglement spectrum, which we illustrate for a Fermi Hubbard model in a topological phase.

We study the phase diagram and critical properties of quantum Ising chains with longrange ferromagnetic interactions decaying in a power-law fashion with exponent α, in regimes of direct interest for current trapped ion experiments. Using large-scale path integral Monte Carlo simulations, we investigate both the ground-state and the nonzerotemperature regimes. We identify the phase boundary of the ferromagnetic phase and obtain accurate estimates for the ferromagnetic-paramagnetic transition temperatures. We further determine the critical exponents of the respective transitions. Our results are in agreement with existing predictions for interaction exponents α > 1 up to small deviations in some critical exponents. We also address the elusive regime α < 1, where we find that the universality class of both the ground-state and nonzero-temperature transition is consistent with the mean-field limit at α = 0. Our work not only contributes to the understanding of the equilibrium properties of long-range interacting quantum Ising models, but can also be important for addressing fundamental dynamical aspects, such as issues concerning the open question of thermalization in such models.

We study the impact of quenched disorder on the dynamics of locally constrained quantum spin chains, that describe 1D arrays of Rydberg atoms in both the frozen (Ising-type) and dressed (XY-type) regime. Performing large-scale numerical experiments, we observe no trace of many-body localization even at large disorder. Analyzing the role of quenched disorder terms in constrained systems we show that they act in two, distinct and competing ways: as an on-site disorder term for the basic excitations of the system, and as an interaction between excitations. The two contributions are of the same order, and as they compete (one towards localization, the other against it), one does never enter a truly strong disorder, weak interaction limit, where many-body localization occurs. Such a mechanism is further clarified in the case of XY-type constrained models: there, a term which would represent a bona fide local quenched disorder term acting on the excitations of the clean model must be written as a series of nonlocal terms in the unconstrained variables. Our observations provide a simple picture to interpret the role of quenched disorder that could be immediately extended to other constrained models or quenched gauge theories.

Quantum many-body systems are characterized by patterns of correlations defining highly nontrivial manifolds when interpreted as data structures. Physical properties of phases and phase transitions are typically retrieved via correlation functions, that are related to observable response functions. Recent experiments have demonstrated capabilities to fully characterize quantum many-body systems via wave-function snapshots, opening new possibilities to analyze quantum phenomena. Here, we introduce a method to data mine the correlation structure of quantum partition functions via their path integral (or equivalently, stochastic series expansion) manifold. We characterize path-integral manifolds generated via state-of-the-art quantum Monte Carlo methods utilizing the intrinsic dimension (ID) and the variance of distances between nearest-neighbor (NN) configurations: the former is related to data-set complexity, while the latter is able to diagnose connectivity features of points in configuration space. We show how these properties feature universal patterns in the vicinity of quantum criticality, that reveal how data structures simplify systematically at quantum phase transitions. This is further reflected by the fact that both ID and variance of NN distances exhibit universal scaling behavior in the vicinity of second-order and Berezinskii-Kosterlitz-Thouless critical points. Finally, we show how non-Abelian symmetries dramatically influence quantum data sets, due to the nature of (noncommuting) conserved charges in the quantum case. Complementary to neural-network representations, our approach represents a first elementary step towards a systematic characterization of path-integral manifolds before any dimensional reduction is taken, that is informative about universal behavior and complexity, and can find immediate application to both experiments and Monte Carlo simulations.

We investigate measurement-induced phase transitions in the quantum Ising chain coupled to a monitoring environment. We compare two different limits of the measurement problem: the stochastic quantum-state diffusion protocol corresponding to infinite small jumps per unit of time and the no-click limit, corresponding to postselection and described by a non-Hermitian Hamiltonian. In both cases we find a remarkably similar phenomenology as the measurement strength γ is increased, namely, a sharp transition from a critical phase with logarithmic scaling of the entanglement to an area-law phase, which occurs at the same value of the measurement rate in the two protocols. An effective central charge, extracted from the logarithmic scaling of the entanglement, vanishes continuously at the common transition point, although with different critical behavior possibly suggesting different universality classes for the two protocols. We interpret the central charge mismatch near the transition in terms of noise-induced disentanglement, as suggested by the entanglement statistics which displays emergent bimodality upon approaching the critical point. The non-Hermitian Hamiltonian and its associated subradiance spectral transition provide a natural framework to understand both the extended critical phase, emerging here for a model which lacks any continuous symmetry, and the entanglement transition into the area law.

Quantum scars are nonthermal eigenstates characterized by low entanglement entropy, initially detected in systems subject to nearest-neighbor Rydberg blockade, the so-called PXP model. While most of these special eigenstates elude an analytical description and seem to hybridize with nearby thermal eigenstates for large systems, some of them can be written as matrix product states with size-independent bond dimension. We study the response of these exact quantum scars to perturbations by analyzing the scaling of the fidelity susceptibility with system size. We find that some of them are anomalously stable at first order in perturbation theory, in sharp contrast to the eigenstate thermalization hypothesis. However, this stability seems to break down when all orders are taken into account. We further investigate models with larger blockade radius and find a set of exact quantum scars that we write down analytically and compare with the PXP exact eigenstates. We show that they exhibit the same robustness against perturbations at first order.

The identification of universal properties from minimally processed data sets is one goal of machine learning techniques applied to statistical physics. Here, we study how the minimum number of variables needed to accurately describe the important features of a data set - the intrinsic dimension (Id) - behaves in the vicinity of phase transitions. We employ state-of-the-art nearest-neighbors-based Id estimators to compute the Id of raw Monte Carlo thermal configurations across different phase transitions: first-order, second-order, and Berezinskii-Kosterlitz-Thouless. For all the considered cases, we find that the Id uniquely characterizes the transition regime. The finite-size analysis of the Id allows us to not only identify critical points with an accuracy comparable to methods that rely on a priori identification of order parameters but also to determine the corresponding (critical) exponent ν in the case of continuous transitions. For the case of topological transitions, this analysis overcomes the reported limitations affecting other unsupervised learning methods. Our work reveals how raw data sets display unique signatures of universal behavior in the absence of any dimensional reduction scheme and suggest direct parallelism between conventional order parameters in real space and the intrinsic dimension in the data space.

We propose a method for detecting bipartite entanglement in a many-body mixed state based on estimating moments of the partially transposed density matrix. The estimates are obtained by performing local random measurements on the state, followed by postprocessing using the classical shadows framework. Our method can be applied to any quantum system with single-qubit control. We provide a detailed analysis of the required number of experimental runs, and demonstrate the protocol using existing experimental data [Brydges et al., Science 364, 260 (2019)SCIEAS0036-807510.1126/science.aau4963].

In our paper, we incorrectly reported that singular value decomposition (SVD) directly evaluates the entanglement entropy of a stabilizer state generated by the dynamics we are interested in. In fact, SVD only provides a rigorous upper bound to the entropy. 1 We are grateful to Y. Li and M. P. A. Fisher for correspondence that elucidated this aspect. Specifically, in the published version of the paper we compute entanglement of the stabilizer state by means of the Hamma-Ionicioiu-Zanardi theorem. This requires the computation of the matrix rank. SVD computes the value of the rank in the field of real numbers R, while due to the algebraic structure of the stabilizer group, we should have computed the rank for the field F2.We observe that any matrix A with binary elements aij=0,1 satisfies.

We study the spectral properties of D-dimensional N = 2 supersymmetric lattice models. We find systematic departures from the eigenstate thermalization hypothesis (ETH) in the form of a degenerate set of ETH-violating supersymmetric (SUSY) doublets, also referred to as many-body scars, that we construct analytically. These states are stable against arbitrary SUSY-preserving perturbations, including inhomogeneous couplings. For the specific case of two-leg ladders, we provide extensive numerical evidence that shows how those states are the only ones violating the ETH, and discuss their robustness to SUSY-violating perturbations. Our work suggests a generic mechanism to stabilize quantum many-body scars in lattice models in arbitrary dimensions.

We study the disconnected entanglement entropy, SD, of the Su-Schrieffer-Heeger model. SD is a combination of both connected and disconnected bipartite entanglement entropies that removes all area and volume law contributions and is thus only sensitive to the non-local entanglement stored within the ground state manifold. Using analytical and numerical computations, we show that SD behaves like a topological invariant, i.e., it is quantized to either 0 or 2 log(2) in the topologically trivial and non-trivial phases, respectively. These results also hold in the presence of symmetry-preserving disorder. At the second-order phase transition separating the two phases, SD displays a finite-size scaling behavior akin to those of conventional order parameters, that allows us to compute entanglement critical exponents. To corroborate the topological origin of the quantized values of SD, we show how the latter remain quantized after applying unitary time evolution in the form of a quantum quench, a characteristic feature of topological invariants associated with particle-hole symmetry.

We show how U(1) lattice gauge theories display key signatures of ergodicity breaking in the presence of a random charge background. We argue that, in such gauge theories, there is a cooperative effect of disorder and interactions in favoring ergodicity breaking: This is due to the confining nature of the Coulomb potential, which suppresses the number of available energy resonances at all distances. Such a cooperative mechanism reflects into very modest finite-volume effects: This allows us to draw a sharp boundary for the ergodic regime, and thus the breakdown of quantum chaos for sufficiently strong gauge couplings, at system sizes accessible via exact diagonalization. Our conclusions are independent on the value of a background topological angle, and are contrasted with a gauge theory with truncated Hilbert space, where instead we observe very strong finite-volume effects akin to those observed in spin chains.

Abstract: Lattice gauge theories, which originated from particle physics in the context of Quantum Chromodynamics (QCD), provide an important intellectual stimulus to further develop quantum information technologies. While one long-term goal is the reliable quantum simulation of currently intractable aspects of QCD itself, lattice gauge theories also play an important role in condensed matter physics and in quantum information science. In this way, lattice gauge theories provide both motivation and a framework for interdisciplinary research towards the development of special purpose digital and analog quantum simulators, and ultimately of scalable universal quantum computers. In this manuscript, recent results and new tools from a quantum science approach to study lattice gauge theories are reviewed. Two new complementary approaches are discussed: first, tensor network methods are presented – a classical simulation approach – applied to the study of lattice gauge theories together with some results on Abelian and non-Abelian lattice gauge theories. Then, recent proposals for the implementation of lattice gauge theory quantum simulators in different quantum hardware are reported, e.g., trapped ions, Rydberg atoms, and superconducting circuits. Finally, the first proof-of-principle trapped ions experimental quantum simulations of the Schwinger model are reviewed. Graphical abstract: [Figure not available: see fulltext.].

We investigate the dynamics of two-dimensional quantum spin systems under the combined effect of random unitary gates and local projective measurements. When considering steady states, a measurement-induced transition occurs between two distinct dynamical phases, one characterized by a volume-law scaling of entanglement entropy, the other by an area law. Employing stabilizer states and Clifford random unitary gates, we numerically investigate square lattices of linear dimension up to L=48 for two distinct measurement protocols. For both protocols, we observe a transition point where the dominant contribution in the entanglement entropy displays multiplicative logarithmic violations to the area law. We obtain estimates of the correlation length critical exponent at the percent level; these estimates suggest universal behavior and are incompatible with the universality class of 3D percolation.