The set of separable states has no finite semidefinite representation except in dimension $3\times 2$

TL;DR

This paper proves that the set of separable quantum states cannot be described by finite semidefinite programming representations in most dimensions, providing a new counterexample to a longstanding conjecture.

Contribution

It demonstrates that for larger dimensions, the set of separable states lacks finite semidefinite representations, offering an elementary proof and a new counterexample to the Helton-Nie conjecture.

Findings

Sep(n,m) has no finite semidefinite representation for most (n,m)

Provides an elementary proof contrasting Scheiderer's approach

Offers a new counterexample to the Helton-Nie conjecture

Abstract

Given integers n $\geq$ m, let Sep(n,m) be the set of separable states on the Hilbert space $\mathbb{C}^n \otimes \mathbb{C}^m$. It is well-known that for (n,m)=(3,2) the set of separable states has a simple description using semidefinite programming: it is given by the set of states that have a positive partial transpose. In this paper we show that for larger values of n and m the set Sep(n,m) has no semidefinite programming description of finite size. As Sep(n,m) is a semialgebraic set this provides a new counterexample to the Helton-Nie conjecture, which was recently disproved by Scheiderer in a breakthrough result. Compared to Scheiderer's approach, our proof is elementary and relies only on basic results about semialgebraic sets and functions.

Equations

$$\text{Sep}(n,m)=\mathrm{conv}\left\{x{x}^{†}\otimes y{y}^{†}:x\in {\mathbb{C}}^n,|x|=1,y\in {\mathbb{C}}^m,|y|=1\right\}.$$

$$\text{PPT}(n,m)=\Bigl{\{}\rho\in\mathbf{H}^{nm}:\rho\geq 0,\;\;(I\otimes\mathsf{T})(\rho)\geq 0,\text{ and }\operatorname{Tr}[\rho]=1\Bigr{\}}.$$

$$C=\pi (S)$$

$$S=\{w\in {\mathbb{R}}^D:M_0+w_1{M}_1+\cdots +w_D{M}_D\ge 0\}$$

$$\text{Sep}(n,m)=\left\{\rho \in {\mathbf{H}}^{nm}:\mathrm{Tr}[\rho ]=1\text{ and }(I\otimes \Phi )(\rho )\ge 0\forall \Phi :{\text{M}}_m\to {\text{M}}_n\text{ positive}\right\}.$$

$$\text{Sep}\subseteq \cdots \subseteq {\text{DPS}}_k\subseteq {\text{DPS}}_{k-1}\subseteq \cdots \subseteq {\text{DPS}}_1=\text{PPT}.$$

$$\text{Sep}(n,m)\simeq \left\{\rho \in \text{Sep}(N,m):({\mathrm{Tr}}_2\rho {)}_{ii}=0\forall i=n+1,\dots ,N\right\}$$

$${\mathcal{C}}_A=\mathrm{cl}\mathrm{conv}\left\{{\mathbf{m}}_A(z):z\in {\mathbb{C}}^n\right\}$$

$$\mathcal{S}\mathcal{E}\mathcal{P}(n,m)$$

$$=\mathrm{conv}\left\{{\left[x^{\alpha }{\bar{x}}^{\beta }{y}^{\gamma }{\bar{y}}^{\delta }\right]}_{|\alpha |=|\beta |=|\gamma |=|\delta |=1}:(x,y)\in {\mathbb{C}}^n\times {\mathbb{C}}^m\right\}$$

$$A=\left\{(\alpha ,\beta ,\gamma ,\delta )\in ({\mathbb{N}}^n\times {\mathbb{N}}^n)\times ({\mathbb{N}}^m\times {\mathbb{N}}^m):|\alpha |=|\beta |=|\gamma |=|\delta |=1\right\}.$$

$$\mathrm{cl}\mathrm{conv}\left\{{\mathbf{m}}_A(x):x\in {\mathbb{R}}^n\right\}$$

$$p(z)=\sum _{(u,v)\in A}{p}_{uv}{z}^u{\bar{z}}^v(A\subset {\mathbb{N}}^n\times {\mathbb{N}}^n)$$

$$p=\sum _k{q}_k^2$$

$$\mathcal{S}\mathcal{E}\mathcal{P}(n,m)=\mathrm{conv}\left\{x{x}^{†}\otimes y{y}^{†}:x\in {\mathbb{C}}^n,y\in {\mathbb{C}}^m\right\}.$$

$$\mathcal{PP}\mathcal{T}(n,m)=\left\{\rho \in {\mathbf{H}}^{nm}:\rho \ge 0\text{ and }(I\otimes \mathsf{T})(\rho )\ge 0\right\}$$

$$p(x,y)=\sum _{\begin{array}{c}1\le i,j\le n\\ 1\le k,l\le m\end{array}}{p}_{ijkl}{x}_i{\bar{x}}_j{y}_k{\bar{y}}_l(x\in {\mathbb{C}}^n,y\in {\mathbb{C}}^m)$$

$$q_t(x,y)=\underset{g_t}{\underbrace{x^{\mathsf{T}}{M}_{t}y}}+\underset{\overline{g_t}}{\underbrace{{\bar{x}}^{\mathsf{T}}\overline{M_t}\bar{y}}}+\underset{h_t}{\underbrace{x^{\mathsf{T}}{N}_t\bar{y}}}+\underset{\overline{h_t}}{\underbrace{{\bar{x}}^{\mathsf{T}}{N}_{t}y}}$$

$$q_t^2=g_t^2+2|g_t{|}^2+2g_t{h}_t+2g_t\overline{h_t}+{\overline{g_t}}^2+2\overline{g_t}{h}_t+2\overline{g_t}\overline{h_t}+h_t^2+2|h_t{|}^2+{\overline{h_t}}^2.$$

$$p=\sum _k|x^{\mathsf{T}}{M}_{t}y{|}^2+|x^{\mathsf{T}}{N}_t\bar{y}{|}^2.$$

$$\begin{array}[]{ccc}\mathcal{SEP}&\subset&\mathcal{PPT}\\ \text{(duality)}\Bigg{\updownarrow}&&\hskip 42.67912pt\Bigg{\updownarrow}\text{(duality)}\\ \begin{array}[]{c}\text{Nonnegative}\\ \text{Hermitian polynomials \eqref{eq:pdual}}\end{array}&\supset&\begin{array}[]{c}\text{Sum-of-squares}\\ \text{Hermitian polynomials \eqref{eq:pdual}}\end{array}\end{array}$$

$$\displaystyle\Bigl{\{}$$

$$\displaystyle\qquad\lambda_{1},\ldots,\lambda_{n+1}\geq 0,s_{1},\ldots,s_{n+1}\in S,\;x=\sum_{i=1}^{n+1}\lambda_{i}s_{i}\text{ and }\sum_{i=1}^{n+1}\lambda_{i}=1\Bigr{\}}.$$

$$\left\{(\rho ,x,y)\in {\mathbf{H}}^{nm}\times {\mathbb{C}}^n\times {\mathbb{C}}^m\text{ s.t. }\rho =x{x}^{†}\otimes y{y}^{†}\text{ and }|x{|}^2=|y{|}^2=1\right\}$$

$${\mathbf{m}}_A(z)={\left[z^u{\bar{z}}^v\right]}_{(u,v)\in A}.$$

$${\mathcal{C}}_A=\mathrm{cl}\mathrm{conv}\left\{{\mathbf{m}}_A(z):z\in {\mathbb{C}}^n\right\}$$

$$S=\left\{w\in {\mathbb{R}}^N:M(w):=M_0+M_1{w}_1+\cdots +M_N{w}_N\ge 0\right\}.$$

