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:1212.5269 · Physical Review Letters

Entanglement at a Two-Dimensional Quantum Critical Point: A Numerical Linked-Cluster Expansion Study

Ann B. Kallin, Katharine Hyatt, Rajiv R. P. Singh, Roger G. Melko

Published Mar 26, 2013• 112 citations
arXivHTMLPDF

Abstract

We develop a method to calculate the bipartite entanglement entropy of quantum models, in the thermodynamic limit, using a numerical linked-cluster expansion (NLCE) involving only rectangular clusters. It is based on exact diagonalization of all n×m rectangular clusters at the interface between entangled subsystems A and B. We use it to obtain the Renyi entanglement entropy of the two-dimensional transverse field Ising model, for arbitrary real Renyi index α. Extrapolating these results as a function of the order of the calculation, we obtain universal pieces of the entanglement entropy associated with lines and corners at the quantum critical point. They show NLCE to be one of the few methods capable of accurately calculating universal properties of arbitrary Renyi entropies at higher dimensional critical points.

DOI
10.1103/physrevlett.110.135702

Related papers

Thermodynamic singularities in the entanglement entropy at a two-dimensional quantum critical point

Rajiv R. P. Singh, Roger G. Melko, J. Oitmaa

We study the bipartite entanglement entropy of the 2D transverse-field Ising model in the thermodynamic limit. Series expansions are developed for the Renyi entropy around both the small-field and large-field limits, allowing the separate calculation of the entanglement associated with lines and corners at the boundary between subsystems. Series extrapolations are used to extract power laws and logarithmic singularities as the quantum critical point is approached, giving access to new universal quantities. In 1D, we find excellent agreement with exact results as well as quantum Monte Carlo simulations. In 2D, we find compelling evidence that the entanglement at a corner is significantly different from a free boson field theory. These results demonstrate the power of the series expansion method for calculating entanglement entropy in interacting systems, a fact that will be particularly useful in searches for exotic quantum criticality in models with and without the sign problem.

arXiv:1204.1340•Physical Review B•Aug 6, 2012•54 citations

Corner contribution to the entanglement entropy of strongly interacting O(2) quantum critical systems in 2+1 dimensions

Edwin Miles Stoudenmire, Peter Gustainis, Ravi Johal +2 more

In a $D=2+1$ quantum critical system, the entanglement entropy across a boundary with a corner contains a subleading logarithmic scaling term with a universal coefficient. It has been conjectured that this coefficient is, to leading order, proportional to the number of field components $N$ in the associated $\text{O}(N)$ continuum ${\ensuremath{\phi}}^{4}$ field theory. Using density matrix renormalization group calculations combined with the powerful numerical linked cluster expansion technique, we confirm this scenario for the $\text{O}(2)$ Wilson-Fisher fixed point in a striking way, through direct calculation at the quantum critical points of two very different microscopic models. The value of this corner coefficient is, to within our numerical precision, twice the coefficient of the Ising fixed point. Our results add to the growing body of evidence that this universal term in the R\'enyi entanglement entropy reflects the number of low-energy degrees of freedom in a system, even for strongly interacting theories.

arXiv:1409.6327•Physical Review B•Dec 2, 2014•62 citations

Precise computation of universal corner entanglement entropy at 2+1 dimension: From Ising to Gaussian quantum critical points

Ben Lee-Yeung Ngai, Justin Tim-Lok Chau, Junchen Rong +3 more

Computing the subleading logarithmic term in the entanglement entropy (EE) of (2+1)d quantum many-body systems remains a significant challenge, despite its central role in revealing universal information about quantum states and quantum critical points (QCPs). Building on recent algorithmic advances that enable the stable calculation of EE as an exponential observable~\cite{zhouIncremental2024,zhangIntegral2024,liaoExtracting2024}, we develop a {\it bubble basis} projector quantum Monte Carlo (QMC) algorithm to precisely and efficiently compute the universal corner of EE at QCPs in a (2+1)d square-lattice transverse-field Ising model augmented with a four-body interaction. Turning on this interaction allows us to trace an Ising critical line, reaching the tricritical point, and then a line of first-order phase transition. In (2+1)d, the tricritical point is described by the Gaussian theory, where a theoretical calculation of the corner logarithmic term in the 2nd Rényi entropy term is available~\cite{UniversalCasini2007}. Our QMC results are in quantitative agreement with this theoretical value, providing a highly nontrivial benchmark of the algorithm. Furthermore, we also study the Rényi EE at the Ising critical line and on the first-order transition line, obtaining results consistent with theoretical expectations. These findings establish the long-sought connection between the universal values of an exactly solvable limit and those of a strongly correlated regime at (2+1)d.

cond-mat.str-elCross-listedarXiv:2512.00382•Reports on Progress in Physics•Dec 2, 2025•0 citations

Unusual corrections to scaling and convergence of universal Renyi properties at quantum critical points

Sharmistha Sahoo, Edwin Miles Stoudenmire, Jean-Marie Stéphan +3 more

