Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1601.07492 · Physical review. B./Physical review. B

Entanglement negativity in free-fermion systems: An overlap matrix approach

Po-Yao Chang, Xueda Wen

Published May 20, 2016• 31 citations
arXivHTMLPDF

Abstract

In this paper, we calculate the entanglement negativity in free-fermion systems by use of the overlap matrices. For a tripartite system, if the ground state can be factored into triples of modes, we show that the partially transposed reduced density matrix can be factorized and the entanglement negativity has a simple form. However, the factorability of the ground state in a tripartite system does not hold in general. In this situation, the partially transposed reduced density matrix can be expressed in terms of the Kronecker product of matrices. We explicitly compute the entanglement negativity for the Su-Schrieffer-Heeger model, the integer Quantum Hall state, and a homogeneous one-dimensional chain. We find that the entanglement negativity for the integer quantum Hall states shows an area law behavior. For the entanglement negativity of two adjacent intervals in a homogeneous one-dimensional gas, we find agreement with the conformal field theory. Our method provides a numerically feasible way to study the entanglement negativity in free-fermion systems.

DOI
10.1103/physrevb.93.195140

Related papers

Fermionic partial transpose in the overlap matrix framework for entanglement negativity

Jun Qi Fang, Xiao Yan Xu

Over the past two decades, the overlap matrix approach has been developed to compute quantum entanglement in free-fermion systems, particularly to calculate entanglement entropy and entanglement negativity. This method involves the use of partial trace and partial transpose operations within the overlap matrix framework. However, in previous studies, only the conventional partial transpose in fermionic systems has been considered, which does not account for fermionic anticommutation relations. Although the concept of a fermionic partial transpose was introduced by Shapourian et al. [Phys. Rev. B 95, 165101 (2017)], it has not yet been systematically incorporated into the overlap matrix framework. In this paper, we introduce the fermionic partial transpose into the overlap matrix approach, provide a systematic analysis of the validity of partial trace and partial transpose operations, and derive an explicit formula for calculating entanglement negativity in bipartite systems. Additionally, we numerically compute the logarithmic negativity of two lattice models to verify the Gioev-Klich-Widom scaling law. For tripartite geometries, we uncover limitations of the overlap matrix method and demonstrate that the previously reported logarithmic negativity result for a homogeneous one-dimensional chain in a disjoint interval geometry exceeds its theoretical upper bound.

arXiv:2503.07742•Physical review. B./Physical review. B•Aug 13, 2025•2 citations

Entanglement negativity in extended systems: a field theoretical approach

Pasquale Calabrese, John Cardy, Erik Tonni C

We report on a systematic approach for the calculation of the negativity in the ground state of a one-dimensional quantum field theory. The partial transpose of the reduced density matrix of a subsystem A = A 1 ∪ A 2 is explicitly constructed as an imaginary-time path integral and from this the replicated traces are obtained. The logarithmic negativity is then the continuation to n → 1 of the traces of the even powers. For pure states, this procedure reproduces the known results. We then apply this method to conformally invariant field theories (CFTs) in several different physical situations for infinite and finite systems and without or with boundaries. In particular, in the case of two adjacent intervals of lengths ℓ 1 ,ℓ 2 in an infinite system, we derive the result ∼ ( c /4)ln(ℓ 1 ℓ 2 /(ℓ 1 + ℓ 2 )), where c is the central charge. For the more complicated case of two disjoint intervals, we show that the negativity depends only on the harmonic ratio of the four end points and so is manifestly scale invariant. We explicitly calculate the scale invariant functions for the replicated traces in the case of the CFT for the free compactified boson, but we have not so far been able to obtain the n → 1 continuation for the negativity even in the limit of large compactification radius. We have checked all our findings against exact numerical results for the harmonic chain which is described by a non-compactified free boson.

arXiv:1210.5359•Journal of Statistical Mechanics Theory and Experiment•Feb 1, 2013•311 citations

Aspects of Entanglement Negativity in One Dimensional Critical Systems

Andrea Coser

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.

•Oct 12, 2015•0 citations

Entanglement Rényi negativity of interacting fermions from quantum Monte Carlo simulations

Fo-Hong Wang, Xiao Yan Xu

Many-body entanglement unveils additional aspects of quantum matter and offers insights into strongly correlated physics. While ground-state entanglement has received much attention in the past decade, the study of mixed-state quantum entanglement using negativity in interacting fermionic systems remains largely unexplored. We demonstrate that the partially transposed density matrix of interacting fermions, similar to their reduced density matrix, can be expressed as a weighted sum of Gaussian states describing free fermions, enabling the calculation of rank-n Rényi negativity within the determinant quantum Monte Carlo framework. We calculate the rank-two Rényi negativity for the half-filled Hubbard model and the spinless t-V model. Our calculation reveals that the area law coefficient of the Rényi negativity for the spinless t-V model has a logarithmic finite-size scaling at the finite-temperature transition point. Our work contributes to the calculation of entanglement and sets the stage for future studies on quantum entanglement in various fermionic many-body mixed states.

arXiv:2312.14155•Nature Communications•Mar 17, 2025•14 citations

Entanglement and negativity Hamiltonians for the massless Dirac field on the half line

Federico Rottoli, Sara Murciano, Erik Tonni C +1 more

Abstract We study the ground-state entanglement Hamiltonian of several disjoint intervals for the massless Dirac fermion on the half-line. Its structure consists of a local part and a bi-local term that couples each point to another one in each other interval. The bi-local operator can be either diagonal or mixed in the fermionic chiralities and it is sensitive to the boundary conditions. The knowledge of such entanglement Hamiltonian is the starting point to evaluate the negativity Hamiltonian, i.e. the logarithm of the partially transposed reduced density matrix, which is an operatorial characterisation of entanglement of subsystems in mixed states. We find that the negativity Hamiltonian inherits the structure of the corresponding entanglement Hamiltonian. We finally show how the continuum expressions for both these operators can be recovered from exact numerical computations in free-fermion chains.