$$⟨{\pi }^{\ast }(\tilde{p}),w⟩+c\ge 0\forall w\in {\mathbb{R}}^N\text{ s.t. }M(w)\ge 0.$$

$$⟨{\pi }^{\ast }(\tilde{p}),w⟩+c=⟨B,M(w)⟩+b\forall w\in {\mathbb{R}}^N.$$

$$⟨\tilde{p},\pi (w(z))⟩+c=⟨B,M(w(z))⟩+b=⟨B,F(z)F(z{)}^{†}⟩+b.$$

$$p(z)+c=\mathrm{Tr}\left[D{D}^{†}F(z)F(z{)}^{†}\right]+b=\sum _{ij=1}^k|D^{†}F(z){|}_{ij}^2+b=\sum _{ij=1}^k\text{Re}[(D^{†}F(z){)}_{ij}{]}^2+\text{Im}[(D^{†}F(z){)}_{ij}{]}^2+b.$$

Full text

1 Introduction

Entanglement is a fundamental aspect of quantum mechanics. The set of separable states (i.e., nonentangled states) on the Hilbert space $\mathbb{C}^{n}\otimes\mathbb{C}^{m}$ is defined as:

[TABLE]

Here $x^{\dagger}=\bar{x}^{\mathsf{T}}$ indicates conjugate transpose, $|x|^{2}=x^{\dagger}x=\sum_{i=1}^{n}|x_{i}|^{2}$ and $\operatorname{{\bf conv}}$ denotes the convex hull. The set $\text{Sep}(n,m)$ lives in the space $\mathbf{H}^{nm}$ of Hermitian matrices of size $nm\times nm$, and it is full-dimensional in the subspace of matrices of trace equal to one.

A fundamental computational task in quantum information is to decide membership in the convex set $\text{Sep}(n,m)$. One of the first tests designed to check whether a state $\rho\in\mathbf{H}^{nm}$ is separable is the Peres-Horodecki criterion [Per96, HHH96] (also known as the Positive Partial Transpose (PPT) criterion). It is based on the observation that for any $\rho\in\text{Sep}(n,m)$, $(I\otimes\mathsf{T})(\rho)$ is positive semidefinite where $I$ is the identity map, and $\mathsf{T}$ the transpose map. Indeed one can easily verify that if $\rho=xx^{\dagger}\otimes yy^{\dagger}$ then $(I\otimes\mathsf{T})(\rho)=xx^{\dagger}\otimes(yy^{\dagger})^{\mathsf{T}}=xx^{\dagger}\otimes\bar{y}\bar{y}^{\dagger}\geq 0$. In other words we have the inclusion $\text{Sep}(n,m)\subseteq\text{PPT}(n,m)$ where

[TABLE]

It is known, from earlier work of Woronowicz [Wor76], that we have equality $\text{Sep}(n,m)=\text{PPT}(n,m)$ if, and only if $n+m\leq 5$. Thus, the smallest cases where $\text{Sep}(n,m)\neq\text{PPT}(n,m)$ are $(n,m)=(4,2)$ and $(n,m)=(3,3)$.

Semidefinite programming

The description of the set $\text{PPT}(n,m)$ in Equation (1) allows us to decide membership, and optimize linear functions on $\text{PPT}(n,m)$, via semidefinite programming. Semidefinite programming is a fundamental tool in optimization that has played a crucial role in recent developments in quantum information theory. We say that a convex set $C$ has a semidefinite representation (also called a semidefinite lift) of size $r$ if it can be expressed as

[TABLE]

where $\pi:\mathbb{R}^{D}\rightarrow\mathbb{R}^{d}$ is a linear map and $S\subset\mathbb{R}^{D}$ is a convex set defined using a linear matrix inequality

[TABLE]

where $M_{0},\ldots,M_{D}$ are Hermitian matrices of size $r\times r$. A set $S$ of the form (3) is known as a spectrahedron. In this paper we are most interested in when a semidefinite representation of finite size exists, and call this simply a semidefinite representation throughout. If a convex set $C$ admits a semidefinite representation, then optimizing a linear function on $C$ can be cast as a semidefinite program. Equation (1) gives a semidefinite representation of $\text{PPT}(n,m)$.

Horodecki’s criterion

The set of separable states has the following well-known description due to the Horodeckis [HHH96]:

[TABLE]

Here $\text{M}_{k}=\mathbb{C}^{k\times k}$ and a $\mathbb{C}$-linear map $\Phi:\text{M}_{m}\rightarrow\text{M}_{n}$ is positive if it is Hermitian preserving and if $\Phi(X)\geq 0$ for all $X\geq 0$. When $n=m$, the relaxation $\text{PPT}(n,m)$ corresponds to having only the identity and transpose maps in Equation (4), which are both positive. A recent result of Skowronek [Sko16] shows that when $n=m=3$, there is no finite family of positive maps $\Phi_{1},\ldots,\Phi_{k}:\text{M}_{3}\rightarrow\text{M}_{3}$ such that $\text{Sep}(3,3)=\left\{\rho\in\mathbf{H}^{9}:(I\otimes\Phi_{i})(\rho)\geq 0\;\forall i=1,\ldots,k\right\}$. Note that the right-hand side of the previous equation is a specific semidefinite representable set. Thus Skowronek’s result rules out certain specific semidefinite representations for $\text{Sep}(3,3)$.111The result of Skrownek is in fact more general than this, and rules any formulation of the form $\left\{\rho:(I\otimes\Phi_{i})((I\otimes B)\rho(I\otimes B^{\dagger}))\geq 0\;\forall i=1,\ldots,k,\;\forall B\in\text{M}_{3}\right\}$ where $\Phi_{1},\ldots,\Phi_{k}$ is a finite set of positive maps.

DPS hierarchy

In [DPS04], Doherty, Parrilo and Spedalieri proposed a complete hierarchy of approximations to the set of separable states based on semidefinite programming. The first level of the hierarchy coincides with the PPT test, and subsequent levels form tighter and tighter convex relaxations of the set of separable states. If we denote the convex relaxation at level $k$ by $\text{DPS}_{k}(n,m)$ we have (dropping the $(n,m)$):

[TABLE]

The key property of the DPS hierarchy is that each set $\text{DPS}_{k}$ has a semidefinite representation of size $\min(n,m)^{O(k)}$. The hierarchy is known to be complete, meaning that if $\rho\notin\text{Sep}$, then there exists a finite $k$ such that $\rho\notin\text{DPS}_{k}$. The integer $k$ however depends on the state $\rho$ and it is known that, unless $n+m\leq 5$, there is no finite $k$ such that $\text{DPS}_{k}(n,m)=\text{Sep}(n,m)$ [DPS04, Section VIII.B].

Contributions

The main result of this paper is

Theorem 1**.**

If $\text{Sep}(n,m)\neq\text{PPT}(n,m)$ then $\text{Sep}(n,m)$ has no (finite) semidefinite representation. In other words, $\text{Sep}(n,m)$ has no semidefinite representation when $n+m>5$.

Remark 1*.*

Some remarks concerning Theorem 1:

• •

Note that Theorem 1 contains as a special case the fact that whenever $\text{Sep}(n,m)\neq\text{PPT}(n,m)$, then there is no representation of $\text{Sep}(n,m)$ as $\{\rho\in\mathbf{H}^{nm}:\operatorname{Tr}[\rho]=1\text{ and }(I\otimes\Phi_{i})(\rho)\geq 0,\;\forall i=1,\ldots,k\}$ where $\Phi_{1},\ldots,\Phi_{k}:\text{M}_{m}\rightarrow\text{M}_{n}$ is a finite family of positive maps. This is because the latter set is semidefinite representable (in fact it is a spectrahedron). Let us mention that if one is only interested in approximating the set $\text{Sep}(n,m)$, Aubrun and Szarek [AS17b] gave a lower bound on the number $k$ of positive maps needed.

