All Publications

In this paper, we apply the form factor bootstrap approach to branch point twist fields in the q-state Potts model for q ≤ 3. For q = 3 this is an integrable interacting quantum field theory with an internal discrete ℤ3 symmetry and therefore provides an ideal starting point for the investigation of the symmetry resolved entanglement entropies. However, more generally, for q ≤ 3 the standard Rényi and entanglement entropies are also accessible through the bootstrap programme. In our work we present form factor solutions both for the standard branch point twist field with q ≤ 3 and for the composite (or symmetry resolved) branch point twist field with q = 3. In both cases, the form factor equations are solved for two particles and the solutions are carefully checked via the ∆-sum rule. Using our analytic predictions, we compute the leading finite-size corrections to the entanglement entropy and entanglement equipartition for a single interval in the ground state.

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 consider the quantum quench dynamics of a Heisenberg-Ising spin ladder which is an archetypal model in which confinement of elementary excitations is triggered by internal interactions rather than an external field. We show that the confinement strongly affects the light cone structure of correlation functions providing signatures of the velocities of the mesons of the model. We also show that the meson masses can be measured from the real time analysis of the evolution of the order parameter.

The confinement of elementary excitations induces distinctive features in the non-equilibrium quench dynamics. One of the most remarkable is the suppression of entanglement entropy, which in several instances turns out to oscillate rather than grow indefinitely. While the qualitative physical origin of this behavior is clear, till now no quantitative understanding away from the field theory limit was available. Here we investigate this problem in the weak quench limit, when mesons are excited at rest, hindering entropy growth and exhibiting persistent oscillations. We provide analytical predictions of the entire entanglement dynamics based on a Gaussian approximation of the many-body state, which captures numerical data with great accuracy and is further simplified to a semiclassical quasiparticle picture in the regime of weak confinement. Our methods are valid in general and we apply explicitly to two prototypical models: the Ising chain in a tilted field and the experimentally relevant long-range Ising model.

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 this paper we continue the program 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 Rényi 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.

We apply the theory of quantum generalized hydrodynamics (QGHD) introduced in (2020 Phys. Rev. Lett. 124 140603) to derive asymptotically exact results for the density fluctuations and the entanglement entropy of a one-dimensional trapped Bose gas in the Tonks-Girardeau (TG) or hard-core limit, after a trap quench from a double well to a single well. On the analytical side, the quadratic nature of the theory of QGHD is complemented with the emerging conformal invariance at the TG point to fix the universal part of those quantities. Moreover, the well-known mapping of hard-core bosons to free fermions, allows to use a generalized form of the Fisher-Hartwig conjecture to fix the non-trivial spacetime dependence of the ultraviolet cutoff in the entanglement entropy. The free nature of the TG gas also allows for more accurate results on the numerical side, where a higher number of particles as compared to the interacting case can be simulated. The agreement between analytical and numerical predictions is extremely good. For the density fluctuations, however, one has to average out large Friedel oscillations present in the numerics to recover such agreement.

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.

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.

The false vacuum decay has been a central theme in physics for half a century with applications to cosmology and to the theory of fundamental interactions. This fascinating phenomenon is even more intriguing when combined with the confinement of elementary particles. Due to the astronomical timescales involved, the research has so far focused on theoretical aspects of this decay. The purpose of this Letter is to show that the false vacuum decay is accessible to current optical experiments as quantum analog simulators of spin chains with confinement of the elementary excitations, which mimic the high energy phenomenology but in one spatial dimension. We study the nonequilibrium dynamics of the false vacuum in a quantum Ising chain and in an XXZ ladder. The false vacuum is the metastable state that arises in the ferromagnetic phase of the model when the symmetry is explicitly broken by a longitudinal field. This state decays through the formation of "bubbles"of true vacuum. Using infinite volume time evolving block decimation (iTEBD) simulations, we are able to study the real-time evolution in the thermodynamic limit and measure the decay rate of local observables. We find that the numerical results agree with the theoretical prediction that the decay rate is exponentially small in the inverse of the longitudinal field.

We consider the non-equilibrium dynamics of the entanglement entropy of a one-dimensional quantum gas of hard-core particles, initially confined in a box potential at zero temperature. At t = 0 the right edge of the box is suddenly released and the system is let free to expand. During this expansion, the initially correlated region propagates with a non-homogeneous profile, leading to the growth of entanglement entropy. This setting is investigated in the hydrodynamic regime, with tools stemming from semi-classical Wigner function approach andwith recent developments of quantum fluctuating hydrodynamics. Within this framework, the entanglement entropy can be associated to a correlation function of chiral twist-fields of the conformal field theory that lives along the Fermi contour and it can be exactly determined. Our predictions for the entanglement evolution are found in agreement with and generalize previous results in literature based on numerical calculations and heuristic arguments.

