Contributor info: Prof. Pasquale Calabrese

Filiberto Ares, Sara Murciano, Eric Vernier, Pasquale Calabrese

SciPost Phys. 15, 089 (2023) · published 12 September 2023 | Toggle abstract · pdf

We consider the quantum quench in the XX spin chain starting from a tilted Néel state which explicitly breaks the $U(1)$ symmetry of the post-quench Hamiltonian. Very surprisingly, the $U(1)$ symmetry is not restored at large time because of the activation of a non-Abelian set of charges which all break it. The breaking of the symmetry can be effectively and quantitatively characterised by the recently introduced entanglement asymmetry. By a combination of exact calculations and quasi-particle picture arguments, we are able to exactly describe the behaviour of the asymmetry at any time after the quench. Furthermore we show that the stationary behaviour is completely captured by a non-Abelian generalised Gibbs ensemble. While our computations have been performed for a non-interacting spin chain, we expect similar results to hold for the integrable interacting case as well because of the presence of non-Abelian charges also in that case.

Paola Ruggiero, Pasquale Calabrese, Thierry Giamarchi, Laura Foini

SciPost Phys. 13, 111 (2022) · published 18 November 2022 | Toggle abstract · pdf

The study of non-equilibrium dynamics of many-body systems after a quantum quench received a considerable boost and a deep theoretical understanding from the path integral formulation in imaginary time. However, the celebrated problem of a quench in the Luttinger parameter of a one dimensional quantum critical system (massless quench) has so far only been solved in the real-time Heisenberg picture. In order to bridge this theoretical gap and to understand on the same ground massive and massless quenches, we study the problem of a gaussian field characterized by a coupling parameter $K$ within a strip and a different one $K_0$ in the remaining two semi-infinite planes. We give a fully analytical solution using the electrostatic analogy with the problem of a dielectric material within a strip surrounded by an infinite medium of different dielectric constant, and exploiting the method of charge images. After analytic continuation, this solution allows us to obtain all the correlation functions after the quench within a path integral approach in imaginary time, thus recovering and generalizing the results in real time. Furthermore, this imaginary-time approach establishes a remarkable connection between the quench and the famous problem of the conductivity of a Tomonaga-Luttinger liquid coupled to two semi-infinite leads: the two are in fact related by a rotation of the spacetime coordinates.

Stefano Scopa , Pasquale Calabrese, Jérôme Dubail

SciPost Phys. 12, 207 (2022) · published 28 June 2022 | Toggle abstract · pdf

We investigate the non-equilibrium dynamics of a one-dimensional spin-1/2 XXZ model at zero-temperature in the regime $|\Delta|< 1$, initially prepared in a product state with two domain walls i.e, $|\downarrow\dots\downarrow\uparrow\dots\uparrow\downarrow\dots\downarrow\rangle$. At early times, the two domain walls evolve independently and only after a calculable time a non-trivial interplay between the two emerges and results in the occurrence of a split Fermi sea. For $\Delta=0$, we derive exact asymptotic results for the magnetization and the spin current by using a semi-classical Wigner function approach, and we exactly determine the spreading of entanglement entropy exploiting the recently developed tools of quantum fluctuating hydrodynamics. In the interacting case, we analytically solve the Generalized Hydrodynamics equation providing exact expressions for the conserved quantities. We display some numerical results for the entanglement entropy also in the interacting case and we propose a conjecture for its asymptotic value.

Vittorio Vitale, Andreas Elben, Richard Kueng, Antoine Neven, Jose Carrasco, Barbara Kraus, Peter Zoller, Pasquale Calabrese, Benoit Vermersch, Marcello Dalmonte

SciPost Phys. 12, 106 (2022) · published 25 March 2022 | Toggle abstract · pdf

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.

David X. Horvath, Pasquale Calabrese, Olalla A. Castro-Alvaredo

SciPost Phys. 12, 088 (2022) · published 10 March 2022 | Toggle abstract · pdf

