Detecting the Largest Correlations using the Correlation Density Matrix: a Quantum Monte Carlo Approach
Pith reviewed 2026-05-25 08:15 UTC · model grok-4.3
The pith
The correlation density matrix between two small subsystems reveals the dominant correlations in large many-body systems.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The correlation density matrix between two small subsystems embedded in the full sample allows detection and computation of the most dominant correlations for many-body systems without prior knowledge, as shown through benchmarks on quantum Ising and Heisenberg models.
What carries the argument
The correlation density matrix measured between two small subsystems, which captures and ranks the dominant correlations of the large system.
Load-bearing premise
That the correlation density matrix between two small subsystems captures the dominant correlations of the entire large system.
What would settle it
Observing that in a system with known dominant long-range order, the correlation density matrix from small subsystems fails to identify it as the largest correlation.
Figures
read the original abstract
We present a quantum Monte Carlo-based approach to detect and compute the most dominant correlations for many-body systems without prior knowledge. It is based on the measurement and analysis of the correlation density matrix between two (small) subsystems embedded in the full (large) sample. In order to benchmark this procedure, we investigate zero-temperature quantum phase transitions in one- and two-dimensional quantum Ising model as well as the two-dimensional bilayer Heisenberg antiferromagnet. The method paves the way for a systematic identification of unknown or exotic order parameters in unexplored phases on large systems accessible to quantum Monte Carlo methods.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a quantum Monte Carlo method to detect and rank the dominant correlations in many-body systems without prior knowledge of the order parameter. The approach relies on constructing and analyzing the correlation density matrix between two small subsystems embedded in a larger sample, with benchmarks on zero-temperature quantum phase transitions in the 1D and 2D quantum Ising models as well as the 2D bilayer Heisenberg antiferromagnet.
Significance. If the local two-subsystem correlation density matrix measurement reliably identifies the globally dominant correlations, the method would enable systematic exploration of unknown or exotic ordered phases in large systems accessible to QMC, providing a practical tool beyond conventional order-parameter searches.
major comments (1)
- [Abstract] Abstract and benchmarking description: the claim that measurements on any two small embedded subsystems suffice to rank the largest correlations present anywhere in the system is supported only by benchmarks on models whose order parameters produce clear local signatures (Ising and bilayer Heisenberg); no demonstration or counter-example analysis is given for cases where the dominant correlations are non-local, multi-point, or otherwise not reducible to two-subsystem operator correlations.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive feedback on our manuscript. We address the major comment below.
read point-by-point responses
-
Referee: [Abstract] Abstract and benchmarking description: the claim that measurements on any two small embedded subsystems suffice to rank the largest correlations present anywhere in the system is supported only by benchmarks on models whose order parameters produce clear local signatures (Ising and bilayer Heisenberg); no demonstration or counter-example analysis is given for cases where the dominant correlations are non-local, multi-point, or otherwise not reducible to two-subsystem operator correlations.
Authors: We agree with the referee that our benchmarks are restricted to models (1D/2D quantum Ising and bilayer Heisenberg) whose dominant correlations have clear local two-point character that is directly captured by the correlation density matrix (CDM) between two small subsystems. The method as formulated extracts the dominant correlation operators from the spectrum of this two-subsystem CDM; it therefore identifies the largest correlations that can be expressed as operators acting on the chosen pair of subsystems. We have not provided explicit demonstrations or counter-examples for cases in which the globally dominant correlations are inherently non-local, require multi-point operators, or cannot be reduced to two-subsystem correlations. We will revise the abstract to clarify the scope of the claim and add a dedicated paragraph in the discussion section that states the current limitations and the conditions under which the two-subsystem CDM is expected to recover the leading correlations. revision: yes
Circularity Check
No circularity: method is a direct QMC measurement procedure validated on benchmarks
full rationale
The paper introduces a computational procedure to measure the correlation density matrix between two small embedded subsystems and use it to rank dominant correlations. This is framed as an empirical detection method benchmarked on standard models (1D/2D Ising, bilayer Heisenberg) with known order parameters. No derivation chain reduces a claimed prediction or result to a fitted input or self-citation by construction; the approach relies on standard QMC sampling rather than any self-referential definition or ansatz smuggling. The central claim about sufficiency of local subsystem measurements is presented as a testable hypothesis, not a tautology.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
P. W. Anderson, More Is Different, Science 177, 393 (1972)
work page 1972
-
[2]
J. Carrasquilla and R. G. Melko, Machine learning phases of matter, Nature Physics 13, 431 (2017)
work page 2017
-
[3]
E. P. L. van Nieuwenburg, Y.-H. Liu, and S. D. Huber, Learning phase transitions by confusion, Nature Physics 13, 435 (2017)
work page 2017
-
[4]
W. Hu, R. R. P. Singh, and R. T. Scalettar, Discovering phases, phase transitions, and crossovers through unsu- pervised machine learning: A critical examination, Phys. 13 Rev. E 95, 062122 (2017)
work page 2017
- [5]
- [6]
-
[7]
Laflorencie, Quantum entanglement in condensed matter systems, Physics Reports 646, 1 (2016)
N. Laflorencie, Quantum entanglement in condensed matter systems, Physics Reports 646, 1 (2016)
work page 2016
-
[8]
S.-A. Cheong and C. L. Henley, Many-body density ma- trices on a two-dimensional square lattice: Noninteract- ing and strongly interacting spinless fermions, Physical Review B 74, 165121 (2006)
work page 2006
-
[9]
S.-A. Cheong and C. L. Henley, Correlation density ma- trix: An unbiased analysis of exact diagonalizations, Physical Review B 79, 212402 (2009)
work page 2009
-
[10]
S. Furukawa, G. Misguich, and M. Oshikawa, Systematic Derivation of Order Parameters through Reduced Den- sity Matrices, Physical Review Letters96, 047211 (2006)
work page 2006
-
[11]
S.-J. Gu, W. C. Yu, and H.-Q. Lin, Construct order pa- rameters from the reduced density matrix spectra, An- nals of Physics 336, 118 (2013)
work page 2013
-
[12]
J. Sudan, Uncovering the Physics of Frustrated Quantum Magnets Using the Correlation Density Matrix Approach, Ph.D. thesis, EPFL, Lausanne (2010)
work page 2010
-
[13]
W. M¨ under, A. Weichselbaum, A. Holzner, J. von Delft, and C. L. Henley, Correlation density matrices for 1- dimensional quantum chains based on the density ma- trix renormalization group, New Journal of Physics 12, 075027 (2010)
work page 2010
-
[14]
A. W. Sandvik, Stochastic series expansion method with operator-loop update, Physical Review B 59, R14157 (1999)
work page 1999
-
[15]
A. W. Sandvik, Stochastic series expansion method for quantum Ising models with arbitrary interactions, Phys- ical Review E 68, 056701 (2003)
work page 2003
-
[16]
L. Wang, K. S. D. Beach, and A. W. Sandvik, High- precision finite-size scaling analysis of the quantum- critical point of S = 1 / 2 Heisenberg antiferromagnetic bilayers, Physical Review B 73, 014431 (2006)
work page 2006
- [17]
- [18]
-
[19]
For instance, using translation symmetry, the RDM could be computed up to a 24-site spin-1/2 chain in Ref. 17
-
[20]
H. G. Evertz, The loop algorithm, Advances in Physics 52, 1 (2003)
work page 2003
-
[21]
W. Marshall and R. E. Peierls, Antiferromagnetism, Pro- ceedings of the Royal Society of London. Series A. Math- ematical and Physical Sciences 232, 48 (1997)
work page 1997
-
[22]
E. Lieb, T. Schultz, and D. Mattis, Two soluble models of an antiferromagnetic chain, Annals of Physics 16, 407 (1961)
work page 1961
-
[23]
K. Harada, Bayesian inference in the scaling analysis of critical phenomena, Physical Review E 84, 056704 (2011)
work page 2011
-
[24]
P. Di Francesco, P. Mathieu, and D. S´ en´ echal,Conformal Field Theory , Graduate Texts in Contemporary Physics (Springer, New York, NY, 1997)
work page 1997
- [25]
-
[26]
The distinction should be clear from the context
The critical exponent related to the order parameter is commonly denoted as β, not to be confused with the inverse temperature. The distinction should be clear from the context
-
[27]
Kardar, Statistical Physics of Fields (Cambridge Uni- versity Press, 2007)
M. Kardar, Statistical Physics of Fields (Cambridge Uni- versity Press, 2007)
work page 2007
-
[28]
Sachdev, Quantum Phase Transitions, 2nd ed
S. Sachdev, Quantum Phase Transitions, 2nd ed. (Cam- bridge University Press, 2011)
work page 2011
-
[29]
Y. Zou, A. Milsted, and G. Vidal, Conformal fields and operator product expansion in critical quantum spin chains, Physical Review Letters 124, 040604 (2020)
work page 2020
- [30]
-
[31]
D. Simmons-Duffin, The lightcone bootstrap and the spectrum of the 3d Ising CFT, Journal of High Energy Physics 2017, 86 (2017)
work page 2017
- [32]
-
[33]
M. Campostrini, M. Hasenbusch, A. Pelissetto, P. Rossi, and E. Vicari, Critical exponents and equation of state of the three-dimensional Heisenberg universality class, Physical Review B 65, 144520 (2002)
work page 2002
-
[34]
M. den Nijs and K. Rommelse, Preroughening transitions in crystal surfaces and valence-bond phases in quantum spin chains, Phys. Rev. B 40, 4709 (1989)
work page 1989
-
[35]
Wen, Colloquium: Zoo of quantum-topological phases of matter, Rev
X.-G. Wen, Colloquium: Zoo of quantum-topological phases of matter, Rev. Mod. Phys. 89, 041004 (2017)
work page 2017
-
[36]
F. Pauw, F. A. Palm, U. Schollw¨ ock, A. Bohrdt, S. Paeckel, and F. Grusdt, Detecting hidden order in fractional Chern insulators, Phys. Rev. Res. 6, 023180 (2024)
work page 2024
-
[37]
J. A. Mydosh and P. M. Oppeneer, Colloquium: Hidden order, superconductivity, and magnetism: The unsolved case of URu 2Si2, Rev. Mod. Phys. 83, 1301 (2011)
work page 2011
-
[38]
P. Young, Everything You Wanted to Know About Data Analysis and Fitting but Were Afraid to Ask , Springer- Briefs in Physics (Springer International Publishing, Cham, 2015). Appendix A: Appendix
work page 2015
-
[39]
Decomposition of the correlation density matrix We reshape the correlation density matrix to separate the indices on A and B: (ρc AB)i,j = (ρc AB)ab,a′b′ → ˜ρaa′,bb′ (A1) 14 Using a singular value decomposition, we can write: ˜ρ = UΣV † (A2) X(A) i a,a′ = Uaa′,i (A3) Y (B) j † b,b′ = V † j,bb′ (A4) σAB i = Σii (A5) The operators are orthonormal under the ...
-
[40]
Details of simulations and analysis a. Simulation details Since we are interested in the groundstate properties of the models considered, we seek to suppress Monte Carlo moves that correspond to energy changes greater than the finite-size energy gap. It is, therefore, important to choose the value of the inverse temperature β appropri- ately and scale it ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.