0. 链接验证表
以下 URL 均在本次报告制作时实际访问。核心论文在 PMLR 页面可直接取得正式 PDF;OpenReview 页面存在但自动访问遇到验证页,因此只确认 PMLR 页面给出的链接,不把 OpenReview 元数据作为定量依据。精确标题检索未找到核心论文的 arXiv、Hugging Face Daily Papers、Papers with Code 或 Semantic Scholar 独立条目;ACL Anthology 与本主题不对应。
| 工作 | 会议 / 出版页 | arXiv(版本与日期) | 代码 | 数据 / 权重 | 项目页 / 视频 |
|---|---|---|---|---|---|
| PemNN(核心) | ICML 2025 / PMLR 267;正式 PDF;OpenReview | 未找到(截至 2026-09-13) | 官方 GitHub | 同一仓库公开四类实验/仿真数据;仅找到 Titin 的两个 Keras checkpoint,未找到四蛋白完整权重 | 独立项目页、视频:未找到 |
| Waite 等:few-shot AFM | Patterns 2023 / PMC 全文;DOI | 未找到可确认的对应记录 | 代码随 Zenodo | Zenodo 公开曲线、标注和代码 | 视频:未找到 |
| Dudko–Hummer–Szabo 力谱理论 | PNAS 2008;PubMed | 未找到 | 未公开 / 不适用 | 未公开 / 不适用 | 未找到 |
| Dingeldein 等:SMFS 的 SBI | MLST 2023 DOI | arXiv:2209.10392v2,2022-11-14 | SBIsmfs | 代码仓库含模拟与推断流程;预训练权重未找到 | SBI 资源页;视频未找到 |
| InceptionTime | Data Mining and Knowledge Discovery 2020 | arXiv:1909.04939v3,2020-12-05 | 官方 GitHub | UCR 等公开基准由仓库脚本获取;本论文专用权重未找到 | 视频未找到 |
| FUSION | Nano Letters 2023;PubMed | 未找到 | 未找到 | 数据 / 权重未找到 | 视频未找到 |
1. 选题审计:候选池、排重与随机选择
运行日对应“AI4Science 与跨界”。索引排重排除了最近十四天已深度覆盖的 Foundation Molecular Grammar(分子生成)、TIMING(科学/医疗时序解释)与坐标鲁棒回归等主题。本轮以 ICML 2025 正式 proceedings 为主源,并由其 OpenReview 链接、arXiv 精确标题检索及聚合站检索交叉核验;这样候选的 venue 层级一致,随机数不会暗中混入“主会对预印本”的质量偏置。
候选池共二十九篇,均为 ICML 2025 主会论文,因而核心候选均在二十四个月时限内:
| # | 候选论文 | 子领域 |
|---|---|---|
| 1 | From Mechanistic Interpretability to Mechanistic Biology | 蛋白语言模型解释 |
| 2 | Training Flexible Models of Genetic Variant Effects | 基因组 |
| 3 | Raptor: Scalable Train-Free Embeddings for 3D Medical Volumes | 医学影像 |
| 4 | Visual and Domain Knowledge for Professional-level Graph-of-Thought Medical Reasoning | 医疗推理 |
| 5 | LDMol | 文本到分子扩散 |
| 6 | Hierarchical Graph Tokenization for Molecule-Language Alignment | 分子—语言 |
| 7 | Scalable Non-Equivariant 3D Molecule Generation via Rotational Alignment | 三维分子 |
| 8 | A Physics-Augmented Deep Learning Framework for Classifying Single Molecule Force Spectroscopy Data | 单分子生物物理 |
| 9 | Geometric Representation Condition Improves Equivariant Molecule Generation | 分子几何 |
| 10 | Gradient Guided Bayesian Flow Networks for Molecule Optimization | 分子优化 |
| 11 | Compositional Flows for 3D Molecule and Synthesis Pathway Co-design | 分子与合成 |
| 12 | CombiMOTS | 双靶点药物设计 |
| 13 | Elucidating the Design Space of Multimodal Protein Language Models | 蛋白语言模型 |
| 14 | PPDiff | 蛋白复合物设计 |
| 15 | Flexibility-conditioned Protein Structure Design with Flow Matching | 蛋白结构 |
| 16 | Fast and Low-Cost Genomic Foundation Models via Outlier Removal | 基因组基础模型 |
| 17 | SPACE | DNA 基础模型 |
| 18 | Near-optimal Sketchy Natural Gradients for PINNs | PINN 优化 |
| 19 | Adaptive Flow Matching for Resolving Small-Scale Physics | 物理模拟 |
| 20 | Calibrated Physics-Informed Uncertainty Quantification | 物理不确定性 |
| 21 | DISCO | 多物理算子学习 |
| 22 | Physics-Informed Weakly Supervised Learning for Interatomic Potentials | 原子间势 |
| 23 | A Generalizable Physics-Enhanced State Space Model | 动力学预测 |
| 24 | Holistic Physics Solver | PDE 求解 |
| 25 | OmniArch | 科学计算基础模型 |
| 26 | Dictionary Learning from Microscopy Foundation Models | 自动生物发现 |
| 27 | Adapting While Learning: Grounding LLMs for Scientific Problems | 科学 Agent |
| 28 | SCENT | 连续科学时空数据 |
| 29 | LLM-SRBench | 科学方程发现 |
按要求在系统 shell 执行 shuf -i 1-29 -n 1,输出
8;第八篇即
PemNN。正式论文、附录、代码、数据均足以支撑方法与定量审计,故重抽次数为零。Hugging
Face、Papers with Code、Semantic Scholar 与 ACL Anthology
的精确标题检索没有返回可确认的核心论文条目,这些“未找到”不是“论文不存在”的证明,只代表本轮检索无法核实。
2. Motivation:真正昂贵的不是预测,而是获得可信标签
单分子力谱(SMFS)以探针拉伸蛋白,记录力随 extension
变化的轨迹。科学目标常是估计单个蛋白域的解折叠力,但一次试验可能没有抓到分子、抓到一个分子,或同时抓到多个分子;后两者的轨迹形态会混杂。论文把这一步改写为三分类:no molecule / single molecule / multiple molecules。传统流程靠专家逐条看曲线,核心论文称分析时间可从“一天”降到“一小时以内”(核心论文,第
1 页摘要与第 8
页结论);这个说法没有独立计时表、样本量或硬件区间,因此只能视为工作流量级陈述,不能当作严格加速比。
痛点有三层。第一,标注少。四个实验蛋白的数据量按类别分别是 Titin
[181,164,191]、细菌 Utrophin
[181,181,178]、昆虫 Utrophin
[175,166,200]、Dystrophin
[191,185,177];每种蛋白跨五至十一个实验日、三至四种拉伸速度(核心论文,表
1,第 6
页)。默认低数据设置只使用约三十条训练曲线,论文报告四蛋白平均准确率
79.6±5.2%(核心论文,第 1、7 页,图
6)。第二,曲线长度不一,实验轨迹从 736 到
9974 点不等,而网络需要统一输入(核心论文,表 1,第 6
页)。第三,纯数据模型可能把仪器噪声、速度或某天的漂移学成捷径;但手写物理模型又不够表达真实蛋白的异质性。PemNN
的下注是:不要求物理模型独自做判别,而是把物理变换作为一个可学习支路的输入。
这个问题的量化基线来自 Waite 等的 DDRs 数据:Triplet 网络总体准确率
66.70%,三个类别分别为
73.30/63.30/63.30%,总体 F1 为
66.70%(核心论文复现表 3,第 19 页;原工作见 Patterns
全文)。核心论文没有把 DDRs 用于 PemNN,因为该数据缺少
extension,无法计算 WLC
轮廓长度;这恰好说明方法不是仅靠力向量即可落地的即插即用分类器。
3. 论文故事线、假设与权衡
论文的故事线很清晰:先构造物理 Monte Carlo 引擎,生成四个蛋白、每个蛋白六百条三分类曲线;再让网络同时看原始力和 WLC 反演的轮廓长度;最后用模拟数据做结构选择,并在真实实验上验证少样本和“纯模拟预训练”两种用法(核心论文,第 4–8 页;附录 A–E,第 13–24 页)。
它隐含四个假设。
- 低维物理坐标更可分。
如果单分子解折叠事件在轮廓长度空间中比原始力空间更稳定,那么把
Lc显式输入可降低样本复杂度。 - 物理模型可以不精确但方向正确。 WLC、FJC、FRC 三种链模型在实验与模拟上的准确率相差不超过约三个百分点(核心论文,附录 E.1,第 19–20 页,图 10);作者据此主张精确模型不是成败关键。由于论文只给曲线,没有逐项矩阵或置信区间,这里不能判断三个百分点是否统计等价。
- 仿真能提供有用的表征初始化。
仅在对应模拟数据上预训练、完全不看实验标签时,四蛋白平均准确率为
72.0±5.9%、平均 ROC-AUC 为0.88(核心论文,第 7 页与附录 E.4,第 23–24 页)。这证明“非零迁移”,不等于完成 sim-to-real。 - 早期融合比晚期决策融合更容易学习互补关系。 Early-Conv 在关键差异图中最佳;只把两个支路的最终 softmax 做晚融合更弱(核心论文,附录 E.1,第 19–20 页,图 9)。
权衡也同样明确:加入物理支路会引入
kBT、固定持久长度、轮廓长度阈值和 extension
质量等依赖;模拟器需要蛋白的能垒与解折叠参数。换来的好处是低数据区显著更稳,以及可把最终筛选后的解折叠力分布与物理常识对照。
4. 方程级形式化与流程
给定第 i 条轨迹的 extension
X^(i)=(x_1,…,x_T)、force F^(i)=(f_1,…,f_T)
与标签 y^(i)∈{0,1,2},分类器输出
p_θ(y|X,F),训练目标是标准类别交叉熵:
L_CE(θ) = -(1/N) Σ_i Σ_{c=0}^2 1[y^(i)=c] log p_θ(c | X^(i),F^(i)).
这是核心论文式(1),第 3 页。真正的物理增强发生在输入变换。WLC
模型把每个正的 (x,f) 对映射到轮廓长度 Lc:
f = (k_B T / L_p) · [ 1/{4(1-x/L_c)^2} - 1/4 + x/L_c ].
这是核心论文式(2),第 4 页。作者固定
L_p=300 pm,逐点数值求解
L_c,保留满足物理约束且 L_c≤1000 nm 的点,再把
force 与 contour-length 两条序列分别采样到四百点(核心论文,第 4–5
页;附录 D,第 18 页)。注意:Lc 不是只由 F
决定,它由 (X,F) 共同得到;如果通用基线只看
F,那么 PemNN 同时获得了额外传感变量和物理坐标变换。
默认网络可以写成:
h_F = Branch_F(Resample400(Normalize(F)))
L_c = SolveWLC(X,F; L_p, L_c,max)
h_L = Branch_L(Resample400(Normalize(L_c)))
h = Fuse_early(h_F,h_L)
p = Softmax(Head(h))
两个支路先各过一维卷积、批归一化和 ReLU;Early-Conv
随后拼接特征,再过两个卷积块与全局平均池化。另一条实现分支在每路使用隐藏维度八、dropout
0.8 的 LSTM 后拼接(核心论文,第 5 页与附录 D,第 18
页)。完整训练/部署流程如下:
┌────────── Monte Carlo simulator ──────────┐
蛋白物理参数 ────>│ WLC/FJC/FRC + DHS off-rate + 随机解折叠 │──> 模拟三类轨迹
└────────────────────────────────────────────┘ │
真实 AFM: (X,F) ───────────────────────────────────────────────────────────┤
v
┌──────────── force branch ─────────────┐ 监督交叉熵
(X,F) ─────────>│ F → 归一化/400 点 → Conv/LSTM │──┐ │
└───────────────────────────────────────┘ ├─融合─>三分类
┌──────────── physics branch ───────────┐ │
│ (X,F)→WLC 反演 Lc→过滤/400点→Conv/LSTM│──┘
└───────────────────────────────────────┘
│
└─ single 类轨迹 → 解折叠力分布
模拟器的关键不是画出“像真的”曲线,而是把链伸长、cantilever
平衡、Dudko–Hummer–Szabo(DHS)力依赖脱离率和随机解折叠串起来。附录算法
1 每个时间步增加 base extension,求力平衡,计算 DHS off-rate
与事件概率,再随机决定 domain unfolding 或 molecule
detachment(核心论文,附录 A,第 13–14 页)。表 2 给出的 Titin、细菌
Utrophin、Dystrophin、昆虫 Utrophin 的 ln k0 分别为
-9.21/-6.50/-4.50/-2.50,Δx‡ 为
0.40/0.85/0.60/0.41 nm,ΔG‡ 为
20/14.5/11.5/9.8 kBT(核心论文,表 2,第 15
页)。这些数直接决定仿真域,因而也是迁移误差的来源。
5. 关键机制与工程细节
| 部件 | 论文实现 | 科学作用 | 工程风险 / 复现备注 |
|---|---|---|---|
| Force 支路 | force 归一化并重采样至 400 点 | 保留仪器观测的形态 | 不等长轨迹的重采样可能抹去窄峰;论文未给插值核消融 |
| Physics 支路 | WLC 逐点反演
Lc;Lp=300 pm;Lc≤1000 nm |
把伸长—受力关系投影到聚合物坐标 | 需要 extension;阈值和固定 Lp
带来模型偏差(核心论文,附录 D,第 18 页) |
| 卷积骨干 | filter 128/256/128,kernel 8/5/3 |
多尺度局部模式 | 参数量、FLOPs、推理延迟均未报告(核心论文,附录 D,第 18 页) |
| LSTM 变体 | 每支路 hidden 8,dropout 0.8 |
建模顺序依赖 | 高 dropout 与 tiny hidden 的敏感性未报告(同上) |
| 优化 | Adam,lr 0.001,β₁ 0.9,β₂
0.999,ε 1e-8;batch 16;200
epochs |
统一训练协议 | baseline 与 PemNN 未做等预算调参;best checkpoint 监控训练 loss(同上) |
| 选模 | validation_split=0.001,checkpoint 监控
loss |
实现上近似以训练损失选 epoch | 公开 API 中 split 的 seed 参数被注释;示例脚本另行设
seed,路径不一致(官方仓库代码审计,commit 7875181) |
| 硬件 | Apple M1 Pro,10-core CPU / 16-core GPU | 表明可在消费硬件训练 | 没有 wall-clock、能耗、方差或并行设置(核心论文,附录 D,第 18 页) |
公开仓库可复现到什么程度?数据目录约含四种实验/模拟数据,且提供训练调用脚本;CallScripts/ML_models/Titin
下可见两个约 8.39 MB 的 Keras
checkpoint,但未找到四种蛋白的完整权重集合。本报告读取的最新仓库提交为
7875181(2025-06-05)。因此“数据与代码公开”成立,“全部报告结果可直接用预训练权重复核”不成立。
6. 数据集、划分、指标
| 数据集 | 三类曲线数 [无/单/多] |
原始长度 | 天数 | 拉伸速度 nm/s | 用途 |
|---|---|---|---|---|---|
| DDRs | [102,102,136] |
400 | 1 | 2000 | 仅通用 baseline;缺 extension,PemNN 不适用 |
| Titin I27O | [181,164,191] |
736–4859 | 5 | 500/1000/2000 | 实验验证 |
| Bact UtrN-R3 | [181,181,178] |
957–9974 | 10 | 500/1000/2000/5000 | 实验验证 |
| Insect UtrN-R3 | [175,166,200] |
1777–8890 | 11 | 500/1000/2000 | 实验验证 |
| DysN-R3 | [191,185,177] |
1659–4814 | 6 | 500/1000/2000 | 实验验证 |
来源均为核心论文表 1,第 6 页。另有四个模拟数据集,每个六百条、三类均衡,用于结构比较与预训练(核心论文,第 6 页)。默认实验采用分层随机划分并重复五次;低数据试验固定百分之八十为测试集,再逐步增加训练子集(核心论文,第 6–7 页与附录 E.2,第 21 页)。论文没有报告按实验日留一、按拉伸速度留一或跨蛋白留一的结果。
主指标是 accuracy;类别层面还给
ROC-AUC。科学下游用“single”预测筛出轨迹,然后比较解折叠力分布的 mode 与
IQR。Accuracy 在三类大致平衡的表 1 数据中可读,但部署时 single
可能是稀有事件,故还需要 precision、recall、PR-AUC
与每百条漏检数。论文给出某阈值规则可达平均 precision
0.93,却没有给对应 recall(核心论文,附录 E.4,第 24 页,式
19 后正文)。
7. 主结果:完整矩阵、复算与边界
7.1 实验数据不同标注预算:六方法 × 四设置
这是论文唯一同时给出六方法逐项均值与标准差的主矩阵。括号为五次运行标准差,accuracy
越高越好;粗体为逐列最高。表头的 5/80
表示固定八成测试、只用约百分之五训练,而后续 40/60
等是常规互补划分,二者不是连续且同定义的比例轴(核心论文,表 4,第 21
页)。
| 方法 | 5/80 | 40/60 | 60/40 | 80/20 |
|---|---|---|---|---|
| PemNN | 79.61 (5.24) | 86.82 (4.86) | 88.50 (5.09) | 88.54 (5.68) |
| LSTMFCN | 68.26 (12.12) | 85.12 (7.49) | 85.51 (6.87) | 87.66 (7.66) |
| ResNet | 37.47 (6.33) | 81.01 (9.13) | 84.62 (7.13) | 84.47 (9.26) |
| FCN | 35.88 (7.37) | 75.50 (8.92) | 75.48 (13.15) | 79.43 (10.18) |
| Triplet | 55.58 (10.17) | 66.78 (12.06) | 69.15 (12.26) | 66.92 (11.45) |
| InceptionTime | 35.32 (3.97) | 80.78 (12.78) | 81.35 (9.07) | 84.76 (11.58) |
来源:核心论文表 4,第 21 页。相对逐列最强非 PemNN 方法,本报告复算如下:
| 设置 | 最强竞品 | PemNN | 绝对提升 | 相对提升 |
|---|---|---|---|---|
| 5/80 | 68.26 | 79.61 | +11.35 pp | +16.63% |
| 40/60 | 85.12 | 86.82 | +1.70 pp | +2.00% |
| 60/40 | 85.51 | 88.50 | +2.99 pp | +3.50% |
| 80/20 | 87.66 | 88.54 | +0.88 pp | +1.00% |
复算公式是 (PemNN−竞品)/竞品;原数来自表
4。四列简单平均为 PemNN 85.87%、LSTMFCN
81.64%,差
4.23 pp。但四列共享数据且训练规模不同,简单平均不是独立样本上的总体效应量。
这组结果最值得看的不是“全部第一”,而是增益随标注量快速收缩:约三十条训练曲线时为
11.35 pp,充分标注时只剩
0.88 pp。它更支持“物理坐标降低低数据样本复杂度”,不支持“在所有预算下都有很大上限优势”。
7.2 DDRs:PemNN 不适用时的 baseline 矩阵
| 方法 | Overall accuracy | No molecule | Single | Multiple | F1 |
|---|---|---|---|---|---|
| Triplet | 66.70 | 73.30 | 63.30 | 63.30 | 66.70 |
| FCN | 77.10 (5.15) | 81.22 (4.96) | 77.49 (13.59) | 73.48 (10.06) | 77.30 (5.32) |
| ResNet | 79.03 (4.84) | 81.90 (10.52) | 83.97 (4.74) | 70.92 (7.44) | 79.02 (4.88) |
| LSTMFCN | 75.81 (3.42) | 70.91 (2.49) | 88.24 (4.16) | 71.30 (9.02) | 76.17 (3.32) |
| InceptionTime | 85.48 (1.14) | 97.14 (2.61) | 84.00 (8.94) | 72.50 (12.96) | 85.26 (1.02) |
| PemNN | 不适用 | 不适用 | 不适用 | 不适用 | 不适用 |
来源:核心论文表 3,第 19 页。这里不能补出 PemNN:DDRs 没有 extension。矩阵说明通用 InceptionTime 在可用信息相同时很强,也提醒部署方先审计采集格式。
7.3 科学下游:筛选后解折叠力是否更“像物理”
| 蛋白 | Raw mode (IQR) pN | Heuristic mode (IQR) pN | PemNN mode (IQR) pN | IQR 相对 Raw 缩减 | IQR 相对 Heuristic 缩减 |
|---|---|---|---|---|---|
| Titin | 192.21 (234.52) | 217.41 (270.46) | 206.68 (52.63) | 77.56% | 80.54% |
| Bact UtrN | 60.50 (63.22) | 72.84 (66.55) | 70.40 (36.97) | 41.52% | 44.45% |
| Insect UtrN | 81.41 (85.29) | 90.58 (80.57) | 84.66 (39.98) | 53.12% | 50.38% |
| DysN | 63.76 (53.18) | 74.38 (53.10) | 72.54 (26.79) | 49.62% | 49.55% |
前三列来自核心论文表 5,第 21 页;后两列为本报告按
(baseline IQR−PemNN IQR)/baseline IQR 复算。四种蛋白的 IQR
都明显收窄,这是比 mode
接近更强的证据。但“更窄”不自动等于“更真”:筛选器也可能过度删除真实异质构象。
正文称 PemNN 与 heuristic 的 mode 差“不超过 10 pN”,但表 5 的 Titin
是 |206.68−217.41|=10.73 pN;其余三种为
2.44/5.92/1.84 pN。因此表格事实与文字概括存在轻微不一致。Titin
的文献期望区间为 204±26 pN,Raw、Heuristic、PemNN 三个 mode
都在该范围内(核心论文,第 7 页);主要增量来自 IQR,而不是把 mode
从错误位置纠正到正确位置。
7.4 sim-to-real:纯模拟训练与实验训练的类别 AUC
| 蛋白 | 模拟训练:No / Single / Multiple | 实验训练:No / Single / Multiple | Single 类绝对差 |
|---|---|---|---|
| Titin | .98 / .86 / .90 | 1.00 / .95 / .96 | +.09 |
| Insect UtrN | .92 / .61 / .75 | .99 / .81 / .84 | +.20 |
| Bact UtrN | .97 / .67 / .82 | 1.00 / .92 / .93 | +.25 |
| DysN | .97 / .61 / .77 | 1.00 / .88 / .90 | +.27 |
来源:核心论文表 7,第 24 页。十二个类别 AUC 的简单平均从模拟训练
0.819 提至实验训练 0.932,本报告复算差
0.113。最难且最科学相关的 Single 类在三个异质蛋白上仍差
0.20–0.27;因此 72.0±5.9%
的纯模拟准确率应理解为一个有用冷启动,而不是可替代专家标签的终点。
论文还比较了模拟与实验筛选后的 mode/IQR:Titin
216.35(50.17) 对 217.75(36.17),Bact
81.58(64.07) 对 85.38(22.28),Insect
89.45(70.86) 对 96.65(40.20),Dys
91.25(65.64) 对 79.75(25.09)(核心论文表 6,第
23 页)。四组 mode 大致接近,但模拟训练的 IQR
每组更宽,继续表明分类边界未完全跨域。
物理支路的优势集中在极低标注预算
纵轴为实验数据平均准确率(%);蓝色为每列最强非 PemNN 方法,橙色为 PemNN。比例定义以论文表 4 为准。
8. 逐层消融
论文的消融分成三层,但只有部分精确数字。
支路层。 论文把双支路 PemNN 与只使用 force、只使用
physics 输入,以及让 baseline 接收 physics 输入的变体比较。图 5a
报告四个模拟数据集平均 accuracy 为 98.9±1.9%;相对最佳
baseline,Titin 的最大差约十五个百分点、Insect UtrN
的最小差约六个百分点(核心论文,第 6 页,图 5a
正文)。图中没有逐数据集×逐方法数表,所以不能构造五基准×六方法的完整模拟矩阵。
融合层。 Early-Conv、Early-LSTM
及晚期融合被比较;Early-Conv 排名最佳,双 LSTM 变体的平均 rank 为
2.4167(核心论文,附录 E.1,第 19–20 页,图
9)。论文没有给每种融合的准确率均值、标准差或 pairwise p
值,所以只能说“排序支持早融合”,不能量化每个部件贡献。
物理模型层。 WLC、FJC、FRC 在模拟与实验曲线上相差不超过约三个百分点(核心论文,附录 E.1,第 19–20 页,图 10)。这提示“合适的物理坐标族”可能比“唯一正确的链模型”重要;但无等价性检验,不能声称三者等价。
缺少的关键消融是:F-only、F+raw X、F+shuffled Lc、F+learned monotonic transform(X,F)、F+WLC Lc
五路 head-to-head。因为当前 physics 分支引入了基线没有的
X,所以现有消融无法回答增益来自额外观测、非线性变换,还是物理正确性。
9. 超参敏感性与“相变”
轮廓长度阈值从 1,000 到 1,024,000 nm
按四倍递增,并加入“不过滤”。1,000–64,000 nm
区间接近满分;到 1,024,000 nm 时 PemNN 仍至少领先 baseline
七个百分点;完全取消阈值时差距约二十五个百分点(核心论文,第 6 页,图
5b)。这不是数学相变证明,但呈现了明显机制:阈值放宽后 physics
支路纳入大量非物理解,通用模型退化更快,而双支路仍能依赖 force
支路兜底。
持久长度 Lp 在 0.3–30,000 pm
上按对数尺度扫描。30–30,000 pm 区间双方接近满分;在
0.3 pm 和 3 pm 的极端失配下,PemNN
分别至少领先约三和十个百分点(核心论文,第 6 页,图
5c)。这支持“支路冗余带来容错”,而非“WLC 参数必须准确”。
论文没有报告学习率、batch、dropout、LSTM hidden size、重采样长度、训练 epoch 或模拟量的敏感性;也没有性能—算力曲线。因此能推出 physics 输入对两个物理参数较鲁棒,不能推出整体训练配方稳健,更不能推出其他仪器或蛋白无需调参。
10. 统计显著性与不确定性
论文以 Friedman 检验后接 Wilcoxon signed-rank 与 Holm 校正,并在
critical-difference 图中给 PemNN 平均 rank
1.4167,称其显著优于五个 baseline(核心论文,第 6 页,图
4)。这是合理的多数据集比较框架;但论文没有报告 Friedman
统计量、自由度、原始/校正 p
值、效应量或置信区间。只有八个任务(四模拟加四实验)的 rank
统计也容易受任务相关性影响。
表 4 给五次运行的标准差,却没有说明同一 split
还是每次重采样,也没有给 paired seed 结果。5/80 下 PemNN
的标准差 5.24 小于 LSTMFCN 的
12.12(核心论文表 4,第 21
页),说明低数据时不仅均值更高且波动更小;但没有每次成绩,无法做配对检验或计算可信区间。表
5 的 mode/IQR 是分布描述,不是分类性能的显著性检验。
结论边界因此是:“均值与 rank 证据一致地支持 PemNN”成立;“优势在每个新实验域都有统计保证”不成立。
11. 真正贡献、可复用设计与迁移潜力
真正贡献不是发明 WLC,也不是发明双支路,而是把三件事闭环:可生成标签的物理模拟器、物理坐标与原始观测的容错融合、以及以解折叠力分布验证科学下游。对 AI4Science 最可复用的设计是“可错的 physics branch + 可兜底的 data branch”:不要强迫网络满足一个可能失配的方程,而是让方程产生另一种表示,再由融合层决定信任度。
第二个可复用点是把 sim-to-real
拆成两个部署档位:零实验标签时用模拟预训练得到 72.0±5.9%
的冷启动;有少量标签时用约三十条曲线把平均准确率推到
79.6±5.2%(核心论文,第 1、7 页)。这比笼统地说“synthetic
data works”更有决策价值。
第三个可复用点是用物理下游量做外部检查。IQR 相对 raw 缩减
41.52%–77.56%(本报告据表 5
复算),说明分类器确实改变了科学统计量;然而这种检查必须同时看
mode、IQR、已知区间和独立实验,避免把“分布更窄”误当成“更准确”。
迁移到其他领域时,可考虑:质谱中的同位素/碎裂规则支路,电化学中的等效电路参数支路,医学波形中的生理方程支路,遥感中的辐射传输坐标支路。必要条件是能从原始观测稳定计算物理表征,并设置 raw branch 防止模型失配;如果物理表征使用了额外传感器,比较基线必须获得同等信息。
12. 相关工作链:问题从哪里来,设计如何传递
12.1 Dudko–Hummer–Szabo:把受力下的逃逸率写进模拟器
解决什么。 DHS
理论把外力作用下的能垒越迁与断裂/解折叠速率连接,使动态力谱可以反演本征速率、能垒位置与高度(PNAS
2008)。 遗留缺口。
它解释单事件动力学,不直接判断一条复杂 AFM
轨迹含零、单还是多个分子;真实仪器还有
cantilever、噪声与异质路径。
核心如何回应。 PemNN 的 Monte Carlo 引擎在每个
extension 步用 DHS off-rate
采样解折叠/脱离事件,再产生三类带标签曲线(核心论文附录 A,第 13–15
页)。
设计如何传递。
传递的不是一个固定分类阈值,而是可生成数据的机制模型:动力学理论成为训练分布的生成器。
12.2 Waite 等:把人工看曲线改成 few-shot 分类
解决什么。 Waite 等把 AFM
曲线字符化为三分类,用少样本 Triplet 学习,公开 DDRs
数据与标注;核心论文复现的总体 accuracy 为
66.70%(核心论文表 3,第 19 页;Patterns
2023)。 遗留缺口。 数据只有 force
序列、单日单速度,且学习器没有显式聚合物物理;少样本性能仍有限。
核心如何回应。 PemNN
收集四蛋白、多日、多速度数据,并增加 WLC 支路和模拟预训练。
设计如何传递。 三分类任务定义、公开数据和 Triplet
baseline
被保留;“few-shot”从纯表征学习转向“物理坐标降低样本复杂度”。
12.3 InceptionTime:强通用时序骨干成为压力测试
解决什么。 InceptionTime 以多尺度卷积和 ensemble
兼顾时序分类准确性与可扩展性;原论文举例,HIVE-COTE 在
N=1500,T=46 上超过八天,而 InceptionTime
约一小时,并可在十三小时处理八百万条序列(arXiv:1909.04939v3,摘要)。
遗留缺口。 它不利用 SMFS 的 extension—force
关系,也不提供 sim-to-real 接口。
核心如何回应。 PemNN
让通用卷积负责观测形态,让物理支路提供任务特定坐标;DDRs 上
InceptionTime 仍以 85.48±1.14% 第一,说明通用骨干不是弱
strawman(核心论文表 3,第 19 页)。
设计如何传递。
多尺度局部卷积成为每条支路的基础元件,但不再独自承担物理不变量学习。
12.4 Dingeldein 等:从分类走向模拟器后验
解决什么。 Simulation-based inference
把机制模型编码进
simulator,用神经密度估计对隐藏分子参数做后验推断,绕开非 Markov
轨迹的不可解 likelihood(arXiv:2209.10392v2)。
遗留缺口。
目标是参数推断而非零/单/多分子筛选;模拟器偏差仍会进入 posterior。
核心如何回应。 PemNN 选择更直接的分类目标,并用 raw
branch 减少对模拟器完全正确的依赖。
设计如何传递。 “物理模拟器 + 神经网络”的范式保留,但从
amortized posterior 变成 synthetic pretraining 与显式 physics
feature。
12.5 FUSION:把人类留在罕见路径发现的闭环里
解决什么。 FUSION
将轨迹编码为图像,通过模型与专家迭代标注来自动识别单分子路径,并可发现初始标签中没有的路径(Nano
Letters 2023)。 遗留缺口。
人在环降低冷启动风险,但扩展到多蛋白仍需要专家迭代,物理坐标未显式进入分类器。
核心如何回应。 PemNN 用仿真生成三类标签,并以 WLC
支路减少少样本负担。
设计如何传递。 两者形成互补:PemNN
提供物理冷启动,FUSION
提供发现新路径后的人工纠偏;最有潜力的系统应把二者串联,而不是二选一。
12.6 非线性技术演进树
┌─ WLC/FJC/FRC:链弹性坐标 ───────────┐
SMFS 物理 ─ DHS 动力学 ─┤ ├─> 物理 Monte Carlo 标签
└─ cantilever / loading protocol ─────┘ │
├─> PemNN
时序学习 ─ FCN/ResNet ─ InceptionTime ───────── raw-force branch ───────────┤
│ │
└─ metric few-shot ─ Waite/Triplet ─ 三分类任务与公开 benchmark ────────┘
机制生成 ─ simulation-based inference ─ 参数 posterior
╲ ╱
╲──────────── sim-to-real 校准空白 ───╱
│
人在环发现 ─────────── FUSION ──> 未知路径 / 异常类 ──> 下一代开放集 PemNN
演进动力有三个转折:从手工阈值转到少样本分类;从纯数据序列转到显式物理坐标;从“模拟器必须拟合真实世界”转到“模拟器提供可错但有用的辅助视角”。下一个转折应是把 simulator discrepancy、实验日漂移和未知类别一起建模。
12.7 跨论文对标
| 工作 | 目标 | 物理进入位置 | 标签需求 | 不确定性 | 开放集能力 | PemNN 吸收 / 未吸收 |
|---|---|---|---|---|---|---|
| DHS | 动力学解释 | 解析 off-rate | 不需要分类标签 | 参数拟合不确定性 | 否 | 吸收进模拟器;未传播参数不确定性 |
| Waite/Triplet | 三类筛选 | 无 | 少量实验标签 | 多次运行 | 否 | 吸收任务定义;扩为多蛋白/多速度 |
| InceptionTime | 通用时序分类 | 无 | 监督标签 | ensemble | 否 | 作为强 baseline;未采用 ensemble 校准 |
| SBI-smFS | 隐参数推断 | simulator + posterior | 模拟数据与观测 | 显式 posterior | 非分类 | 吸收 simulator 思路;未吸收后验校准 |
| FUSION | 路径识别与发现 | 隐式专家知识 | 人在环 | 迭代审阅 | 部分支持 | 未吸收未知路径发现 |
| PemNN | 零/单/多分类 | simulator + WLC feature | 模拟或少量实验标签 | run std / ROC | 否 | 双支路容错,但域外与未知类仍空缺 |
13. 论文自述限制
- 模拟参数难获得。 作者明确指出 Monte Carlo 引擎需要能量景观参数,而这些参数并非所有蛋白都可得;未来需要在多拉伸速度与恒力实验上更完整地约束模拟器(核心论文,Limitations,第 25 页)。
- 同质域假设。 当前分类把同一蛋白/设置的训练与测试曲线视为同域;作者把异质域与 domain adaptation 列为未来方向(核心论文,Limitations,第 25 页)。
- sim-to-real 对异质蛋白不足。 附录指出 Utrophin 与
Dystrophin 的构象/实验域异质性更强,实验标签让 Single 类 AUC 至少增加
0.20(核心论文,附录 E.4,第 24 页,表 7)。这不是独立新增限制,但给第二点提供了定量证据。
14. 独立批判:因为 X,所以 Y
- [因果归因] 因为 physics 支路的
Lc由(X,F)共同计算,而通用 baseline 默认只接收F,所以 PemNN 的+11.35 pp低数据优势不能唯一归因于“物理”;它也可能来自额外 extension 通道。必须加入F+raw X与F+shuffled Lc对照。 - [数据泄漏 / 外推] 因为实验数据跨五至十一个日期、三至四种速度,却只报告曲线级分层随机划分(核心论文表 1,第 6 页),所以同一天的仪器漂移、同一速度的形态可能同时进入训练和测试;随机 split 可能高估跨日部署性能。
- [模型选择] 因为附录称 best model 按最低 training
loss 选择,代码 checkpoint 也监控
loss,且 API 使用validation_split=0.001(核心论文附录 D,第 18 页;官方代码),所以两百 epoch 结果可能偏向训练集最优,而非泛化最优。 - [复现性] 因为公开 API 的
train_test_split中random_state被注释,而示例训练脚本另行固定 seed,所以 API 用户与论文脚本的划分协议不同;即使同一数据和超参,也不能保证复现表 4 的五次均值。 - [统计] 因为论文宣称 Holm 校正后显著,却没有报告 p 值、统计量或效应量,且 rank 只来自四模拟与四实验任务(核心论文图 4,第 6 页),所以读者无法审计显著性是否由相关任务重复计数驱动。
- [sim-to-real] 因为十二个类别 AUC 的均值从模拟训练
0.819到实验训练0.932,且三个异质蛋白 Single 类差0.20–0.27(本报告据表 7 复算),所以模拟器适合冷启动但尚不足以取代实验标签;“无需专家输入”应限定在平均 accuracy72.0±5.9%的档位。 - [部署指标] 因为阈值
α=10时论文只给平均 precision0.93、未给 recall(核心论文式 19 后,第 24 页),所以无法计算每百条真实单分子曲线会漏掉多少;对稀有事件筛选,accuracy 与 precision 不能单独决定可用性。 - [科学有效性] 因为表 5 的主要改善是 IQR 收窄,而三个
Titin mode 都落在
204±26 pN期望区间,所以模型可能是在提高纯度,也可能丢弃了真实异质路径;没有独立金标准或盲审,不能把窄分布等同于更真实。 - [基线公平性] 因为 baseline 基本沿用 sktime 配方、统一两百 epoch,但 PemNN 的结构和输入专为 SMFS 设计,论文未报告等预算超参搜索,所以“架构优于五个 baseline”混合了任务先验与调参投入差异。
- [内部一致性] 因为正文称 mode 差不超过
10 pN,表 5 的 Titin 差却是10.73 pN,所以叙述存在舍入越界;虽然不推翻总体结论,却说明科学量化应以表格为准。
15. 开放挑战:可直接执行的 head-to-head 实验
理论
- 信息增益还是物理归纳偏置? 在完全相同
backbone、参数量和训练预算下比较
F、F+X、F+Lc、F+随机置换Lc、F+可学习单调变换(X,F);用按日留一准确率、Single PR-AUC 与校准误差判断。若F+X≈F+Lc,贡献是额外信息;若只有真实Lc稳定领先,才支持物理坐标解释。 - 模拟参数不确定性如何传播? 比较固定 DHS/WLC
参数、从先验采样的 domain randomization、SBI
后验采样三种预训练;在四蛋白上报告 sim-only AUC
与加入五、十五、三十条实验标签后的学习曲线。目标是检验后验覆盖能否缩小表
7 的
0.113平均 AUC 差。 - 双支路何时会负迁移? 人为把
Lp、能垒和 cantilever stiffness 从真实值偏移多个数量级,对比 hard-physics、gated fusion、uncertainty-weighted fusion;定位 physics branch 从帮助到伤害的临界域,而不是只看图 5 的两个参数。
工程与应用
- 真实部署划分。 对每个蛋白做 leave-one-day-out、leave-one-speed-out、leave-one-instrument-out,与论文随机 split 同预算比较;同时报告 wall-clock、人工复核分钟数、每百条漏检 single 数和 p95 延迟。这个实验能直接检验“一天到一小时以内”是否端到端成立。
- 公平 baseline。 给 InceptionTime、LSTMFCN 与
Transformer 同样的
(F,X)两通道、相同超参搜索次数和参数/FLOPs上限,与 PemNN head-to-head;若 PemNN 仍领先,才可排除输入与预算不公平。 - 未知类别与主动复核。 加入粘附伪影、部分折叠、断裂和新蛋白作为 unknown,比较 closed-set softmax、energy-based OOD、deep ensemble 与 FUSION 式人在环;指标用 risk–coverage、unknown AUROC 和专家复核量。
- 完整复现包。 冻结数据 split manifest、全部四蛋白 checkpoint、容器与随机种子;由独立团队在 Linux GPU 和 Apple Silicon 上复现表 4。成功标准不是单次相近,而是均值落入原报一个标准差内。
新兴方向
- 物理基础模型预训练。 比较“每蛋白单独模拟器预训练”与“跨蛋白 masked modeling + physics tokens”,再对未见蛋白做五样本适配;检验共享链物理能否转化为跨蛋白 foundation representation。
- 分类与参数后验联合。 把三分类头与 SBI 后验头联合训练,比较仅分类、仅参数推断、联合多任务;除 accuracy 外评 posterior coverage 与解折叠力 credible interval,避免只有点预测。
- 闭环实验设计。 用不确定性主动选择下一条要人工标注或下一种拉伸速度,和随机标注、entropy sampling、纯 FUSION 人在环比较;横轴用专家分钟数,纵轴用 leave-day-out Single PR-AUC,评价真正的数据效率。
16. 近期工作雷达
- Synthetic data-driven deep learning for label-free autonomous atomic force microscopy(Nature Communications 2026):用合成 AFM 数据支持无标签、自治成像工作流,代表从“离线筛曲线”走向“模型控制下一次测量”;应重点看跨样品泛化是否比 PemNN 的同蛋白随机划分更严格。
- Automating the Expert
Eye: A System-Agnostic Deep Learning Framework for Rare Event Discovery
in Imbalanced Force
Spectroscopy(arXiv:2606.09541,2026-06-08):在目标事件占
`13/970=1.34%` 时报告 accuracy `0.9196`、recall
0.9231,并可丢弃八百八十条非目标轨迹(论文摘要/结果);它正面处理了 PemNN 平衡三分类没有覆盖的稀有事件部署。 - Atomic Force Microscopy Based Single Molecule Recognition: Deciphering Biological Interactions at the Nanoscale(Journal of Molecular Recognition 2026):从探针功能化、力谱到识别任务梳理误差来源,适合作为判断“分类更准是否等于生物结论更准”的实验背景。
- Automated atomic force microscopy analysis using convolutional neural networks(2025):把 CNN 用于 AFM 分析自动化,体现通用视觉/信号骨干仍是强路线;与 PemNN 对比时应统一输入传感量,而不是只比较模型名字。
- Machine learning approaches for improving atomic force microscopy instrumentation and data analytics(Frontiers in Physics 2024):覆盖仪器控制、图像和力曲线分析,强调真正收益来自采集—分析闭环;PemNN 目前只完成闭环中的离线过滤段。
- Deep Learning Regressors of Surface Properties from Atomic Force Microscopy(Applied Sciences 2024):把 AFM 从分类推进到连续表面属性回归;对 PemNN 的启发是将解折叠参数与类别联合预测,而不是只输出 hard label。
- Machine Learning for Analyzing AFM Images Generated from Polymer Blends(arXiv:2409.11438v2,2024-09-20):在一百四十四张图像上比较频域特征与 ResNet,提醒 AI4Science 中手工变换有时比大 backbone 更数据高效;这与 PemNN 的物理坐标支路形成方法论呼应(该文摘要与版本记录)。
17. 结论与证据边界
PemNN 最扎实的发现是低标签区的样本效率:表 4 中约三十条训练曲线时为
79.61±5.24%,较最强 LSTMFCN 高
11.35 pp;当训练比例到八成,差距缩为
0.88 pp。这使它值得被理解为一个“物理先验在少样本时最值钱”的案例,而非新的通用时序
SOTA。
最关键的未闭环问题是归因与域外性。Physics branch 同时带入
extension,随机曲线 split 又没有隔离实验日和速度;sim-only 的 Single AUC
在三个异质蛋白上仍落后实验训练
0.20–0.27。所以当前证据足以支持“物理双视角是有效工程设计”,不足以支持“WLC
物理本身导致增益”或“可无专家部署”。
报告脚注。 候选池、随机命令、URL 核验、正式 PDF
与公开代码均在本轮实际检查;精确数字只来自已读论文表/正文或标明的算术复算。ml-report-index.md
已追加本报告条目。HTML 由 Warm Editorial
模板生成,包含浅/深主题与移动端响应规则,并通过模板静态验证器(零警告);当前环境有
Playwright 包,但
/usr/bin/chromium、/usr/bin/chromium-browser
与 /opt/pw-browsers 均无可执行
Chromium,故无法完成真实浏览器的浅色、深色和移动端截图检查,不以静态校验冒充该项。