We investigate mixed-state entanglement measures, namely reflected entropy, mutual information and logarithmic negativity, for two disjoint intervals in one-dimensional systems excited by a finite number of quasiparticles. While whole system is in a pure state, the two disjoint intervals are in a generically mixed state. To address the problem that natural subsystem bases are generically non-orthonormal in such excited states, we use a general and efficient algorithm that computes these measures directly from the density matrix expressed in an arbitrary non-orthonormal basis. Applying this method to classical, bosonic, and fermionic quasiparticle excitations on a circle, we discover a universal additivity property: in the limit of large momentum differences, the mixed-state entanglement of a multi-quasiparticle state decomposes exactly into the sum of independent contributions. This additivity unifies the entanglement behavior across classical and quantum statistics, with the classical result emerging naturally as a special case. Our findings establish a robust computational framework for mixed-state entanglement in excited many-body systems and reveal a generic decoupling mechanism that governs entanglement distribution beyond the ground state.
We investigate mixed-state entanglement measures, namely reflected entropy, mutual information and logarithmic negativity, for two disjoint intervals in one-dimensional systems excited by a finite number of quasiparticles. While whole system is in a pure state, the two disjoint intervals are in a generically mixed state. To address the problem that natural subsystem bases are generically non-orthonormal in such excited states, we use a general and efficient algorithm that computes these measures directly from the density matrix expressed in an arbitrary non-orthonormal basis. Applying this method to classical, bosonic, and fermionic quasiparticle excitations on a circle, we discover a universal additivity property: in the limit of large momentum differences, the mixed-state entanglement of a multi-quasiparticle state decomposes exactly into the sum of independent contributions. This additivity unifies the entanglement behavior across classical and quantum statistics, with the classical result emerging naturally as a special case. Our findings establish a robust computational framework for mixed-state entanglement in excited many-body systems and reveal a generic decoupling mechanism that governs entanglement distribution beyond the ground state.
The Rényi entanglement entropies of quasiparticle excitations in the many-body gapped systems show a remarkable universal picture which is independent of the model, the quasiparticle momenta, and the connectedness of the subsystem. In this letter we calculate exactly the single-interval and double-interval Rényi entanglement entropies of quasiparticle excitations in the many-body gapped fermions, bosons, and XY chains and find additional contributions to the universal Rényi entanglement entropy in the excited states with quasiparticles of different momenta. The additional terms are different in different models, depend on the momentum differences of the quasiparticles, and are different for the single interval and the double interval. We derive the analytical Rényi entanglement entropy in the extremely gapped limit, matching perfectly the numerical results as long as either the intrinsic correlation length of the model or all the de Broglie wavelengths of the quasiparticles are small. When the momentum difference of any pair of distinct quasiparticles is small, the additional terms are non-negligible. We argue that the derived formulas can be applied to a wider range of models than those discussed here.
Abstract The statistical mechanics characterization of finite subsystems embedded in an infinite system is a fundamental question of quantum physics. Nevertheless, a full closed form for all required entropic measures does not exist in the general case, even for free systems, when the finite system in question is composed of several disjoint intervals. Here we develop a mathematical framework based on the Riemann–Hilbert approach to treat this problem in the one-dimensional case where the finite system is composed of two disjoint intervals and in the thermodynamic limit (both intervals and the space between them contains an infinite number of lattice sites and the result is given as a thermodynamic expansion). To demonstrate the usefulness of our method, we compute the change in the entanglement and negativity spectra namely the spectrum of eigenvalues of the reduced density matrix with or without time reversal of one of the intervals. We do this in the case that the distance between the intervals is much larger than their size. The method we use can be easily applied to compute any power in an expansion in the ratio of the distance between the intervals to their size. We expect these results to provide the necessary mathematical apparatus to address relevant questions in concrete physical scenarios, namely the structure and extent of quantum correlations in fermionic systems subject to local environments.
In a non-equilibrium many-body system, the quantum information dynamics between non-complementary regions is a crucial feature to understand the local relaxation towards statistical ensembles. Unfortunately, its characterization is a formidable task, as non-complementary parts are generally in a mixed state. We show that for integrable systems, this quantum information dynamics can be quantitatively understood within the quasiparticle picture for the entanglement spreading. Precisely, we provide an exact prediction for the time evolution of the logarithmic negativity after a quench. In the space-time scaling limit of long times and large subsystems, the negativity becomes proportional to the Rényi mutual information with Rényi index α = 1 / 2 . We provide robust numerical evidence for the validity of our results for free-fermion and free-boson models, but our framework applies to any interacting integrable system.
We investigate the quantum entanglement content of quasiparticle excitations in extended many-body systems. We show that such excitations give an additive contribution to the bipartite von Neumann and Rényi entanglement entropies that takes a simple, universal form. It is largely independent of the momenta and masses of the excitations and of the geometry, dimension, and connectedness of the entanglement region. The result has a natural quantum information theoretic interpretation as the entanglement of a state where each quasiparticle is associated with two qubits representing their presence within and without the entanglement region, taking into account quantum (in)distinguishability. This applies to any excited state composed of finite numbers of quasiparticles with finite de Broglie wavelengths or finite intrinsic correlation length. This includes particle excitations in massive quantum field theory and gapped lattice systems, and certain highly excited states in conformal field theory and gapless models. We derive this result analytically in one-dimensional massive bosonic and fermionic free field theories and for simple setups in higher dimensions. We provide numerical evidence for the harmonic chain and the two-dimensional harmonic lattice in all regimes where the conditions above apply. Finally, we provide supporting calculations for integrable spin chain models and other interacting cases without particle production. Our results point to new possibilities for creating entangled states using many-body quantum systems.
We investigate the quasiparticle excited state R\'enyi entropy of single, double and more intervals in the fermionic chain, bosonic chain and XY chain. When the gap of the theory is large or all the momenta of the quasiparticles are large, the R\'enyi Entropy takes a super universal form, which is independent of, among others, the model, the quasiparticle momenta, and the subsystem connectedness. We calculate analytically the R\'enyi Entropy in the extremely gapped limit by writing the excited states in terms of the local excitations and find different additional contributions to the super universal R\'enyi entropy in various models. The additional contributions to the super universal R\'enyi entropy cannot be neglected when the momentum differences of the excited quasiparticles are small, however, they are negligible when all the momentum differences are large. The R\'enyi entropy with additional terms derived in the extremely gapped limit is universal in the sense that it is valid in the slightly gapped and critical models as long as all the momenta of the excited quasiparticles are large. In the case of double interval in the XY chain we find new exact super universal as well as the universal results. In the case of the bosonic chain in the extremely massive limit we find analytically a novel formula, written as a permanent of a certain matrix, for the R\'enyi entropy. We support all of our analytical results with numerical calculations.
In this thesis we study entanglement in one-dimensional critical quantum many-body systems and in particular we will focus on disconnected regions. Given a system in a pure state, to quantify the amount of entanglement between a multicomponent subsystem and the rest of the full system we can use as an entanglement measure the renowned entanglement entropy. However, if we are interested in the entanglement shared among the disconnected regions, the entanglement entropy fails to be a good quantifier. The reason is that the state of the subsystem is in general mixed once the rest is traced out, and the entanglement entropy is a good measure of entanglement only for pure states. A good measure of entanglement in mixed states is the logarithmic negativity, which is the quantitative version of Peres' criterion of separability. The main advantage of the negativity with respect to other entanglement measures is the simplicity of its definition in terms of the density matrix describing the quantum state. Since the definition does not require any variational calculus, it is much more easily computable than any other entanglement measure, and therefore we can obtain some results also in complicated settings such as many-body systems. We will be mostly interested in the configuration where the full system is divided into two parts, A and B. If we do not have access to the degrees of freedom in B, we can describe subsystem A through its reduced density matrix, where the degrees of freedom in B have been traced out. The residual subsystem in A will in general be left in a mixed state. Subsystem A is then divided again into two parts, A_1 and A_2, and we will be interested in the entanglement shared by these two components. Knowing the full density matrix, the logarithmic negativity can be easily computed from the eigenvalues of its partial transpose with respect to the degrees of freedom living in one of the two subsystems A_1 or A_2. However, the full density matrix of a many-body state is in general unaccessible, even numerically, since the size of the matrix grows exponentially with the size of the system. If we concentrate on critical systems whose low-energy physics can be described by a quantum field theory, we can resort to its powerful tools to compute the entanglement properties. In particular the main tool is the replica trick. The entanglement entropy is obtained from the moments of the reduced density matrix, while the negativity can be computed from the moments of the partial transpose. However, even the computation of the integer-order moments is not an easy task, and exact analytical results are known only for the simplest quantum field theories. Hence, throughout the thesis we will usually consider conformal field theories, which have an enhanced set of symmetries and therefore allow for some exact computations. In the Introduction, we will try to frame this work by briefly reviewing some of the main topics in the literature of quantum many-body systems and quantum field theories where entanglement plays a crucial role. We will also describe what are the main features that are requested to a `good' measure of entanglement, and we will review some of the main measures that have been considered so far in the context of quantum information. We will define the entanglement negativity and stress its main advantages, as well as its drawbacks. After reviewing some basic facts of conformal field theories, we will give some technical details on the computation of entanglement entropy and logarithmic negativity in general quantum field theories, which will be needed in the rest of the thesis. In Chap. 2 we will start the study of the entanglement of several disjoint disconnected regions by considering the entanglement that the union of these regions share with the rest of the system. We will compute the integer-order Rényi entropies of the subsystem from which the entanglement entropy can be obtained through the replica limit. Unfortunately, the exact analytic continuation to real order is still out of reach and therefore the entanglement entropy cannot be computed exactly. However numerical extrapolations allow to obtain some accurate estimates. In the remaining chapters we will focus on the entanglement shared between two non complementary disjoint regions, through the computations of the entanglement negativity in different settings. In Chap. 3, we will consider a global quantum quench starting from a conformal boundary state and evolving through a conformal evolution. A general formula for the mutual information and the logarithmic negativity for two adjacent and disjoint intervals is given in the spacetime limit. In Chap. 4 we focus on the XY chain and we recover a formula for the partial transpose as a sum of four fermionic auxiliary Gaussian density matrices, by generalizing some previous results for the reduced density matrix of spin systems and for the partial transpose of pure fermionc Gaussian states. Even if the computation of the negativity is still out of reach, we can obtain formulas for the integer-order moments of the partial transpose in terms of the correlation matrices relative to the two components and to the region in between them. In Chap. 5 we will study the entanglement negativity for a free Dirac fermion field. Starting from the lattice results, we will obtain a path integral representation of the partial transpose which can be easily generalized to the continuum. With this representation of the partial transpose, we can construct all its integer-order moments and compute them for the simple conformal field theory of the free fermion. Again, the computation of the negativity could not be accomplished due to technical difficulties in the analytic continuation to real order. The same computation has been extended also to the modular invariant Dirac fermion and the Ising model, and the obtained formulas coincide with the ones already present in the literature. This analysis draws an interesting connection between some terms appearing in the formulas for the moments of the partial transpose (as well as of the R\'enyi entropies) found for the lattice models and the ones found for the corresponding theories describing their scaling limit. Whenever possible, all the analytical results will be checked against numerical calculations performed on simple free chain models. In Chaps. 2 and 3 we will consider bosonic Gaussian states, specifically the harmonic chain in its ground state and out of equilibrium, while in Chaps. 4 and 5 we will consider fermionic systems, specifically the XX and Ising spin chains, and the tight-binding model. Finally, we will draw some conclusions and discuss some open problems.