• •

Our result also includes as a special case the fact that there is no finite $k$ such that $\text{DPS}_{k}(n,m)=\text{Sep}(n,m)$, when $n+m>5$. We note that our result is a strict generalization of this fact. Indeed, the failure of the DPS hierarchy to converge in a finite number of levels does not preclude by itself the existence of another semidefinite program that represents $\text{Sep}(n,m)$ exactly. There are well-known examples of convex sets where the sum-of-squares hierarchy (of which DPS can be seen as a particular instance) is never exact and yet a finite semidefinite representation does exist, see e.g., [NPS10, Example 3.7].

• •

Observe that if $\text{Sep}(n,m)$ has no semidefinite representation, then the same is true for $\text{Sep}(N,m)$ for $N\geq n$. This is because $\text{Sep}(n,m)$ can be realized as a linear section of $\text{Sep}(N,m)$ as follows:

[TABLE]

where $\operatorname{Tr}_{2}\rho$ is the result of tracing out the second subsystem from $\rho$. Indeed, setting $(\operatorname{Tr}_{2}\rho)_{ii}=0$ implies that in any representation of $\rho$ as $\rho=\sum_{k}p_{k}x_{k}x_{k}^{\dagger}\otimes y_{k}y_{k}^{\dagger}$, the vectors $x_{k}$ must satisfy $(x_{k})_{i}=0$ for all $i=n+1,\ldots,N$, i.e., that $x_{k}\in\mathbb{C}^{n}\times\{0\}^{N-n}\simeq\mathbb{C}^{n}$. To prove our theorem it thus suffices to prove that $\text{Sep}(3,3)$ and $\text{Sep}(4,2)$ have no semidefinite representations.

Helton-Nie conjecture

The question of finding semidefinite representations for convex sets has attracted a lot of attention in the optimization community [Nem06, GPT13]. Helton and Nie [HN09] gave sufficient conditions for a set to have a semidefinite representation, and conjectured that any convex semialgebraic set has a semidefinite representation. (A set is semialgebraic set if it can be described using a finite boolean combinations of polynomial equations and inequalities. One can verify that $\text{Sep}(n,m)$ is a semialgebraic set, see Section 2.3.) In his breakthrough paper, Scheiderer [Sch18] disproved this conjecture and exhibited convex semialgebraic sets that have no semidefinite representations.

Our proof of Theorem 1 is inspired from the arguments of Scheiderer. Compared to the paper of Scheiderer the present paper has two main contributions. First, the proof we give simplifies the arguments of Scheiderer and does not rely on any specialized results from algebraic geometry. We only use basic results from analysis (Taylor expansions), and some standard facts about semialgebraic sets and functions which are elementary to state. The proof should thus be accessible to readers in quantum information and optimization. The second contribution is the application of the method of proof for $\text{Sep}(n,m)$ which is defined in terms of complex numbers. This turns out to cause certain difficulties as certain standard facts about real polynomials are not true about Hermitian polynomials, particularly on the relation between homogeneous polynomials and their dehomogenizations (see Section 5 and Appendix A for more details).

Main technical result

Our main technical result, Theorem 2 below and of which Theorem 1 is a corollary, gives a general way to construct a convex set with no semidefinite representation from a nonnegative Hermitian polynomial that is not a sum of squares. We recall that a Hermitian polynomial $p(z)$ is a polynomial with complex coefficients in the indeterminates $(z,\bar{z})=(z_{1},\ldots,z_{n},\bar{z}_{1},\ldots,\bar{z}_{n})$ such that $p(z)\in\mathbb{R}$ for all $z\in\mathbb{C}^{n}$. A Hermitian polynomial is a sum of squares if it can be written as a sum of squares of Hermitian polynomials. (More details about Hermitian polynomials are given in Section 2.) For the statement of the theorem, we use the monomial notation $z^{u}=\prod_{i=1}^{n}z_{i}^{u_{i}}$ for $u\in\mathbb{N}^{n}$.

Theorem 2** (General theorem).**

Let $p(z)=\sum_{(u,v)\in A}p_{uv}z^{u}\bar{z}^{v}$ be a Hermitian polynomial supported on $A\subset\mathbb{N}^{n}\times\mathbb{N}^{n}$, and assume that $p$ is nonnegative on $\mathbb{C}^{n}$ but not a sum of squares. Assume furthermore that $A$ is downward closed, i.e., if $(u,v)\in A$ then all $(u^{\prime},v^{\prime})\in\mathbb{N}^{n}\times\mathbb{N}^{n}$ with $0\leq u^{\prime}\leq u$ and $0\leq v^{\prime}\leq v$ are in $A$. Define the monomial map $\mathbf{m}_{A}:\mathbb{C}^{n}\rightarrow\mathbb{C}^{|A|}$, $z\mapsto\left[z^{u}\bar{z}^{v}\right]_{(u,v)\in A}$ for $z\in\mathbb{C}^{n}$. Then the convex set

[TABLE]

is not semidefinite representable, where $\operatorname{{\bf cl}}$ denotes topological closure.

The set of separable states is of the form (5) for well-chosen set $A$. Indeed, dropping the normalization condition and letting $\mathcal{SEP}(n,m)$ be the convex cone of separable states, we have:

[TABLE]

where for $\zeta\in\mathbb{N}^{k}$ we let $|\zeta|=\sum_{i=1}^{k}\zeta_{i}$. This shows that $\mathcal{SEP}(n,m)=\mathcal{C}_{A}$ where

[TABLE]

The attentive reader will notice that this set $A$ is not downward closed, and so does not satisfy the condition of Theorem 2. As a matter of fact, to prove Theorem 1 we apply Theorem 2 with a dehomogenization of $A$ which satisfies the downward closed condition, and then homogenize back to get the desired convex cone. The details are explained in Section 5.

Overview of proof

We briefly sketch the main ideas for the proof of Theorem 2.

• •

We first show that if the set $\mathcal{C}_{A}$ has a semidefinite representation, then there exists a finite number of functions $f_{1},\ldots,f_{r}:\mathbb{R}^{2n}\simeq\mathbb{C}^{n}\rightarrow\mathbb{R}$ such that any nonnegative Hermitian polynomial supported on $A\cup\{(\mathbf{0},\mathbf{0})\}$ can be written as a sum of squares from $\operatorname{span}_{\mathbb{R}}(f_{1},\ldots,f_{r})$. This characterization of semidefinite representations via sums of squares is not new: it follows from the factorization theorem of Gouveia, Parrilo and Thomas [GPT13] and its sum-of-squares interpretation see e.g., [Faw16]. We note that a similar characterization is also used in Scheiderer’s paper, see [Sch18, Theorem 3.4].

• •

One of the main observations needed to prove Theorem 2 is to note that the functions $f_{1},\ldots,f_{r}$ can be chosen to be semialgebraic. (We recall the precise definition of semialgebraic functions in Section 2.) One key property of such functions that turns out to be particularly important is that they are smooth almost everywhere. Combining this property with a simple observation regarding smooth sum of squares decompositions of homogeneous polynomials allows us to prove Theorem 2 already in the special case where $p$ is a homogeneous polynomial. This allows us to prove that $\text{Sep}(n,m)$ is not semidefinite representable when $(n,m)=(5,3)$ or $(4,4)$. The complete proof of Theorem 2 which allows us to cover the cases $(n,m)=(4,2)$ and $(3,3)$ for separable states, requires an additional technical argument using Puiseux expansions for univariate continuous semialgebraic functions.

