Mixed state entanglement measures can act as a versatile probes of many-body systems. However, they are generally hard to compute, often relying on tricky optimizations. One measure that is straightforward to compute is the logarithmic negativity, yet done naively even this is still limited to small system sizes. Here, we introduce a method to compute the logarithmic negativity for arbitrary subsystems of a densely represented state, as well as block subsystems of matrix product states. The method combines lazily evaluated, tensor network representations of the partially transposed density matrix with stochastic Lanczos quadrature, and is easily extendible to other quantities and classes of many-body states. As examples, we compute the entanglement within random pure states for density matrices of up to 30 qubits, explore scrambling in a many-body quench, and match the results of conformal field theory in the ground-state of the Heisenberg model for density matrices of up to 1000 spins. An implementation of the algorithm has been made available in the open-source library \textit{quimb}.
Mixed state entanglement measures can act as a versatile probes of many-body systems. However, they are generally hard to compute, often relying on tricky optimizations. One measure that is straightforward to compute is the logarithmic negativity, yet done naively even this is still limited to small system sizes. Here, we introduce a method to compute the logarithmic negativity for arbitrary subsystems of a densely represented state, as well as block subsystems of matrix product states. The method combines lazily evaluated, tensor network representations of the partially transposed density matrix with stochastic Lanczos quadrature, and is easily extendible to other quantities and classes of many-body states. As examples, we compute the entanglement within random pure states for density matrices of up to 30 qubits, explore scrambling in a many-body quench, and match the results of conformal field theory in the ground-state of the Heisenberg model for density matrices of up to 1000 spins. An implementation of the algorithm has been made available in the open-source library \textit{quimb}.
Entanglement is not only the key resource for many quantum technologies, but essential in understanding the structure of many-body quantum matter. At the interface of these two crucial areas are simulators, controlled systems capable of mimicking physical models that might escape analytical tractability. Traditionally, these simulations have been performed classically, where recent advancements such as tensor-networks have made explicit the limitation entanglement places on scalability. Increasingly however, analog quantum simulators are expected to yield deep insight into complex systems. This thesis advances the field in across various interconnected fronts. Firstly, we introduce schemes for verifying and distributing entanglement in a quantum dot simulator, tailored to specific experimental constraints. We then confirm that quantum dot simulators would be natural candidates for simulating many-body localization (MBL) - a recently emerged phenomenon that seems to evade traditional statistical mechanics. Following on from that, we investigate MBL from an entanglement perspective, shedding new light on the nature of the transition to it from a ergodic regime. As part of that investigation we make use of the logarithmic negativity, an entanglement measure applicable to many-body mixed states. In order to tie back into quantum simulators, we then propose an experimental scheme to measure the logarithmic negativity in realistic many-body settings. This method uses choice measurements on three or more copies of a mixed state along with machine learning techniques. We also introduce a fast method for computing many-body entanglement in classical simulations that significantly increases the size of system addressable. Finally, we introduce quimb, an open-source library for interactive but efficient quantum information and many-body calculations. It contains general purpose tensor-network support alongside other novel algorithms.
Can the entanglement of an arbitrary many-body state be quantified without reconstructing the state? We answer the question affirmatively by showing how the measure of logarithmic negativity can be estimated using a finite number of even partial transpose (PT) moments of a many-body state. Considering the squared PT spectrum as a probability measure, we solve exactly the task of a truncated Hausdorff moment problem with Gauss and Gauss-Radau quadrature. We obtain a hierarchy of lower bounds on the logarithmic negativity, which improves with increasing order of moments. Convergence to the exact value of the measure can be achieved for every state of finite dimension. We illustrate the bounds for ground states of three different interacting spin models, as well as for a thermal state. An error analysis demonstrates the efficacy of our scheme towards handling the limitation of finite statistics encountered in actual experiments with randomized measurements. Moreover, an example of density matrix renormalization group based numerical estimation is provided for handling system sizes beyond the scope of exact diagonalization.
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.
Entanglement plays a crucial role in quantum information science and many-body physics, yet quantifying it in mixed quantum many-body systems has remained a notoriously difficult problem. Here, we introduce families of quantitative entanglement witnesses, constructed from partial transpose and realignment moments, which provide rigorous bounds on entanglement monotones as well as entanglement dimensionality. Our witnesses can be efficiently measured using SWAP tests or variants of Bell measurements, thus making them directly implementable on current hardware. Leveraging our witnesses, we present several novel results on entanglement properties of mixed states, both in quantum information and many-body physics. We develop efficient algorithms to test whether mixed states with bounded entropy have low or high entanglement, which previously was only possible for pure states. We also provide an efficient algorithm to test the Schmidt rank using only two-copy measurements, and the operator Schmidt rank using four-copy measurements. Further, our witnesses robustly certify the quantum circuit depth in the presence of noise, as well as the Schmidt rank of mixed states. Finally, we show that the entanglement phase diagram of Haar random states, quantified by the partial transpose negativity, can be fully established solely by computing our witness, a result that also applies to any state 4 -design. Our witnesses can also be efficiently computed for matrix product states, thus enabling the characterization of entanglement in extensive many-body systems. Finally, we make progress on the entanglement required for quantum cryptography, establishing rigorous limits on pseudoentanglement and pseudorandom density matrices with bounded entropy. Our work opens new avenues for quantifying entanglement in large and noisy quantum systems.
Quantum entanglement plays a leading role in modern understanding of physical systems, from quantum phases of matter to quantum gravity. In quantum information theory, one seeks operationally meaningful quantifiers of entanglement, which for large quantum systems are notoriously difficult to evaluate due to the lack of computationally efficient algorithms. In this Letter, we show that for large random induced mixed states the logarithmic negativity, an efficiently computable entanglement measure, generically coincides with the exact entanglement cost under positive-partial-transpose-preserving operations, thereby acquiring a precise operational interpretation. Our results establish logarithmic negativity as an exact characterization of entanglement in generic many-body states and provide a tractable route for quantifying entanglement in complex quantum systems.
We study the entanglement between disjoint subregions in quantum critical systems through the lens of the logarithmic negativity. We work with systems in arbitrary dimensions, including conformal field theories and their corresponding lattice Hamiltonians, as well as resonating valence-bond states. At small separations, the logarithmic negativity is big and displays universal behavior, but we show nonperturbatively that it decays faster than any power at large separations. This can already be seen in the minimal setting of single-spin subregions. The corresponding absence of distillable entanglement at large separations generalizes the one-dimensional result, and indicates that quantum critical ground states do not possess long-range bipartite entanglement, at least for bosons. For systems with fermions, a more suitable definition of the logarithmic negativity exists that takes into account fermion parity, and we show that it decays algebraically. Along the way we obtain general results for the moments of the partially transposed density matrix. Published by the American Physical Society 2024