At a quantum critical point, bipartite entanglement entropies have universal quantities which are subleading to the ubiquitous area law. For Renyi entropies, these terms are known to be similar to the von Neumann entropy, while being much more amenable to numerical and even experimental measurement. We show here that when calculating universal properties of Renyi entropies, it is important to account for unusual corrections to scaling that arise from relevant local operators present at the conical singularity in the multisheeted Riemann surface. These corrections grow in importance with increasing Renyi index. We present studies of Renyi correlation functions in the $1+1$ transverse-field Ising model (TFIM) using conformal field theory, mapping to free fermions, and series expansions, and the logarithmic entropy singularity at a corner in $2+1$ for both free bosonic field theory and the TFIM, using numerical linked cluster expansions. In all numerical studies, accurate results are only obtained when unusual corrections to scaling are taken into account. In the worst case, an analysis ignoring these corrections can get qualitatively incorrect answers, such as predicting a decrease in critical exponents with the Renyi index, when they are actually increasing. We discuss a two-step extrapolation procedure that can be used to account for the unusual corrections to scaling.

arXiv:1509.00468•Physical review. B./Physical review. B•Feb 12, 2016•20 citations

Quantum quench dynamics in the transverse-field Ising model: A numerical expansion in linked rectangular clusters

Jonas Richter, Tjark Heitmann, Robin Steinigeweg

We study quantum quenches in the transverse-field Ising model defined on different lattice geometries such as chains, two- and three-leg ladders, and two-dimensional square lattices. Starting from fully polarized initial states, we consider the dynamics of the transverse and the longitudinal magnetization for quenches to weak, strong, and critical values of the transverse field. To this end, we rely on an efficient combination of numerical linked cluster expansions (NLCEs) and a forward propagation of pure states in real time. As a main result, we demonstrate that NLCEs comprising solely rectangular clusters provide a promising approach to study the real-time dynamics of two-dimensional quantum many-body systems directly in the thermodynamic limit. By comparing to existing data from the literature, we unveil that NLCEs yield converged results on time scales which are competitive to other state-of-the-art numerical methods.

arXiv:2005.03104•SciPost Physics•Sep 3, 2020•28 citations

Entanglement at a two-dimensional quantum critical point: aT= 0 projector quantum Monte Carlo study

Stephen Inglis, Roger G. Melko

Although the leading-order scaling of entanglement entropy is non-universal at a quantum critical point (QCP), sub-leading scaling can contain universal behaviour. Such universal quantities are commonly studied in non-interacting field theories, however it typically requires numerical calculation to access them in interacting theories. In this paper, we use large-scale T = 0 quantum Monte Carlo simulations to examine in detail the second Rényi entropy of entangled regions at the QCP in the transverse-field Ising model in 2 + 1 space–time dimensions—a fixed point for which there is no exact result for the scaling of entanglement entropy. We calculate a universal coefficient of a vertex-induced logarithmic scaling for a polygonal entangled subregion, and compare the result to interacting and non-interacting theories. We also examine the shape-dependence of the Rényi entropy for finite-size toroidal lattices divided into two entangled cylinders by smooth boundaries. Remarkably, we find that the dependence on cylinder length follows a shape-dependent function calculated previously by Stephan et al (2013 New J. Phys. 15 015004) at the QCP corresponding to the 2 + 1 dimensional quantum Lifshitz free scalar field theory. The quality of the fit of our data to this scaling function, as well as the apparent cutoff-independent coefficient that results, presents tantalizing evidence that this function may reflect universal behaviour across these and other very disparate QCPs in 2 + 1 dimensional systems.