In this paper we continue the programme initiated in Part I, that is the study of entanglement measures in the sine-Gordon model. In both parts, we have focussed on one specific technique, that is the well-known connection between branch point twist field correlators and measures of entanglement in 1+1D integrable quantum field theory. Our papers apply this technique for the first time to a non-diagonal theory with an involved particle spectrum, the sine-Gordon model. In this Part II we focus on a different entanglement measure, the symmetry resolved entanglement, and develop its associated twist field description, exploiting the underlying U(1) symmetry of the theory. In this context, conventional branch point twist fields are no longer the fields required, but instead we must work with one of their composite generalisations, which can be understood as the field resulting from the fusion of a standard branch point twist field and the sine-Gordon exponential field associated with U(1) symmetry. The resulting composite twist field has correlators which as usual admit a form factor expansion. In this paper we write the associated form factor equations and solve them for various examples in the breather sector by using the method of angular quantisation. We show that, in the attractive regime, this is the sector which provides the leading contribution to the symmetry resolved entropies, both Renyi and von Neumann. We compute the latter in the limit of a large region size and show that they satisfy the property of equipartition, that is the leading contribution to the symmetry resolved entanglement is independent of the symmetry sector.

Elena Tartaglia, Pasquale Calabrese, Bruno Bertini

SciPost Phys. 12, 028 (2022) · published 18 January 2022 | Toggle abstract · pdf

We consider the real-time evolution of the Hubbard model in the limit of infinite coupling. In this limit the Hamiltonian of the system is mapped into a number-conserving quadratic form of spinless fermions, i.e. the tight binding model. The relevant local observables, however, do not transform well under this mapping and take very complicated expressions in terms of the spinless fermions. Here we show that for two classes of interesting observables the quench dynamics from product states in the occupation basis can be determined exactly in terms of correlations in the tight-binding model. In particular, we show that the time evolution of any function of the total density of particles is mapped directly into that of the same function of the density of spinless fermions in the tight-binding model. Moreover, we express the two-point functions of the spin-full fermions at any time after the quench in terms of correlations of the tight binding model. This sum is generically very complicated but we show that it leads to simple explicit expressions for the time evolution of the densities of the two separate species and the correlations between a point at the boundary and one in the bulk when evolving from the so called generalised nested Néel states.

Paola Ruggiero, Pasquale Calabrese, Laura Foini, Thierry Giamarchi

SciPost Phys. 11, 055 (2021) · published 13 September 2021 | Toggle abstract · pdf

We study the quantum quench in two coupled Tomonaga-Luttinger Liquids (TLLs), from the off-critical to the critical regime, relying on the conformal field theory approach and the known solutions for single TLLs. We consider a squeezed form of the initial state, whose low energy limit is fixed in a way to describe a massive and a massless mode, and we encode the non-equilibrium dynamics in a proper rescaling of the time. In this way, we compute several correlation functions, which at leading order factorize into multipoint functions evaluated at different times for the two modes. Depending on the observable, the contribution from the massive or from the massless mode can be the dominant one, giving rise to exponential or power-law decay in time, respectively. Our results find a direct application in all the quench problems where, in the scaling limit, there are two independent massless fields: these include the Hubbard model, the Gaudin-Yang gas, and tunnel-coupled tubes in cold atoms experiments.

Sara Murciano, Riccarda Bonsignori, Pasquale Calabrese

SciPost Phys. 10, 111 (2021) · published 19 May 2021 | Toggle abstract · pdf

We consider the problem of symmetry decomposition of the entanglement negativity in free fermionic systems. Rather than performing the standard partial transpose, we use the partial time-reversal transformation which naturally encodes the fermionic statistics. The negativity admits a resolution in terms of the charge imbalance between the two subsystems. We introduce a normalised version of the imbalance resolved negativity which has the advantage to be an entanglement proxy for each symmetry sector, but may diverge in the limit of pure states for some sectors. Our main focus is then the resolution of the negativity for a free Dirac field at finite temperature and size. We consider both bipartite and tripartite geometries and exploit conformal field theory to derive universal results for the charge imbalance resolved negativity. To this end, we use a geometrical construction in terms of an Aharonov-Bohm-like flux inserted in the Riemann surface defining the entanglement. We interestingly find that the entanglement negativity is always equally distributed among the different imbalance sectors at leading order. Our analytical findings are tested against exact numerical calculations for free fermions on a lattice.

Pasquale Calabrese

SciPost Phys. Lect. Notes 20 (2020) · published 1 December 2020 | Toggle abstract · pdf

These are lecture notes for a short course given at the Les Houches Summer School on ``Integrability in Atomic and Condensed Matter Physics'', in summer 2018. Here, I pedagogically discuss recent advances in the study of the entanglement spreading during the non-equilibrium dynamics of isolated integrable quantum systems. I first introduce the idea that the stationary thermodynamic entropy is the entanglement accumulated during the non-equilibrium dynamics and then join such an idea with the quasiparticle picture for the entanglement spreading to provide quantitive predictions for the time evolution of the entanglement entropy in arbitrary integrable models, regardless of the interaction strength.

