pith. sign in

arxiv: 2412.05202 · v4 · pith:KW7INAGUnew · submitted 2024-12-06 · 🪐 quant-ph

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

classification 🪐 quant-ph
keywords matrix product statesentanglement scalingsmooth functionsquantum circuit approximationtensor cross interpolationLevy distributionsshallow circuits
0
0 comments X

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.

The paper establishes rigorous asymptotic expansions that link the rate at which entanglement decays across successive bonds in a matrix-product-state encoding to the differentiability or analyticity order of the input function. Lower entanglement for smoother functions directly permits construction of shallower quantum circuits that still approximate the encoded distribution to high accuracy. The authors incorporate localization and support properties into the same scaling analysis. From these expansions they derive an improved MPS construction algorithm that employs tensor cross interpolation to build the required circuits in a memory-efficient manner. The resulting circuits are validated by sampling heavy-tailed distributions, including Levy distributions, on IBM quantum hardware with up to 156 qubits.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2412.05202 by Georgios Korpas, Illia Lukin, Maciej Koch-Janusz, Mykola Luhanko, Mykola Maksymenko, Philippe J.S. De Brouwer, Vladyslav Bohun.

Figure 1
Figure 1. Figure 1: a) A real or complex function f(x) on [0, 1] is discretized on a binary grid x = P i σi2 −i . Successive digits σi of the binary expansion describe decreasing spatial scales: parts of the interval [0, 1] sharing the same value of σi are shown in green and blue. b) Discretized function can be represented as a rank-N tensor fσ1σ2...σN . c) A matrix product state (MPS) can approximate fσ1σ2...σN , but the bon… view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of the analytical (see Theorem [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Entanglement decay for exponential and polynomial localization: [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: An overview of our MPS based-circuit construction algorithm. First, the standard canonical form MPS is transformed [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The dependence of the infidelity on the number and [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: a) The log-log plot of infidelity of a 1-layer quantum circuit as a function of the size of the truncated support interval L. We fix the discretization step size, with the qubit number given by N = 10 + ⌈log2 L⌉. Here we consider the truncated normal distribution with σ = 2, show the numerically computed infidelity, and compare it to the analytical estimates. The label “Analytical” refers to Eq.(17) and “A… view at source ↗
Figure 7
Figure 7. Figure 7: The encoded L´evy distribution (c = 1) for different truncated support and discretization (number of qubits) sizes on an ideal noiseless simulator. Amplitudes of the quantum state created using a 2-layer encoding circuit are compared to the exact distribution. In a) for N = 6 qubits and L = 1; in b) for N = 10 qubits and L = 5. a general unitary synthesis requires, however, 3 CNOT gates and is not optimal.… view at source ↗
Figure 8
Figure 8. Figure 8: Encoding circuits execution on the ibm torino QPU. a) L´evy distribution with c = 50 truncated to L = 29 with N = 20 qubits b) L´evy distribution with c = 1310 on L = 214 with N = 25 qubits c) Log-normal distribution with on N = 25 qubits d) Gamma distribution with shape parameter k = 1 on N = 25 qubits. In each case 5000 shots were used. All executions pass the KS test with 200 samples. The better of the … view at source ↗
Figure 9
Figure 9. Figure 9: Tensor Cross Interpolation (TCI) constructed 2-layer encoding circuits for the [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: (left) The amplitudes of the quantum state encoding the L´evy distribution on 40 qubits obtained using Tensor Cross Interpolation (TCI), compared to the exact Matrix Product State (MPS) amplitudes (right) The results of sampling from the 40-qubit quantum circuit constructed using TCI approach on an ideal simulator [PITH_FULL_IMAGE:figures/full_fig_p020_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Tensor Cross Interpolation (TCI) constructed 2-layer encoding circuits for the L´evy distributions on 50 [PITH_FULL_IMAGE:figures/full_fig_p020_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The L´evy distribution encoded with 1- and [PITH_FULL_IMAGE:figures/full_fig_p021_12.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

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)
  1. §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.
  2. 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.
  3. 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

0 responses · 0 unresolved

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

0 steps flagged

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

0 free parameters · 0 axioms · 0 invented entities

The central claim rests on the existence of an MPS representation whose entanglement is controlled by function smoothness; no free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

pith-pipeline@v0.9.0 · 5728 in / 1158 out tokens · 15274 ms · 2026-05-23T07:37:03.648038+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. State preparation with parallel-sequential circuits

    quant-ph 2025-03 unverdicted novelty 7.0

    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.

  2. Minimizing entanglement entropy for enhanced quantum state preparation

    quant-ph 2025-07 unverdicted novelty 5.0

    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.

  3. Solving the Gross-Pitaevskii equation on multiple different scales using the quantics tensor train representation

    quant-ph 2025-07 unverdicted novelty 5.0

    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

60 extracted references · 60 canonical work pages · cited by 3 Pith papers · 1 internal anchor

  1. [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)

  2. [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)

  3. [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)

  4. [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)

  5. [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

  6. [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)

  7. [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)

  8. [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)

  9. [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)

  10. [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)

  11. [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)

  12. [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)

  13. [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

  14. [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

  15. [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

  16. [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

  17. [18]

    Quantum-inspired algorithms for multivariate analysis: from interpolation to partial differential equations

    Juan Jos´ e Garc´ ıa-Ripoll. “Quantum-inspired algorithms for multivariate analysis: from interpolation to partial differential equations”. Quantum 5, 431 (2021)

  18. [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)

  19. [20]

    Multiscale interpolative construction of quantized tensor trains

    Michael Lindsey. “Multiscale interpolative construction of quantized tensor trains” (2024). arXiv:2311.12554

  20. [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)

  21. [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)

  22. [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)

  23. [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

  24. [25]

    Quantum state prepa- ration for probability distributions with mirror symmetry using matrix product states

    Yuichi Sano and Ikko Hamamura. “Quantum state prepa- ration for probability distributions with mirror symmetry using matrix product states” (2024). arXiv:2403.16729

  25. [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)

  26. [27]

    Quantum risk analysis

    Stefan Woerner and Daniel J. Egger. “Quantum risk analysis”. npj Quantum Information 5 (2019)

  27. [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)

  28. [29]

    Th´ eorie des risques financiers

    Jean-Philippe Bouchaud and Marc Potters. “Th´ eorie des risques financiers”. Al´ ea-Saclay: Eyrolles. Paris (1997)

  29. [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)

  30. [31]

    Random walks and diffusion

    M. Z. Bazant. “Random walks and diffusion”. MIT OpenCourseWare (2003). url: https://ocw.mit.edu/

  31. [32]

    Fractals and scaling in finance

    Benoit Mandelbrot. “Fractals and scaling in finance”. Springer. New York (1997)

  32. [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)

  33. [34]

    Target-oriented investment advice

    Philippe J. S. De Brouwer. “Target-oriented investment advice”. Journal of Asset Management (2011)

  34. [35]

    L´ evy processes in finance

    Aleksander Weron and Rafa l Weron. “L´ evy processes in finance”. Physica A 299, 38–46 (2001)

  35. [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)

  36. [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)

  37. [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

  38. [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)

  39. [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)

  40. [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

  41. [42]

    Matrix product states

    Tao Xiang. “Matrix product states”. Page 154–165. Cambridge University Press. (2023)

  42. [43]

    Cheby- shev approximation and composition of functions in ma- trix product states for quantum-inspired numerical anal- ysis

    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

  43. [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)

  44. [45]

    Quantics tensor cross interpolation for high- resolution parsimonious representations of multivariate functions

    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)

  45. [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)

  46. [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)

  47. [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)

  48. [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)

  49. [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)

  50. [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)

  51. [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)

  52. [53]

    “torchTT”

    Ion Gabriel Ion. “torchTT”. https://github.com/ ion-g-ion/torchTT. 15 Appendix A: Proofs of theorems

  53. [54]

    σk, and vk = σk+1σk+2

    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...

  54. [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...

  55. [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...

  56. [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...

  57. [58]

    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 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...

  58. [59]

    The leading term 21L2/8c2 defines the characteristic scale k = log L/c, where the entanglement spectra become small

    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...

  59. [60]

    − Γ( 1 2 , L2 8σ2 ) − − L6 32σ6 exp h − L2 4σ2 i [Γ( 1

  60. [61]

    function values

    − Γ( 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...