REVIEW 3 major objections 3 minor 2 cited by
The sum-of-squares hierarchy on the sphere, and applications in quantum information theory
T0 review · 3 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper proves that the sum-of-squares hierarchy on the sphere converges at rate O(d^2/ℓ^2) and transfers that rate to quantum entanglement bounds.
desk verdict A valuable paper with a genuinely new quadratic rate and a clean DPS/SOS duality, but Section 3 contains two false constant bounds that must be fixed before it is citable as is. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The engine is the polynomial kernel $K(x,y) = q(\langle x,y\rangle)^2$ on $\mathbb{S}^{d-1} \times \mathbb{S}^{d-1}$. The Funk-Hecke formula diagonalizes this integral operator: its eigenvalues are the Gegenbauer coefficients $\lambda_{2k}$ of the univariate polynomial $\varphi(t)=q(t)^2$, and if those coefficients are close enough to $1$, then $K^{-1}(F+\delta I)$ is nonnegative and $F+\delta I = K K^{-1}(F+\delta I)$ is an explicit sum of squares. Choosing $q$ is recast as a generalized Toeplitz eigenvalue problem: $\rho_{2n}(d,\ell) = \min_{\|e\|=1} \sum_{k=1}^n |(e^T T[C_{2k}/C_{2k}(1)]e)^{-1}-1|$, with $T[h]$ the truncated multiplication matrices for normalized Gegenbauer polynomials. The asymptotic rate comes from lower-bounding the largest eigenvalue of $T[h]$ for $h = \frac{1}{n}\sum_{k=1}^n C_{2k}/C_{2k}(1)$ via the largest root of $C_{\ell+1}$, and the matrix-valued case follows by applying the same harmonic-component bounds entrywise with spectral norms.
What would settle it
Compute the largest zero of $C_{\ell+1}$ for $d=2$, $\ell=10$: it is about $0.9898$, below the lower bound $1 - d^2/(4\ell^2) = 0.99$, so checking the quoted zero bound over the paper's full parameter range would settle whether the proof's constants hold. One can also evaluate $h'(1) = (n+1)(3d+4n-4)/(3(d-1))$ at $n=d$ and compare it with $7n/3$.
Extended reading notes
Core claim
The central claim is Theorem 2: if $F(x)$ is a homogeneous symmetric matrix-valued polynomial of degree $2n$ in $d$ variables with $n \le d$ and $0 \le F(x) \le I$ on $\mathbb{S}^{d-1}$, then for every $\ell \ge C_n d$ the polynomial $F + C'_n (d/\ell)^2 I$ is a sum of squares of polynomials of degree at most $\ell$ on the sphere. For the scalar polynomial $F = (p_{\max}-p)/(p_{\max}-p_{\min})$, this gives $(p_\ell-p_{\min})/(p_{\max}-p_{\min}) \le 1 + (C_n d/\ell)^2$, a quadratic improvement over the previously known rate. The more general Theorem 6 bounds every level $\ell$ through a quantity $\rho_{2n}(d,\ell)$ that can be computed for degree $2$ and $4$, and the $O(d^2/\ell^2)$ rate is shown tight for quadratic polynomials. A duality theorem identifies the DPS hierarchy of separable quantum states with a real-sum-of-squares condition on Hermitian polynomials, so the same rate transfers to the Best Separable State problem: $h_{\mathrm{Sep}}(M) \le h_{\mathrm{DPS}_\ell}(M) \le (1 + C d_B^2/\ell^2) h_{\mathrm{Sep}}(M)$ for $\ell \ge C' d_B$.
Load-bearing premise
The proof's explicit constants in the $O(d^2/\ell^2)$ estimate rest on two literature bounds holding throughout the asserted range: the largest zero of the Gegenbauer polynomial $C_{\ell+1}$ is at least $1 - d^2/(4\ell^2)$, and the derivative $h'(1)$ of $h = \frac{1}{n}\sum_{k=1}^n C_{2k}/C_{2k}(1)$ is at most $7n/3$; if either fails, those constants are not established.
Editorial extensions
If this is right
- For a homogeneous polynomial of fixed degree $2n$, reaching relative accuracy $\varepsilon$ requires level $\ell = O(d/\sqrt{\varepsilon})$ rather than $O(d/\varepsilon)$.
- The certificate level is independent of the size of the matrix-valued polynomial, so operator-valued or multi-constraint problems do not require deeper levels.
- The DPS hierarchy approximates the best separable state value with relative error $O(d_B^2/\ell^2)$ at level $\ell \ge C d_B$, with constants independent of the other subsystem dimension.
- The duality theorem makes the quantum convergence rate a special case of the SOS rate, so the same bound applies to nonquadratic Hermitian polynomials.
Reading between the lines
- A direct numerical computation of $\rho_{2n}(d,\ell)$ for small $d$ would test whether the $O(d^2/\ell^2)$ scaling persists outside the proven regime $\ell \ge C_n d$.
- The same kernel construction transfers naturally to products of spheres, since a matrix-valued polynomial in $x$ is a bihomogeneous polynomial in $(x,y)$; this would give SOS-based bounds for multilinear optimization.
- If the quadratic-case tightness is representative, the $d^2/\ell^2$ rate is likely the true order of the SOS upper-bound hierarchy on the sphere, and further speedups would require a different relaxation family rather than a sharper kernel.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the convergence rate of the Lasserre/sum-of-squares hierarchy for maximizing a homogeneous polynomial of degree 2n on the unit sphere S^{d-1}. Its main result, Theorem 2, states that for n ≤ d and a homogeneous matrix-valued polynomial F of degree 2n with 0 ≤ F ≤ I on the sphere, F + C'_n (d/ℓ)^2 I is ℓ-sos whenever ℓ ≥ C_n d, with constants depending only on n; Theorem 1 is the scalar corollary asserting (p_ℓ − p_min)/(p_max − p_min) ≤ 1 + (C_n d/ℓ)^2. The proof uses positive kernels K(x,y) = q(⟨x,y⟩)^2 and analyzes their Gegenbauer coefficients through generalized Toeplitz matrices. The paper also proves a duality theorem (Theorem 12) identifying the dual of the DPS hierarchy in quantum information with real sum-of-squares certificates of the associated Hermitian polynomial, and uses it to transfer the rate result to the Best Separable State problem, recovering and generalizing the quadratic convergence rate of Navascués, Owari and Plenio. A tightness result (Theorem 10) shows the quadratic rate cannot be improved in the quadratic case.
Significance. If the proof is repaired, the paper settles an open question raised by de Klerk and Laurent by improving the known O(d/ℓ) rate for the upper-bound SOS hierarchy on the sphere to O(d^2/ℓ^2). The matrix-valued formulation with constants independent of the matrix size is a real strength, as is the self-contained proof of the DPS/SOS duality, which fills a gap in the literature. The recovery of the known NOP09 rate provides a useful external check, and the tightness analysis for quadratic polynomials is a good addition. However, the proof of the general-degree rate rests on Proposition 7, and two inequalities used there are false as stated. Because Proposition 7 is the only step producing the O(d^2/ℓ^2) estimate on ρ_{2n}, the constants in Theorems 6, 2, and 1 are not established in the present manuscript. The errors concern constants and thresholds rather than the overall proof architecture, so the central claim remains plausible and likely repairable.
major comments (3)
- [Section 3, Proposition 7] The proof invokes the bound x_{ℓ+1,ℓ+1} ≥ 1 − d^2/(4ℓ^2) from [DJ12, Section 2.3]. This bound is false for d = 2, which is in the admissible range n ≤ d. For d = 2 the Gegenbauer polynomials are (up to scaling) the Chebyshev polynomials, so the largest zero of C_{ℓ+1} is cos(π/(2(ℓ+1))) = 1 − π^2/(8ℓ^2) + O(ℓ^{-3}), which is strictly smaller than 1 − 1/ℓ^2 for ℓ = 10 and asymptotically for all larger ℓ. The authors need either a valid universal root bound with a larger constant, or an explicit dimension restriction together with a separate argument for the excluded small dimensions.
- [Section 3, Proposition 7] The assertion h'(1) ≤ 7n/3, derived from h'(1) = (n+1)(3d+4n−4)/(3(d−1)) and n ≤ d, is false. For n = d = 2 one has h'(1) = 10, while 7n/3 = 14/3. In fact, for d = n the exact value is h'(1) = 7n/3 + (10n−4)/(3(n−1)), which exceeds 7n/3 for every n ≥ 2. Consequently the final lower bound λ_max(T[h]) ≥ 1 − (7n/12)d^2/ℓ^2 is not proven. A correct uniform bound such as h'(1) ≤ 5n for n ≤ d would repair the argument with a different constant.
- [Section 3, Theorems 6, 2, and 1] Because Proposition 7 is the only step that bounds ρ̃_{2n}(d,ℓ) by O(d^2/ℓ^2), and because Proposition 9 transfers that bound to ρ_{2n}(d,ℓ), the false inequalities in Proposition 7 leave the constants C_n and C'_n in Theorems 6, 2, and 1 unproven as written. The O(d^2/ℓ^2) rate may well survive with corrected constants, but the current proof is incomplete at this load-bearing point.
minor comments (3)
- [Section 4.4, Theorem 15] The phrase “where C, C′ > 0 is some absolute constant” should be “where C, C′ > 0 are absolute constants.”
- [Section 3, Eq. (13)] The inverse λ_{2k}^{-1} in the definition of ρ_{2n}(d,ℓ) presupposes λ_{2k} ≠ 0; the authors could state explicitly that the minimization can be restricted to polynomials for which the even Gegenbauer coefficients are positive, since otherwise the objective is not well defined.
- [Section 4.1, Eq. (22)] The notation in Eq. (22) uses a bracehtipup/down artefact in the displayed formula; this should be cleaned up to a standard partial-transpose notation such as (I_A ⊗ T_{B_1} ⊗ ⋯ ⊗ T_{B_s} ⊗ I_{B_{s+1}} ⊗ ⋯ ⊗ I_{B_ℓ}).
Circularity Check
No circularity: derivation is self-contained; Proposition 7's constant gap is a correctness issue, not circularity.
full rationale
I found no circular step in the paper. The central result, Theorem 2, is proved from an explicit kernel construction: for a nonnegative normalized polynomial F, the paper writes F + δ I = K K^{-1}(F + δ I) with K defined by φ(⟨x,y⟩) = q(⟨x,y⟩)^2, and the error ‖K^{-1}(F + δ I) - (F + δ I)‖∞ is controlled by the quantity ρ_{2n}(d,ℓ). The proof of the O(d²/ℓ²) rate then reduces to the analysis of generalized Toeplitz matrices associated with Gegenbauer polynomials (Proposition 7 and surrounding text). This analysis uses external, explicitly cited results on extreme zeros (DJ12, ADGR04) and standard orthogonal-polynomial facts (Par65, DLMF); no parameter is fitted to the target rate, no prior work by the same authors is invoked as load-bearing, and the recovered Navascués-Owari-Plenio DPS rate is presented as an external consistency check rather than an input. Theorem 1 follows from Theorem 2 by the legitimate normalization F = (p_max - p)/(p_max - p_min), which is an application of the theorem, not a disguised restatement of the conclusion. I note a genuine correctness gap in the written constants: the cited bound x_{ℓ+1,ℓ+1} ≥ 1 - d²/(4ℓ²) fails for d=2 at the relevant Chebyshev zeros, and the asserted inequality h'(1) ≤ 7n/3 from the displayed formula h'(1) = (n+1)(3d+4n-4)/(3(d-1)) and n ≤ d is false at n=d=2, where h'(1)=10 > 14/3. That is an arithmetic/root-bound error affecting the stated constants in Proposition 7, but it is not circularity: the conclusion is not equivalent to an input by construction, and a corrected constant would preserve the O(d²/ℓ²) rate. Since circularity is about the reduction of the claimed derivation to its own inputs, and no such reduction occurs, the appropriate circularity score is 0.
Assumptions & free parameters
assumptions (5)
- standard math Funk-Hecke formula and Gegenbauer basis diagonalize SO(d)-invariant kernels on the sphere.
- standard math Eigenvalues of T[f] for linear f are the roots of p_{ℓ+1} (Parter).
- domain assumption Gegenbauer C_i lies above its tangent at t=1 (Proposition 18).
- domain assumption Largest zero bound x_{ℓ+1,ℓ+1} ≥ 1 - d^2/(4ℓ^2) from Driver-Jordaan.
- standard math Integral SOS representations can be truncated to finite SOS by convexity.
Cite this review
Pith. "Pith review of The sum-of-squares hierarchy on the sphere, and applications in quantum information theory." pith.science (2026). https://pith.science/paper/JE3MNTER
@misc{pith2026190805155,
author = {Pith},
title = {Pith review of: The sum-of-squares hierarchy on the sphere, and applications in quantum information theory},
year = {2026},
howpublished = {\url{https://pith.science/paper/JE3MNTER}},
note = {Machine review of arXiv:1908.05155}
}
abstract
We consider the problem of maximizing a homogeneous polynomial on the unit sphere and its hierarchy of Sum-of-Squares (SOS) relaxations. Exploiting the polynomial kernel technique, we obtain a quadratic improvement of the known convergence rate by Reznick and Doherty & Wehner. Specifically, we show that the rate of convergence is no worse than $O(d^2/\ell^2)$ in the regime $\ell \geq \Omega(d)$ where $\ell$ is the level of the hierarchy and $d$ the dimension, solving a problem left open in the recent paper by de Klerk & Laurent (arXiv:1904.08828). Importantly, our analysis also works for matrix-valued polynomials on the sphere which has applications in quantum information for the Best Separable State problem. By exploiting the duality relation between sums of squares and the DPS hierarchy in quantum information theory, we show that our result generalizes to nonquadratic polynomials the convergence rates of Navascu\'es, Owari & Plenio.
Figures
Forward citations
Cited by 2 Pith papers
-
The power of unentanglement without destructive interference
StoqMA(2) contains NP with Õ(√n)-qubit proofs and completeness error 2^{-polylog(n)}, is contained in EXP, and satisfies StoqMA(k)=StoqMA(2) for k≥2 when completeness error is negligible.
-
A refinement of Reznick's Positivstellensatz with applications to quantum information theory
The complex Reznick Positivstellensatz is proved with improved bounds via an explicit inversion of the Chiribella identity, yielding better exponential de Finetti error estimates; the real-case version is not yet esta...
Reference graph
Works this paper leans on
-
[1]
I. Area, D. Dimitrov, E. Godoy, and A. Ronveaux. Zeros of G egenbauer and H ermite polynomials and connection coefficients. Mathematics of Computation , 73(248):1937--1951, 2004
work page 1937
-
[2]
J. V. Baxley. Extreme eigenvalues of T oeplitz matrices associated with certain orthogonal polynomials. SIAM Journal on Mathematical Analysis , 2(3):470--482, 1971
work page 1971
- [3]
-
[4]
V. Bhattiprolu, M. Ghosh, V. Guruswami, E. Lee, and M. Tulsiani. Weak decoupling, polynomial folds and approximate optimization over the sphere. In 2017 IEEE 58th Annual Symposium on Foundations of Computer Science (FOCS) , pages 1008--1019. IEEE, 2017
work page 2017
- [5]
-
[6]
G. Blekherman. Convexity properties of the cone of nonnegative polynomials. Discrete & Computational Geometry , 32(3):345--371, 2004
work page 2004
- [7]
-
[8]
M. Christandl, R. K \"o nig, G. Mitchison, and R. Renner. One-and-a-half quantum de F inetti theorems. Communications in Mathematical Physics , 273(2):473--498, 2007
work page 2007
Show all 31 references
-
[9]
Driver and K
K. Driver and K. Jordaan. Bounds for extreme zeros of some classical orthogonal polynomials. Journal of Approximation Theory , 164(9):1200--1204, 2012
2012
-
[10]
De Klerk
E. De Klerk. The complexity of optimizing over a simplex, hypercube or sphere: a short survey. Central European Journal of Operations Research , 16(2):111--125, 2008
2008
-
[11]
de Klerk and M
E. de Klerk and M. Laurent. Convergence analysis of a L asserre hierarchy of upper bounds for polynomial minimization on the sphere. arXiv preprint arXiv:1904.08828 , 2019
1904 arXiv
-
[12]
de Klerk, M
E. de Klerk, M. Laurent, and P. Parrilo. On the equivalence of algebraic approaches to the minimization of forms on the simplex. In Positive Polynomials in Control , pages 121--132. Springer, 2005
2005
-
[13]
http://dlmf.nist.gov/, Release 1.0.22 of 2019-03-15
NIST Digital Library of Mathematical Functions . http://dlmf.nist.gov/, Release 1.0.22 of 2019-03-15. F. W. J. Olver, A. B. Olde Daalhuis , D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller and B. V. Saunders, eds
2019
-
[14]
J. P. D'Angelo and M. Putinar. Polynomial optimization on odd-dimensional spheres. In Emerging applications of algebraic geometry , pages 1--15. Springer, 2009
2009
-
[15]
A. C. Doherty, P. A. Parrilo, and F. M. Spedalieri. Distinguishing separable and entangled states. Physical Review Letters , 88(18):187904, 2002
2002
-
[16]
A. C. Doherty, P. A. Parrilo, and F. M. Spedalieri. Complete family of separability criteria . Physical Review A , 69(2):022308, 2004
2004
-
[17]
A. C. Doherty and S. Wehner. Convergence of SDP hierarchies for polynomial optimization on the hypersphere. arXiv:1210.5048 , 2012
2012 arXiv
-
[18]
L. Gurvits. Classical deterministic complexity of E dmonds' problem and quantum entanglement. In Proceedings of the thirty-fifth annual ACM symposium on Theory of computing , pages 10--19. ACM, 2003
2003
-
[19]
Horodecki, P
M. Horodecki, P. Horodecki, and R. Horodecki. Mixed-state entanglement and quantum communication. In Quantum information , pages 151--195. Springer, 2001
2001
-
[20]
H.-y. Hsu. Certain integrals and infinite series involving ultra-spherical polynomials and Bessel functions . Duke Mathematical Journal , 4(2):374--383, 1938
1938
-
[21]
Koenig and G
R. Koenig and G. Mitchison. A most compendious and facile quantum de F inetti theorem. Journal of Mathematical Physics , 50(1):012105, 2009
2009
-
[22]
J. B. Lasserre. Global optimization with polynomials and the problem of moments. SIAM Journal on optimization , 11(3):796--817, 2001
2001
-
[23]
Lewenstein, D
M. Lewenstein, D. Bru , J. I. Cirac, B. Kraus, M. Ku \' s , J. Samsonowicz, A. Sanpera, and & . R. Tarrach. Separability and distillability in composite quantum systems-a primer . Journal of Modern Optics , 47(14):2481--2499, 2000
2000
-
[24]
Nesterov
Y. Nesterov. Random walk in a simplex and quadratic optimization over convex polytopes. Technical report, CORE, 2003
2003
-
[25]
Navascues, M
M. Navascues, M. Owari, and M. B. Plenio. The power of symmetric extensions for entanglement detection . Physical Review A - Atomic, Molecular, and Optical Physics , 80(5):1--16, 2009
2009
-
[26]
S. V. Parter. Remarks on the extreme eigenvalues of T oeplitz forms associated with orthogonal polynomials . Journal of Mathematical Analysis and Applications , 12(3):456--470, 1965
1965
-
[27]
P. A. Parrilo. Structured semidefinite programs and semialgebraic geometry methods in robustness and optimization . PhD thesis, California Institute of Technology, 2000
2000
-
[28]
P. A. Parrilo. Approximation quality of SOS relaxations, 2013. Talk at ICCOPT 2013
2013
-
[29]
P \'o lik and T
I. P \'o lik and T. Terlaky. A survey of the s-lemma. SIAM Review , 49(3):371--418, 2007
2007
-
[30]
B. Reznick. Uniform denominators in Hilbert's seventeenth problem . Mathematische Zeitschrift , 220(1):75--97, 1995
1995
-
[31]
S. L. Woronowicz. Positive maps of low dimensional matrix algebras. Reports on Mathematical Physics , 10(2):165--183, 1976
1976
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.