arXiv:1305.1069•New Journal of Physics•Jul 26, 2013•47 citations
We develop a method to calculate the bipartite entanglement entropy of quantum models, in the thermodynamic limit, using a numerical linked-cluster expansion (NLCE) involving only rectangular clusters. It is based on exact diagonalization of all n×m rectangular clusters at the interface between entangled subsystems A and B. We use it to obtain the Renyi entanglement entropy of the two-dimensional transverse field Ising model, for arbitrary real Renyi index α. Extrapolating these results as a function of the order of the calculation, we obtain universal pieces of the entanglement entropy associated with lines and corners at the quantum critical point. They show NLCE to be one of the few methods capable of accurately calculating universal properties of arbitrary Renyi entropies at higher dimensional critical points.
We study the bipartite entanglement entropy of the 2D transverse-field Ising model in the thermodynamic limit. Series expansions are developed for the Renyi entropy around both the small-field and large-field limits, allowing the separate calculation of the entanglement associated with lines and corners at the boundary between subsystems. Series extrapolations are used to extract power laws and logarithmic singularities as the quantum critical point is approached, giving access to new universal quantities. In 1D, we find excellent agreement with exact results as well as quantum Monte Carlo simulations. In 2D, we find compelling evidence that the entanglement at a corner is significantly different from a free boson field theory. These results demonstrate the power of the series expansion method for calculating entanglement entropy in interacting systems, a fact that will be particularly useful in searches for exotic quantum criticality in models with and without the sign problem.
In a D=2+1D=2+1D=2+1 quantum critical system, the entanglement entropy across a boundary with a corner contains a subleading logarithmic scaling term with a universal coefficient. It has been conjectured that this coefficient is, to leading order, proportional to the number of field components NNN in the associated O(N)\text{O}(N)O(N) continuum ${\ensuremath{\phi}}^{4}$ field theory. Using density matrix renormalization group calculations combined with the powerful numerical linked cluster expansion technique, we confirm this scenario for the O(2)\text{O}(2)O(2) Wilson-Fisher fixed point in a striking way, through direct calculation at the quantum critical points of two very different microscopic models. The value of this corner coefficient is, to within our numerical precision, twice the coefficient of the Ising fixed point. Our results add to the growing body of evidence that this universal term in the R\'enyi entanglement entropy reflects the number of low-energy degrees of freedom in a system, even for strongly interacting theories.
Computing the subleading logarithmic term in the entanglement entropy (EE) of (2+1)d quantum many-body systems remains a significant challenge, despite its central role in revealing universal information about quantum states and quantum critical points (QCPs). Building on recent algorithmic advances that enable the stable calculation of EE as an exponential observable~\cite{zhouIncremental2024,zhangIntegral2024,liaoExtracting2024}, we develop a {\it bubble basis} projector quantum Monte Carlo (QMC) algorithm to precisely and efficiently compute the universal corner of EE at QCPs in a (2+1)d square-lattice transverse-field Ising model augmented with a four-body interaction. Turning on this interaction allows us to trace an Ising critical line, reaching the tricritical point, and then a line of first-order phase transition. In (2+1)d, the tricritical point is described by the Gaussian theory, where a theoretical calculation of the corner logarithmic term in the 2nd Rényi entropy term is available~\cite{UniversalCasini2007}. Our QMC results are in quantitative agreement with this theoretical value, providing a highly nontrivial benchmark of the algorithm. Furthermore, we also study the Rényi EE at the Ising critical line and on the first-order transition line, obtaining results consistent with theoretical expectations. These findings establish the long-sought connection between the universal values of an exactly solvable limit and those of a strongly correlated regime at (2+1)d.
At a quantum critical point, bipartite entanglement entropies have universal quantities which are subleading to the ubiquitous area law. For Renyi entropies, these terms are known to be similar to the von Neumann entropy, while being much more amenable to numerical and even experimental measurement. We show here that when calculating universal properties of Renyi entropies, it is important to account for unusual corrections to scaling that arise from relevant local operators present at the conical singularity in the multisheeted Riemann surface. These corrections grow in importance with increasing Renyi index. We present studies of Renyi correlation functions in the 1+11+11+1 transverse-field Ising model (TFIM) using conformal field theory, mapping to free fermions, and series expansions, and the logarithmic entropy singularity at a corner in 2+12+12+1 for both free bosonic field theory and the TFIM, using numerical linked cluster expansions. In all numerical studies, accurate results are only obtained when unusual corrections to scaling are taken into account. In the worst case, an analysis ignoring these corrections can get qualitatively incorrect answers, such as predicting a decrease in critical exponents with the Renyi index, when they are actually increasing. We discuss a two-step extrapolation procedure that can be used to account for the unusual corrections to scaling.
We study quantum quenches in the transverse-field Ising model defined on different lattice geometries such as chains, two- and three-leg ladders, and two-dimensional square lattices. Starting from fully polarized initial states, we consider the dynamics of the transverse and the longitudinal magnetization for quenches to weak, strong, and critical values of the transverse field. To this end, we rely on an efficient combination of numerical linked cluster expansions (NLCEs) and a forward propagation of pure states in real time. As a main result, we demonstrate that NLCEs comprising solely rectangular clusters provide a promising approach to study the real-time dynamics of two-dimensional quantum many-body systems directly in the thermodynamic limit. By comparing to existing data from the literature, we unveil that NLCEs yield converged results on time scales which are competitive to other state-of-the-art numerical methods.
Although the leading-order scaling of entanglement entropy is non-universal at a quantum critical point (QCP), sub-leading scaling can contain universal behaviour. Such universal quantities are commonly studied in non-interacting field theories, however it typically requires numerical calculation to access them in interacting theories. In this paper, we use large-scale T = 0 quantum Monte Carlo simulations to examine in detail the second Rényi entropy of entangled regions at the QCP in the transverse-field Ising model in 2 + 1 space–time dimensions—a fixed point for which there is no exact result for the scaling of entanglement entropy. We calculate a universal coefficient of a vertex-induced logarithmic scaling for a polygonal entangled subregion, and compare the result to interacting and non-interacting theories. We also examine the shape-dependence of the Rényi entropy for finite-size toroidal lattices divided into two entangled cylinders by smooth boundaries. Remarkably, we find that the dependence on cylinder length follows a shape-dependent function calculated previously by Stephan et al (2013 New J. Phys. 15 015004) at the QCP corresponding to the 2 + 1 dimensional quantum Lifshitz free scalar field theory. The quality of the fit of our data to this scaling function, as well as the apparent cutoff-independent coefficient that results, presents tantalizing evidence that this function may reflect universal behaviour across these and other very disparate QCPs in 2 + 1 dimensional systems.