精读笔记
Problem Setting
论文实际处理的是超立方体上非负多项式的“effective SOS approximation”问题：不是证明 SOS cone 稠密，而是给出随 SOS degree r 变化的显式误差率，并把它转化为 Lasserre 层级的 convergence rate。
真正困难点在于三者要同时满足：输出必须是 SOS；算子必须接近 identity；误差要在系数 ℓ1 范数中可控。已有 Fejér/Christoffel–Darboux 风格 kernel 能给正性和近似性，但在超立方体 Chebyshev 系数范数下主要得到 O(1/r) 级别。卡点本质上是：一阶平滑核的 bias 由权重的一阶差分控制，无法自然产生二阶误差。
这篇论文绕开 monomial basis，转到 Chebyshev basis。这个选择不是 cosmetic：Chebyshev 正交性和乘积公式使 kernel 对基函数的作用可以被精确展开，从而能看清误差项到底是一阶还是二阶。代价是最终结论主要是 Chebyshev ℓ1，而非最初 Lasserre projection 中最自然的 monomial ℓ1。
Motivation
已有路线不够的地方在于：它们知道如何用正核/SOS 核近似非负多项式，但核的归一化和截断误差通常只给 O(1/r)。另一方面，Lasserre measure-based upper bounds 和一些正交多项式极值理论中已经出现 O(1/r^2)，暗示超立方体 lower hierarchy 的 O(1/r) 可能只是分析或核设计不够好，而非根本障碍。
作者的关键缺口判断是：需要一个仍可转成 SOS 的正核，但其 identity approximation error 应该由“权重差分平方”而非“权重差分”控制。平方核正好提供这个结构。Fejér kernel 本身非负；平方后自动正，并且在 Chebyshev product formula 下产生二阶误差项。这是本文最重要的观察。
为什么会想到这个方向：upper-bound spectral hierarchy 的最优 density 是平方多项式形式，形如 σ(y)=(∑γ_α T_α(y))^2。作者实际上把这种 squared density 的结构反向用于 lower-bound SOS approximation：先构造一个平方核，再用它作为正线性算子逼近 identity。
Core Idea
核心思想可以压缩成一句：用 squared Fejér-type Chebyshev kernel 构造一个正的近似恒等算子，并通过 SOS 归一化把它变成对非负多项式输出 SOS 的算子。不同于普通 kernel smoothing，这里的平方操作改变了误差代数结构：当核作用在 T_k 上时，误差不再主要由 |a_{j+k}-a_j| 累积，而由 (a_{j+k}-a_j)^2 累积，因此在选择 a_j=j/r 的 Fejér 权重时得到近似二阶的行为。
这个建模方式的本质变化是：不再把核看成一个直接归一化的 averaging operator，而是先允许 squared kernel 产生位置相关质量 M_r(x)，再用可控的 SOS 多项式近似 1/M_r(x) 做归一化。也就是说，论文把“正性”和“归一化”拆开处理：平方核负责正性和二阶误差，截断几何级数负责把归一化重新塞回 SOS cone。
相对 prior 的本质区别不是用了 Chebyshev，也不是用了 Fejér，而是“平方后再归一化”的误差机制。Chebyshev basis 和 Fejér 权重是让该机制可分析、可实现的载体。
Method
1. Chebyshev basis 下的 kernel calculus：解决的是如何精确追踪 kernel 对低次多项式的扰动。需要它是因为 monomial basis 中缺乏这种干净的正交卷积结构。核心变化是将 SOS 逼近问题转化为对每个 T_k 的 operator bias 控制。
2. Squared univariate kernel：构造 S_r^{(a)} 后使用 K=(S)^2。它解决的是普通正核误差只有一阶的问题。关键推导显示 K(T_k)-M T_k 的 Chebyshev 系数由 a_k^2 和权重差分平方控制。这个平方差分结构是 rate 改善的来源。
3. SOS 归一化：由于 K(1)=M_r^{(a)}(x) 不是常数，必须除以 M 才能近似 identity。直接除法会离开多项式/SOS 世界，因此用截断几何级数 p_N≈1/M，并利用偶数 N 时几何和是 SOS。它解决的是“归一化 vs SOS 保持”的冲突。代价是 degree 增加一个 log r 因子。
4. 多变量 product kernel：通过逐坐标相乘扩展到 [-1,1]^n。它解决的是多变量构造问题，但没有捕捉跨坐标结构；本质上是 tensorization。误差通过 product perturbation lemma 累积，带来 n 相关因子。
5. 从基函数误差到一般 f：对 f=∑f_αT_α 线性叠加，用 Chebyshev ℓ1 范数吸收所有系数。这一步简单但决定了最终 bound 依赖 ||f||_{1,cheb}，也解释了为何该结果对系数稀疏/低 Chebyshev complexity 的多项式更友好。
Key Insight / Why It Works
最核心的贡献是 squared kernel 使误差二阶化。普通 Fejér-type kernel 的近似误差本质上是 spectral window 的一阶平滑误差；本文通过平方 kernel，在 Chebyshev product algebra 中让误差出现为权重差分平方。对 a_j=j/r，有 |a_{j+k}-a_j|≈k/r，因此平方后累积为 roughly k^2/r，再除以归一化质量 M≈r，得到 k^2/r^2 量级。这个机制非常干净，也很可能是本文最可迁移的 insight。
归一化步骤是必要辅助，但不是核心创新。p_N≈1/M 的截断几何级数保证 SOS，是技术上优雅的构造；不过它引入 log r degree overhead，并最终把单变量近似中的 log r 推成主定理中的 log^3 r。这里的 log 损失看起来更像 proof/degree bookkeeping 的副产物，而不是已证明的结构性下界。文中未充分说明 log^3 是否不可避免。
多变量部分基本是 tensorization，不是新的结构发现。它证明了机制可扩展，但也暴露了上限：product kernel 对变量交互没有特殊利用，n^3 因子和 r/log r ≥ Θ(n^2d^2) 的门槛很自然地出现。这里谈不上 hidden scalability breakthrough。
从 Lasserre lower hierarchy 的角度看，本文真正做的是把一个 SOS approximation theorem plug into 已有的误差转移框架。hierarchy 收敛率的改善不是来自新的 SDP 分析，而来自更好的 Positivstellensatz/SOS approximation kernel。
对 upper-bound spectral hierarchy，本文的新证明更像 conceptual unification：用同一个 squared kernel 构造 feasible density，从而重证 O(1/r^2)。这部分不是改善 rate，但解释了 measure-based hierarchy 中 O(1/r^2) 的核方法来源。
Relation To Prior Work
最接近的是 Baldi–Slot 的超立方体 Putinar/Lasserre 收敛率分析，以及 Gribling–de Klerk–Vera 之前用 polynomial kernel method 得到 O(1/r) 的路线。本文沿用同一技术谱系：Chebyshev basis、kernel polynomial method、SOS approximation、Lasserre hierarchy error transfer。真正不同点是 kernel 设计：从普通 Fejér/positive kernel 变成 squared Fejér kernel，并引入 SOS 归一化。
与 Lasserre 早期 SOS approximation 的关系：Lasserre 给出密度和 monomial ℓ1 projection 结构，但没有显式 rate；本文给出 Chebyshev ℓ1 下的显式 rate。因此它不是替代 Lasserre projection theorem，而是在更可分析的 norm/basis 中给 effective version。
与 De Klerk–Laurent 的 measure-based upper bounds 的关系：后者通过 Jacobi polynomial extremal roots 等非构造分析得到 O(1/r^2)；本文给出 constructive kernel proof，并显式展示 n,d,||f||_{1,cheb} 依赖。这里的新增信息是机制解释，而非更强 rate。
看似新的部分中，product kernel、多变量误差叠加、Chebyshev norm equivalence 都是已有思想重组；实质创新集中在 squared kernel 的误差展开与 SOS 归一化组合。
Dataset / Evaluation
这不是实验论文，没有 dataset，也没有 empirical benchmark。evaluation 完全是理论性质：主指标是 SOS degree r 下的 Chebyshev coefficient ℓ1 approximation error，以及由此推出的 Lasserre lower/upper hierarchy 收敛率。
这些理论结果确实验证了核心 claim：相对于已有 O(1/r)，在超立方体上获得 O(log^3 r/r^2) 的显式 SOS approximation 和 lower hierarchy 收敛率。但验证范围也很明确：它验证的是 worst-case analytic upper bound，不说明实际 SDP 求解中会有同等收益，也不说明 monomial ℓ1 projection 的真实误差达到该 rate。
对 upper hierarchy 的部分，结果支持“squared kernel 可构造 O(1/r^2) density”的 claim，但不是新的 asymptotic rate，只是给出更构造化、参数依赖更清楚的证明。
Limitation
第一，核心结果绑定 Chebyshev ℓ1。虽然文中给出 Chebyshev ℓ1 与 monomial ℓ1 的等价关系，但常数含 (5/2)^d，不能视为无损转换。结论不能直接回答 introduction 中最自然的 monomial projection rate 问题。
第二，log^3 r 是明显的未闭合部分。单变量 squared Fejér 机制本身指向 O(k^2 log r/r^2)，主定理中 log^3 主要来自归一化 degree 与参数重标定。是否可以通过更好的 reciprocal approximation、更直接的 normalized SOS kernel 或 sharper bookkeeping 去掉这些 log，文中未充分说明。
第三，多变量扩展是 product construction，scalability 上限清楚：维度依赖较重，且没有利用 f 的低秩、稀疏交互、变量图结构等。对高维但结构化多项式，这个 bound 可能非常松。
第四，非负性假设是全局的，构造 p=Kf 输出 SOS 依赖 f≥0 以及 kernel 正性。若 f 仅在约束集上非负但表达复杂，或者换成一般 semialgebraic set，需要额外转移定理；论文最后提到可通过 Łojasiewicz 常数得到一般集合上的弱 rate，但那会显著恶化。
第五，Lasserre lower-bound 改善的常数依赖 ||f||_{1,cheb}，而不是 f_max-f_min。对 Chebyshev 系数膨胀的多项式，新 bound 可能在有限 r 区间不优于 Baldi–Slot。这里的增益是 asymptotic rate，而非 uniformly better practical certificate。
Takeaway
- 1. 最值得记住的是 squared kernel trick：把正核平方后，Chebyshev product algebra 会把 identity approximation error 变成权重差分平方，这是从 O(1/r) 推向近 O(1/r^2) 的结构性原因。
- 2. SOS approximation 的一个有效设计原则是拆分“正性”和“归一化”：先用平方结构保证正性和二阶 bias，再用 SOS reciprocal approximation 处理 normalization。
- 这一范式可能可迁移到其他正交多项式体系。
- 3. 论文推动的是 effective Positivstellensatz / Lasserre rate 的核方法，而不是新的 SDP 算法。
一句话总结
这篇论文在 Chebyshev-kernel/SOS approximation 路线中引入 squared Fejér kernel，把超立方体上 Lasserre lower hierarchy 的已知理论收敛率从一阶推进到近二阶，核心贡献是误差二阶化的核设计而非新的 SDP 层级。