精读笔记
Problem Setting
《Optimal control of symmetry-breaking dynamics near criticality》（arXiv preprint / 2026）研究的是近 pitchfork 分岔的非线性系统在有限时域二次代价下的最优控制结构。它实际解决的问题不是“如何数值求一个最优控制”，而是：当系统接近临界点、中心模态变慢、稳定模态仍快衰减时，最优控制律是否像动力学 normal form 一样具有低维普适形式。
关键矛盾在于：pitchfork 的自然幅度是 c~ε^{1/2}，中心流形上的线性漂移和三次非线性漂移同为 O(ε^{3/2})。因此传统按状态多项式次数排序的非线性最优控制展开在这里失去主导项层级；LQR 只在控制足够强、能压过中心流形非线性时成立。另一方面，常规奇异摄动也不直接适用，因为哪些模态被 slaved 不只由 plant 决定，还由 tracking penalty 的尺度决定。论文的核心问题就是在这个尺度冲突中重新建立 optimality 的 reduced description。
Motivation
已有路线缺的是“最优性条件的 normal form”。经典 bifurcation control 通常通过 feedback 移动分岔点或稳定化不稳定平衡，关注 controllability / stabilizability 退化；本文处理的是相反情形：临界模态失稳，但 critical mode 仍可直接控制。障碍不是控制不存在，而是通常构造控制律的 asymptotic ordering 不再可靠。
作者的关键观察是：near criticality 下，控制强度不能被当成普通小扰动；它相对 ε^{3/2} 的大小决定了 reduced optimality problem 的类型。也就是说，optimal control near bifurcation 需要把 cost parameters 的 scaling 当成模型的一部分。这个方向自然来自一个缺口：动力学 normal form 已经告诉我们 state dynamics 的普适性，但 external cue / control 是否继承这种普适性，此前没有被系统化分析。
Core Idea
论文真正的核心思想是：不要先把 plant 降维再套标准控制器，而是直接对 PMP 给出的 state-costate Hamiltonian 系统做中心流形式的渐近展开。这样 reduced object 不是单纯的 amplitude equation，而是 amplitude-costate 的 control normal form。这个建模方式把“最优性”本身纳入 normal form，而不是把控制看作 normal form 之后的外接模块。
本质区别在于它引入了一个新的 inductive bias：控制律的结构由临界漂移、控制耦合和代价尺度之间的 distinguished balance 决定。strong regime 回到全系统 LQR；intermediate regime 变成 driftless critical-mode Riccati；weak regime 则保留 pitchfork 的非线性 amplitude dynamics，得到真正的 nonlinear Hamiltonian tracking problem。这个分类比单一 LQR 或单一中心流形控制更 general，因为它承认 near criticality 下没有唯一的 leading-order control law。
Method
第一，PMP reduction 解决的是最优控制的必要条件如何与分岔展开兼容的问题。通过 state、costate 和控制的一阶条件，控制被写成 costate 投影 u=-r_u^T p，从而控制强度直接由 costate scaling 和 penalty scaling 决定。这一步的作用是把“代价如何驱动系统”显式放进 Hamiltonian flow。
第二，三种 penalty scaling 解决的是 leading-order balance 的选择问题。χ,ζ=O(1) 时控制为 O(ε^{1/2})，稳定模态被 leading-order 激发，得到全系统 LQR；χ=O(ε^{1/2}), ζ=O(ε) 时控制为 O(ε)，稳定模态准静态 slaving，中心模态上得到 driftless scalar Riccati；χ=O(ε), ζ=O(ε^2) 时控制为 O(ε^{3/2})，刚好与 μA+νA^3 平衡，得到非线性 amplitude-costate 系统。
第三，weak regime 的 Hamiltonian / bifurcation / boundary-layer analysis 解决的是非线性 reduced BVP 仍然难解的问题。作者把 reduced Hamiltonian 变换成 canonical potential V(X;κ,Xc)，用相图刻画 branch structure，再在小 κ、长时间下把轨迹分解为 heteroclinic layers，用 front-interaction force balance 定位 transition layers，并给出 leading-order cost。这不是数值技巧，而是把最优轨迹的全局几何结构显式化。
Key Insight / Why It Works
最重要的 insight 是：near pitchfork 的最优控制不是“LQR 加非线性修正”，而是由控制强度相对中心流形内禀漂移的尺度决定的 regime problem。传统 LQR 在 strong regime 有效，是因为控制强到足以让内禀 pitchfork 非线性降为高阶；但在 weak regime，三次项与控制同阶，删掉它会改变问题类型。这一点是论文最实质的贡献。
方法成立的原因来自两个结构性事实。第一，pitchfork 有单一临界模态，稳定模态有谱隙，所以在控制不太强时稳定模态可被 slaved；这给了低维化。第二，PMP Hamiltonian 保留了可投影的 symplectic / energy structure，所以 weak-regime 的非线性最优轨迹可以用 level set 和 separatrix geometry 组织。这里不是 scaling trick alone，而是 scaling 加 latent dynamical structure 的结合。
我认为最核心的是三 regime 分类和 weak-regime amplitude-costate normal form。boundary-layer 的长时解析很漂亮，但更像对 weak regime 的深入 exploitation；它依赖小 κ、长时、constant target 和层状轨迹结构，迁移范围较窄。strong regime 的 LQR 部分则主要是 sanity check：说明在强控制下会退回已有理论，并不是主要创新。
需要注意，weak regime 的“closed-form feedback”并不完全闭合，因为能量 E 仍依赖初值和终端时间，文中未给出闭式 E(A0,T)。数值图中通过 fitting / shooting 选择 E，这意味着实际控制律仍带有全局 BVP 信息。若把它当作完全显式 feedback，会高估贡献。
Relation To Prior Work
最接近的谱系有三条：bifurcation control、Al'brekht / HJB 非线性最优控制展开、singular perturbation / turnpike optimal control。本文和 bifurcation control 的差异在于，它不主要移动分岔点或处理 controllability 退化，而是在 stabilizability retained 的情况下分析 optimality conditions 的尺度结构。和 Al'brekht 方法的差异在于，它不按状态多项式次数展开 value function，因为 near pitchfork 下线性项和三次项同阶；按次数排序会误判 leading dynamics。
和中心流形控制的关系更微妙。本文看似是“中心流形 + PMP”，但实质新增信息是 cost scaling 会决定 reduced optimality 的维度和形式；这不是普通 plant reduction 可以给出的。intermediate regime 的 scalar Riccati 是已有 LQR 思想在临界模态上的重组；weak regime 的 nonlinear Hamiltonian normal form、κ-Xc 相图和多 layer force balance 才是更实质的创新。
large deviation / Freidlin-Wentzell 的联系也很准确：ζ=0 时 reduced weak Hamiltonian 类似最小作用路径；tracking penalty 打破异宿轨道平移不变性并决定 transition timing。这里的新增不是借用大偏差形式，而是把 tracking tilt 如何组织最优控制相图讲清楚。
Dataset / Evaluation
evaluation 是单一 canonical bistable biochemical switch，二维、对称、critical mode 直接受控，并且满足作者假设。它覆盖了两类转移情形：unstable-to-stable 和 stable-to-stable，并在三种 scaling regime 下与 full numerical optimal control 解比较。这个设计适合验证渐近推导是否在一个干净模型中工作。
但它没有真正验证“general n-dimensional system”的经验 claim。高维系统、多个稳定模态、控制只通过稳定模态间接耦合、非对称扰动、非局部 target、时间变化 target 都没有测试。benchmark 支撑的是局部 normal form 的正确性，而不是工程上可泛化的控制算法。
数值解还使用 asymptotic feedback 作为 warm start，再优化 full discretized objective。作为验证可以接受，但如果要证明该方法独立改善求解效率或全局最优性，证据不够。增益来源主要是理论降维和 near-critical scaling，而不是 benchmark 上展示出的 numerical superiority。
Limitation
第一，成立前提很强：单一零特征值、稳定谱隙、pitchfork 非退化、critical mode 直接可控、初值和目标都在 O(ε^{1/2}) 邻域内。任何一个条件破坏，三 regime 结构都可能改变。尤其是间接控制 critical mode 的情形被明确排除，这在实际系统中并不少见。
第二，weak regime 的最优性没有被充分保证。Appendix A 显示二阶变分包含 6νAQ 项，符号不受控，因此一阶 PMP extremal 不一定是 minimizer。文中未充分说明如何排除 conjugate points；这对 weak-control nonlinear feedback 的控制合成意义是实质风险。
第三，所谓普适性是局部渐近普适性，不是全局控制普适性。只要轨迹离开临界邻域，normal form 系数不足以描述控制效果。scalability 上限也由此决定：高维系统可以在理论上通过少数系数降维，但前提是目标任务确实被单一中心模态支配。
第四，weak-regime 的闭式解依赖 constant target、小 κ、长 horizon 和可分解 layer sequence。对于一般 finite-time BVP，energy E 仍需数值 shooting，branch choice 也依赖全局相图。这里有一部分问题只是从 full Hamiltonian BVP 转移到了 reduced Hamiltonian BVP。
Takeaway
- 最值得记住的不是某个公式，而是 regime selection：near criticality 下，最优控制律的 universality class 由 cost scaling 相对 ε 决定。
- 以后做临界系统控制时，不能默认 LQR 是 leading-order baseline；必须先判断控制是否强到压过 center-manifold drift。
- 可迁移的 insight 是：对含临界慢模态的问题，应该 reduction optimality conditions，而不只是 reduction dynamics。
- costate 的尺度包含了任务和代价信息，它会改变 slow-fast decomposition。
一句话总结
这篇论文把 pitchfork normal form 从 plant dynamics 推进到 PMP optimality dynamics，核心贡献是用 cost scaling 划分近临界最优控制的三类普适 regime，其中 weak-control nonlinear Hamiltonian normal form 是最有价值的新增结构。