Real version of Theorem 2

We note that one can state an analogue of Theorem 2 dealing with real polynomials instead of Hermitian polynomials. The proof is similar, and we state it below just for convenience and for future reference.

Theorem 3** (Main theorem for real polynomials).**

Let $p(x)=\sum_{u}p_{u}x^{u}\in\mathbb{R}[x]$ where $A\subset\mathbb{N}^{n}$ finite, be a real polynomial that is nonnegative on $\mathbb{R}^{n}$ but not a sum of squares. Assume furthermore that $A$ is downward closed, i.e., if $u\in A$ then all $u^{\prime}\in\mathbb{N}^{n}$ with $0\leq u^{\prime}\leq u$ are in $A$. Define the monomial map $\mathbf{m}_{A}(x)=\left[x^{u}\right]_{u\in A}$ for $x\in\mathbb{R}^{n}$. Then the convex set

[TABLE]

is not semidefinite representable.

The theorem above can be used to recover the result of Scheiderer [Sch18, Corollary 4.25], that the cone $P_{n,2d}$ of nonnegative (real) forms in $n$ variables of degree $2d$ is not semidefinite representable when it is distinct from $\Sigma_{n,2d}$, the cone of sums of squares. Indeed, it suffices to take $p$ in Theorem 3 to be a dehomogenization of a nonnegative form that is not a sum of squares, and to use the well-known fact that a convex set has a semidefinite representation if and only if its dual has one.

Organization

The paper is organized as follows. In Section 2 we set some of the notations and present some background material on Hermitian polynomials, sums of squares, and semialgebraic sets and functions that are useful for the proof of the main theorem. In Section 3 we review the connection between the existence of semidefinite programming representations, and sums of squares. The proof of Theorem 2 is in Section 4 and the proof of Theorem 1 in Section 5.

2 Preliminaries

We recall in this section some results on Hermitian polynomials, the duality Sep/nonnegative polynomials and PPT/sums of squares and semialgebraic sets and functions.

2.1 Hermitian polynomials

For $z\in\mathbb{C}^{n}$, we denote the elementwise complex conjugate of $z$ by $\bar{z}=(\bar{z}_{1},\ldots,\bar{z}_{n})$. If $u\in\mathbb{N}^{n}$ we define the monomial $z^{u}=z_{1}^{u_{1}}\dots z_{n}^{u_{n}}$. A Hermitian polynomial $p(z)$ is a polynomial in $z$ and $\bar{z}$ of the form

[TABLE]

such that $p(z)\in\mathbb{R}$ for all $z\in\mathbb{C}^{n}$. This is equivalent to saying that $p_{uv}=\overline{p_{vu}}$ for all $u,v$. The support of $p$ is $\operatorname{supp}(p)=\{(u,v):p_{uv}\neq 0\}\subset\mathbb{N}^{n}\times\mathbb{N}^{n}$. The Hermitian polynomial $p$ is nonnegative if $p(z)\geq 0$ for all $z\in\mathbb{C}^{n}$. Further, we say that $p$ is a sum of squares if we can write

[TABLE]

for Hermitian polynomials $q_{k}$. If $p(z)$ is a Hermitian polynomial we will often consider the real polynomial $P(a,b)=p(a+ib)$ in $\mathbb{R}[a_{1},\ldots,a_{n},b_{1},\ldots,b_{n}]$. One can check that $p$ is a sum-of-squares if and only if $P$ is a sum-of-squares of real polynomials.

Remark 2* (Sums of squares for Hermitian polynomials).*

Another common definition of a Hermitian polynomial $p(z)$ being a sum-of-squares is that $p$ can be written as $p(z)=\sum_{k}|g_{k}(z)|^{2}$ where $g_{k}$ are (holomorphic) polynomials in $z$ only (and not in $\bar{z}$). Clearly if $p$ has such a representation then it is a sum-of-squares in the sense (8) since then $p=\sum_{k}\text{Re}[g_{k}]^{2}+\text{Im}[g_{k}]^{2}$ and $\text{Re}[g_{k}]$ and $\text{Im}[g_{k}]$ are both Hermitian polynomials. The converse however is not true. It is possible that a polynomial $p$ has a representation (8) and cannot be written as a sum of modulus squares of holomorphic polynomial mappings. See e.g., [DP09] for more on this distinction. In this paper we only work with the definition (8) of sums of squares.

2.2 Sep, PPT, nonnegative polynomials, and sums of squares

For convenience, we will work in this paper with the cone of separable states, where we drop the normalization condition:

[TABLE]

One can verify that $\text{Sep}(n,m)$ is the compact slice $\text{Sep}(n,m)=\mathcal{SEP}(n,m)\cap\{\rho:\operatorname{Tr}\rho=1\}$.222Indeed if $\rho=\sum_{i}p_{k}x_{k}x_{k}^{\dagger}\otimes y_{k}y_{k}^{\dagger}$ with $\operatorname{Tr}\rho=1$ and $p_{k}\geq 0$, then by redefining $p_{k}\leftarrow p_{k}|x_{k}|^{2}|y_{k}|^{2}$ we can assume without loss of generality that $|x_{k}|=|y_{k}|=1$. Taking the trace on both sides of $\rho=\sum_{i}p_{k}x_{k}x_{k}^{\dagger}\otimes y_{k}y_{k}^{\dagger}$ tells us that $1=\sum_{k}p_{k}$ since $\operatorname{Tr}\rho=1$, i.e., $\rho\in\text{Sep}(n,m)$. Let also $\mathcal{PPT}$ be the cone of states that have positive partial transpose, i.e.,

[TABLE]

so that $\text{PPT}(n,m)=\mathcal{PPT}(n,m)\cap\left\{\rho:\operatorname{Tr}\rho=1\right\}$.

Dual of Sep

For any integer $k$, let $\text{M}_{k}=\mathbb{C}^{k\times k}$. A $\mathbb{C}$-linear map $\Phi:\text{M}_{n}\rightarrow\text{M}_{m}$ that is Hermitian preserving is positive if $\Phi(\rho)\geq 0$ whenever $\rho\geq 0$. Equivalently, $\Phi$ is positive if the degree-four Hermitian polynomial $p(x,y)=y^{\dagger}\Phi\left(xx^{\dagger}\right)y$ is nonnegative on $\mathbb{C}^{n+m}\simeq\mathbb{C}^{n}\times\mathbb{C}^{m}$. It is well-known that the dual of $\mathcal{SEP}(n,m)$ can be identified, via the Choi isomorphism, with the cone of positive maps $\text{M}_{n}\rightarrow\text{M}_{m}$ (see e.g., [AS17a, Table 2.2]). Equivalently, the dual of $\mathcal{SEP}(n,m)$ can be identified with nonnegative degree-four Hermitian polynomials of the form

[TABLE]

where $p_{ijkl}=\Phi(E_{ij})_{lk}$. Polynomials of the form (9) have a biquadratic structure: they are quadratic independently in each block of variables $x$ and $y$. The duality between $\mathcal{SEP}$ and nonnegative polynomials of the form (9) is in fact immediate from the definition of $\mathcal{SEP}$.

Dual of PPT