Jiaju Zhang, Pasquale Calabrese, Marcello Dalmonte, M. A. Rajabpour

SciPost Phys. Core 2, 007 (2020) · published 7 May 2020 | Toggle abstract · pdf

We carry out a comprehensive comparison between the exact modular Hamiltonian and the lattice version of the Bisognano-Wichmann (BW) one in one-dimensional critical quantum spin chains. As a warm-up, we first illustrate how the trace distance provides a more informative mean of comparison between reduced density matrices when compared to any other Schatten $n$-distance, normalized or not. In particular, as noticed in earlier works, it provides a way to bound other correlation functions in a precise manner, i.e., providing both lower and upper bounds. Additionally, we show that two close reduced density matrices, i.e. with zero trace distance for large sizes, can have very different modular Hamiltonians. This means that, in terms of describing how two states are close to each other, it is more informative to compare their reduced density matrices rather than the corresponding modular Hamiltonians. After setting this framework, we consider the ground states for infinite and periodic XX spin chain and critical Ising chain. We provide robust numerical evidence that the trace distance between the lattice BW reduced density matrix and the exact one goes to zero as $\ell^{-2}$ for large length of the interval $\ell$. This provides strong constraints on the difference between the corresponding entanglement entropies and correlation functions. Our results indicate that discretized BW reduced density matrices reproduce exact entanglement entropies and correlation functions of local operators in the limit of large subsystem sizes. Finally, we show that the BW reduced density matrices fall short of reproducing the exact behavior of the logarithmic emptiness formation probability in the ground state of the XX spin chain.

Sara Murciano, Giuseppe Di Giulio, Pasquale Calabrese

SciPost Phys. 8, 046 (2020) · published 26 March 2020 | Toggle abstract · pdf

We study the symmetry resolved entanglement entropies in gapped integrable lattice models. We use the corner transfer matrix to investigate two prototypical gapped systems with a U(1) symmetry: the complex harmonic chain and the XXZ spin-chain. While the former is a free bosonic system, the latter is genuinely interacting. We focus on a subsystem being half of an infinitely long chain. In both models, we obtain exact expressions for the charged moments and for the symmetry resolved entropies. While for the spin chain we found exact equipartition of entanglement (i.e. all the symmetry resolved entropies are the same), this is not the case for the harmonic system where equipartition is effectively recovered only in some limits. Exploiting the gaussianity of the harmonic chain, we also develop an exact correlation matrix approach to the symmetry resolved entanglement that allows us to test numerically our analytic results.

Hassan Shapourian, Paola Ruggiero, Shinsei Ryu, Pasquale Calabrese

SciPost Phys. 7, 037 (2019) · published 25 September 2019 | Toggle abstract · pdf

A basic diagnostic of entanglement in mixed quantum states is known as the positive partial transpose (PT) criterion. Such criterion is based on the observation that the spectrum of the partially transposed density matrix of an entangled state contains negative eigenvalues, in turn, used to define an entanglement measure called the logarithmic negativity. Despite the great success of logarithmic negativity in characterizing bosonic many-body systems, generalizing the operation of PT to fermionic systems remained a technical challenge until recently when a more natural definition of PT for fermions that accounts for the Fermi statistics has been put forward. In this paper, we study the many-body spectrum of the reduced density matrix of two adjacent intervals for one-dimensional free fermions after applying the fermionic PT. We show that in general there is a freedom in the definition of such operation which leads to two different definitions of PT: the resulting density matrix is Hermitian in one case, while it becomes pseudo-Hermitian in the other case. Using the path-integral formalism, we analytically compute the leading order term of the moments in both cases and derive the distribution of the corresponding eigenvalues over the complex plane. We further verify our analytical findings by checking them against numerical lattice calculations.

Alvise Bastianello, Pasquale Calabrese

SciPost Phys. 5, 033 (2018) · published 15 October 2018 | Toggle abstract · pdf

