We study the dynamics of the entanglement in one-dimensional critical quantum systems after a local quench in which two independently thermalized semi-infinite halves are joined to form a homogeneous infinite system and left to evolve unitarily. We show that under certain conditions a nonequilibrium steady state (NESS) is reached instantaneously as soon as the entanglement interval is within the light cone emanating from the contact point. In this steady state, the exact expressions for the entanglement entropy and the logarithmic negativity are in agreement with the steady state density matrix being a boosted thermal state, as expected. We derive various general identities: relating the negativity after the quench with unequal left and right initial temperatures with that where the left and right temperatures are equal; and relating these with the negativity in equilibrium thermal states. In certain regimes the resulting expressions can be analytically evaluated. Immediately after the interval intersects the light cone, we find logarithmic growth. For a very long interval, we find that the negativity approaches a plateau after sufficiently long times, different from its NESS value. The NESS value is reached instantly as soon as the entire interval is contained in the light cone. This provides a theoretical framework explaining recently obtained numerical results.
We study the dynamics of the entanglement in one-dimensional critical quantum systems after a local quench in which two independently thermalized semi-infinite halves are joined to form a homogeneous infinite system and left to evolve unitarily. We show that under certain conditions a nonequilibrium steady state (NESS) is reached instantaneously as soon as the entanglement interval is within the light cone emanating from the contact point. In this steady state, the exact expressions for the entanglement entropy and the logarithmic negativity are in agreement with the steady state density matrix being a boosted thermal state, as expected. We derive various general identities: relating the negativity after the quench with unequal left and right initial temperatures with that where the left and right temperatures are equal; and relating these with the negativity in equilibrium thermal states. In certain regimes the resulting expressions can be analytically evaluated. Immediately after the interval intersects the light cone, we find logarithmic growth. For a very long interval, we find that the negativity approaches a plateau after sufficiently long times, different from its NESS value. The NESS value is reached instantly as soon as the entire interval is contained in the light cone. This provides a theoretical framework explaining recently obtained numerical results.
We study the time evolution of the entanglement negativity after a local quantum quench in (1 + 1)-dimensional conformal field theories (CFTs), which we introduce by suddenly joining two initially decoupled CFTs at their end points. We calculate the negativity evolution for both adjacent intervals and disjoint intervals explicitly. For two adjacent intervals, the entanglement negativity grows logarithmically in time right after the quench. After developing a plateau-like feature, the entanglement negativity drops to the ground-state value. For the case of two spatially separated intervals, a light-cone behavior is observed in the negativity evolution; in addition, a long-range entanglement, which is independent of the distance between two intervals, can be created. Our results agree with the heuristic picture that quasiparticles, which carry entanglement, are emitted from the joining point and propagate freely through the system. Our analytical results are confirmed by numerical calculations based on a critical harmonic chain.
We study the entanglement in a chain of harmonic oscillators driven out of equilibrium by preparing the two sides of the system at different temperatures, and subsequently joining them together. The steady state is constructed explicitly and the logarithmic negativity is calculated between two adjacent segments of the chain. We find that, for low temperatures, the steady-state entanglement is a sum of contributions pertaining to left- and right-moving excitations emitted from the two reservoirs. In turn, the steady-state entanglement is a simple average of the Gibbs-state values and thus its scaling can be obtained from conformal field theory. A similar averaging behaviour is observed during the entire time evolution. As a particular case, we also discuss a local quench where both sides of the chain are initialized in their respective ground states.
We study the entanglement in a chain of harmonic oscillators driven out of equilibrium by preparing the two sides of the system at different temperatures, and subsequently joining them together. The steady state is constructed explicitly and the logarithmic negativity is calculated between two adjacent segments of the chain. We find that, for low temperatures, the steady-state entanglement is a sum of contributions pertaining to left-and right-moving excitations emitted from the two reservoirs. In turn, the steady-state entanglement is a simple average of the Gibbs-state values and thus its scaling can be obtained from conformal field theory. A similar averaging behaviour is observed during the entire time evolution. As a particular case, we also discuss a local quench where both sides of the chain are initialized in their respective ground states.
We show that the dynamics resulting from preparing a one-dimensional quantum system in the ground state of two decoupled parts, then joined together and left to evolve unitarily with a translational invariant Hamiltonian (a local quench), can be described by means of quantum field theory. In the case when the corresponding theory is conformal, we study the evolution of the entanglement entropy for different bipartitions of the line. We also consider the behavior of one- and two-point correlation functions. All our findings may be explained in terms of a picture, that we believe to be valid more generally, whereby quasiparticles emitted from the joining point at the initial time propagate semiclassically through the system.
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.
We study entanglement properties in a nonequilibrium steady state of a free-fermion chain, that emerges after connecting two half-chains prepared at different temperatures. The entanglement negativity and the R\'enyi mutual information between two adjacent intervals scale logarithmically in the system size, with prefactors that we calculate analytically as a function of the bath temperatures. In particular, we show that the negativity and the R\'enyi mutual information with index $\ensuremath{\alpha}=1/2$ are described by different prefactors, and thus the two quantities provide inequivalent information about the state. Furthermore, we show that the logarithmic growth of the negativity during time evolution is also governed by the steady-state prefactor.