Using the identification above, it turns out that the dual of $\mathcal{PPT}(n,m)$ corresponds to polynomials $p(x,y)$ that are sums of squares. Indeed, it is well-known (see again [AS17a, Table 2.2]) that the dual cone of $\mathcal{PPT}(n,m)$ can be identified, via the Choi isomorphism, with the cone of maps $\Phi:\text{M}_{n}\rightarrow\text{M}_{m}$ that are decomposable, i.e., that can be written $\Phi=S_{1}+S_{2}\circ\mathsf{T}$ where $S_{1}$ and $S_{2}$ are two completely positive maps, and $\mathsf{T}$ is the transpose map. Recall that a map $S:\text{M}_{n}\rightarrow\text{M}_{m}$ is completely positive if there exist matrices $V_{t}$ such that $S(X)=\sum_{t}V_{t}^{*}XV_{t}$. One can verify that a map $\Phi$ is decomposable if and only if, the associated Hermitian polynomial (9) is a sum of squares. We did not find any reference for this equivalence, so we include a proof here. (The proofs we found in the literature are only for the direction $\Rightarrow$ in Proposition 1. The proof of Proposition 1 is a special case of a more general result in [FF19], joint with Kun Fang, where it is shown that the dual of $\text{DPS}_{k}$ can be identified with a sum-of-squares condition of degree $k$.)

Proposition 1**.**

A map $\Phi:\text{M}_{n}\rightarrow\text{M}_{m}$ is decomposable if, and only if, the Hermitian polynomial $p(x,y)=y^{\dagger}\Phi\left(xx^{\dagger}\right)y$ is a sum of squares.

Proof.

If $S(\rho)=\sum_{t}V_{t}\rho V_{t}^{\dagger}$ is a completely positive map then $y^{\dagger}S(xx^{\dagger})y=\sum_{t}y^{\dagger}V_{t}xx^{\dagger}V_{t}^{\dagger}y=\sum_{t}|\bar{y}^{\mathsf{T}}V_{t}x|^{2}$ is a sum-of-squares. Also for the transpose map $\mathsf{T}$, we have $y^{\dagger}(S\circ\mathsf{T})(xx^{\dagger})y=y^{\dagger}S(\bar{x}\bar{x}^{\dagger})y=\sum_{t}|y^{\mathsf{T}}\bar{V_{t}}x|^{2}$ is also a sum-of-squares. It follows that if $\Phi$ is decomposable then $p(x,y)=y^{\dagger}\Phi\left(xx^{\dagger}\right)y$ is a sum-of-squares.

We now prove the converse. Assume $p(x,y)=y^{\dagger}\Phi(xx^{\dagger})y$ is a sum-of-squares, i.e., $p(x,y)=\sum_{t}q_{t}(x,y)^{2}$ for some Hermitian polynomials $q_{t}$. We need to show that $\Phi$ is decomposable. Since the coefficient of the monomial $x_{i}^{2}\bar{x_{i}}^{2}$ in $p$ is 0, we see that $q_{t}$ cannot have monomials $x_{i}^{2},\bar{x_{i}}^{2}$ or $x_{i}\bar{x_{i}}$. To be sure, let $\alpha_{t},\bar{\alpha_{t}},\beta_{t}$ be the coefficients in $q_{t}$ of these monomials (note that $\beta_{t}\in\mathbb{R}$ since $x_{i}\bar{x_{i}}$ is real). The coefficient of $x_{i}^{2}\bar{x_{i}}^{2}$ in $\sum_{t}q_{t}^{2}$ is $\sum_{t}2|\alpha_{t}|^{2}+\beta_{t}^{2}=0$ which implies that $\alpha_{t}=\beta_{t}=0$ for all $t$. Similarly, by looking at the coefficient of $x_{i}\bar{x_{i}}x_{j}\bar{x_{j}}$ in $p$, we see that the $q_{t}$ cannot have monomials of the form $x_{i}x_{j},\bar{x_{i}}\bar{x_{j}},x_{i}\bar{x_{j}}$ or $\bar{x_{i}}x_{j}$. The same of course is true for the $y$’s. Thus this means that $q_{t}$ must have the form

[TABLE]

where $M\in\mathbb{C}^{n\times n}$ and $N\in\mathbb{C}^{m\times m}$. Squaring $q_{t}$ we get

[TABLE]

When summing $\sum_{k}q_{t}^{2}$ we see that the only terms that can produce monomials of the form $x_{i}\bar{x_{j}}y_{k}\bar{y_{l}}$ (the monomials that appear in $p$) are the terms $2|g_{t}|^{2}$ and $2|h_{t}|^{2}$. The sum of all the other terms must thus be equal to 0. At the end we get (including the constant 2 in $M_{t}$ and $N_{t}$):

[TABLE]

From here it easily follows that $\Phi=S_{1}+S_{2}\circ\mathsf{T}$ where $S_{1}(\rho)=\sum_{t}N_{t}\rho N_{t}^{\dagger}$ and $S_{2}(\rho)=\sum_{t}\bar{M_{t}}\rho\bar{M_{t}}^{\dagger}$. ∎

The following diagram summarizes the discussion above.

[TABLE]

2.3 Semialgebraic sets and functions

Semialgebraic sets

A semialgebraic subset of $\mathbb{R}^{n}$ is a subset that can be defined by a finite boolean combination of polynomial equations ($P=0$) and inequalities ($P>0$) where $P\in\mathbb{R}[x_{1},\ldots,x_{n}]$. For example a set of the form $\{w\in\mathbb{R}^{D}:M_{0}+w_{1}M_{1}+\dots+w_{D}M_{D}\geq 0\}$ is semialgebraic since the condition that a matrix is positive semidefinite can be expressed by a finite number of polynomial inequalities. The set of separable states can be shown to be semialgebraic. One can prove this using the celebrated and powerful result of Tarski stating that the projection of a semialgebraic set is semialgebraic. A consequence of Tarski’s theorem is that the convex hull of a semialgebraic set $S\subset\mathbb{R}^{n}$ is semialgebraic. Indeed this is because we can write $\operatorname{{\bf conv}}(S)$ as the projection on the $x$ component of the following semialgebraic set

[TABLE]

(Note that, by Carathéodory theorem any element in $\operatorname{{\bf conv}}(S)$ is a convex combination of at most $n+1$ points in $S$.) To see why the set $\text{Sep}(n,m)$ is a semialgebraic subset of $\mathbf{H}^{nm}\simeq\mathbb{R}^{2(nm)^{2}}$ first note that the following set

[TABLE]

is a semalgebraic subset of $\mathbf{H}^{nm}\times\mathbb{C}^{n}\times\mathbb{C}^{m}\simeq\mathbb{R}^{2(nm)^{2}}\times\mathbb{R}^{2n}\times\mathbb{R}^{2m}$ since the equations can all be written as real polynomial equations in the real and imaginary components. By Tarski’s theorem it follows that the projection of (10) on the $\mathbf{H}^{nm}$ component, which is precisely the set of pure product states, is semialgebraic. Thus Sep is semialgebraic as the convex hull of a semialgebraic set.

Semialgebraic functions

A function $f:\mathbb{R}^{n}\rightarrow\mathbb{R}^{m}$ is called semialgebraic if its graph $\{(x,f(x)):x\in\mathbb{R}^{n}\}\subseteq\mathbb{R}^{n}\times\mathbb{R}^{m}$ is a semialgebraic set. Even though semialgebraic functions form a very broad class of functions, they are tame and possess nice regularity properties. Examples of semialgebraic functions are polynomials, rational functions, or power functions (with rational exponent). Functions that are not semialgebraic are e.g., $\exp(x)$, or the indicator function of the rationals in $\mathbb{R}$. We state two basic results about semialgebraic functions that will be crucial for us.

Theorem 4** (Almost everywhere smoothness of semialgebraic functions).**

