On the set of reduced states of translation invariant, infinite quantum systems

The set of two-body reduced states of translation invariant, infinite quantum spin chains can be approximated from inside and outside using matrix product states and marginals of finite systems, respectively. These lead to hierarchies of algebraic approximations that become tight only in the limit of infinitely many auxiliary variables. We show that this is necessarily so for any algebraic ansatz by proving that the set of reduced states is not semialgebraic. We also provide evidence that additional elementary transcendental functions cannot lead to a finitary description.


Introduction
Quantum correlations become non-trivially constrained when considering parts of a larger system that obeys symmetry rules.This fact is exhibited most prominently by the monogamy of entanglement.More generally, it is manifestly inherent in all quantum marginal problems that concern the consistency of reduced density matrices with global, often symmetry-based, constraints.A central motivation for studying these problems is the fact that ground-state energies in quantum chemistry [CY00] and condensed matter physics [VC06] depend solely on the reduced density matrices under global symmetry constraints.Also, phase transitions and symmetry breaking are reflected in the geometry of the set of reduced density matrices [ZDV + 16].
In this paper, we consider the set of reduced two-body density matrices that arise in translationally invariant, infinite quantum spin chains.In these systems, the ground state energy density of a nearest-neighbor interaction Hamiltonian could in principle be obtained by optimizing over all admissible reduced states.However, even in the case of two-dimensional local Hilbert spaces, this set is not known exactly.There are inner approximations based on mean-field and matrix-product techniques [VC06] and outer approximations based on relaxations of constraints [KSDN22,HKR20,BP12].When unraveling these approximations, which are usually expressed directly as algorithms for bounding ground state energy densities, they all correspond to semialgebraic sets of reduced density matrices.That is, sets that can be described by a finite number of algebraic (in-)equalities.Our aim is to prove that the set of interest, i.e., the set of admissible two-body reduced density matrices, is itself not semialgebraic.In the last section, we show that the situation is most likely even worse (or more interesting) by providing evidence that additional elementary transcendental functions are still not sufficient to allow an explicit description by finite means.
Our work is inspired by an argument of Fannes, Nachtergaele, and Werner.In [FNW92] they show that the ground state of the antiferromagnetic spin 1/2 Heisenberg chain cannot be a C * -finitely correlated state: its energy is transcendental while the algebraic ansatz can only lead to algebraic energies.This argument holds for any parameter-free algebraic ansatz but fails as soon as transcendental parameters are allowed.For this reason, we invoke transcendental sets and functions rather than numbers.

(Semi-)algebraic preliminaries
We begin with reviewing some mathematical concepts and introducing terms and notations.
A function f : I → R on an interval I ⊆ R is called algebraic over a subfield F ⊆ R if there is a polynomial p ∈ F[y, x] and an interval J ⊆ R s.t.y = f (x) ⇔ p(y, x) = 0 holds for all (x, y) ∈ I × J.A function that is not piecewise algebraic over F is called transcendental over F. We emphasize the field F since these notions are often understood with F = Q in mind, but we will need F = R.
A set A ⊆ R n is called semialgebraic over a subring R ⊆ R if it can be defined by a finite number of polynomial equalities and inequalities with coefficients in R. A powerful tool for showing that a set is semialgebraic is the Tarski-Seidenberg theorem, which enables the use of quantifiers without leaving the semialgebraic world: Theorem 1 (Tarski-Seidenberg quantifier elimination (cf.[BPCR07])).Let R ⊆ R be a subring and {p j (v, x)} l j=1 a finite set of polynomial equalities and inequalities with coefficients in R and variables ) is a Boolean combination (using ∧, ∨, ¬) of the p j 's and then there exists a quantifier-free formula ψ(x) consisting of a Boolean combination of finitely many polynomial (in-)equalities with coefficients in R, s.t.
Moreover, there exists an effective algorithm that constructs ψ from Ψ.
In our case, R n will be a real vector space of complex Hermitian matrices, represented in some Hermitian operator basis.The choice of the basis does not matter as long as R = R is considered.
A function is called semialgebraic over R if its graph is a semialgebraic set.

