Entanglement scaling in matrix product state representation of smooth functions and their shallow quantum circuit approximations
Pith reviewed 2026-05-23 07:37 UTC · model grok-4.3
The pith
Entanglement across bonds in the MPS representation of a function decays asymptotically according to the function's smoothness class, whether real or complex.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For input functions that are smooth (real or complex), the bond entanglement in their matrix-product-state representation admits rigorous asymptotic expansions governed by the smoothness class; these expansions, together with localization and support considerations, yield an improved tensor-cross-interpolation algorithm that produces shallow, accurate quantum circuits for encoding the functions.
What carries the argument
Matrix product state representation of a function, whose bond entanglement scaling is controlled by the smoothness class and is used to construct shallow quantum circuits via tensor cross interpolation.
If this is right
- The derived expansions give explicit decay rates for entanglement as a function of differentiability order or analyticity.
- Tensor cross interpolation combined with the scaling results produces circuits whose gate count scales linearly with system size while maintaining accuracy.
- Heavy-tailed distributions such as Levy distributions can be loaded and sampled on current quantum hardware up to at least 156 qubits.
- The same scaling analysis applies to both real-valued and complex-valued input functions.
Where Pith is reading between the lines
- The same smoothness-based scaling could be used to pre-select bond dimensions before circuit construction, reducing classical preprocessing cost.
- If the asymptotic expansions hold for functions on higher-dimensional domains, the method may extend to multivariate probability distributions without exponential growth in circuit depth.
- The approach offers a concrete route to test whether smoothness-induced entanglement reduction improves performance on other quantum algorithms that rely on function encoding.
Load-bearing premise
The entanglement decay rate in the MPS representation is governed by the smoothness class of the function in a manner that permits a rigorous asymptotic expansion.
What would settle it
A concrete counter-example would be a C^infty function whose measured bond entanglement in an MPS truncation fails to follow the predicted asymptotic decay rate (for instance, exponential rather than the expected polynomial or faster).
Figures
read the original abstract
Encoding classical data in a quantum state is a key prerequisite of many quantum algorithms. Recently matrix product state (MPS) methods emerged as the most promising approach for constructing shallow quantum circuits approximating input functions, including probability distributions, with only linear number of gates. We derive rigorous asymptotic expansions for the decay of entanglement across bonds in the MPS representation depending on the smoothness of the input function, real or complex. We also consider the dependence of the entanglement on localization properties and function support. Based on these analytical results we construct an improved MPS-based algorithm yielding shallow and accurate encoding quantum circuits. By using Tensor Cross Interpolation we are able to construct utility-scale quantum circuits in a compute- and memory-efficient way. We validate our methods by loading heavy-tailed distributions, including Levy, important in finance, but they apply to any smooth function inputs. We test the performance of the resulting quantum circuits by executing and sampling from them on IBM quantum devices, for up to 156 qubits.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives rigorous asymptotic expansions relating the decay of entanglement across bonds in the MPS representation of input functions (real or complex) to the smoothness class of the function, including dependence on localization and support properties. These analytic results are used to construct an improved MPS-based algorithm for shallow quantum circuit approximations of smooth functions, implemented efficiently via Tensor Cross Interpolation. The approach is validated by encoding heavy-tailed distributions (e.g., Lévy) and executing the resulting circuits on IBM quantum hardware for systems up to 156 qubits.
Significance. If the claimed expansions are rigorously established, the work supplies a parameter-free theoretical link between function smoothness and MPS entanglement scaling that directly informs circuit depth and accuracy for quantum data encoding. The hardware demonstration on utility-scale systems and the focus on finance-relevant distributions (Lévy) add practical value; the combination of analytic derivations with reproducible numerical validation on real devices is a clear strength.
minor comments (3)
- §3 (or wherever the expansions are stated): the asymptotic forms should be accompanied by explicit remainder bounds or convergence rates to make the 'rigorous' claim fully verifiable without external references.
- Figure 4 (hardware sampling results): error bars or shot statistics are not visible in the caption; adding them would clarify the statistical significance of the reported fidelities for 156-qubit circuits.
- Notation: the definition of the smoothness class (e.g., C^k vs. analytic) and the precise bond index for the entanglement cut should be restated once in the main text before the first expansion to avoid ambiguity for readers.
Simulated Author's Rebuttal
We thank the referee for their careful reading and positive evaluation of our manuscript, including the accurate summary of our analytic results on entanglement scaling and the practical demonstration on IBM hardware up to 156 qubits. The recommendation for minor revision is noted; however, the report contains no specific major comments requiring point-by-point rebuttal.
Circularity Check
No significant circularity in derivation chain
full rationale
The central claim consists of analytic derivations of asymptotic expansions relating entanglement decay rates in MPS representations directly to the smoothness class (analyticity or differentiability order) of the input function. These steps are presented as independent mathematical results from function class definitions, with no reduction to fitted parameters, self-definitional constructions, or load-bearing self-citations that would make the expansions equivalent to their inputs by construction. Numerical tests and circuit constructions are downstream applications and do not retroactively define the claimed expansions. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
Forward citations
Cited by 3 Pith papers
-
State preparation with parallel-sequential circuits
Parallel-sequential circuits provide a tunable family of quantum circuit layouts that numerically outperform brickwall, sequential, and log-depth circuits for 1D ground-state preparation under realistic noise models.
-
Minimizing entanglement entropy for enhanced quantum state preparation
A two-step method minimizes entanglement entropy of target states before using matrix product state representations to achieve high-accuracy quantum state preparation on NISQ devices.
-
Solving the Gross-Pitaevskii equation on multiple different scales using the quantics tensor train representation
A quantics tensor train solver resolves the Gross-Pitaevskii equation across seven orders of magnitude in length scale in one dimension and on grids larger than a trillion points in two dimensions.
Reference graph
Works this paper leans on
-
[1]
Synthe- sis of quantum-logic circuits
V.V. Shende, S.S. Bullock, and I.L. Markov. “Synthe- sis of quantum-logic circuits”. IEEE Transactions on Computer-Aided Design of Integrated Circuits and Sys- tems 25, 1000–1010 (2006)
work page 2006
-
[2]
Quantum generative adversarial networks for learning and loading random distributions
Christa Zoufal, Aur´ elien Lucchi, and Stefan Woerner. “Quantum generative adversarial networks for learning and loading random distributions”. npj Quantum Infor- mation 5, 103 (2019)
work page 2019
-
[4]
Quantum self-learning monte carlo and quantum-inspired fourier transform sampler
Katsuhiro Endo, Taichi Nakamura, Keisuke Fujii, and Naoki Yamamoto. “Quantum self-learning monte carlo and quantum-inspired fourier transform sampler”. Phys. Rev. Res. 2, 043442 (2020)
work page 2020
-
[5]
Linear-depth quan- tum circuits for loading fourier approximations of arbi- trary functions
Mudassir Moosa, Thomas W Watts, Yiyou Chen, Abhi- jat Sarma, and Peter L McMahon. “Linear-depth quan- tum circuits for loading fourier approximations of arbi- trary functions”. Quantum Science and Technology 9, 015002 (2023)
work page 2023
-
[6]
Loading prob- ability distributions in a quantum circuit
Kalyan Dasgupta and Binoy Paine. “Loading prob- ability distributions in a quantum circuit” (2022). arXiv:2208.13372
-
[7]
The density-matrix renormalization group in the age of matrix product states
Ulrich Schollw¨ ock. “The density-matrix renormalization group in the age of matrix product states”. Annals of Physics 326, 96–192 (2011)
work page 2011
-
[8]
Matrix product states and pro- jected entangled pair states: Concepts, symmetries, the- orems
J. Ignacio Cirac, David P´ erez-Garc´ ıa, Norbert Schuch, and Frank Verstraete. “Matrix product states and pro- jected entangled pair states: Concepts, symmetries, the- orems”. Rev. Mod. Phys. 93, 045003 (2021)
work page 2021
-
[9]
Encoding of matrix product states into quantum circuits of one- and two-qubit gates
Shi-Ju Ran. “Encoding of matrix product states into quantum circuits of one- and two-qubit gates”. Phys. Rev. A 101, 032310 (2020)
work page 2020
-
[10]
Decomposition of matrix product states into shallow quantum circuits
Manuel S Rudolph, Jing Chen, Jacob Miller, Atithi Acharya, and Alejandro Perdomo-Ortiz. “Decomposition of matrix product states into shallow quantum circuits”. Quantum Science and Technology 9, 015012 (2023)
work page 2023
-
[11]
Automati- cally differentiable quantum circuit for many-qubit state preparation
Peng-Fei Zhou, Rui Hong, and Shi-Ju Ran. “Automati- cally differentiable quantum circuit for many-qubit state preparation”. Phys. Rev. A 104, 042601 (2021)
work page 2021
-
[12]
Approximate encoding of quantum states using shallow circuits
Matan Ben-Dov, David Shnaiderov, Adi Makmal, and Emanuele G. Dalla Torre. “Approximate encoding of quantum states using shallow circuits”. npj Quantum Information 10, 65 (2024)
work page 2024
-
[13]
Preparation of matrix product states with log-depth quantum circuits
Daniel Malz, Georgios Styliaris, Zhi-Yuan Wei, and J. Ig- nacio Cirac. “Preparation of matrix product states with log-depth quantum circuits”. Phys. Rev. Lett. 132, 040404 (2024)
work page 2024
-
[14]
Constant-depth prepara- tion of matrix product states with adaptive quantum cir- cuits
Kevin C. Smith, Abid Khan, Bryan K. Clark, S. M. Girvin, and Tzu-Chieh Wei. “Constant-depth prepara- tion of matrix product states with adaptive quantum cir- cuits” (2024). arXiv:2404.16083
-
[15]
Constructive representation of functions in low-rank tensor formats
I. Oseledets. “Constructive representation of functions in low-rank tensor formats”. Constructive Approximation 37, 1 – 18 (2012). url: https://api.semanticscholar. org/CorpusID:124954724
work page 2012
-
[16]
Efficient quantum circuits for accurate state preparation of smooth, differentiable functions
Adam Holmes and Anne Y. Matsuura. “Efficient quantum circuits for accurate state preparation of smooth, differentiable functions”. 2020 IEEE Interna- tional Conference on Quantum Computing and Engi- neering (QCE)Pages 169–179 (2020). url: https://api. semanticscholar.org/CorpusID:218581252
work page 2020
-
[17]
Entanglement prop- erties of quantum superpositions of smooth, differentiable functions
Adam Holmes and A. Y. Matsuura. “Entanglement prop- erties of quantum superpositions of smooth, differentiable functions” (2020). arXiv:2009.09096
-
[18]
Juan Jos´ e Garc´ ıa-Ripoll. “Quantum-inspired algorithms for multivariate analysis: from interpolation to partial differential equations”. Quantum 5, 431 (2021)
work page 2021
-
[19]
Quantum algorithms for approximate func- 14 tion loading
Gabriel Marin-Sanchez, Javier Gonzalez-Conde, and Mikel Sanz. “Quantum algorithms for approximate func- 14 tion loading”. Phys. Rev. Res. 5, 033114 (2023)
work page 2023
-
[20]
Multiscale interpolative construction of quantized tensor trains
Michael Lindsey. “Multiscale interpolative construction of quantized tensor trains” (2024). arXiv:2311.12554
-
[21]
Deterministic and entanglement-efficient preparation of amplitude-encoded quantum registers
Prithvi Gundlapalli and Junyi Lee. “Deterministic and entanglement-efficient preparation of amplitude-encoded quantum registers”. Phys. Rev. Appl. 18, 024013 (2022)
work page 2022
-
[22]
Quantum state preparation of normal distributions us- ing matrix product states
Jason Iaconis, Sonika Johri, and Elton Yechao Zhu. “Quantum state preparation of normal distributions us- ing matrix product states”. npj Quantum Information 10, 15 (2024)
work page 2024
-
[23]
Quantum state preparation using tensor networks
Ar A Melnikov, A A Termanova, S V Dolgov, F Neukart, and M R Perelshtein. “Quantum state preparation using tensor networks”. Quantum Science and Technology 8, 035027 (2023)
work page 2023
-
[24]
A modu- lar engine for quantum monte carlo integration
Ismail Yunus Akhalwaya, Adam Connolly, Roland Guichard, Steven Herbert, Cahit Kargi, Alexandre Kra- jenbrink, Michael Lubasch, Conor Mc Keever, Julien Sorci, Michael Spranger, and Ifan Williams. “A modu- lar engine for quantum monte carlo integration” (2023). arXiv:2308.06081
-
[25]
Yuichi Sano and Ikko Hamamura. “Quantum state prepa- ration for probability distributions with mirror symmetry using matrix product states” (2024). arXiv:2403.16729
-
[26]
Efficient quantum amplitude encoding of polynomial functions
Javier Gonzalez-Conde, Thomas W. Watts, Pablo Rodriguez-Grasa, and Mikel Sanz. “Efficient quantum amplitude encoding of polynomial functions”. Quantum 8, 1297 (2024)
work page 2024
-
[27]
Stefan Woerner and Daniel J. Egger. “Quantum risk analysis”. npj Quantum Information 5 (2019)
work page 2019
-
[28]
Option pricing using quantum computers
Nikitas Stamatopoulos, Daniel J. Egger, Yue Sun, Christa Zoufal, Raban Iten, Ning Shen, and Stefan Woerner. “Option pricing using quantum computers”. Quantum 4, 291 (2020)
work page 2020
-
[29]
Th´ eorie des risques financiers
Jean-Philippe Bouchaud and Marc Potters. “Th´ eorie des risques financiers”. Al´ ea-Saclay: Eyrolles. Paris (1997)
work page 1997
-
[30]
Financial modeling with l´ evy processes and stochastic volatility
Rama Cont and Peter Tankov. “Financial modeling with l´ evy processes and stochastic volatility”. Journal of Com- putational Finance 8, 1–29 (2004)
work page 2004
-
[31]
M. Z. Bazant. “Random walks and diffusion”. MIT OpenCourseWare (2003). url: https://ocw.mit.edu/
work page 2003
-
[32]
Fractals and scaling in finance
Benoit Mandelbrot. “Fractals and scaling in finance”. Springer. New York (1997)
work page 1997
-
[33]
Maslowian portfolio the- ory: An alternative formulation of the behavioural port- folio theory
Philippe J. S. De Brouwer. “Maslowian portfolio the- ory: An alternative formulation of the behavioural port- folio theory”. Journal of Asset Management 9, 359– 365 (2009)
work page 2009
-
[34]
Target-oriented investment advice
Philippe J. S. De Brouwer. “Target-oriented investment advice”. Journal of Asset Management (2011)
work page 2011
-
[35]
Aleksander Weron and Rafa l Weron. “L´ evy processes in finance”. Physica A 299, 38–46 (2001)
work page 2001
-
[36]
Fitting l´ evy stable distributions to financial market data
Ioannis Koutrouvelis. “Fitting l´ evy stable distributions to financial market data”. Journal of Financial Econo- metrics 8, 150–170 (1980)
work page 1980
-
[37]
Scaling behavior in the dynamics of an economic index
Rosario N. Mantegna and H. Eugene Stanley. “Scaling behavior in the dynamics of an economic index”. Nature 376, 46–49 (1995)
work page 1995
-
[38]
Learning tensor networks with tensor cross interpolation: new algorithms and libraries
Yuriel N´ u˜ nez Fern´ andez, Marc K. Ritter, Matthieu Jean- nin, Jheng-Wei Li, Thomas Kloss, Thibaud Louvet, Satoshi Terasaki, Olivier Parcollet, Jan von Delft, Hi- roshi Shinaoka, and Xavier Waintal. “Learning tensor networks with tensor cross interpolation: new algorithms and libraries” (2024). arXiv:2407.02454
-
[39]
O(dlog n)-quantics approximation of n-d tensors in high-dimensional numerical modeling
Boris N. Khoromskij. “O(dlog n)-quantics approximation of n-d tensors in high-dimensional numerical modeling”. Constructive Approximation 34, 257–280 (2011)
work page 2011
-
[40]
Efficient classical simulation of slightly entangled quantum computations
Guifr´ e Vidal. “Efficient classical simulation of slightly entangled quantum computations”. Phys. Rev. Lett. 91, 147902 (2003)
work page 2003
-
[41]
Matrix Product State Representations
D. Perez-Garcia, F. Verstraete, M. M. Wolf, and J. I. Cirac. “Matrix product state representations” (2007). arXiv:quant-ph/0608197
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[42]
Tao Xiang. “Matrix product states”. Page 154–165. Cambridge University Press. (2023)
work page 2023
-
[43]
Juan Jos´ e Rodr´ ıguez-Aldavero, Paula Garc´ ıa-Molina, Luca Tagliacozzo, and Juan Jos´ e Garc´ ıa-Ripoll. “Cheby- shev approximation and composition of functions in ma- trix product states for quantum-inspired numerical anal- ysis” (2024). arXiv:2407.09609
-
[44]
Tt-cross ap- proximation for multidimensional arrays
Ivan Oseledets and Eugene Tyrtyshnikov. “Tt-cross ap- proximation for multidimensional arrays”. Linear Alge- bra and its Applications 432, 70–88 (2010)
work page 2010
-
[45]
Marc K. Ritter, Yuriel N´ u˜ nez Fern´ andez, Markus Waller- berger, Jan von Delft, Hiroshi Shinaoka, and Xavier Waintal. “Quantics tensor cross interpolation for high- resolution parsimonious representations of multivariate functions”. Phys. Rev. Lett. 132, 056501 (2024)
work page 2024
-
[46]
Fast adaptive interpolation of multi-dimensional arrays in tensor train format
Dmitry Savostyanov and Ivan Oseledets. “Fast adaptive interpolation of multi-dimensional arrays in tensor train format”. In The 2011 International Workshop on Multi- dimensional (nD) Systems. Pages 1–8. (2011)
work page 2011
-
[47]
Variational quantum al- gorithms for nonlinear problems
Michael Lubasch, Jaewoo Joo, Pierre Moinier, Martin Kiffner, and Dieter Jaksch. “Variational quantum al- gorithms for nonlinear problems”. Phys. Rev. A 101, 010301 (2020)
work page 2020
-
[48]
Cgmy and the fine structure of asset returns
Dilip Madan and Marc Yor. “Cgmy and the fine structure of asset returns”. Mathematical Finance 18, 1–20 (2008)
work page 2008
-
[49]
The fine structure of asset returns: An empirical investigation
Peter Carr, H´ elyette Geman, Dilip Madan, and Marc Yor. “The fine structure of asset returns: An empirical investigation”. Journal of Business 75, 305–332 (2002)
work page 2002
-
[50]
Quantum circuits for isometries
Raban Iten, Roger Colbeck, Ivan Kukuljan, Jonathan Home, and Matthias Christandl. “Quantum circuits for isometries”. Phys. Rev. A 93, 032318 (2016)
work page 2016
-
[51]
Finding purifications with minimal entangle- ment
Johannes Hauschild, Eyal Leviatan, Jens H. Bardar- son, Ehud Altman, Michael P. Zaletel, and Frank Poll- mann. “Finding purifications with minimal entangle- ment”. Phys. Rev. B 98, 235163 (2018)
work page 2018
-
[52]
Decomposition of or- thogonal matrix and synthesis of two-qubit and three- qubit orthogonal gates
Hai-Rui Wei and Yao-Min Di. “Decomposition of or- thogonal matrix and synthesis of two-qubit and three- qubit orthogonal gates”. Quant. Inf. Comput. 12, 0262– 0270 (2012)
work page 2012
- [53]
-
[54]
Proof of Theorem 1 To obtain the purities we first construct the reduced density matrix ρk in terms of the tensor representation: ρk = 1 2N X vk fuk,vk f ∗ u′ k,vk , (A1) where we introduced collective indices uk = σ1σ2 . . . σk, and vk = σk+1σk+2 . . . σN for the bipartition of the sys- tem (across spatial scales). Here f ∗ denotes a complex conjugate of...
-
[55]
Let the subleading part of the entanglement spec- trum be of order s
Proof of Corollaries 1 and 2 We prove both results simultaneously. Let the subleading part of the entanglement spec- trum be of order s. Then from the state normalizationP i Λ2 k,i = 1, it follows that Λ k,0 = 1 − s2/2 + O(s4). Using Eq.7 the purity has now the form: pk = (1−s2/2+O(s4))4+O(s4) = 1−2s2+O(s4). (A14) Comparing with the expansion of purity in...
-
[56]
Proof of Theorem 2 Let us consider the reduced density matrix of the first k qubits ρk and compute its particular element: ρk(uk, u′ k) = 1 2k Z 1 0 f(uk + v/2k)f ∗(u′ k + v/2k)dv, (A16) where uk, u′ k are the combined indices of the reduced density matrix. Expanding the functions f, f ∗ in v/2k we can write the reduced density matrix as ρk = 1 2k X m,n≥0...
-
[57]
Con- sider the case where the function f has only l abso- lutely continuous derivatives
Breaking the smoothness assumption Here we briefly comment on how our results change when the smoothness assumption is not satisfied. Con- sider the case where the function f has only l abso- lutely continuous derivatives. Then, by Taylor’s the- orem it can be expanded in series: f(u + v/2k) = f(u) + f ′(u)v/2k + ... + f(l)(u)vl/l!2kl + Rl(u, v), where Rl...
-
[58]
The third eigenvalue and their sum For the error analysis we compute the expression for the third entanglement coefficient Λ k,2. To the leading order for large enough k we have: Λk,2 ≈ s ⟨f(2)(u)|P1|f(2)(u)⟩ 720 × 16k = r g2(f) 720 × 16k , (B1) where we defined g2(f) = ⟨f(2)(u)|P1|f(2)(u)⟩. The sum of the remaining eigenvalues appear in com- puting fidel...
-
[59]
Exact functional coefficients for probability distributions Here we give explicit expressions for the g1(f) func- tions for the symmetrical normal ( µ = L/2), log-normal and L´ evy distributions, obtained by performing the inte- grals in Eq.A13 for these examples: g1( p p(x)normal) (B3) = L2 4σ2 1 − r L2 2σ2 exp − L2 8σ2 (Γ(1/2) − Γ(1/2, L2/8σ2)) g...
-
[60]
− Γ( 1 2 , L2 8σ2 ) − − L6 32σ6 exp h − L2 4σ2 i [Γ( 1
-
[61]
− Γ( 1 2 , L2 8σ2 )]2 (B6) If we drop the finitness of the interval, that is, consider the large L limit, the expression simplifies to: g2( p p(x)normal) = 1 8σ4 . (B7) Appendix C: Entanglement and function localization In Sec.III A we showed that the entanglement on the smallest scales exhibits universal exponential decay. We can similarly ask how the en...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.