Let $f:\mathbb{R}^{n}\rightarrow\mathbb{R}$ be a semialgebraic function. Then $f$ is smooth ($C^{\infty}$) everywhere except possibly on the zero set of a polynomial $P\in\mathbb{R}[x_{1},\ldots,x_{n}]\setminus\{0\}$.

Proof.

See e.g., [HP16, Theorem 1.7]. ∎

Clearly Theorem 4 is not true for general functions, cf. the indicator function of the rationals in $\mathbb{R}$. The second result that we will need concerns semialgebraic functions in one variable.

Theorem 5** (Puiseux expansion for one-dimensional semialgebraic functions).**

Assume $f:(0,\eta)\rightarrow\mathbb{R}$ where $\eta>0$, is a semialgebraic continuous function that is bounded. Then $f$ can be extended by continuity to the interval $[0,\eta)$. Additionally, there exists an integer $m$ such that the map $t\mapsto f(t^{m})$ is $C^{\infty}$ on $[0,\epsilon)$ for some $0<\epsilon<\eta$.

Proof.

For the first part, see [BPR06, Proposition 3.18]. For the second part, see [Cos05, page 10]. ∎

We note that the theorem above is not true for arbitrary functions. For example the function $x\mapsto\sin(1/x)$ is bounded and continuous on any interval $(0,\eta)$ but cannot be extended by continuity at [math]. We finally record the following result which will also be needed for our proof. It simply says that any linear map, restricted to a semialgebraic set always admits a semialgebraic inverse.

Theorem 6**.**

Let $S\subset\mathbb{R}^{N}$ be a semialgebraic set and let $\pi:\mathbb{R}^{N}\rightarrow\mathbb{R}^{n}$ be a linear map. Then there exists a semialgebraic function $F:\pi(S)\rightarrow S$ that satisfies $\pi(F(x))=x$ for all $x\in\pi(S)$.

Proof.

See [HP16, Lemma 1.5]. ∎

3 Semidefinite programming lifts

We are now ready to start the proof of Theorem 2 (and thus of Theorem 1 too). The first thing we need is a necessary condition for the existence of a semidefinite representation for a given convex set $C$. The condition we state in Theorem 7 below is very similar to [GPT13, Theorem 1]333The given condition is also sufficient, but we will only need necessity here. Recall that, for $A\subset\mathbb{N}^{n}\times\mathbb{N}^{n}$ we denote by $\mathbf{m}_{A}(z)$ the monomial map:

[TABLE]

We also denote by $\operatorname{{\bf cl}}S$ the (topological) closure of a set $S$.

Theorem 7**.**

Let $A\subset\mathbb{N}^{n}\times\mathbb{N}^{n}$. Assume that

[TABLE]

has a semidefinite representation of size $k$. Then there exists $2k^{2}+1$ semialgebraic functions $f_{j}:\mathbb{C}^{n}\rightarrow\mathbb{R}$ ($j=1,\ldots,2k^{2}+1$) such that any nonnegative Hermitian polynomial $p$ supported on $A\cup\{(\mathbf{0},\mathbf{0})\}$ is a sum-of-squares from $V=\operatorname{span}_{\mathbb{R}}(f_{1},\ldots,f_{2k^{2}+1})$, i.e., $p=\sum_{j}h_{j}^{2}$ for some $h_{j}\in V$. Furthermore, the magnitude of the coefficients expressing the $h_{j}$ in terms of the basis $(f_{1},\ldots,f_{2k^{2}+1})$ are all bounded by $\phi(\|p\|)$ where $\|p\|$ is the largest magnitude of the coefficients of $p$, and $\phi$ is some polynomial that only depends on the semidefinite representation of $\mathcal{C}_{A}$.

The main difference between the statement above and the one in [GPT13] (see also [Faw16, Theorem 5, Chapter 2]) is that here the functions $f_{1},\ldots,f_{2k^{2}+1}$ are semialgebraic. (We say that a function $f:\mathbb{C}^{n}\rightarrow\mathbb{R}$ is semialgebraic if the function $F:\mathbb{R}^{n}\times\mathbb{R}^{n}\rightarrow\mathbb{R}$ defined by $F(a,b)=f(a+ib)$ is semialgebraic.) This observation will be crucial to us. We note that a statement similar to the theorem above appears as Theorem 3.4 in [Sch18]. Instead of working with semialgebraic functions, Scheiderer works with polynomial functions on an algebraic variety $X$.

Proof.

Assume that $\mathcal{C}_{A}=\pi(S)$ where $S$ is a spectrahedron defined using a linear matrix inequality of size $k\times k$:

[TABLE]

We can assume without loss generality that $S$ has nonempty interior in $\mathbb{R}^{N}$. This in turn implies, using standard results about spectrahedra, that there exists $\widetilde{w}\in S$ such that $M(\widetilde{w})$ is positive definite (possibly after changing $M$), see e.g., [RG95, Section 2.4].

For any $z\in\mathbb{C}^{n}$ there exists $w(z)\in S$ such that $\pi(w(z))=\mathbf{m}_{A}(z)$. Since $M(w(z))\geq 0$ we can find $F(z)\in\mathbb{C}^{k\times k}$ such that $M(w(z))=F(z)F(z)^{\dagger}$. Furthermore, by Theorem 6, the function $z\in\mathbb{C}^{n}\mapsto F(z)\in\mathbb{C}^{k\times k}$ can be taken to be semialgebraic.

Let $p(z)$ be a Hermitian polynomial supported on $A\cup\{(\mathbf{0},\mathbf{0})\}$, i.e., $p(z)=\langle\widetilde{p},\mathbf{m}_{A}(z)\rangle+c$ for some $c\in\mathbb{R}$, where $\widetilde{p}$ denotes the coefficients of the polynomial $p(z)$ in the monomial basis. Since $p\geq 0$ we get $\langle\widetilde{p},\mathbf{m}_{A}(z)\rangle+c\geq 0$ for all $z\in\mathbb{C}^{n}$. This implies that $\langle\widetilde{p},\sigma\rangle+c\geq 0$ for all $\sigma\in\mathcal{C}_{A}$. We can lift this linear inequality to an inequality on the spectrahedron $S$, i.e., we have $\langle\widetilde{p},\pi(w)\rangle+c\geq 0$ for all $w\in S$, in other words

[TABLE]

By Farkas’ lemma/duality for SDPs, this means that there exists $B\geq 0$ and $b\geq 0$ such that

[TABLE]

Plugging $w=w(z)$ we get

[TABLE]

Since $\pi(w(z))=\mathbf{m}_{A}(z)$ and $\langle\widetilde{p},\mathbf{m}_{A}(z)\rangle=p(z)$ we get finally that $p(z)+c=\langle B,F(z)F(z)^{\dagger}\rangle+b$ for all $z\in\mathbb{C}^{n}$. Factorizing $B=DD^{\dagger}$ we get

[TABLE]

If we define the $2k^{2}+1$ semialgebraic functions to be the constant function and the $z\mapsto\text{Re}[F_{ij}(z)]$ and $z\mapsto\text{Im}[F_{ij}(z)]$, we get the desired claim.

For the last statement of the theorem, we show that the coefficients of $B$ (and thus of $D$) are bounded by a polynomial in the coefficients of $p$. To get this, we can simply plug the value $w=\widetilde{w}$ that makes $M(w)$ positive definite in (11). If we denote by $\lambda>0$ the smallest eigenvalue of $M(\widetilde{w})$ we get $\langle\pi^{*}(\widetilde{p}),\widetilde{w}\rangle+c=\langle B,M(\widetilde{w})\rangle+b\geq\lambda\operatorname{Tr}[B]+b\geq\lambda\|B\|+b$, thus $\max(\|B\|,b)\leq(\langle\pi^{*}(\widetilde{p}),\widetilde{w}\rangle+c)/\min(\lambda,1)\leq O(\max\{\|p\|,|c|\})$. ∎

