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:2602.18099 · SciPost Physics

A contour for the entanglement negativity of bosonic Gaussian states

Gioele Zambotti, Erik Tonni

cond-mat.stat-mechhep-thquant-phSubmitted Feb 20, 2026• Announced Feb 23, 2026• 0 citations
arXivHTMLPDF

Abstract

We construct a contour function for the logarithmic negativity and the logarithm of the moments of the partial transpose of the reduced density matrix for multimode bosonic Gaussian states of a free lattice model. In one spatial dimension, numerical results are obtained for harmonic chains either in the ground state or at finite temperature, by considering, respectively, either a subsystem made by two adjacent or disjoint blocks on the line or a bipartition of the circle. The contour function of the logarithmic negativity diverges only at the entangling points, while the contour function for the logarithm of the moments of the partial transpose is divergent also at the boundary of the bipartite subsystem, as functions of the position. In a two-dimensional conformal field theory, analytic expressions that describe these divergencies are discussed. In one spatial dimension, we explore the partial derivative of the logarithmic negativity of two adjacent intervals with respect to the logarithm of the harmonic ratio of their lengths while their ratio and the other parameters are kept fixed. Considering the ground state of the harmonic chain on the line and in the massive regime, we report numerical results showing that this quantity displays a monotonically decreasing behaviour.

Journal reference
SciPost Phys. 21, 027 (2026)
DOI
10.21468/SciPostPhys.21.1.027
Comments
62 pages, 14 figures

Related papers

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

Entanglement Negativity in Quantum Field Theory

Pasquale Calabrese, John Cardy, Erik Tonni C

We develop a systematic method to extract the negativity in the ground state of a $1+1$ dimensional relativistic quantum field theory, using a path integral formalism to construct the partial transpose ${\ensuremath{\rho}}_{A}^{{T}_{2}}$ of the reduced density matrix of a subsystem $A={A}_{1}\ensuremath{\cup}{A}_{2}$, and introducing a replica approach to obtain its trace norm which gives the logarithmic negativity $\mathcal{E}=\mathrm{ln}\ensuremath{\Vert}{\ensuremath{\rho}}_{A}^{{T}_{2}}\ensuremath{\Vert}$. This is shown to reproduce standard results for a pure state. We then apply this method to conformal field theories, deriving the result $\mathcal{E}\ensuremath{\sim}(c/4)\mathrm{ln}[{\ensuremath{\ell}}_{1}{\ensuremath{\ell}}_{2}/({\ensuremath{\ell}}_{1}+{\ensuremath{\ell}}_{2})]$ for the case of two adjacent intervals of lengths ${\ensuremath{\ell}}_{1}$, ${\ensuremath{\ell}}_{2}$ in an infinite system, where $c$ is the central charge. For two disjoint intervals it depends only on the harmonic ratio of the four end points and so is manifestly scale invariant. We check our findings against exact numerical results in the harmonic chain.

arXiv:1206.3092•Physical Review Letters•Sep 28, 2012•463 citations

Entanglement negativity in a two dimensional harmonic lattice: area law and corner contributions

Cristiano De Nobili, Andrea Coser, Erik Tonni C

Abstract We study the logarithmic negativity and the moments of the partial transpose in the ground state of a two dimensional massless harmonic square lattice with nearest neighbour interactions for various configurations of adjacent domains. At leading order for large domains, the logarithmic negativity and the logarithm of the ratio between the generic moment of the partial transpose and the moment of the reduced density matrix at the same order satisfy an area law in terms of the length of the curve shared by the adjacent regions. We give numerical evidence that the coefficient of the area law term in these quantities is related to the coefficient of the area law term in the Rényi entropies. Whenever the curve shared by the adjacent domains contains vertices, a subleading logarithmic term occurs in these quantities and the numerical values of the corner function for some pairs of angles are obtained. In the special case of vertices corresponding to explementary angles, we provide numerical evidence that the corner function of the logarithmic negativity is given by the corner function of the Rényi entropy of order 1/2.

arXiv:1604.02609•Journal of Statistical Mechanics Theory and Experiment•Aug 1, 2016•67 citations

A contour for the entanglement entropies in harmonic lattices

Andrea Coser, Cristiano De Nobili, Erik Tonni C

