Ensembles of random quantum states tunable from volume law to area law
Pith reviewed 2026-05-10 11:07 UTC · model grok-4.3
The pith
A single control parameter tunes random quantum states from volume-law to area-law entanglement.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We introduce the σ-ensembles – a family of random quantum states with only a single control parameter. By imposing a probability distribution on the eigenvalues of the successive subsystems and subsequently reconstructing a compatible global state using the matrix product state formalism, the ensemble can be tuned between volume-law and area-law entanglement behavior.
What carries the argument
The σ-ensemble constructed by choosing probability distributions for the eigenvalues of successive reduced subsystems and rebuilding the full pure state via the MPS representation.
If this is right
- Area-law states from the ensemble can be simulated efficiently on classical computers using tensor-network methods.
- The single parameter supplies a continuous dial to explore the crossover between volume-law and area-law regimes.
- The resulting states resemble the entanglement properties of typical Hamiltonian ground states more closely than Haar-random states.
- The construction sidesteps the exponential classical cost of simulating volume-law entangled states.
Where Pith is reading between the lines
- The ensembles could serve as controllable test cases for benchmarking classical simulation algorithms across different entanglement strengths.
- The same subsystem-eigenvalue approach might be adapted to enforce other scaling behaviors or to generate states with prescribed correlation functions.
- Varying the parameter could let researchers isolate how entanglement scaling alone influences the computational hardness of quantum many-body problems.
Load-bearing premise
Assigning probability distributions to the eigenvalues of successive subsystems and reconstructing via the MPS formalism produces valid global pure states whose entanglement scaling can be continuously tuned by the single parameter without introducing inconsistencies.
What would settle it
Generate states at several values of the control parameter and measure the scaling of bipartite entanglement entropy with subsystem size to test whether it changes from volume law to area law.
Figures
read the original abstract
A standard approach to generate random pure quantum states relies on sampling from the Haar measure. However, the entanglement properties of such states present a fundamental challenge for their general applicability. Here, we introduce the $\sigma$-ensembles $\unicode{x2013}$ a family of random quantum states with only a single control parameter. Crucially, these states are designed such that they can be tuned between volume-law and area-law behavior, which has been a major obstacle thus far. We construct representatives of this ensemble by imposing a probability distribution on the eigenvalues of the successive subsystems, and subsequently reconstructing a compatible global state using the matrix product state (MPS) formalism. Due to their area-law entanglement, our approach circumvents the intractability of Haar-random pure states in classical simulations of quantum systems and is more representative of typical Hamiltonian ground states.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces the σ-ensembles, a family of random pure quantum states controlled by a single parameter σ. These states are constructed by assigning probability distributions to the eigenvalues of successive subsystems and then reconstructing compatible global states via the matrix product state (MPS) formalism. The central claim is that this procedure yields an ensemble whose entanglement entropy can be tuned continuously from volume-law to area-law scaling, providing states that are both more representative of typical Hamiltonian ground states and classically simulable.
Significance. If the construction is shown to be consistent, the σ-ensembles would supply a practical, single-parameter family of random states bridging Haar-random volume-law states and area-law states. This could enable systematic numerical studies of entanglement transitions and typical properties in many-body systems without the exponential cost of full Haar sampling.
major comments (1)
- [Abstract and construction method] Abstract and construction procedure: the central claim requires that independent sampling of Schmidt eigenvalues from σ-parameterized distributions on successive subsystems, followed by MPS reconstruction, produces valid global pure states whose entanglement scaling varies continuously with σ. However, the Schmidt spectrum on subsystem 1..k must be exactly the marginal of the spectrum on 1..k+1 after contraction with the next tensor; independent draws generically violate this marginalization. The manuscript supplies no explicit consistency check, rejection/renormalization procedure, or proof that the chosen distributions remain realizable for arbitrary σ without introducing σ-dependent biases that would break the single-parameter tunability.
minor comments (1)
- [Abstract] The abstract states that the states are 'designed such that they can be tuned' but does not specify the functional form of the eigenvalue distribution or the precise manner in which σ enters the sampling.
Simulated Author's Rebuttal
We thank the referee for their positive overall assessment of the manuscript and for identifying an important technical point regarding the consistency of the σ-ensemble construction. We address the major comment below and will make corresponding revisions to improve clarity and rigor.
read point-by-point responses
-
Referee: [Abstract and construction method] Abstract and construction procedure: the central claim requires that independent sampling of Schmidt eigenvalues from σ-parameterized distributions on successive subsystems, followed by MPS reconstruction, produces valid global pure states whose entanglement scaling varies continuously with σ. However, the Schmidt spectrum on subsystem 1..k must be exactly the marginal of the spectrum on 1..k+1 after contraction with the next tensor; independent draws generically violate this marginalization. The manuscript supplies no explicit consistency check, rejection/renormalization procedure, or proof that the chosen distributions remain realizable for arbitrary σ without introducing σ-dependent biases that would break the single-parameter tunability.
Authors: We thank the referee for this observation, which correctly notes that our construction samples the eigenvalue distributions independently for each cut. The MPS reconstruction builds a valid global pure state by sequentially fixing the tensors to realize the prescribed singular values at each bond; the resulting state is always normalized and pure by construction of the MPS. While the sampled spectra are not guaranteed to be exact marginals of one another, the parameter σ directly controls the typical decay rate of the singular values at every cut, thereby tuning the entanglement entropy scaling from volume-law (small σ) to area-law (large σ) as claimed. Numerical sampling of the ensemble for a range of σ and system sizes confirms that the average entanglement entropy follows the expected scaling without introducing additional σ-dependent biases beyond those encoded in the target distributions. We will revise the manuscript to add an explicit description of the sequential reconstruction algorithm, a discussion of the marginalization property, and numerical validation of consistency and tunability in a new appendix. revision: yes
Circularity Check
No circularity: explicit construction of tunable ensembles
full rationale
The paper presents an explicit construction of the σ-ensembles by choosing probability distributions on subsystem eigenvalues (parameterized by σ) and reconstructing global states via MPS. This is framed as a design choice to achieve continuous tuning of entanglement scaling, not as a derivation or prediction from independent principles that later reduces to the inputs. No equations or steps equate a claimed result to a fitted parameter or self-referential definition by construction; the single parameter σ is an input to the distributions, and the resulting volume-to-area-law behavior is a direct consequence of that choice. The approach is self-contained with no load-bearing self-citations or hidden constraints invoked to force the outcome.
Axiom & Free-Parameter Ledger
free parameters (1)
- σ
axioms (2)
- standard math Matrix product states can represent the global pure state consistent with given subsystem eigenvalue distributions.
- domain assumption The chosen probability distributions on subsystem eigenvalues produce valid quantum states.
invented entities (1)
-
σ-ensembles
no independent evidence
Reference graph
Works this paper leans on
- [1]
-
[2]
Here,U [l+1] is left-unitary and W [l+1] is right-unitary and square, hence unitary
0W [l] (18) whereW [l] was defined in the previous iteration, or W [l] = 1 ifl= 0. Here,U [l+1] is left-unitary and W [l+1] is right-unitary and square, hence unitary. We expressM [l+1] via its reduced QR decomposition: M [l+1] =Q [l+1]R[l+1] (19) whereQ [l+1]†Q[l+1] = ˆ1dml×dml andR [l+1] is upper tri- angular. We then set X [l+1] = Σ[l]Ω[l]Q[l...
-
[3]
Guifré Vidal. Efficient simulation of one-dimensional quantum many-body systems.Physical Review Letters, 93(4), July 2004
work page 2004
-
[4]
J. Eisert. Computational difficulty of global variations in the density matrix renormalization group.Physical Review Letters, 97(26), December 2006
work page 2006
-
[5]
M B Hastings. An area law for one-dimensional quantum systems.Journal of Statistical Mechanics: Theory and Experiment, 2007(08):P08024–P08024, August 2007
work page 2007
-
[6]
Norbert Schuch, Michael M. Wolf, Frank Verstraete, and J. Ignacio Cirac. Entropy scaling and simulability by matrix product states.Physical Review Letters, 100(3), January 2008
work page 2008
- [7]
-
[8]
Ulrich Schollwöck. The density-matrix renormalization group in the age of matrix product states.Annals of Physics, 326(1):96–192, January 2011
work page 2011
-
[9]
Fisher, Vedika Khemani, Adam Nahum, and Sagar Vijay
Matthew P.A. Fisher, Vedika Khemani, Adam Nahum, and Sagar Vijay. Random quantum circuits.Annual Re- view of Condensed Matter Physics, 14(1):335–379, March 2023
work page 2023
-
[10]
Random unitaries in extremely low depth, 2025
Thomas Schuster, Jonas Haferkamp, and Hsin-Yuan Huang. Random unitaries in extremely low depth, 2025
work page 2025
-
[11]
Don N. Page. Average entropy of a subsystem.Phys. Rev. Lett., 71:1291–1294, Aug 1993
work page 1993
-
[12]
S. K. Foong and S. Kanno. Proof of page’s conjecture on the average entropy of a subsystem.Phys. Rev. Lett., 72:1148–1151, Feb 1994
work page 1994
-
[13]
Average entropy of a quantum subsys- tem.Physical Review Letters, 77(1):1–3, July 1996
Siddhartha Sen. Average entropy of a quantum subsys- tem.Physical Review Letters, 77(1):1–3, July 1996
work page 1996
- [14]
-
[15]
Michael Wolf, Frank Verstraete, Matthew Hastings, and J Cirac. Area laws in quantum systems: Mutual informa- tion and correlations.Physical review letters, 100:070502, 03 2008
work page 2008
-
[16]
The1d arealawandthe complexity ofquantum states: A combinatorial approach
Dorit Aharonov, Itai Arad, Zeph Landau, and Umesh Vazirani. The1d arealawandthe complexity ofquantum states: A combinatorial approach. pages 324–333, 2011
work page 2011
-
[17]
A spacetime area law bound on quantum correlations
Ilya Kull, Philippe Allard Guérin, and Časlav Brukner. A spacetime area law bound on quantum correlations. npj Quantum Information, 5(1), July 2019
work page 2019
-
[18]
Peter W. Shor. Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum com- puter.SIAM Journal on Computing, 26(5):1484–1509, October 1997
work page 1997
-
[19]
Lov K. Grover. A fast quantum mechanical algorithm for database search. page 212–219, 1996
work page 1996
-
[20]
Sergio Boixo, Troels Rønnow, Sergei Isakov, Zhihui Wang, David Wecker, Daniel Lidar, John Martinis, Mathew Ababei, Mohammad Amin, Alireza Carignan- Dugas, et al. Evidence for quantum annealing with more than one hundred qubits.Nature Physics, 10(3):218–224, Mar 2014
work page 2014
-
[21]
Entanglement entropy scal- ing of noisy random quantum circuits in two dimensions
Meng Zhang, Chao Wang, Shaojun Dong, Hao Zhang, Yongjian Han, and Lixin He. Entanglement entropy scal- ing of noisy random quantum circuits in two dimensions. Physical Review A, 106(5), November 2022
work page 2022
-
[22]
Zhi Li, Shengqi Sang, and Timothy H. Hsieh. Entangle- ment dynamics of noisy random circuits.Physical Review B, 107(1), January 2023
work page 2023
-
[23]
Volume-law entanglement en- tropy of typical pure quantum states.PRX Quantum, 3(3), July 2022
Eugenio Bianchi, Lucas Hackl, Mario Kieburg, Marcos Rigol, and Lev Vidmar. Volume-law entanglement en- tropy of typical pure quantum states.PRX Quantum, 3(3), July 2022
work page 2022
-
[24]
Der massbegriff in der theorie der kon- tinuierlichen gruppen.Annals of Mathematics, 34:147, 1933
Alfred Haar. Der massbegriff in der theorie der kon- tinuierlichen gruppen.Annals of Mathematics, 34:147, 1933
work page 1933
-
[25]
Publications de l’Institut de mathématique de l’Université de Strasbourg
André Weil.L’intégration dans les groupes topologiques et ses applications / par André Weil. Publications de l’Institut de mathématique de l’Université de Strasbourg. Hermann & Cie, éditeurs, Paris, 2e édition edition, 1953
work page 1953
-
[26]
Sur la mesure de haar.Comptes Rendus de l’Académie des Sciences de Paris, 211:759–762, 1940
Henri Cartan. Sur la mesure de haar.Comptes Rendus de l’Académie des Sciences de Paris, 211:759–762, 1940
work page 1940
-
[27]
Elihu Lubkin. Entropy of an n-system from its correla- tion with a k-reservoir.Journal of Mathematical Physics, 19(5):1028–1031, 05 1978
work page 1978
- [28]
- [29]
-
[30]
Cambridge University Press, 2006
Ingemar Bengtsson and Karol Zyczkowski.Geometry of Quantum States: An Introduction to Quantum Entangle- ment. Cambridge University Press, 2006
work page 2006
-
[31]
Distribution ofgconcurrence of random pure states.Phys
Valerio Cappellini, Hans-Jürgen Sommers, and Karol Ży- czkowski. Distribution ofgconcurrence of random pure states.Phys. Rev. A, 74:062322, Dec 2006
work page 2006
-
[32]
Majumdar, Oriol Bohigas, and Arul Laksh- minarayan
Satya N. Majumdar, Oriol Bohigas, and Arul Laksh- minarayan. Exact minimum eigenvalue distribution of an entangled random pure state.Journal of Statistical Physics, 131(1):33–49, February 2008
work page 2008
-
[33]
In- duced measures in the space of mixed quantum states
Karol Zyczkowski and Hans-Jürgen Sommers. In- duced measures in the space of mixed quantum states. Journal of Physics A: Mathematical and General, 34(35):7111–7125, August 2001
work page 2001
-
[34]
Olivier Giraud. Purity distribution for bipartite random pure states.Journal of Physics A: Mathematical and Theoretical, 40(49):F1053, nov 2007
work page 2007
-
[35]
Entanglement of random vectors
Marko Žnidarič. Entanglement of random vectors. Journal of Physics A: Mathematical and Theoretical, 40(3):F105, dec 2006
work page 2006
-
[36]
YangChen, Dang-ZhengLiu, andDa-ShengZhou. Small- est eigenvalue distribution of the fixed-trace laguerre beta-ensemble.Journal of Physics A: Mathematical and Theoretical, 43(31):315303, July 2010. 11
work page 2010
-
[37]
Correlations, spectra and entanglement transitions in ensembles of ma- trix product states, 2025
Hugo Lóio, Guillaume Cecile, Sarang Gopalakrishnan, Guglielmo Lami, and Jacopo De Nardis. Correlations, spectra and entanglement transitions in ensembles of ma- trix product states, 2025
work page 2025
-
[38]
Approach to typicality in many-body quantum systems.Phys
Shawn Dubey, Luciano Silvestri, Justin Finn, Sai Vin- janampathy, and Kurt Jacobs. Approach to typicality in many-body quantum systems.Phys. Rev. E, 85:011141, Jan 2012
work page 2012
-
[39]
Tao Xin, Dawei Lu, Joel Klassen, Nengkun Yu, Zhengfeng Ji, Jianxin Chen, Xian Ma, Guilu Long, Bei Zeng, and Raymond Laflamme. Quantum state tomog- raphy via reduced density matrices.Physical Review Let- ters, 118(2), January 2017
work page 2017
-
[40]
Yi-Kai Liu, Matthias Christandl, and F. Ver- straete. Quantum computational complexity of the n-representability problem: Qma complete.Phys. Rev. Lett., 98:110503, Mar 2007
work page 2007
-
[41]
E. H. Lieb and M. B. Ruskai. Proof of the strong subad- ditivity of quantum-mechanical entropy.J. Math. Phys., 14:1938–1941, 1973
work page 1938
-
[42]
A. Higuchi, A. Sudbery, and J. Szulc. One-qubit reduced states of a pure many-qubit state: Polygon inequalities. Physical Review Letters, 90(10), March 2003
work page 2003
-
[43]
Requirements for compatibility between local and multipartite quantum states, 2003
Sergey Bravyi. Requirements for compatibility between local and multipartite quantum states, 2003
work page 2003
-
[44]
Quantum marginal problem and representations of the symmetric group
Alexander Klyachko. Quantum marginal problem and representations of the symmetric group. 9 2004
work page 2004
-
[45]
Alexander A Klyachko. Quantum marginal problem and n-representability.Journal of Physics: Conference Se- ries, 36:72–86, April 2006
work page 2006
-
[46]
Burak Şahinoğlu, and Michael Walter
Matthias Christandl, M. Burak Şahinoğlu, and Michael Walter. Recoupling coefficients and quantum entropies. Annales Henri Poincaré, 19(2):385–410, December 2017
work page 2017
-
[47]
Morgan Makhina. The quantum marginal problem. Preprint, University of California at San Diego, La Jolla, CA, 2023. Available at UCSD Physics Department web- site
work page 2023
-
[48]
U. Schollwöck. The density-matrix renormalization group.Rev. Mod. Phys., 77:259–315, Apr 2005
work page 2005
- [49]
-
[50]
Benoît Collins, Ion Nechita, and Karol Życzkowski. Ran- dom graph states, maximal flow and fuss–catalan distri- butions.Journal of Physics A: Mathematical and Theo- retical, 43(27):275303, June 2010
work page 2010
-
[51]
Benoît Collins, Ion Nechita, and Karol Życzkowski. Area lawforrandomgraphstates.Journal of Physics A: Math- ematical and Theoretical, 46(30):305302, July 2013
work page 2013
-
[52]
L. E. Blumenson. A derivation of n-dimensional spher- ical coordinates.The American Mathematical Monthly, 67(1):63–66, 1960
work page 1960
-
[53]
J.L. Devore and K.N. Berk.Modern Mathematical Statis- tics with Applications. Thomson Brooks/Cole, 2007
work page 2007
-
[54]
G. M. Fikhtengol’ts.Kurs differentsial’nogo i in- tegral’nogo ischisleniya, volume I–III. Gostekhizdat, Moscow and Leningrad, 1947–1949. Also published in German asDifferential- und Integralrechnung I–III, VEB Deutscher Verlag der Wissenschaften, Berlin, 1986–1987
work page 1947
-
[55]
I. S. Gradshteyn and I. M. Ryzhik.Table of integrals, series, and products. Elsevier/Academic Press, Amster- dam, seventh edition, 2007. Translated from the Russian, Translation edited and with a preface by Alan Jeffrey and Daniel Zwillinger, With one CD-ROM (Windows, Mac- intosh and UNIX)
work page 2007
-
[56]
Ronald L. Graham, Donald Ervin Knuth, and Oren Patashnik.Concrete Mathematics: A Foundation for Computer Science. Addison-Wesley, Reading, MA, sec- ond edition, 1994
work page 1994
-
[57]
Maximally entangled multipartite states: A brief survey
Marco Enriquez, I Wintrowicz, and Karol Zyczkowski. Maximally entangled multipartite states: A brief survey. Journal of Physics: Conference Series, 698:012003, 03 2016
work page 2016
-
[58]
Michael A. Nielsen and Isaac L. Chuang.Quantum Com- putation and Quantum Information: 10th Anniversary Edition. Cambridge University Press, 2010. 12 Appendix A: Generating a set of eigenvalues for a uniform sampling on the unit sphere To generate random eigenvalues, we sample a random point on the positive orthant of the unitn-sphere and use the square of...
work page 2010
-
[59]
PDF of the eigenvalues The Cartesian coordinates(x 1, x2, . . . , xn)of a point on the unitn-sphere are related to its spherical coordinates as follows [50]: x1 = cos(φ1), x2 = sin(φ1) cos(φ2), x3 = sin(φ1) sin(φ2) cos(φ3), ... xn−1 = sin(φ1). . .sin(φ n−2) cos(φn−1), xn = sin(φ1). . .sin(φ n−2) sin(φn−1). 1 PDF of the eigenvalues 14 whereφ 1, φ2, . . . ,...
-
[60]
Expressingρin its eigenbasis,ρ=P j λj |j⟩ ⟨j|, the von Neumann entropy reduces to S=− X j λj lnλ j
Expectation value of the von Neumann entropy For a quantum state with density matrixρ, the von Neumann entropy is given by [28, p.273]: S=−Tr(ρlnρ), 2 Expectation value of the von Neumann entropy 17 wherelndenotes the natural logarithm of a matrix. Expressingρin its eigenbasis,ρ=P j λj |j⟩ ⟨j|, the von Neumann entropy reduces to S=− X j λj lnλ j. For a ra...
-
[61]
Therefore, A(s) =π(1−s) −a+ 1 2 ln(1−s) . Substituting this expression back intoF2, we obtain: F2 = 1 π Z 1 0 (1−λ 1)p λ1(1−λ 1) −a+ ln(1−λ 1)1 2 dλ1 = 1 π −a Z 1 0 r 1−λ 1 λ1 dλ1 + 1 2 Z 1 0 (1−λ 1)p λ1(1−λ 1) ln(1−λ 1) dλ1 = 1 π − aπ 2 + 1 2 Z 1 0 (1−λ 1)p λ1(1−λ 1) ln(1−λ 1) dλ1 . Using the substitutionλ= 1−λ 1 in the remaining integral, we havedλ=−dλ1...
-
[62]
We refer to this point as themaximally entangled point
Point on then-sphere corresponding to the maximally entangled state In this appendix, we determine the coordinates of the point on then-sphere that leads to the maximally entangled state. We refer to this point as themaximally entangled point. A bipartite pure stateρ∈ H A ⊗ HB is maximally entangled if the reduced density matrices of both subsystems are g...
-
[63]
Probability distribution of the eigenvalues Sampling eigenvalues in the vicinity of the maximally entangled point ensures that they remain close to those of a maximally entangled state and therefore exhibit a high degree of entanglement, characteristic of a volume law. In this appendix, we describe how sets of eigenvalues are sampled around the maximally ...
-
[64]
M σl=d−1 ,(E1) wheredagain denotes the local Hilbert space dimension
Stack the site matricesMσl row wise to obtain M [l] = M σl=0 ... M σl=d−1 ,(E1) wheredagain denotes the local Hilbert space dimension
-
[65]
Perform a SVDM[l] =U ˜S[l]V [l]
-
[66]
Compute the distanceϵl =∥S [l] − ˜S[l]∥measuring the deviation between the actual Schmidt values of the MPS and the targeted Schmidt values
-
[67]
V σl=d−1 ,(E2) and replace the current site tensorsVσl ←[M σl
Vertically split V [l] = V σl=0 ... V σl=d−1 ,(E2) and replace the current site tensorsVσl ←[M σl
-
[68]
ContractU S [l] into the next site tensorsAσl−1 U S[l] ←[A σl−1. Once the site tensorsAσ1 have been updated, the total distance to the target Schmidt value distributionsϵ=PL−1 l=1 ϵl is evaluated. If the convergence threshold is reached, i.e.,ϵ < δ, the procedure can be terminated. Otherwise, continue by sweeping reverting the sweep direction and repeat t...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.