We extend the semiclassical picture for the spreading of entanglement and correlations to quantum quenches with several species of quasiparticles that have non-trivial pair correlations in momentum space. These pair correlations are, for example, relevant in inhomogeneous lattice models with a periodically-modulated Hamiltonian parameter. We provide explicit predictions for the spreading of the entanglement entropy in the space-time scaling limit. We also predict the time evolution of one- and two-point functions of the order parameter for quenches within the ordered phase. We test all our predictions against exact numerical results for quenches in the Ising chain with a modulated transverse field and we find perfect agreement.

Stefan Groha, Fabian H. L. Essler, Pasquale Calabrese

SciPost Phys. 4, 043 (2018) · published 28 June 2018 | Toggle abstract · pdf

We consider the full probability distribution for the transverse magnetization of a finite subsystem in the transverse field Ising chain. We derive a determinant representation of the corresponding characteristic function for general Gaussian states. We consider applications to the full counting statistics in the ground state, finite temperature equilibrium states, non-equilibrium steady states and time evolution after global quantum quenches. We derive an analytical expression for the time and subsystem size dependence of the characteristic function at sufficiently late times after a quantum quench. This expression features an interesting multiple light-cone structure.

Vincenzo Alba, Pasquale Calabrese

SciPost Phys. 4, 017 (2018) · published 27 March 2018 | Toggle abstract · pdf

The time evolution of the entanglement entropy in non-equilibrium quantum systems provides crucial information about the structure of the time-dependent state. For quantum quench protocols, by combining a quasiparticle picture for the entanglement spreading with the exact knowledge of the stationary state provided by Bethe ansatz, it is possible to obtain an exact and analytic description of the evolution of the entanglement entropy. Here we discuss the application of these ideas to several integrable models. First we show that for non-interacting systems, both bosonic and fermionic, the exact time-dependence of the entanglement entropy can be derived by elementary techniques and without solving the dynamics. We then provide exact results for interacting spin chains that are carefully tested against numerical simulations. Finally, we apply this method to integrable one-dimensional Bose gases (Lieb-Liniger model) both in the attractive and repulsive regimes. We highlight a peculiar behaviour of the entanglement entropy due to the absence of a maximum velocity of excitations.

Jérôme Dubail, Jean-Marie Stéphan, Pasquale Calabrese

SciPost Phys. 3, 019 (2017) · published 6 September 2017 | Toggle abstract · pdf

The light-cone spreading of entanglement and correlation is a fundamental and ubiquitous feature of homogeneous extended quantum systems. Here we point out that a class of inhomogenous Luttinger liquids (those with a uniform Luttinger parameter $K$) at low energy display the universal phenomenon of curved light cones: gapless excitations propagate along the geodesics of the metric $ds^2=dx^2+v(x)^2 d\tau^2$, with $v(x)$ being the calculable spatial dependent velocity induced by the inhomogeneity. We confirm our findings with explicit analytic and numerical calculations both in- and out-of-equilibrium for a Tonks-Girardeau gas in a harmonic potential and in lattice systems with artificially tuned hamiltonian density.

Jérôme Dubail, Jean-Marie Stéphan, Jacopo Viti, Pasquale Calabrese

SciPost Phys. 2, 002 (2017) · published 13 February 2017 | Toggle abstract · pdf

Conformal field theory (CFT) has been extremely successful in describing large-scale universal effects in one-dimensional (1D) systems at quantum critical points. Unfortunately, its applicability in condensed matter physics has been limited to situations in which the bulk is uniform because CFT describes low-energy excitations around some energy scale, taken to be constant throughout the system. However, in many experimental contexts, such as quantum gases in trapping potentials and in several out-of-equilibrium situations, systems are strongly inhomogeneous. We show here that the powerful CFT methods can be extended to deal with such 1D situations, providing a few concrete examples for non-interacting Fermi gases. The system's inhomogeneity enters the field theory action through parameters that vary with position; in particular, the metric itself varies, resulting in a CFT in curved space. This approach allows us to derive exact formulas for entanglement entropies which were not known by other means.

Lorenzo Piroli, Pasquale Calabrese, Fabian H. L. Essler

SciPost Phys. 1, 001 (2016) · published 14 September 2016 | Toggle abstract · pdf

We study quantum quenches to the one-dimensional Bose gas with attractive interactions in the case when the initial state is an ideal one-dimensional Bose condensate. We focus on properties of the stationary state reached at late times after the quench. This displays a finite density of multi-particle bound states, whose rapidity distribution is determined exactly by means of the quench action method. We discuss the relevance of the multi-particle bound states for the physical properties of the system, computing in particular the stationary value of the local pair correlation function $g_2$.