4 Proof of Theorem 2

We are now ready to prove our main theorem, Theorem 2. We first recall a piece of notation that we will use throughout the proof: for any function $f:\mathbb{C}^{n}\rightarrow\mathbb{R}$ we associate the function of real variables $F:\mathbb{R}^{n}\times\mathbb{R}^{n}\rightarrow\mathbb{R}$, denoted by a capital letter, defined by $F(a,b)=f(a+ib)$. We will also sometimes think of a vector $z\in\mathbb{C}^{n}$ as $z\in\mathbb{R}^{2n}$ and write for instance $F(z)$.

Assume that $\mathcal{C}_{A}$ has an SDP representation. Then, from Theorem 7 there exist semialgebraic functions $f_{1},\ldots,f_{r}:\mathbb{C}^{n}\rightarrow\mathbb{R}$ such that the following is true:

[TABLE]

We will now prove that the functions $f_{i}$ can be taken to be smooth at the origin. (By this we mean that the associated functions $F_{i}:\mathbb{R}^{n}\times\mathbb{R}^{n}\rightarrow\mathbb{R}$ are smooth at $(0,0)$.) This will follow from our assumption that $A$ is downward closed. Since the $F_{i}$ are semialgebraic we know from Theorem 4 that each $F_{i}$ is smooth almost everywhere. Thus we can find a common point $z_{0}=(a_{0},b_{0})\in\mathbb{R}^{n}\times\mathbb{R}^{n}$ such that all functions $F_{i}$ are smooth at $z_{0}$. Now let $\widetilde{f_{i}}(z)=f_{i}(z-z_{0})$ for all $i\in\{1,\ldots,r\}$. We claim that these semialgebraic functions still satisfy the property ($\ast$ ‣ 4). Indeed if $q$ is a nonnegative Hermitian polynomial supported on $A\cup\{(\mathbf{0},\mathbf{0})\}$ then $q(z+z_{0})$ is nonnegative and is also supported on $A$, since $A$ is downward closed. It follows that $q(z+z_{0})$ is a sum-of-squares from $\operatorname{span}_{\mathbb{R}}(f_{1},\ldots,f_{r})$. But this implies that $q(z)$ is a sum-of-squares from $\operatorname{span}_{\mathbb{R}}(\widetilde{f}_{1},\ldots,\widetilde{f}_{r})$.

In the rest of the proof we will thus assume that the $F_{1},\ldots,F_{r}$ are smooth at the origin. If $p$ is homogeneous we are almost done by the following simple observation: if $P$ is a real homogeneous polynomial and $P=\sum_{j}H_{j}^{2}$ for some functions $H_{j}:\mathbb{R}^{m}\rightarrow\mathbb{R}$ that are smooth at the origin, then $P$ is a sum-of-squares of polynomials. This can be proved by a simple Taylor expansion; for example by applying the following proposition to the identity $t^{2k}P(z)=P(tz)=\sum_{j}H_{j}(tz)^{2}$ and observing that $\left.\frac{d^{k}}{dt^{k}}H_{j}(tz)\right|_{t=0}$ is a degree $k$ polynomial in $z\in\mathbb{R}^{2n}$ (it is the $k$’th term in the Taylor expansion of $H_{j}$ at $z$).

Proposition 2**.**

Assume that $g_{j}:[0,\eta)\rightarrow\mathbb{R}$ are smooth functions444Smoothness at [math] is smoothness on the right and that there exists $a\in\mathbb{R}$ such that $at^{2k}=\sum_{j}g_{j}(t)^{2}$ for all $t\in[0,\eta)$. Then $a=\sum_{j}\left(\frac{g_{j}^{(k)}(0)}{k!}\right)^{2}$.

Proof.

If we Taylor expand the $g_{j}$ at [math] we get $at^{2k}=\sum_{j}(g_{j}(0)+tg_{j}^{\prime}(0)+\dots+t^{k}g_{j}^{(k)}(0)/k!+o(t^{k}))^{2}$. By equating powers of $t$ we get that $g_{j}(0)=\dots=g_{j}^{(k-1)}(0)=0$ and that $a=\sum_{j}(g_{j}^{(k)}(0)/k!)^{2}$ as desired. ∎

The case where $p$ is not necessarily homogeneous requires an additional argument. The following argument is inspired from [Sch18, Proposition 4.18]. Let $2d=\deg p$, and for any $t\in\mathbb{R}$ consider the Hermitian polynomial $p_{t}(z)=t^{2d}p(z/t)$. This Hermitian polynomial is nonnegative and is also supported on $A$. Thus we know from property ($\ast$ ‣ 4) that there exist real coefficients $a_{j}(t)\in\mathbb{R}^{r}$ s.t.

[TABLE]

where we let $F(z)=(F_{1}(z),\ldots,F_{r}(z))$. The functions $a_{j}(t)$ are defined by a semialgebraic relation and so can be taken to be semialgebraic. As such the $a_{j}$ must be continuous on some $(0,\eta)$ for $\eta>0$. From the last part of the statement of Theorem 7 we also know that the $a_{j}$ must be bounded on $(0,\eta)$. Thus, by Theorem 5 we know that the $a_{j}$ can be extended by continuity to $[0,\eta)$ and that for large enough $m$, $a_{j}(t^{m})$ is smooth on $[0,\eta^{\prime})$ for some $0<\eta^{\prime}<\eta$. From (12) we get:

[TABLE]

But note that $P_{t^{m}}(t^{m}z)=t^{2dm}P(z)$. Thus

[TABLE]

If we let $g_{j}(t)=a_{j}(t^{m})^{\mathsf{T}}F(t^{m}z)$ we know that the $g_{j}$ are smooth on $[0,\eta^{\prime})$, since the $F$ are smooth at the origin. We can apply the observation of Proposition 2 to get that

[TABLE]

But, from the definition of $g_{j}$, $g_{j}^{(dm)}(0)$ is a polynomial (of degree $d$) in $z$. This contradicts the assumption that $p(z)$ is not a sum of squares of polynomials.

5 Proof of Theorem 1

The case $\text{Sep}(3,3)$: We prove that $\mathcal{SEP}(3,3)$ has no semidefinite representation. Define the Choi polynomial [Cho75] by

[TABLE]

It was shown in [Cho75] (see also [Cho80, Appendix B]) that $p(x,y)\geq 0$ for all $(x,y)\in\mathbb{C}^{3}\times\mathbb{C}^{3}$, and yet $p(x,y)$ is not a sum of squares. In fact the real polynomial $p(x,y)$ when $(x,y)\in\mathbb{R}^{3}\times\mathbb{R}^{3}$ is not a sum of squares. It follows, by a simple homogenization argument, that the Hermitian polynomial $\hat{p}(x_{1},x_{2},y_{1},y_{2})=p(x_{1},x_{2},1,y_{1},y_{2},1)$ is not a sum of squares. Note that the support of $\hat{p}$ satisfies

[TABLE]

Since $\hat{A}$ is downward closed it follows from Theorem 2 that

[TABLE]

does not have a semidefinite representation. To conclude that $\mathcal{SEP}(3,3)$ has no semidefinite representation, it remains to note that $\mathcal{C}_{\hat{A}}$ is a hyperplane section of $\mathcal{SEP}(3,3)$. Indeed, first recall that $\mathcal{SEP}(3,3)$ can be written as

[TABLE]

where

[TABLE]