arXiv:2210.12109•Journal of Statistical Mechanics Theory and Experiment•Jan 1, 2023•19 citations

Estimation for entanglement negativity of free fermions

Christopher P. Herzog, Yihong Wang

Abstract In this letter we study the negativity of one dimensional free fermions. We derive the general form of the Z N symmetric term in moments of the partial transposed (reduced) density matrix, which is an algebraic function of the end points of the system. Such a path integral turns out to be a convenient tool for making estimations for the negativity.

arXiv:1601.00678•Journal of Statistical Mechanics Theory and Experiment•Jul 1, 2016•43 citations
In this paper, we calculate the entanglement negativity in free-fermion systems by use of the overlap matrices. For a tripartite system, if the ground state can be factored into triples of modes, we show that the partially transposed reduced density matrix can be factorized and the entanglement negativity has a simple form. However, the factorability of the ground state in a tripartite system does not hold in general. In this situation, the partially transposed reduced density matrix can be expressed in terms of the Kronecker product of matrices. We explicitly compute the entanglement negativity for the Su-Schrieffer-Heeger model, the integer Quantum Hall state, and a homogeneous one-dimensional chain. We find that the entanglement negativity for the integer quantum Hall states shows an area law behavior. For the entanglement negativity of two adjacent intervals in a homogeneous one-dimensional gas, we find agreement with the conformal field theory. Our method provides a numerically feasible way to study the entanglement negativity in free-fermion systems.
Over the past two decades, the overlap matrix approach has been developed to compute quantum entanglement in free-fermion systems, particularly to calculate entanglement entropy and entanglement negativity. This method involves the use of partial trace and partial transpose operations within the overlap matrix framework. However, in previous studies, only the conventional partial transpose in fermionic systems has been considered, which does not account for fermionic anticommutation relations. Although the concept of a fermionic partial transpose was introduced by Shapourian et al. [Phys. Rev. B 95, 165101 (2017)], it has not yet been systematically incorporated into the overlap matrix framework. In this paper, we introduce the fermionic partial transpose into the overlap matrix approach, provide a systematic analysis of the validity of partial trace and partial transpose operations, and derive an explicit formula for calculating entanglement negativity in bipartite systems. Additionally, we numerically compute the logarithmic negativity of two lattice models to verify the Gioev-Klich-Widom scaling law. For tripartite geometries, we uncover limitations of the overlap matrix method and demonstrate that the previously reported logarithmic negativity result for a homogeneous one-dimensional chain in a disjoint interval geometry exceeds its theoretical upper bound.
We report on a systematic approach for the calculation of the negativity in the ground state of a one-dimensional quantum field theory. The partial transpose of the reduced density matrix of a subsystem A = A 1 ∪ A 2 is explicitly constructed as an imaginary-time path integral and from this the replicated traces are obtained. The logarithmic negativity is then the continuation to n → 1 of the traces of the even powers. For pure states, this procedure reproduces the known results. We then apply this method to conformally invariant field theories (CFTs) in several different physical situations for infinite and finite systems and without or with boundaries. In particular, in the case of two adjacent intervals of lengths ℓ 1 ,ℓ 2 in an infinite system, we derive the result ∼ ( c /4)ln(ℓ 1 ℓ 2 /(ℓ 1 + ℓ 2 )), where c is the central charge. For the more complicated case of two disjoint intervals, we show that the negativity depends only on the harmonic ratio of the four end points and so is manifestly scale invariant. We explicitly calculate the scale invariant functions for the replicated traces in the case of the CFT for the free compactified boson, but we have not so far been able to obtain the n → 1 continuation for the negativity even in the limit of large compactification radius. We have checked all our findings against exact numerical results for the harmonic chain which is described by a non-compactified free boson.
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.
Many-body entanglement unveils additional aspects of quantum matter and offers insights into strongly correlated physics. While ground-state entanglement has received much attention in the past decade, the study of mixed-state quantum entanglement using negativity in interacting fermionic systems remains largely unexplored. We demonstrate that the partially transposed density matrix of interacting fermions, similar to their reduced density matrix, can be expressed as a weighted sum of Gaussian states describing free fermions, enabling the calculation of rank-n Rényi negativity within the determinant quantum Monte Carlo framework. We calculate the rank-two Rényi negativity for the half-filled Hubbard model and the spinless t-V model. Our calculation reveals that the area law coefficient of the Rényi negativity for the spinless t-V model has a logarithmic finite-size scaling at the finite-temperature transition point. Our work contributes to the calculation of entanglement and sets the stage for future studies on quantum entanglement in various fermionic many-body mixed states.
Abstract We study the ground-state entanglement Hamiltonian of several disjoint intervals for the massless Dirac fermion on the half-line. Its structure consists of a local part and a bi-local term that couples each point to another one in each other interval. The bi-local operator can be either diagonal or mixed in the fermionic chiralities and it is sensitive to the boundary conditions. The knowledge of such entanglement Hamiltonian is the starting point to evaluate the negativity Hamiltonian, i.e. the logarithm of the partially transposed reduced density matrix, which is an operatorial characterisation of entanglement of subsystems in mixed states. We find that the negativity Hamiltonian inherits the structure of the corresponding entanglement Hamiltonian. We finally show how the continuum expressions for both these operators can be recovered from exact numerical computations in free-fermion chains.