Abstract We construct a contour function for the entanglement entropies in generic harmonic lattices. In one spatial dimension, numerical analysis are performed by considering harmonic chains with either periodic or Dirichlet boundary conditions. In the massless regime and for some configurations where the subsystem is a single interval, the numerical results for the contour function are compared to the inverse of the local weight function which multiplies the energy-momentum tensor in the corresponding entanglement hamiltonian, found through conformal field theory methods, and a good agreement is observed. A numerical analysis of the contour function for the entanglement entropy is performed also in a massless harmonic chain for a subsystem made by two disjoint intervals.

arXiv:1701.08427•Journal of Physics A Mathematical and Theoretical•Jun 13, 2017•68 citations

Negativity spectrum of one-dimensional conformal field theories

Paola Ruggiero, Vincenzo Alba, Pasquale Calabrese

Logarithmic negativity is a proper measure of quantum entanglement in mixed states. Originally proposed in the context of quantum information, in recent years it is attracting interest in quantum field theory, condensed matter physics, and statistical physics. The key object to construct the negativity is the partial transpose of the reduced density matrix. Here, focusing on two adjacent intervals, the authors provide a complete study of the distribution of the eigenvalues of the partial transpose (negativity spectrum) for gapless one-dimensional systems described by a conformal field theory (CFT). Similar to the entanglement spectrum, the distribution of the negativity spectrum is universal, and it depends only on the central charge of the underlying CFT. The precise form of the spectrum depends on whether the two intervals are in a pure or mixed state, and in both cases, a dependence on the sign of the eigenvalues is found. This dependence is weak for bulk eigenvalues, whereas it is strong at the spectrum edges. The authors also investigate the scaling of the smallest (negative) and largest (positive) eigenvalues of the partial transpose. The analytical results are thoroughly checked against numerical DMRG simulations as well as in exactly solvable models.

arXiv:1607.02992•Physical review. B./Physical review. B•Nov 10, 2016•74 citations

Entanglement negativity and conformal field theory: a Monte Carlo study

Vincenzo Alba

We investigate the behavior of the moments of the partially transposed reduced density matrix ρ^{T_2}_A in critical quantum spin chains. Given subsystem A as union of two blocks, this is the (matrix) transposed of ρ_A with respect to the degrees of freedom of one of the two. This is also the main ingredient for constructing the logarithmic negativity. We provide a new numerical scheme for calculating efficiently all the moments of ρ_A^{T_2} using classical Monte Carlo simulations. In particular we study several combinations of the moments which are scale invariant at a critical point. Their behavior is fully characterized in both the critical Ising and the anisotropic Heisenberg XXZ chains. For two adjacent blocks we find, in both models, full agreement with recent CFT calculations. For disjoint ones, in the Ising chain finite size corrections are non negligible. We demonstrate that their exponent is the same governing the unusual scaling corrections of the mutual information between the two blocks. Monte Carlo data fully match the theoretical CFT prediction only in the asymptotic limit of infinite intervals. Oppositely, in the Heisenberg chain scaling corrections are smaller and, already at finite (moderately large) block sizes, Monte Carlo data are in excellent agreement with the asymptotic CFT result.

