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:1312.1168 · Physical Review B

Entanglement negativity via the replica trick: A quantum Monte Carlo approach

Chia-Min Chung, Vincenzo Alba, Lars Bonnes, Pochung Chen, Andreas Martin Läuchli

Published Aug 4, 2014• 76 citations
arXivHTMLPDF

Abstract

Motivated by recent developments in conformal field theory (CFT), we devise a quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. These quantities can serve as witnesses of criticality. In particular, we study several scale invariant combinations of the moments for the one-dimensional (1D) hard-core boson model. For two adjacent intervals unusual finite-size corrections are present, showing parity effects that oscillate with a filling dependent period. These are more pronounced in the presence of boundaries. For large chains we find perfect agreement with CFT. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. Remarkably, the finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderately large sizes.

DOI
10.1103/physrevb.90.064401

Related papers

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

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 the critical Ising chain

Pasquale Calabrese, Luca Tagliacozzo, Erik Tonni C

We study the scaling of the traces of the integer powers of the partially transposed reduced density matrix and of the entanglement negativity for two spin blocks as function of their length and separation in the critical Ising chain. For two adjacent blocks, we show that tensor network calculations agree with universal conformal field theory (CFT) predictions. In the case of two disjoint blocks the CFT predictions are recovered only after taking into account the finite size corrections induced by the finite length of the blocks.

arXiv:1302.1113•Journal of Statistical Mechanics Theory and Experiment•May 7, 2013•115 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

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

Quantum Monte Carlo Study of Two Component Bosons and Entanglement

鍾佳民

This thesis is divided into two parts. In the first part we study quantum phases in bosonic systems in an optical lattice by using numerical unbiased Quantum Monte Carlo (QMC) method. Ultracold atoms experiments in optical lattices provide a possibility to simulate the materials of condensed matter in a clean and well-controled way. We focus on two components bosonic systems where exotic phases such as pair-superfluid and pair-supersolid appear. In the second part we provide numerical methods based on QMC simulations to study the entanglement properties of strong correlated systems. By employing the replica trick, we reconstruct the entanglement spectrum from the Renyi entropies. Furthermore, we study the trace of the power of the partial transposed reduced density matrix, which is related to the entanglement measurement negativity for mixed states. In the following we briefly summarize each chapter. In Chapter 1 we very briefly introduce the cold atom systems and the connection with QMC simulations. In Chapter 2 we review the QMC method based on the path integral (world line) representation. We focus on the directed worm algorithm which is an efficient algorithm with global updtae for bosonic and spin systems. All the following methods in this thesis – the method of two component bosonic systems, the reconstruction of entanglement spectrum and the measuring of partial transposed quantities, are based on directed worm algorithm introduced in this chapter. In Chapter 3 we discuss two component bosons in a square lattice. We show that the interspecies attraction and nearest-neighbor intraspecies repulsion result in the pair-supersolid phase, where a diagonal solid order coexists with an off-diagonal pair-superfluid order. The quantum and thermal transitions out of the pair-supersolid phase are characterized. It is found that there is a direct first-order transition from the pair-supersolid phase to the double-superfluid phase without an intermediate region. Furthermore, the melting of the pair-supersolid occurs in two steps. Upon heating, first the pair-superfluid is destroyed via a Kosterlitz-Thouless transition, then the solid order melts via an Ising transition. In Chapter 4 we represent a new method to reconstruct a subset of the entanglement spectrum of quantum many body systems by QMC, where the method can in principle be applied to two or higher dimension. The approach builds on the replica trick to evaluate particle number resolved traces of the first n of powers of a reduced density matrix. From this information we reconstruct first n entanglement spectrum levels using a polynomial root solver. We illustrate the power and limitations of the method by an application to the extended Bose-Hubbard model in one dimension where we are able to resolve the quasidegeneracy of the entanglement spectrum in the Haldane-insulator phase. In general, the method is able to reconstruct the largest few eigenvalues in each symmetry sector. In Chapter 5 we devise a Quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. In particular, we study several scale invariant combinations of the moments for the 1D hard-core boson model. For two adjacent intervals unusual finite size corrections are present, showing parity effects that oscillate with a filling dependent period. For large chains we find perfect agreement with conformal field theory (CFT) calculations. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. The finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderate large sizes.

