The entanglement between two arbitrary subsystems of random pure states is studied via properties of the density matrix's partial transpose ${\ensuremath{\rho}}_{12}^{{T}_{2}}$. The density of states of ${\ensuremath{\rho}}_{12}^{{T}_{2}}$ is close to the semicircle law when both subsystems have dimensions which are not too small and are of the same order. A simple random matrix model for the partial transpose is found to capture the entanglement properties well, including a transition across a critical dimension. Log negativity is used to quantify entanglement between subsystems and analytic formulas for this are derived based on the simple model. The skewness of the eigenvalue density of ${\ensuremath{\rho}}_{12}^{{T}_{2}}$ is derived analytically, using the average of the third moment over the ensemble of random pure states. The third moment after partial transpose is also shown to be related to a generalization of the Kempe invariant. The smallest eigenvalue after partial transpose is found to follow the extreme value statistics of random matrices; namely, the Tracy-Widom distribution. This distribution, with relevant parameters obtained from the model, is found to be useful in calculating the fraction of entangled states at critical dimensions. These results are tested in a quantum dynamical system of three coupled standard maps, where one finds that if the parameters represent a strongly chaotic system, the results are close to those of random states, although there are some systematic deviations at critical dimensions.
The entanglement between two arbitrary subsystems of random pure states is studied via properties of the density matrix's partial transpose \ensuremathρ12T2. The density of states of \ensuremathρ12T2 is close to the semicircle law when both subsystems have dimensions which are not too small and are of the same order. A simple random matrix model for the partial transpose is found to capture the entanglement properties well, including a transition across a critical dimension. Log negativity is used to quantify entanglement between subsystems and analytic formulas for this are derived based on the simple model. The skewness of the eigenvalue density of \ensuremathρ12T2 is derived analytically, using the average of the third moment over the ensemble of random pure states. The third moment after partial transpose is also shown to be related to a generalization of the Kempe invariant. The smallest eigenvalue after partial transpose is found to follow the extreme value statistics of random matrices; namely, the Tracy-Widom distribution. This distribution, with relevant parameters obtained from the model, is found to be useful in calculating the fraction of entangled states at critical dimensions. These results are tested in a quantum dynamical system of three coupled standard maps, where one finds that if the parameters represent a strongly chaotic system, the results are close to those of random states, although there are some systematic deviations at critical dimensions.
The probability of large deviations of the smallest Schmidt eigenvalue for random pure states of bipartite systems, denoted as A and B, is computed analytically using a Coulomb gas method. It is shown that this probability, for large N, goes as exp[-βN^{2}Φ(ζ)], where the parameter β is the Dyson index of the ensemble, ζ is the large deviation parameter, while the rate function Φ(ζ) is calculated exactly. Corresponding equilibrium Coulomb charge density is derived for its large deviations. Effects of the large deviations of the extreme (largest and smallest) Schmidt eigenvalues on the bipartite entanglement are studied using the von Neumann entropy. Effect of these deviations is also studied on the entanglement between subsystems 1 and 2, obtained by further partitioning the subsystem A, using the properties of the density matrix's partial transpose ρ_{12}^{Γ}. The density of states of ρ_{12}^{Γ} is found to be close to the Wigner's semicircle law with these large deviations. The entanglement properties are captured very well by a simple random matrix model for the partial transpose. The model predicts the entanglement transition across a critical large deviation parameter ζ. Log negativity is used to quantify the entanglement between subsystems 1 and 2. Analytical formulas for it are derived using the simple model. Numerical simulations are in excellent agreement with the analytical results.
The entanglement properties of random pure states are relevant to a variety of problems ranging from chaotic quantum dynamics to black-hole physics.The averaged bipartite entanglement entropy of such states admits a volume law and upon increasing the subregion size follows the Page curve.In this paper, we generalize this setup to random mixed states by coupling the system to a bath and use the partial transpose to study their entanglement properties.We develop a diagrammatic method to incorporate partial transpose within random matrix theory and formulate a perturbation theory in 1/L, the inverse of the Hilbert-space dimension.We compute several quantities including the spectral density of partial transpose (or entanglement negativity spectrum), two-point correlator of eigenvalues, and the logarithmic negativity.As long as the bath is smaller than the system, we find that upon sweeping the subregion size, the logarithmic negativity shows an initial increase and a final decrease similar to the Page curve, while it admits a plateau in the intermediate regime where the logarithmic negativity depends only on the size of the system and of the bath but not on how the system is partitioned.This intermediate phase has no analog in random pure states, and is separated from the two other regimes by a critical point.We further show that when the bath is larger than the system by at least two extra qubits the logarithmic negativity is identically zero, which implies that there is no distillable entanglement.Using the diagrammatic approach, we provide a simple derivation of the semicircle law of the entanglement negativity spectrum in the latter two regimes.We show that despite the appearance of a semicircle distribution, reminiscent of Gaussian unitary ensemble (GUE), the higher-order corrections to the negativity spectrum and two-point correlator deviate from those of GUE.
The entanglement properties of random pure states are relevant to a variety of problems ranging from chaotic quantum dynamics to black-hole physics. The averaged bipartite entanglement entropy of such states admits a volume law and upon increasing the subregion size follows the Page curve. In this paper, we generalize this setup to random mixed states by coupling the system to a bath and use the partial transpose to study their entanglement properties. We develop a diagrammatic method to incorporate partial transpose within random matrix theory and formulate a perturbation theory in 1/L, the inverse of the Hilbert-space dimension. We compute several quantities including the spectral density of partial transpose (or entanglement negativity spectrum), two-point correlator of eigenvalues, and the logarithmic negativity. As long as the bath is smaller than the system, we find that upon sweeping the subregion size, the logarithmic negativity shows an initial increase and a final decrease similar to the Page curve, while it admits a plateau in the intermediate regime where the logarithmic negativity depends only on the size of the system and of the bath but not on how the system is partitioned. This intermediate phase has no analog in random pure states, and is separated from the two other regimes by a critical point. We further show that when the bath is larger than the system by at least two extra qubits the logarithmic negativity is identically zero, which implies that there is no distillable entanglement. Using the diagrammatic approach, we provide a simple derivation of the semicircle law of the entanglement negativity spectrum in the latter two regimes. We show that despite the appearance of a semicircle distribution, reminiscent of Gaussian unitary ensemble (GUE), the higher-order corrections to the negativity spectrum and two-point correlator deviate from those of GUE.
We analyze the properties of entangled random pure states of a quantum system partitioned into two smaller subsystems of dimensions N and M. Framing the problem in terms of random matrices with a fixed-trace constraint, we establish, for arbitrary N≤M, a general relation between the n-point densities and the cross moments of the eigenvalues of the reduced density matrix, i.e., the so-called Schmidt eigenvalues, and the analogous functionals of the eigenvalues of the Wishart-Laguerre ensemble of the random matrix theory. This allows us to derive explicit expressions for two-level densities, and also an exact expression for the variance of von Neumann entropy at finite N,M. Then, we focus on the moments E{K^{a}} of the Schmidt number K, the reciprocal of the purity. This is a random variable supported on [1,N], which quantifies the number of degrees of freedom effectively contributing to the entanglement. We derive a wealth of analytical results for E{K^{a}} for N=2 and 3 and arbitrary M, and also for square N=M systems by spotting for the latter a connection with the probability P(x_{min}^{GUE}≥sqrt[2N]ξ) that the smallest eigenvalue x_{min}^{GUE} of an N×N matrix belonging to the Gaussian unitary ensemble is larger than sqrt[2N]ξ. As a by-product, we present an exact asymptotic expansion for P(x_{min}^{GUE}≥sqrt[2N]ξ) for finite N as ξ→∞. Our results are corroborated by numerical simulations whenever possible, with excellent agreement.
For a system of N identical particles in a random pure state, there is a threshold ${k}_{0}={k}_{0}(N)\ensuremath{\sim}N/5$ such that two subsystems of k particles each typically share entanglement if k>k0, and typically do not share entanglement if k<k0. By ``random'' we mean here uniformly distributed on the sphere of the corresponding Hilbert space. The analogous phase transition for the positive partial transpose (PPT) property can be described even more precisely. For example, for N qubits the two subsystems of size k are typically in a PPT state if $k<{k}_{1}:=N/4\ensuremath{-}1/2$ and typically in a non-PPT state if k>k1. Since, for a given state of the entire system, the induced state of a subsystem is given by the partial trace, the above facts can be rephrased as properties of random induced states. An important step in the analysis depends on identifying the asymptotic spectral density of the partial transposes of such random induced states, a result which is interesting in its own right.
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.