Continuous matrix product operators for quantum fields
Pith reviewed 2026-05-18 00:55 UTC · model grok-4.3
The pith
Continuous matrix product operators for quantum fields admit a closed-form expression in finite matrix-valued functions and preserve the entanglement area law in the continuum.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors present an ansatz for continuous matrix product operators that admits a closed-form expression in terms of a finite number of matrix-valued functions without any lattice parameter. The ansatz is obtained as a suitable continuum limit of matrix product operators, preserves the entanglement area law directly in the continuum, and maps continuous matrix product states to other continuous matrix product states. As an application, several families of continuous matrix product unitaries are constructed that go beyond quantum cellular automata.
What carries the argument
The ansatz for continuous matrix product operators, which are operators built from a finite collection of matrix-valued functions on the real line that act on continuous matrix product states while preserving the area law.
Load-bearing premise
A suitable continuum limit of discrete matrix product operators exists that yields operators expressible with a finite number of matrix-valued functions while exactly preserving the entanglement area law without residual lattice artifacts or divergences.
What would settle it
A explicit calculation for a free scalar field or Ising chain in which the continuum limit either requires an infinite number of matrix functions or produces an O(1) violation of the area law would falsify the claim.
read the original abstract
In this work we introduce an ansatz for continuous matrix product operators for quantum field theory. We show that (i) they admit a closed-form expression in terms of finite number of matrix-valued functions without reference to any lattice parameter; (ii) they are obtained as a suitable continuum limit of matrix product operators; (iii) they preserve the entanglement area law directly in the continuum, and in particular they map continuous matrix product states (cMPS) to another cMPS. As an application, we use this ansatz to construct several families of continuous matrix product unitaries beyond quantum cellular automata.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces an ansatz for continuous matrix product operators (cMPOs) suitable for quantum field theory. It claims that these operators admit a closed-form expression in terms of a finite number of matrix-valued functions with no reference to lattice spacing, are obtained via a suitable continuum limit of discrete matrix product operators, and preserve the entanglement area law directly in the continuum while mapping continuous matrix product states (cMPS) to other cMPS. The ansatz is applied to construct families of continuous matrix product unitaries beyond quantum cellular automata.
Significance. If the central constructions are rigorous, the work would provide a useful extension of tensor-network methods to continuum quantum fields, enabling operator representations that respect area-law entanglement without lattice artifacts. This could support variational approaches in QFT and the study of continuous unitary dynamics.
major comments (1)
- [continuum limit section] The continuum limit construction (detailed in the section following the ansatz definition): the claim that the resulting cMPO is exactly free of lattice artifacts and maps cMPS to cMPS requires an explicit demonstration that the partial trace over auxiliary degrees of freedom commutes with the limit; without a specified scaling of bond dimension with lattice spacing a or a check against a known solvable limit (e.g., free scalar field), residual a-dependent terms in the entanglement entropy cannot be ruled out.
minor comments (2)
- [Abstract] The abstract states that 'several families' of unitaries are constructed but provides no count or distinguishing features; adding a brief enumeration or reference to the relevant subsection would improve clarity.
- [ansatz definition] Notation for the matrix-valued functions (A(x), B(x), …) is introduced without an immediate comparison table to the discrete MPO tensors; a short table relating the continuum objects to their lattice precursors would aid readability.
Simulated Author's Rebuttal
We thank the referee for their careful reading of our manuscript and for the constructive major comment on the continuum limit. We address the point in detail below and have revised the manuscript to include additional explicit demonstrations as requested.
read point-by-point responses
-
Referee: [continuum limit section] The continuum limit construction (detailed in the section following the ansatz definition): the claim that the resulting cMPO is exactly free of lattice artifacts and maps cMPS to cMPS requires an explicit demonstration that the partial trace over auxiliary degrees of freedom commutes with the limit; without a specified scaling of bond dimension with lattice spacing a or a check against a known solvable limit (e.g., free scalar field), residual a-dependent terms in the entanglement entropy cannot be ruled out.
Authors: We thank the referee for this observation, which helps clarify the presentation. In the construction, the continuum limit a to 0 is taken with fixed bond dimension D, as is standard for cMPS. The auxiliary space is discrete and the partial trace is defined on the same auxiliary indices both before and after the limit; because the limit acts only on the physical-space operators while the auxiliary contraction remains a finite matrix product, the operations commute by construction. To make this fully explicit, the revised manuscript now includes a dedicated paragraph deriving the commutation explicitly for the translation-invariant case and adds a short comparison to the free scalar field, confirming that the resulting entanglement entropy exhibits a clean area law with no residual a-dependent contributions. We have not introduced a scaling of D with a, as that would alter the continuum nature of the ansatz. revision: yes
Circularity Check
No circularity: cMPO ansatz constructed explicitly from continuum limit of discrete MPOs
full rationale
The paper defines the continuous matrix product operator ansatz directly as the result of a stated continuum limit applied to known discrete matrix product operators, yielding closed-form expressions in a finite set of matrix-valued functions A(x), B(x), … with no lattice spacing remaining. The area-law preservation and cMPS-to-cMPS mapping are shown to follow from this limit commuting with the partial trace over auxiliary indices. No equation reduces a claimed prediction to a fitted parameter by construction, no uniqueness theorem is imported from self-citations to forbid alternatives, and the central claims rest on the explicit construction rather than on renaming or smuggling prior ansatzes. The derivation is therefore self-contained against the external benchmark of discrete MPO theory.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Existence of a suitable continuum limit of discrete matrix product operators that eliminates lattice dependence while preserving area-law entanglement.
Forward citations
Cited by 1 Pith paper
-
An operational continuum limit of quantum combs
A continuous process tensor is defined by embedding the discrete multi-partite Choi matrix of a quantum comb into bosonic Fock space, closing the gap between discrete and continuum descriptions of multi-time quantum p...
Reference graph
Works this paper leans on
-
[1]
J. I. Cirac, D. P´ erez-Garc ´ ıa, N. Schuch, and F. Ver- straete, Matrix product states and projected entangled pair states: Concepts, symmetries, theorems, Rev. Mod. Phys.93, 045003 (2021)
work page 2021
-
[2]
M. B. Hastings, An area law for one-dimensional quan- tum systems, Journal of Statistical Mechanics: Theory and Experiment2007, P08024 (2007). 5
work page 2007
-
[3]
Matrix Product State Representations
D. Perez-Garcia, F. Verstraete, M. M. Wolf, and J. I. Cirac, Matrix product state representations, arXiv preprint quant-ph/0608197 (2006)
work page internal anchor Pith review Pith/arXiv arXiv 2006
-
[4]
F. Verstraete, J. J. Garc ´ ıa-Ripoll, and J. I. Cirac, Ma- trix product density operators: Simulation of finite- temperature and dissipative systems, Phys. Rev. Lett. 93, 207204 (2004)
work page 2004
-
[5]
M. Zwolak and G. Vidal, Mixed-state dynamics in one- dimensional quantum lattice systems: A time-dependent superoperator renormalization algorithm, Phys. Rev. Lett.93, 207205 (2004)
work page 2004
-
[6]
G. D. l. Cuevas, N. Schuch, D. P´ erez-Garc ´ ıa, and J. Ig- nacio Cirac, Purifications of multipartite states: limita- tions and constructive methods, New Journal of Physics 15, 123021 (2013)
work page 2013
- [7]
- [8]
-
[9]
K. Kato, Exact renormalization group flow for matrix product density operators, arXiv preprint arXiv:2410.22696 (2024)
-
[10]
J. Garre-Rubio, L. Lootens, and A. Moln´ ar, Classifying phases protected by matrix product operator symmetries using matrix product states, Quantum7, 927 (2023)
work page 2023
-
[11]
L. Lootens, J. Fuchs, J. Haegeman, C. Schweigert, and F. Verstraete, Matrix product operator symmetries and intertwiners in string-nets with domain walls, SciPost Phys.10, 053 (2021)
work page 2021
-
[12]
N. Bultinck, M. Mari¨ en, D. J. Williamson, M. B. S ¸ahino˘ glu, J. Haegeman, and F. Verstraete, Anyons and matrix product operator algebras, Annals of physics378, 183 (2017)
work page 2017
- [13]
- [14]
-
[15]
M. P. Zaletel, R. S. K. Mong, C. Karrasch, J. E. Moore, and F. Pollmann, Time-evolving a matrix product state with long-ranged interactions, Phys. Rev. B91, 165112 (2015)
work page 2015
-
[16]
M. M. Wolf, F. Verstraete, M. B. Hastings, and J. I. Cirac, Area laws in quantum systems: Mutual infor- mation and correlations, Phys. Rev. Lett.100, 070502 (2008)
work page 2008
-
[17]
L. Piroli and J. I. Cirac, Quantum cellular automata, tensor networks, and area laws, Phys. Rev. Lett.125, 190402 (2020)
work page 2020
-
[18]
R. Orus, A practical introduction to tensor networks: Matrix product states and projected entangled pair states, Annals of Physics349, 117 (2014)
work page 2014
- [19]
-
[20]
F. G. Brandao and M. Horodecki, Exponential decay of correlations implies area law, Communications in math- ematical physics333, 761 (2015)
work page 2015
-
[21]
F. Verstraete and J. I. Cirac, Matrix product states rep- resent ground states faithfully, Phys. Rev. B73, 094423 (2006)
work page 2006
- [22]
-
[23]
Reversible quantum cellular automata
B. Schumacher and R. F. Werner, Reversible quan- tum cellular automata, arXiv preprint quant-ph/0405174 (2004)
work page internal anchor Pith review Pith/arXiv arXiv 2004
- [24]
-
[25]
Farrelly, A review of Quantum Cellular Automata, Quantum4, 368 (2020)
T. Farrelly, A review of Quantum Cellular Automata, Quantum4, 368 (2020)
work page 2020
-
[26]
J. I. Cirac, D. Perez-Garcia, N. Schuch, and F. Ver- straete, Matrix product unitaries: structure, symmetries, and topological invariants, Journal of Statistical Mechan- ics: Theory and Experiment2017, 083105 (2017)
work page 2017
-
[27]
M. B. S ¸ahino˘ glu, S. K. Shukla, F. Bi, and X. Chen, Matrix product representation of locality preserving uni- taries, Phys. Rev. B98, 245122 (2018)
work page 2018
-
[28]
F. Verstraete and J. I. Cirac, Continuous matrix product states for quantum fields, Phys. Rev. Lett.104, 190405 (2010)
work page 2010
-
[29]
J. Haegeman, J. I. Cirac, T. J. Osborne, and F. Ver- straete, Calculus of continuous matrix product states, Phys. Rev. B88, 085118 (2013)
work page 2013
-
[30]
T. J. Osborne, J. Eisert, and F. Verstraete, Holographic quantum states, Physical review letters105, 260401 (2010)
work page 2010
-
[31]
J. Rinc´ on, M. Ganahl, and G. Vidal, Lieb-liniger model with exponentially decaying interactions: A continuous matrix product state study, Phys. Rev. B92, 115107 (2015)
work page 2015
- [32]
-
[33]
Tilloy, Relativistic continuous matrix product states for quantum fields without cutoff, Phys
A. Tilloy, Relativistic continuous matrix product states for quantum fields without cutoff, Phys. Rev. D104, 096007 (2021)
work page 2021
-
[34]
A. Tilloy, A study of the quantum sinh-gordon model with relativistic continuous matrix product states, arXiv preprint arXiv:2209.05341 (2022)
-
[35]
B. Tuybens, J. De Nardis, J. Haegeman, and F. Ver- straete, Variational optimization of continuous matrix product states, Phys. Rev. Lett.128, 020501 (2022)
work page 2022
-
[36]
G. De las Cuevas, N. Schuch, D. Perez-Garcia, and J. I. Cirac, Continuum limits of matrix product states, Phys. Rev. B98, 174303 (2018)
work page 2018
-
[37]
M. Balanz´ o-Juand´ o and G. De las Cuevas, Generalized ansatz for continuous matrix product states, Phys. Rev. A101, 052312 (2020)
work page 2020
-
[38]
We note that in [47] a notion of cMPO was introduced to capture the continuous (imaginary) time evolution of a lattice system in terms of tensor networks. However, it is defined implicitly as a limit of discrete MPO tensors along the time direction, and our ansatz provides the explicit cMPO directly in the limit by viewing the spatial direction as the tim...
-
[39]
G. Styliaris, R. Trivedi, D. Perez-Garcia, and J. I. Cirac, Matrix-product unitaries: Beyond quantum cellular au- tomata, Quantum9, 1645 (2025). 6
work page 2025
-
[40]
Here we useL(x) rather thanR(x) as is common in the cMPS literature [28, 29], as we will associate L(x),R(x) with left- and right-multiplication maps re- spectively when defining cMPO
-
[41]
In practice, one may also need to impose additional smoothness conditions ong(such ask-differentiability) to ensure that the the cMPO action on a cMPS does not worsen the well-behavedness of the initial cMPS [29]
-
[42]
Here we pick the convention that Θ(0) = 0 (left- continuous) so that [Π(x),ψ †(x)] = 0 for allx∈ [−ℓ/2,ℓ/2]
-
[43]
C. Kruszynska and B. Kraus, Local entanglability and multipartite entanglement, Phys. Rev. A79, 052304 (2009)
work page 2009
- [44]
-
[45]
Q. Mortier, L. Devos, L. Burgelman, B. Vanhecke, N. Bultinck, F. Verstraete, J. Haegeman, and L. Van- derstraeten, Fermionic tensor network methods, SciPost Phys.18, 012 (2025)
work page 2025
-
[46]
A. Tilloy and J. I. Cirac, Continuous tensor network states for quantum fields, Phys. Rev. X9, 021040 (2019)
work page 2019
-
[47]
W. Tang, H.-H. Tu, and L. Wang, Continuous matrix product operator approach to finite temperature quan- tum states, Phys. Rev. Lett.125, 170604 (2020)
work page 2020
-
[48]
R. Gilmore,Lie Groups, Physics, and Geometry: An In- troduction for Physicists, Engineers and Chemists(Cam- bridge University Press, 2008). End matter Product relations.To show that the product of two cMPO is another cMPO with bond matrices given by (10), we rely on the fact that only these products con- tribute in the Dyson series expansion: Ω·rx(Ω) =r x...
work page 2008
-
[49]
Phase MPU Our starting point is the concept of locally maximally entangleable (LME) states first introduced in [43]. Definition 3(LME state [43]).A multipartite state |ΨLME⟩is LME if there exists isometries Vj :H′ j→HA,j⊗HB,j, wheredimH ′ j =d ′ j anddimH α,j=d j forα=A,B, withV † jVj =1 d′ j such that the state |˜Ψ LME⟩AB = ⨂ j Vj|ΨLME⟩(38) is maximally ...
-
[50]
The first step is to take unnormalized phase MPS state in Eq
Phase cMPU We are now ready to construct our first non-trivial family of cMPUs by following the construction of phase MPUs. The first step is to take unnormalized phase MPS state in Eq. (43) and consider its continuum limit, i.e., a phase cMPS. Definition 7(Phase cMPS).We say that|Φθ⟩is a phase cMPS if it is a cMPS|ψ[B,Q,L]⟩defined in Eq.(2)such that the ...
-
[51]
Displaced phase cMPU Using phase cMPU as a basis, we can now go beyond the diagonal cMPUs using the displacement unitaries. To do so, we first recall the following result from the discrete MPUs that is closely related to theLU-equivalence of (2,2)-LME states in Proposition 1. Proposition 4([39]).EveryN-qubit unitaryUadmit- ting an(2,2)-LME state compressi...
-
[52]
Finite-dimensional cMPU We say that a cMPU is finite-dimensional if it acts non-trivially only on a finite-dimensional subspace of the full bosonic Hilbert space. The first such example we had was from the number-controlled phase cMPU in Propo- sition 3. We might expect that more examples should be possible because it acts non-trivially only on finite- di...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.