It is easy to see that there is a one-to-one correspondence between $\hat{A}$ and $A$. In terms of the monomial map $\mathbf{m}$ this simply means that $\mathbf{m}_{A}$ is the homogenization of $\mathbf{m}_{\hat{A}}$. For example under an appropriate ordering of the monomials we have $\mathbf{m}_{A}(x_{1},x_{2},x_{3},y_{1},y_{2},y_{3})=|x_{3}|^{2}|y_{3}|^{2}\mathbf{m}_{\hat{A}}\left(\frac{x}{x_{3}},\frac{y}{y_{3}}\right)$ where we let $x=(x_{1},x_{2})$ and $y=(y_{1},y_{2})$. It thus follows that $\mathcal{SEP}(3,3)$ can be written as

[TABLE]

It can then be readily verified that $\mathcal{C}_{\hat{A}}$ is a hyperplane section of $\mathcal{SEP}(3,3)$ where the appropriate coordinate (corresponding to the monomial $|x_{3}|^{2}|y_{3}|^{2}$) is set to 1.

The case $\text{Sep}(4,2)$: Following the same approach as above, we need to exhibit a Hermitian polynomial $\hat{p}(x_{1},x_{2},x_{3},y)$ supported on

[TABLE]

that is nonnegative but not a sum-of-squares. Since $\text{Sep}(4,2)\neq\text{PPT}(4,2)$ we know that there exists a Hermitian homogeneous polynomial $p(x,y)$ on $(x,y)\in\mathbb{C}^{4}\times\mathbb{C}^{2}$ of the form

[TABLE]

that is nonnegative but not a sum of squares. Note that such a $p$ satisfies $p(\lambda x,\mu y)=|\lambda|^{2}|\mu|^{2}p(x,y)$ for any $(\lambda,\mu)\in\mathbb{C}^{2}$. To get the desired polynomial $\hat{p}$ it would suffice to dehomogenize the polynomial $p$ by setting one of the $x$ variables to 1, and one of the $y$ variables to 1. It turns out, however, that one cannot guarantee in general that this dehomogenized polynomial is not a sum of squares. (We give an explicit example in Appendix A.) The reason we could dehomogenize the Choi polynomial in the $(3,3)$ case was that the Choi polynomial is not a sum of squares when the variables are real. One cannot expect this to be true for our polynomial in (4,2) variables as it is known that any biquadratic real polynomial in $(n,2)$ variables is a sum of squares [Cal73]. Nevertheless we show in Appendix A that by choosing an appropriate polynomial $p$, and an appropriate dehomogenization we can get a polynomial $\hat{p}(x,y)$ supported on $\hat{A}$ of Equation (14) that is not a sum-of-squares. This implies that $\mathcal{C}_{\hat{A}}$ is not semidefinite representable. Using a similar argument as for the $(3,3)$ case we get that $\mathcal{SEP}(4,2)$ is not semidefinite representable.

Appendix A The case $\text{Sep}(4,2)$

In this section we exhibit a nonnegative Hermitian polynomial $\hat{p}(x_{1},x_{2},x_{3},y)$ supported on $\hat{A}$, defined in Equation (14), that is not a sum of squares.

Consider the following map $\Phi:\text{M}_{2}\rightarrow\text{M}_{4}$ studied in [HK16]:

[TABLE]

It is shown in [HK16] that the map $\Phi$ is positive but not decomposable. We associate to $\Phi$ the Hermitian polynomial

[TABLE]

Then we know from Proposition 1 that $\mathbf{W}$ is positive but not a sum of squares. The purpose of this section is to prove the following:

Proposition 3**.**

The (nonhomogeneous) Hermitian polynomial $\mathbf{W}(1,x_{2},x_{3},x_{4},1,y_{2})$ is not a sum of squares.

The proof of this proposition involves some computations, and we use the specific properties of $\mathbf{W}$ (its zeros) which have been studied in [HK16]. We note that there are other dehomogenizations of $\mathbf{W}$ that are sums of squares. For example, we show later that $\mathbf{W}(x_{1},x_{2},-6,x_{4},1,y_{2})$ is a sum of squares.

Proof of Proposition 3.

To lighten the notation we let $y_{2}=\alpha$. In [HK16], the zeros of the polynomial $\mathbf{W}$ were identified. Namely it was shown that

[TABLE]

where

[TABLE]

Let $p(x_{2},x_{3},x_{4},\alpha)=\mathbf{W}(1,x_{2},x_{3},x_{4},1,\alpha)$. The explicit formula of $p$ is

[TABLE]

Assume that $p=\sum_{i}g_{i}^{2}$ where $g_{i}$ are Hermitian polynomials in $x_{2},x_{3},x_{4},\alpha$. Since $p$ has no terms $|x_{3}|^{2}$ we see that the $g_{i}$ cannot contain monomials $x_{3}$ or $\bar{x_{3}}$. Similarly $p$ does not have a term $|\alpha|^{2}|x_{2}|^{2}$ and so $p$ cannot contain monomials $\alpha x_{2}$, $\bar{\alpha}x_{2}$, $\alpha\bar{x_{2}}$ or $\bar{\alpha}\bar{x_{2}}$. It follows that each $g_{i}$ must be a linear combination of the monomials

[TABLE]

and their conjugates. In other words, each $g_{i}$ is of the form:

[TABLE]

where $a_{i},b_{i},\ldots,h_{i}\in\mathbb{C}$. We will now use the information about the zeros of $\mathbf{W}$ (and thus of $p$) to deduce relations about these coefficients and reach a contradiction.

Since $\mathbf{W}$ is bihomogeneous in the first set of variables, we have (dividing by $|x_{1}(\alpha)|^{2}$) from (15) that

[TABLE]

Since $p=\sum_{i}g_{i}^{2}$ we get that for all $i$,

[TABLE]

We can clear denominators in (17) by multiplying the expression by $|x_{1}(\alpha)|^{2}$. As a result we get that

[TABLE]

The left-hand side of (18) is a Hermitian polynomial in $\alpha$ that is identically zero. Hence all its coefficients must be equal to 0. This allows us to derive conditions on the coefficients of $g_{i}$ in (16). More precisely:

• •

The coefficient of $\alpha^{4}$ is $4\bar{h_{i}}$. Setting $4\bar{h_{i}}$ to zero yields $h_{i}=0$.

• •

The coefficient of $\alpha^{4}\bar{\alpha}$ is $4\bar{g_{i}}+2\bar{h_{i}}$. Setting to zero we get $g_{i}=0$.

• •

The coefficient of $\alpha^{2}$ is $-4\bar{d_{i}}-8\bar{g_{i}}$. Setting to zero we get $d_{i}=0$.

This gives a contradiction: indeed the coefficient of $|x_{4}|^{2}$ in $p$ is $4>0$ and yet $\sum_{i}|d_{i}|^{2}=0$. ∎

We conclude this appendix by proving, as promised, that there is another dehomogenization of $\mathbf{W}$ that is a sum of squares. Let:

[TABLE]

Let

[TABLE]

One can verify that $A-B\geq 0$ and $A+B\geq 0$, i.e., that $\left[\begin{smallmatrix}A&B\\ B&A\end{smallmatrix}\right]\geq 0$. Let $m(x_{1},x_{2},x_{4},\alpha)=(\alpha,x_{1},x_{2},x_{4},\bar{\alpha}x_{1},\bar{\alpha}x_{4})$. Then one can check that we have the following sum of squares decomposition of $q$:

[TABLE]

Acknowledgements

I would like to thank Omar Fawzi for his encouragements and for helpful comments on the paper. I would also like to thank Claus Scheiderer for useful discussions and exchanges related to the material of this paper, and James Saunderson for comments that helped improved the exposition.