•Jan 1, 2014•0 citations
Motivated by recent developments in conformal field theory (CFT), we devise a quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. These quantities can serve as witnesses of criticality. In particular, we study several scale invariant combinations of the moments for the one-dimensional (1D) hard-core boson model. For two adjacent intervals unusual finite-size corrections are present, showing parity effects that oscillate with a filling dependent period. These are more pronounced in the presence of boundaries. For large chains we find perfect agreement with CFT. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. Remarkably, the finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderately large sizes.
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.
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 study the scaling of the traces of the integer powers of the partially transposed reduced density matrix and of the entanglement negativity for two spin blocks as function of their length and separation in the critical Ising chain. For two adjacent blocks, we show that tensor network calculations agree with universal conformal field theory (CFT) predictions. In the case of two disjoint blocks the CFT predictions are recovered only after taking into account the finite size corrections induced by the finite length of the blocks.
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.
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.
This thesis is divided into two parts. In the first part we study quantum phases in bosonic systems in an optical lattice by using numerical unbiased Quantum Monte Carlo (QMC) method. Ultracold atoms experiments in optical lattices provide a possibility to simulate the materials of condensed matter in a clean and well-controled way. We focus on two components bosonic systems where exotic phases such as pair-superfluid and pair-supersolid appear. In the second part we provide numerical methods based on QMC simulations to study the entanglement properties of strong correlated systems. By employing the replica trick, we reconstruct the entanglement spectrum from the Renyi entropies. Furthermore, we study the trace of the power of the partial transposed reduced density matrix, which is related to the entanglement measurement negativity for mixed states. In the following we briefly summarize each chapter. In Chapter 1 we very briefly introduce the cold atom systems and the connection with QMC simulations. In Chapter 2 we review the QMC method based on the path integral (world line) representation. We focus on the directed worm algorithm which is an efficient algorithm with global updtae for bosonic and spin systems. All the following methods in this thesis – the method of two component bosonic systems, the reconstruction of entanglement spectrum and the measuring of partial transposed quantities, are based on directed worm algorithm introduced in this chapter. In Chapter 3 we discuss two component bosons in a square lattice. We show that the interspecies attraction and nearest-neighbor intraspecies repulsion result in the pair-supersolid phase, where a diagonal solid order coexists with an off-diagonal pair-superfluid order. The quantum and thermal transitions out of the pair-supersolid phase are characterized. It is found that there is a direct first-order transition from the pair-supersolid phase to the double-superfluid phase without an intermediate region. Furthermore, the melting of the pair-supersolid occurs in two steps. Upon heating, first the pair-superfluid is destroyed via a Kosterlitz-Thouless transition, then the solid order melts via an Ising transition. In Chapter 4 we represent a new method to reconstruct a subset of the entanglement spectrum of quantum many body systems by QMC, where the method can in principle be applied to two or higher dimension. The approach builds on the replica trick to evaluate particle number resolved traces of the first n of powers of a reduced density matrix. From this information we reconstruct first n entanglement spectrum levels using a polynomial root solver. We illustrate the power and limitations of the method by an application to the extended Bose-Hubbard model in one dimension where we are able to resolve the quasidegeneracy of the entanglement spectrum in the Haldane-insulator phase. In general, the method is able to reconstruct the largest few eigenvalues in each symmetry sector. In Chapter 5 we devise a Quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. In particular, we study several scale invariant combinations of the moments for the 1D hard-core boson model. For two adjacent intervals unusual finite size corrections are present, showing parity effects that oscillate with a filling dependent period. For large chains we find perfect agreement with conformal field theory (CFT) calculations. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. The finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderate large sizes.