We carry out a comprehensive comparison between the exact modular Hamiltonian and the lattice version of the Bisognano-Wichmann (BW) one in one-dimensional critical quantum spin chains. As a warm-up, we first illustrate how the trace distance provides a more informative mean of comparison between reduced density matrices when compared to any other Schatten n n -distance, normalized or not. In particular, as noticed in earlier works, it provides a way to bound other correlation functions in a precise manner, i.e., providing both lower and upper bounds. Additionally, we show that two close reduced density matrices, i.e. with zero trace distance for large sizes, can have very different modular Hamiltonians. This means that, in terms of describing how two states are close to each other, it is more informative to compare their reduced density matrices rather than the corresponding modular Hamiltonians. After setting this framework, we consider the ground states for infinite and periodic XX spin chain and critical Ising chain. We provide robust numerical evidence that the trace distance between the lattice BW reduced density matrix and the exact one goes to zero as \ell^{-2} ℓ − 2 for large length of the interval \ell ℓ . This provides strong constraints on the difference between the corresponding entanglement entropies and correlation functions. Our results indicate that discretized BW reduced density matrices reproduce exact entanglement entropies and correlation functions of local operators in the limit of large subsystem sizes. Finally, we show that the BW reduced density matrices fall short of reproducing the exact behavior of the logarithmic emptiness formation probability in the ground state of the XX spin chain.
We carry out a comprehensive comparison between the exact modular Hamiltonian and the lattice version of the Bisognano-Wichmann (BW) one in one-dimensional critical quantum spin chains. As a warm-up, we first illustrate how the trace distance provides a more informative mean of comparison between reduced density matrices when compared to any other Schatten n n -distance, normalized or not. In particular, as noticed in earlier works, it provides a way to bound other correlation functions in a precise manner, i.e., providing both lower and upper bounds. Additionally, we show that two close reduced density matrices, i.e. with zero trace distance for large sizes, can have very different modular Hamiltonians. This means that, in terms of describing how two states are close to each other, it is more informative to compare their reduced density matrices rather than the corresponding modular Hamiltonians. After setting this framework, we consider the ground states for infinite and periodic XX spin chain and critical Ising chain. We provide robust numerical evidence that the trace distance between the lattice BW reduced density matrix and the exact one goes to zero as \ell^{-2} ℓ − 2 for large length of the interval \ell ℓ . This provides strong constraints on the difference between the corresponding entanglement entropies and correlation functions. Our results indicate that discretized BW reduced density matrices reproduce exact entanglement entropies and correlation functions of local operators in the limit of large subsystem sizes. Finally, we show that the BW reduced density matrices fall short of reproducing the exact behavior of the logarithmic emptiness formation probability in the ground state of the XX spin chain.
We consider a lattice version of the Bisognano-Wichmann (BW) modular Hamiltonian as an ansatz for the bipartite entanglement Hamiltonian of the quantum critical chains. Using numerically unbiased methods, we check the accuracy of the BW ansatz by both comparing the BW R\'enyi entropy to the exact results and investigating the size scaling of the norm distance between the exact reduced density matrix and the BW one. Our study encompasses a variety of models, scanning different universality classes, including integrable models such as the transverse field Ising, three-state Potts and XXZ chains, and the nonintegrable bilinear-biquadratic model. We show that the R\'enyi entropies obtained via the BW ansatz properly describe the scaling properties predicted by conformal field theory. Remarkably, the BW R\'enyi entropies also faithfully capture the corrections to the conformal field theory scaling associated with the energy density operator. In addition, we show that the norm distance between the discretized BW density matrix and the exact one asymptotically goes to zero with the system size: this indicates that the BW ansatz also can be employed to predict properties of the eigenvectors of the reduced density matrices and is thus potentially applicable to other entanglement-related quantities such as negativity.
The modular (or entanglement) Hamiltonian correspondent to the half-space bipartition of a quantum state uniquely characterizes its entanglement properties. However, in the context of lattice models, its explicit form is analytically known only for the two spin chains and certain free theories in one dimension. In this work, we provide a thorough investigation of entanglement Hamiltonians in lattice models obtained via the Bisognano-Wichmann theorem, which provides an explicit functional form for the entanglement Hamiltonian itself in quantum field theory. Our study encompasses a variety of one- and two-dimensional models, supporting diverse quantum phases and critical points, and, most importantly, scanning several universality classes, including Ising, Potts, and Luttinger liquids. We carry out extensive numerical simulations based on the density matrix renormalization group method, exact diagonalization, and quantum Monte Carlo. In particular, we compare the exact entanglement properties and correlation functions to those obtained applying the Bisognano-Wichmann theorem on the lattice. We carry out this comparison on both the eigenvalues and eigenvectors of the entanglement Hamiltonian, and expectation values of correlation functions and order parameters. Our results evidence that as long as the low-energy description of the lattice model is well captured by a Lorentz-invariant quantum field theory, the Bisognano-Wichmann theorem provides a qualitatively and quantitatively accurate description of the lattice entanglement Hamiltonian. The resulting framework paves the way to direct studies of entanglement properties utilizing well-established statistical mechanics methods and experiments.
We report on a systematic replica approach to calculate the subsystem trace distance for a generic quantum field theory. This method has been recently introduced in [J. Zhang, P. Ruggiero, P. Calabrese, Phys. Rev. Lett. 122, 141602 (2019)], of which this work is a completion. The trace distance between two reduced density matrices ρA and σA is obtained from the moments tr(ρA−σA)n and taking the limit n→1 of the traces of the even powers. We focus here on the case of a subsystem consisting of a single interval of length ℓ embedded in the low lying eigenstates of a one-dimensional critical system of length L, a situation that can be studied exploiting the path integral form of the reduced density matrices of two-dimensional conformal field theories. The trace distance turns out to be a {\it scale invariant universal function} of ℓ/L. Here we complete our previous work by providing detailed derivations of all results and further new formulas for the distances between several low-lying states in two-dimensional free massless compact boson and fermion theories. Remarkably, for one special case in the bosonic theory and for another in the fermionic one, we obtain the exact trace distance, as well as the Schatten n-distance, for an interval of arbitrary length, while in generic case we have a general form for the first term in the expansion in powers of ℓ/L. The analytical predictions in conformal field theories are tested against exact numerical calculations in XX and Ising spin chains, finding perfect agreement. As a byproduct, new results are also obtained for other entanglement-related quantities, such as the relative entropy and the fidelity.
Abstract We report on a systematic replica approach to calculate the subsystem trace distance for a quantum field theory. This method has been recently introduced in [J. Zhang, P. Ruggiero and P. Calabrese, Phys. Rev. Lett.122(2019) 141602 ], of which this work is a completion. The trace distance between two reduced density matricesρAandσAis obtained from the moments tr(ρA− σA)nand taking the limitn →1 of the traces of the even powers. We focus here on the case of a subsystem consisting of a single interval of lengthℓembedded in the low lying eigenstates of a one-dimensional critical system of lengthL, a situation that can be studied exploiting the path integral form of the reduced density matrices of two-dimensional conformal field theories. The trace distance turns out to be ascale invariant universal functionofℓ/L. Here we complete our previous work by providing detailed derivations of all results and further new formulas for the distances between several low-lying states in two-dimensional free massless compact boson and fermion theories. Remarkably, for one special case in the bosonic theory and for another in the fermionic one, we obtain the exact trace distance, as well as the Schattenn-distance, for an interval of arbitrary length, while in generic case we have a general form for the first term in the expansion in powers ofℓ/L. The analytical predictions in conformal field theories are tested against exact numerical calculations in XX and Ising spin chains, finding perfect agreement. As a byproduct, new results in two-dimensional CFT are also obtained for other entanglement-related quantities, such as the relative entropy and the fidelity.
As a powerful theoretical construct, the entanglement Hamiltonian (EH) encapsulates the essential entanglement properties of a quantum many-body system. From the EH, one can extract a variety of entanglement quantities, such as entanglement entropies, negativity, and the entanglement spectrum. However, its general analytical form remains largely unknown. While the Bisognano-Wichmann theorem gives an exact EH form for Lorentz-invariant field theories, its validity on lattice systems is limited, especially when Lorentz invariance is absent. In this work, we propose a general scheme based on the lattice-Bisognano-Wichmann (LBW) ansatz and multi-replica-trick quantum Monte Carlo methods to numerically reconstruct the entanglement Hamiltonian in two-dimensional systems and systematically explore its applicability to systems without translational invariance, going beyond the original scope of the primordial Bisognano-Wichmann theorem. Various quantum phases--including gapped and gapless phases, critical points, and phases with either discrete or continuous symmetry breaking--are investigated, demonstrating the versatility of our method in reconstructing entanglement Hamiltonians. Furthermore, we find that when the entanglement boundary of a system is ordinary (i.e., free from surface anomalies), the LBW ansatz provides an accurate approximation well beyond Lorentz-invariant cases. Our work thus establishes a general framework for investigating the analytical structure of entanglement in the complex quantum many-body systems.
We develop a systematic method to calculate the trace distance between two reduced density matrices in 1+1 dimensional quantum field theories. The approach exploits the path integral representation of the reduced density matrices and an ad hoc replica trick. We then extensively apply this method to the calculation of the distance between reduced density matrices of one interval of length ℓ in eigenstates of conformal field theories. When the interval is short, using the operator product expansion of twist operators, we obtain a universal form for the leading order in ℓ of the trace distance. We compute the trace distances among the reduced density matrices of several low lying states in two-dimensional free massless boson and fermion theories. We compare our analytic conformal results with numerical calculations in XX and Ising spin chains finding perfect agreement.