arXiv:1302.1110•Journal of Statistical Mechanics Theory and Experiment•May 17, 2013•86 citations
We construct a contour function for the logarithmic negativity and the logarithm of the moments of the partial transpose of the reduced density matrix for multimode bosonic Gaussian states of a free lattice model. In one spatial dimension, numerical results are obtained for harmonic chains either in the ground state or at finite temperature, by considering, respectively, either a subsystem made by two adjacent or disjoint blocks on the line or a bipartition of the circle. The contour function of the logarithmic negativity diverges only at the entangling points, while the contour function for the logarithm of the moments of the partial transpose is divergent also at the boundary of the bipartite subsystem, as functions of the position. In a two-dimensional conformal field theory, analytic expressions that describe these divergencies are discussed. In one spatial dimension, we explore the partial derivative of the logarithmic negativity of two adjacent intervals with respect to the logarithm of the harmonic ratio of their lengths while their ratio and the other parameters are kept fixed. Considering the ground state of the harmonic chain on the line and in the massive regime, we report numerical results showing that this quantity displays a monotonically decreasing behaviour.
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.
We develop a systematic method to extract the negativity in the ground state of a 1+11+11+1 dimensional relativistic quantum field theory, using a path integral formalism to construct the partial transpose ${\ensuremath{\rho}}_{A}^{{T}_{2}}$ of the reduced density matrix of a subsystem $A={A}_{1}\ensuremath{\cup}{A}_{2}$, and introducing a replica approach to obtain its trace norm which gives the logarithmic negativity $\mathcal{E}=\mathrm{ln}\ensuremath{\Vert}{\ensuremath{\rho}}_{A}^{{T}_{2}}\ensuremath{\Vert}$. This is shown to reproduce standard results for a pure state. We then apply this method to conformal field theories, deriving the result $\mathcal{E}\ensuremath{\sim}(c/4)\mathrm{ln}[{\ensuremath{\ell}}_{1}{\ensuremath{\ell}}_{2}/({\ensuremath{\ell}}_{1}+{\ensuremath{\ell}}_{2})]$ for the case of two adjacent intervals of lengths ${\ensuremath{\ell}}_{1}$, ${\ensuremath{\ell}}_{2}$ in an infinite system, where ccc is the central charge. For two disjoint intervals it depends only on the harmonic ratio of the four end points and so is manifestly scale invariant. We check our findings against exact numerical results in the harmonic chain.
Abstract We study the logarithmic negativity and the moments of the partial transpose in the ground state of a two dimensional massless harmonic square lattice with nearest neighbour interactions for various configurations of adjacent domains. At leading order for large domains, the logarithmic negativity and the logarithm of the ratio between the generic moment of the partial transpose and the moment of the reduced density matrix at the same order satisfy an area law in terms of the length of the curve shared by the adjacent regions. We give numerical evidence that the coefficient of the area law term in these quantities is related to the coefficient of the area law term in the Rényi entropies. Whenever the curve shared by the adjacent domains contains vertices, a subleading logarithmic term occurs in these quantities and the numerical values of the corner function for some pairs of angles are obtained. In the special case of vertices corresponding to explementary angles, we provide numerical evidence that the corner function of the logarithmic negativity is given by the corner function of the Rényi entropy of order 1/2.
Abstract We construct a contour function for the entanglement entropies in generic harmonic lattices. In one spatial dimension, numerical analysis are performed by considering harmonic chains with either periodic or Dirichlet boundary conditions. In the massless regime and for some configurations where the subsystem is a single interval, the numerical results for the contour function are compared to the inverse of the local weight function which multiplies the energy-momentum tensor in the corresponding entanglement hamiltonian, found through conformal field theory methods, and a good agreement is observed. A numerical analysis of the contour function for the entanglement entropy is performed also in a massless harmonic chain for a subsystem made by two disjoint intervals.
Logarithmic negativity is a proper measure of quantum entanglement in mixed states. Originally proposed in the context of quantum information, in recent years it is attracting interest in quantum field theory, condensed matter physics, and statistical physics. The key object to construct the negativity is the partial transpose of the reduced density matrix. Here, focusing on two adjacent intervals, the authors provide a complete study of the distribution of the eigenvalues of the partial transpose (negativity spectrum) for gapless one-dimensional systems described by a conformal field theory (CFT). Similar to the entanglement spectrum, the distribution of the negativity spectrum is universal, and it depends only on the central charge of the underlying CFT. The precise form of the spectrum depends on whether the two intervals are in a pure or mixed state, and in both cases, a dependence on the sign of the eigenvalues is found. This dependence is weak for bulk eigenvalues, whereas it is strong at the spectrum edges. The authors also investigate the scaling of the smallest (negative) and largest (positive) eigenvalues of the partial transpose. The analytical results are thoroughly checked against numerical DMRG simulations as well as in exactly solvable models.
We investigate the behavior of the moments of the partially transposed reduced density matrix ρ^{T_2}_A in critical quantum spin chains. Given subsystem A as union of two blocks, this is the (matrix) transposed of ρ_A with respect to the degrees of freedom of one of the two. This is also the main ingredient for constructing the logarithmic negativity. We provide a new numerical scheme for calculating efficiently all the moments of ρ_A^{T_2} using classical Monte Carlo simulations. In particular we study several combinations of the moments which are scale invariant at a critical point. Their behavior is fully characterized in both the critical Ising and the anisotropic Heisenberg XXZ chains. For two adjacent blocks we find, in both models, full agreement with recent CFT calculations. For disjoint ones, in the Ising chain finite size corrections are non negligible. We demonstrate that their exponent is the same governing the unusual scaling corrections of the mutual information between the two blocks. Monte Carlo data fully match the theoretical CFT prediction only in the asymptotic limit of infinite intervals. Oppositely, in the Heisenberg chain scaling corrections are smaller and, already at finite (moderately large) block sizes, Monte Carlo data are in excellent agreement with the asymptotic CFT result.