We develop a systematic approach to compute the subsystem trace distances and relative entropies for subsystem reduced density matrices associated to excited states in different symmetry sectors of a 1+1 dimensional conformal field theory having an internal U(1) symmetry. We provide analytic expressions for the charged moments corresponding to the resolution of both relative entropies and distances for general integer n. For the relative entropies, these formulas are manageable and the analytic continuation to n = 1 can be worked out in most of the cases. Conversely, for the distances the corresponding charged moments become soon untreatable as n increases. A remarkable result is that relative entropies and distances are the same for all symmetry sectors, i.e. they satisfy entanglement equipartition, like the entropies. Moreover, we exploit the OPE expansion of composite twist fields, to provide very general results when the subsystem is a single interval much smaller than the total system. We focus on the massless compact boson and our results are tested against exact numerical calculations in the XX spin chain.

We consider the problem of the decomposition of the Rényi entanglement entropies in theories with a non-abelian symmetry by doing a thorough analysis of Wess-Zumino-Witten (WZW) models. We first consider SU(2)k as a case study and then generalise to an arbitrary non-abelian Lie group. We find that at leading order in the subsystem size L the entanglement is equally distributed among the different sectors labelled by the irreducible representation of the associated algebra. We also identify the leading term that breaks this equipartition: it does not depend on L but only on the dimension of the representation. Moreover, a log log L contribution to the Rényi entropies exhibits a universal prefactor equal to half the dimension of the Lie group.

We revisit early suggestions to observe spin-charge separation (SCS) in cold-atom settings in the time domain by studying one-dimensional repulsive Fermi gases in a harmonic potential, where pulse perturbations are initially created at the center of the trap. We analyze the subsequent evolution using generalized hydrodynamics (GHD), which provide an exact description, at large space-time scales, for arbitrary temperature T, particle density, and interactions. At T=0 and vanishing magnetic field, we find that, after a nontrivial transient regime, spin and charge dynamically decouple up to perturbatively small corrections which we quantify. In this limit, our results can be understood based on a simple phase-space hydrodynamic picture. At finite temperature, we solve numerically the GHD equations, showing that for low T>0 effects of SCS survive and characterize explicitly the value of T for which the two distinguishable excitations melt.

The study of the entanglement dynamics plays a fundamental role in understanding the behaviour of many-body quantum systems out of equilibrium. In the presence of a globally conserved charge, further insights are provided by the knowledge of the resolution of entanglement in the various symmetry sectors. Here, we carry on the program we initiated in Parez et al (2021 Phys. Rev. B 103 L041104), for the study of the time evolution of the symmetry resolved entanglement in free fermion systems. We complete and extend our derivations also by defining and quantifying a symmetry resolved mutual information. The entanglement entropies display a time delay that depends on the charge sector that we characterise exactly. Both entanglement entropies and mutual information show effective equipartition in the scaling limit of large time and subsystem size. Furthermore, we argue that the behaviour of the charged entropies can be quantitatively understood in the framework of the quasiparticle picture for the spreading of entanglement, and hence we expect that a proper adaptation of our results should apply to a large class of integrable systems. We also find that the number entropy grows logarithmically with time before saturating to a value proportional to the logarithm of the subsystem size.

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.

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.

We generalise the form factor bootstrap approach to integrable field theories with U(1) symmetry to derive matrix elements of composite branch-point twist fields associated with symmetry resolved entanglement entropies. The bootstrap equations are solved for the free massive Dirac and complex boson theories, which are the simplest theories with U(1) symmetry. We present the exact and complete solution for the bootstrap, including vacuum expectation values and form factors involving any type and arbitrarily number of particles. The non-trivial solutions are carefully cross-checked by performing various limits and by the application of the ∆-theorem. An alternative and compact determination of the novel form factors is also presented. Based on the form factors of the U(1) composite branch-point twist fields, we re-derive earlier results showing entanglement equipartition for an interval in the ground state of the two models.

We study the symmetry resolved entanglement entropies in one-dimensional systems with boundaries. We provide some general results for conformal invariant theories and then move to a semi-infinite chain of free fermions. We consider both an interval starting from the boundary and away from it. We derive exact formulas for the charged and symmetry resolved entropies based on theorems and conjectures about the spectra of Toeplitz+Hankel matrices. En route to characterise the interval away from the boundary, we prove a general relation between the eigenvalues of Toeplitz+Hankel matrices and block Toeplitz ones. An important aspect is that the saddle-point approximation from charged to symmetry resolved entropies introduces algebraic corrections to the scaling that are much more severe than in systems without boundaries.

The time evolution of the entanglement entropy is a key concept to understand the structure of a nonequilibrium quantum state. In a large class of models, such evolution can be understood in terms of a semiclassical picture of moving quasiparticles spreading the entanglement throughout the system. However, it is not yet known how the entanglement splits between the sectors of an internal local symmetry of a quantum many-body system. Here, guided by the examples of conformal field theories and free-fermion chains, we show that the quasiparticle picture can be adapted to this goal, leading to a general conjecture for the charged entropies whose Fourier transform gives the desired symmetry-resolved entanglement Sn(q). We point out two physically relevant effects that should be easily observed in atomic experiments: a delay time for the onset of Sn(q) which grows linearly with |Δq| (the difference between the charge q and its mean value) and an effective equipartition when |Δq| is much smaller than the subsystem size.