精读笔记
Problem Setting
论文处理的问题是凸正则化最小二乘 min_x 1/2||Ax-b||² + λ r(x) 的解映射 S(A,b,λ) 在给定解 x̄ 附近的局部稳定性。这里的稳定性不是泛泛的 well-posedness，而是更细的性质：局部单值、Lipschitz 连续、方向可微、semismooth/semismooth*。
真正困难点在于：最优性条件可以立即写成 0 ∈ (1/λ)Aᵀ(Ax-b)+∂r(x)， 但要判定这个 generalized equation 是否强正则，标准 Mordukhovich criterion 会落到 ker A ∩ ker ∂²r(x̄|v̄)={0}。 这里 ∂²r 是 coderivative of subdifferential，即 generalized Hessian。对于非光滑凸正则项，这个对象通常比原问题还难算。
因此关键矛盾是：稳定性由二阶变分几何决定，但应用中可获得的往往只是正则项的一阶/面结构信息。论文真正解决的是把“二阶 kernel 判别”替换成“共轭子微分平行子空间判别”。
Motivation
已有路线的不足很明确。LASSO/extended ℓ1 类结果依赖 polyhedral 活动集，表达清楚但结构窄；partial smoothness 路线强调 active manifold，但对一般 cone-reducible 或 PLQ penalty 不够统一；Cui et al. 的对偶强正则性路线给出了漂亮的一阶几何判别，但更偏 dual construction，定量导数和 primal generalized equation 机制没有完全展开。
作者的核心观察是：在隐函数判别中真正进入条件的不是整个 ∂²r，而只是其 kernel。若能刻画 kernel，就能绕开最难的二阶对象。对 conjugate C²-cone reducible regularizer，这个 kernel 恰好等于 par ∂r*(v̄)，后者通常由 face/active constraints/normal cone 直接给出。
所以这篇的动机不是提出新的优化算法，而是补上一个理论接口：从单调广义方程的强正则性，到 regularized LS 的一阶可检验稳定性条件，再到结构化正则项的定量 sensitivity formula。
Core Idea
核心思想可以压缩为一句：把解映射稳定性看成 primal monotone generalized equation 的隐函数问题，并把其中不可计算的二阶正则项 kernel 用共轭侧的一阶平坦子空间替代。
这改变了建模方式。传统上，非光滑正则化问题的稳定性要么沿着显式活动集做 case analysis，要么在对偶问题上构造强正则性；本文直接在 primal inclusion 上工作：smooth least-squares part 提供 AᵀA 的曲率，regularizer 提供 ∂r 的最大单调性，二者的横截性决定解是否稳定。横截性的最终形式是 ker A ∩ par ∂r*(v̄)={0}。 这非常直观：如果 A 在正则项无法提供二阶辨识的平坦方向上也看不见扰动，那么解必然不稳定；反之，只要 A 在这些方向上非退化，局部唯一性和 Lipschitz 稳定性就成立。
和 prior 的本质区别在于，它不是为每类正则项重新做活动集分析，而是抽象出“平坦方向 = 共轭子微分的 parallel subspace”这一可迁移机制。这个机制比显式 polyhedral 计算更 general，也比完整 generalized Hessian 更 scalable。
Method
1. Monotone generalized equation 隐函数定理： 论文先证明 Proposition 3.6。对象是 S(p)={x | 0∈f(p,x)+F(x)}，其中 F maximally monotone，f 对 x 单调且光滑。它解决的问题是：如何从 coderivative criterion 推出解映射的局部 Lipschitz、方向可微和 semismoothness。关键变化是利用单调性把 metric regularity 升级到 strong metric regularity，从而得到局部单值。
2. Mordukhovich criterion 的结构化简： 对 least-squares inclusion，D_x f=(1/λ)AᵀA，criterion 化为 ker A ∩ ker ∂²r(x̄|v̄)={0}。 这一步把数据项和正则项的作用完全分离：A 控制观测方向，∂²r 的 kernel 控制正则项平坦方向。
3. generalized Hessian kernel 表达： Theorem 4.6 是技术中心。它证明在 g* C²-cone reducible 时， ker ∂²g(x̄|v̄)=par ∂g*(v̄)。 证明借助 tilt stability 的等价刻画：如果某方向落入 generalized Hessian kernel，就能构造一个二次光滑项使 tilt stability 失败；反过来利用 C²-cone reducibility 提供的二阶 epi-differentiability/regularity 得到等号。
4. 结构化正则项实例化： 对 r=σ_P∘M，par ∂r*(v̄) 可写为 span M^{-1}(N_P(ȳ))；对 PLQ penalty θ_{P,B}∘M，在 compactness qualification 下仍能得到 kernel=par ∂r*(v̄)。这些结果把抽象条件变成 active normal cone 上的线性代数条件，并进一步给出方向导数公式。
Key Insight / Why It Works
最重要的 insight 是：稳定性的本质不是“正则项是否强凸”，而是“loss 的曲率是否覆盖正则项的平坦方向”。对 least squares，loss 在 ker A 上没有曲率；对 nonsmooth r，某些方向上 ∂r 的图形也可能不提供二阶抵抗。二者若重合，就出现多解或非 Lipschitz；若横截，则稳定。
这也是为什么 par ∂r*(v̄) 是正确对象。共轭子微分 ∂r*(v̄) 描述在当前 dual certificate v̄ 下 primal 解可能移动的 face；其 parallel subspace 就是正则项在当前 certificate 下不能区分的局部平坦自由度。把它和 ker A 相交，正好检测“数据项看不见、正则项也不辨识”的方向。
最核心贡献显然是 Theorem 4.6，而不是后面的公式推导。后面 polyhedral/PLQ 的导数和 Lipschitz bound 更多是把该机制落地，技术上有价值，但本质上是 active-cone 线性代数。Proposition 3.6 也重要，但更像把已有 variational-analysis 工具重新组织成适合 regularized LS 的接口。
这不是 scaling，不是数据覆盖，不是 benchmark engineering；这是一个理论归约：把不可计算的二阶 coderivative kernel 归约到可计算的共轭侧一阶几何。值得迁移的正是这种“只刻画二阶对象中真正进入正则性判别的部分”的思路。
Relation To Prior Work
这篇属于 Mordukhovich/Rockafellar-Wets/Dontchev-Rockafellar 变分分析谱系，靠 graphical derivative、coderivative、metric regularity、tilt stability 推稳定性。和 Bolte/Pauwels/Silveti-Falls 的 conservative Jacobian/path differentiability 路线相比，它更偏 classical variational stability，优势是条件与正则项几何直接绑定；劣势是前提更重、表达更抽象。
和 Berk et al. 的 LASSO/Square-root LASSO 分析最接近。本文把那类 LASSO 特例中的 nullspace-active-set 条件抽象成 ker A ∩ par ∂r*(v̄)=0，并推广到 conjugate C²-cone reducible 与 PLQ penalty。
和 Cui et al. 2026 的工作关系最关键。Cui et al. 已经通过对偶强正则性给出 C²-cone reducible conjugate 下的稳定性条件。本文的新增点是：走 primal monotone generalized equation，而不是 dual strong regularity；并补充 semismoothness、方向导数、定量 Lipschitz bound，以及 generalized Hessian kernel 的显式机制解释。
哪些是重组？metric regularity、Mordukhovich criterion、tilt stability、C²-cone reducibility 都是已有工具。实质创新在于把这些工具拼成一个非常干净的桥：ker ∂²r = par ∂r*，然后用它统一 regularized LS 稳定性。
Dataset / Evaluation
这篇没有数据集或实验评估，evaluation 是纯理论的。验证方式是：在一般 convex regularized least-squares 上给出定理；在 weighted polyhedral support functions、LASSO、PLQ penalties 上实例化条件、方向导数和 Lipschitz bound。
从理论论文角度，覆盖范围是合理的：polyhedral support 覆盖 weighted LASSO/analysis sparsity 类，PLQ penalty 覆盖一批凸分段二次正则项，C²-cone reducible conjugate 覆盖 polyhedral、SOC/PSD/spectral-like 几何。它确实支持“统一并扩展若干稳定性结果”的 claim。
但定量部分的 evaluation 有明显上限：文中没有数值案例展示 bound 的 tightness，也没有说明这些方向导数公式在实际算法中是否易于稳定计算。对于 semismooth Newton 的潜在用途，论文只指出方向，没有真正验证。
Limitation
1. 适用性依赖强结构假设。 最干净的结论依赖 r* 在 v̄ 处 C²-cone reducible。很多重要正则项满足，但这不是泛化到任意 convex regularizer。Example 5.9 也说明 twice epi-differentiable/PLQ 并不自动给出相同 kernel 公式。
2. 结论是局部的。 所有稳定性都是 around (Ā,b̄,λ̄,x̄)。它不说明正则化路径上的全局行为，不处理活动集切换的全局复杂性，也不保证远离该点时的 well-conditioning。
3. 可计算性没有完全解决。 par ∂r*(v̄) 比 ∂²r 好算很多，但对复杂 spectral/矩阵正则项或线性复合 PLQ，实际构造 basis、计算投影、枚举活动锥仍可能昂贵。论文把困难从 coderivative 转移到共轭面结构识别；这是显著简化，但不是消除计算问题。
4. 定量 bound 的实用价值未充分说明。 Corollary 5.6/5.14 给了 Lipschitz bound，但是否 tight、是否数值稳定、是否能作为算法步长或 condition estimator，文中未充分说明。
5. PLQ 复合的理论缺口仍在。 作者明确指出不清楚 θ_{P,B}∘M 的共轭是否 C²-cone reducible；因此需要额外 compactness qualification。这里的增益来自补丁式几何论证，而不是完整统一理论。
Takeaway
- 1. 以后分析非光滑正则化稳定性时，不必执着于完整 generalized Hessian；很多时候只需要其 kernel，而 kernel 可能有简单的共轭几何表达。
- 2. ker A ∩ par ∂r*(v̄)=0 是一个非常可迁移的稳定性模板：data curvature 必须覆盖 regularizer 当前 face 的自由方向。
- 这个模板可以迁移到其他 smooth convex fidelity，只需把 ker A 换成 fidelity Hessian 的 nullspace。
- 3. primal monotone generalized equation 路线值得保留。
一句话总结
这篇论文把正则化最小二乘的局部稳定性从难算的二阶 coderivative 判别，归约为共轭子微分平坦子空间与观测 nullspace 的横截性判别，是一篇以变分分析重构 LASSO 型稳定性理论接口的工作。