Reduced density matrices
Consider a multipartite quantum system of spins aligned on a one-dimensional chain such that each spin is assigned a Hilbert space C d .Our interest lies in the reduced density matrices of neighboring spins that are admissible under the constraint of global translation symmetry.
Definition 1.Let Red d denote the subset of density operators on C d ⊗ C d s.t. for every N ∈ N there is a density operator ρ (N +1) on (C d ) ⊗(N +1) with the property that all its successive bipartite reduced states are equal to ρ, i.e., More informally, ρ ∈ Red d if and only if ρ can be extended in a translational invariant way to an infinite quantum spin chain.
Proposition 1. Red d is convex and compact.
Proof.Convexity follows readily from the convexity of the set of density matrices and the linearity of the partial trace.Moreover, boundedness is inherited from the set of all density operators.In order to see closedness, we write ρ on the set of (N + 1)-partite density operators.Since f is continuous, the set is compact as the image of a compact set under a continuous map.Finally, since R N +1 is the set of bipartite density operators that can be extended to an (N + 1)-partite state with equal successive bipartite marginals, we see that Given a two-body Hamiltonian described by a Hermitian operator h acting The corresponding ground state energy density in the thermodynamic limit (N → ∞) is then where the infimum is taken over all N -partite density operators ρ (N ) .
The following is well known [VC06], but we provide proof for a more selfconsistent and coherent presentation.
Proposition 2. For every two-body Hamiltonian h the corresponding ground state energy density in the thermodynamic limit is given by Proof.First note that compactness of Red d (Prop.1)guarantees the existence of a minimizer.The r.h.s. of Eq.( 6) can be seen to be an upper bound on ϵ by restricting the infimum in Eq.( 5) to density operators with equal bipartite marginals.This leads to In order to obtain an inequality in the other direction, we add an additional interaction between the N 'th and the first particle so that the resulting overall Hamiltonian becomes cyclic on a ring of size N .Since the change in energy is at most ∥h∥, we can bound Here, we have identified the N + 1'st with the first site in Eq.( 7) and used for Eq.( 8) that the energy of a translational invariant Hamiltonian with periodic boundary condition is minimized by a translational invariant state.That is, in particular, by a state with equal bipartite marginals.□ Theorem 2. For d ≥ 2, the set Red d is not semialgebraic over R.
Proof.We will show this by proof of contradiction, starting from the assumption that Red d is semialgebraic.On an interval I ⊆ R consider a mapping I ∋ γ → h(γ) into the set of Hermitian matrices on C d ⊗ C d with the property that each matrix entry of h(γ) is a polynomial function of γ.We may, for instance, consider a simple affine relation of the form h(γ) = h 0 + γh 1 .Regarding h(γ) as a parameter-dependent interaction Hamiltonian, the graph of the corresponding ground state energy density ϵ(γ) on the interval I (a semialgebraic set) can be expressed as Under the assumptions that h is polynomial and Red d semialgebraic over R, the Tarski-Seidenberg theorem guarantees that the set defined in Eq.( 9) is semialgebraic since all equations are polynomial and all quantifiers run over semialgebraic sets.The fact that the quantifiers in Eq.( 9) are not at the beginning is not an obstacle since every first-order formula can be brought to prenex normal form, which is used in Thm.1.Consequently, I ∋ γ → ϵ(γ) is a semialgebraic function and therefore piecewise algebraic over R. Hence, the proof is completed by any example whose ground state energy density is transcendental over R (despite the fact that γ → h(γ) is polynomial).The next section will provide such an example for the case d = 2.This also covers the case of higher dimensions by simply embedding smaller Hilbert spaces into larger ones.□

Transcendental ground state energy density
As a special example for the completion of the proof of Thm.2, we consider the anisotropic XY-model, which is specified by the two-body Hamiltonian Theorem 3. In the thermodynamic limit, the ground state energy density ϵ(γ) of the anisotropic XY-model is an analytic function that is transcendental over

R.
Proof.The ground state energy density in the thermodynamic limit was derived in [LSM61] and shown to be where E(z) is the complete elliptic integral of the second kind.Suppose that ϵ is algebraic over R in some neighborhood.Since the set of algebraic functions is closed under composition, the same would be true for E. The function E is known to be an analytic solution to the hypergeometric differential equation for (a, b, c) = (−1/2, 1/2, 1) (see 17.3.10and 15.5.1 in [AS64]).Due to the analyticity of E, we can extend the view to the complex plane such that a supposed polynomial relation over R would imply one over C.
Consequently, E and therefore also ϵ are transcendental functions over R. □

Extensions and implications
The fact that Red d is not semialgebraic over R implies that no algebraic ansatz (with or without transcendental parameters) can describe the set exactly without invoking a limit of infinitely many variables.The best-known example is the set of matrix product states with increasing bond dimension.This raises the question of whether supplementing algebraic methods with transcendental functions could enable a finitary, explicit description of Red d .In [FFS23], for instance, entropy constraints, which are known to be transcendental [BW23,BM23], are used in addition.
Evidence against such finitary descriptions we can see from two sides: first, functions related to elliptic integrals (as in Eq.( 11)) can often be shown to be non-elementary [Kho14,Kas80] in the sense that they can not be expressed in terms of algebraic and elementary transcendental functions-neither explicitly nor implicitly.At the moment, however, we do not know exactly what the implications are for Red d .We will therefore pursue a second line of thought more closely.
We want to allow the additional use of the exponential function exp (and thereby implicitly also of ln).To this end, denote by R exp := (R, +, •, exp , < , 0, 1) the ordered field of real numbers with exponentiation, and call a subset A ⊆ R n definable in the first-order language of R exp if a ∈ A ⇔ Φ(a) for some first-order formula Φ.That is, there is a k ∈ N and a Boolean combination ϕ of finitely many (in-)equalities of polynomial 1 and exponential functions in n + k real variables, s.t.
where each Q i ∈ {∃, ∀} is a quantifier.Unlike for R there is no quantifier elimination for R exp [vd84].However, in [MW96] Macintyre and Wilkie proved that the first-order theory of R exp is decidable conditioned on the validity of the following: Conjecture 1 (Schanuel's conjecture (cf.Chap.21 in [MR14])).If z 1 , . . ., z n ∈ C are linearly independent over Q, then {z 1 , . . ., z n , exp z 1 , . . ., exp z n } contains at least n algebraically independent numbers. 2 Schanuel's conjecture can be considered the central conjecture of transcendental number theory.The best-known proven special case is the one where all z i 's are algebraic numbers.Then this becomes the content of the Lindemann-Weierstrass theorem.
Returning to the set of reduced density matrices, we obtain the following: Theorem 4. Assuming Schanuel's conjecture, there exists a d ∈ N such that the set Red d is not definable in the first-order language of R exp .
1 with algebraic coefficients 2 A set of numbers is called algebraically independent if there is no non-zero multivariate polynomial with coefficients in Q that has these numbers as roots.is undecidable.Assuming that Red d was definable in the first-order language of R exp , then Eq.( 14) would be expressible as a sentence in the first-order language of R exp .So according to the result of Macintyre and Wilkie [MW96], Eq.( 14) would be decidable if we assume Schanuel's conjecture.□ As shown in [BS04], Schanuel's conjecture in Thm.4 can be replaced by the so-called transfer conjecture, which also implies decidability of the first order theory of R exp .Moreover, there are closely related results that do not rely on any unproven conjecture.For instance, it was proven in [MW12] that sentences in the first-order theory of R trans , for a transcendental function trans ∈ {exp, ln, arctan, . ..}, are decidable if only the outermost quantified variable occurs in the transcendental function.

Proof.
By the results of [BCLPG20, CPGW22] there is a d ∈ N and a family of interaction Hamiltonians h on C d ⊗ C d with matrix entries in Q[ √ 2] (and thus expressible in R exp ) such that the ground state energy density problem ∀ρ ∈ Red d : tr [hρ] > 0 (14)