Review & Revision Report · 有限元分析审查
头条 −42.34% 站不住，收敛值约 −44.5%；高渗组 h 收敛非单调；UQ 全部跑在 710 顶点粗网格。会被直接拒。
141 个点不翻号不是稳健性证据；−42% 里约 3/4 来自一个纯假设的预张力。"epistemic robustness" 论证力接近 0。
没做均质核对照，异质染色质输入可能对头条结论完全不起作用；预测不可证伪，方向检验无区分力。
做得好的地方（要保留，是这份稿子最强的部分）：
但有三个层级的问题，其中第一层是会被直接拒稿的（见页首三卡）。
定位判断：这是一个纯计算、通用几何、零实验、结论方向由构造决定的工作。以现状投 Nature 主刊会在编辑部台面被拒（desk reject），不是因为质量差，而是因为不存在"意外发现"。合适的去处是 Biophysical Journal / Biomech Model Mechanobiol / PLoS Comput Biol / J Mech Behav Biomed Mater。如果坚持冲高影响因子，必须先把 Flipper 实验做完并入正文，且必须有一个"模型预测了非平凡的事、实验验证了"的环节（见 §3.1）。
把 Table S2 的五个离散化各自算了组间变化：
| 离散化 | h (µm) | dofs | Control | Hyper | 变化 |
|---|---|---|---|---|---|
| coarse P1 | 1.25 | 2 130 | 0.060352 | 0.037475 | −37.91% |
| inter P1 | 0.75 | 6 771 | 0.065003 | 0.038412 | −40.91% |
| main P1 ← 现在报的 | 0.45 | 22 593 | 0.066064 | 0.038092 | −42.34% |
| coarse P2 | 1.25 | 14 646 | 0.062100 | 0.034376 | −44.64% |
| inter P2 | 0.75 | 48 909 | 0.067416 | 0.037408 | −44.51% |
两族的行为：
也就是说：coarse_h1p25_P2（14 646 dof）比 main_h0p45_P1（22 593 dof）又便宜又准。现在把最不收敛的 P1 结果当头条，是把结论往小里报了 2 个百分点，还给审稿人送了一把刀。
把 production 结果换成 P2。最低限度用 intermediate_h0p75_P2；有算力就补 main_h0p45_P2。头条数字改为 −44.5%，或更稳妥地写 −44%（跨已测离散化的收敛族范围 −44.5 ~ −44.6%）。同时在 Table S2 里加一列"组间变化 %"，把五行 −37.9 ~ −44.6 全部亮出来——审稿人自己会算，你不如自己先算。
Control：0.060352 → 0.065003 → 0.066064，单调，增量比 4.38，收敛阶 p≈2.89，GCI≈0.59%。这条没问题。
Hyperosmotic：0.037475 → 0.038412 → 0.038092。增量符号翻转（+0.000937 → −0.000320）。这是振荡收敛，Richardson/GCI 无法计算。现在报的"高渗组 h 加密差 0.84%"是振荡的巧合，不是收敛度量。
不要再用"intermediate-to-main P1 差 1.61% / 0.84%"论证收敛。改为：Control 报 GCI（0.59%）；Hyperosmotic 明确写"h 收敛非单调，故不报 GCI，改以 P2 双网格一致性（0.13 pp）作为离散化证据"。Table S3 里 "Mean endpoint passes h- and p-refinement" 这条 reason 必须重写。
Control 的 P1 三点 Richardson 极限是 0.066378，而 P2 在中等网格就给出 0.067416，差 1.54%。同一个连续问题的两族离散化不可能收敛到不同极限。只有两种解释：
第 2 条几乎肯定是主因，理由见下条。
从 dof 反推（P1 向量场，3 dof/节点，椭球体积 604.6 µm³）：
| 网格 | 顶点数 | 平均间距 | 与壳厚 0.12 µm 之比 |
|---|---|---|---|
| coarse | 710 | 0.95 µm | 7.9× |
| intermediate | 2 257 | 0.64 µm | 5.4× |
| main | 7 531 | 0.43 µm | 3.6× |
coarse_h1p25_P1 / Hyperosmotic 的残差是 2.67e-8，比其余所有算例（~1e-10）大 100 倍，而全部 UQ 都跑在这个粗网格上。得说明这是不是普遍情况——建议把 141 对 UQ 的残差分布（直方图/最大值）放进 Figure S5。核的 K/G = 30 → ν = 0.4835（近不可压）。线性四面体是近不可压有限变形弹性里最差的单元，体积锁定是教科书级别的问题。这正好解释了：为什么 P1/P2 同网格差 3.58%/2.69%（p 误差 > h 误差），为什么 P1 族收敛慢而 P2 族一到粗网格就稳。
production 用 P2（同上），或改用混合 u–p / F-bar / MINI 单元。并在 Figure S2 里加一张"P1 锁定诊断"图（例如同网格 P1/P2 的 J 场对比，或体积应变能占比）。这条改完，§1.1–1.3 的问题一起解决。
Newton 收敛 + J>0 不等于平衡是稳定的。你有一个薄壳（0.12 µm），UQ 里还出现了整体平均 resultant 穿越零（→ 净受压）的算例。真实 lamina 在净压缩下会皱曲/屈曲，这个模型（光滑 Newton、对称网格、无缺陷、无接触）既不会捕捉屈曲也不会告诉你它错过了分叉点。
你已经有投影后的一致切线算子和 MINRES，算最小几个特征值几乎免费（LOBPCG）。对所有被接受的平衡态报告投影切线的最小特征值 > 0。并把任何 resultant 均值 ≤ 0 的算例明确划出模型有效域（那个 −164.35% 的点尤其）。
自写的 scikit-fem 工作流，带乘法生长分解 Fe = F·Fg⁻¹ 和初应力残差化——这两块最容易写错，但全文只有"projected MINRES 与 augmented direct 在粗 P2 上机器精度比对"（那只验证了线性解法器，没验证残差/切线）。
没有这三样，"自定义代码 + 22 593 dof + 6 位有效数字"的组合会让审稿人直接质疑代码正确性。
这是最担心的一条，也是审稿人一定会问的。
模型里壳的张力来自两项：初始预张力 0.05 mN/m 和 375 Pa 探针。现在施加各向同性自然收缩 Fg：
在参数盒的任何一点，几何上都不存在张力升高的路径。所以 64/64 和 77/77 不翻号是必然的，跟"稳健"没关系。把它当作 Main Figure 3 的核心内容，等于用一整张主图证明了一件由模型定义决定的事。
更要命的是 §5 里的表述读起来像发现，其实是定义："All 64 paired Latin-hypercube points and all 77 parameter-domain boundary/anchor points predict a lower Hyperosmotic global mean"。
换算一下（t = 0.12 µm）：
| 量 | 值 | 等效膜应力 |
|---|---|---|
| Control 端点 | 0.066064 mN/m | 551 Pa |
| Hyper 端点 | 0.038092 mN/m | 317 Pa |
| 假设的初始预张力 | 0.05 mN/m | 417 Pa |
预张力占 Control 端点的 76%。头条那个 −42%（正确值 −44.5%）主要不是"文献参数化的力学"，而是一个 provenance 标为 assumption、无任何文献来源的 0.05 mN/m 在几何收缩下自然松弛。而且 UQ 的 pretension 范围是 0.025–0.1，从来没有包含 0，所以这个分解从头到尾没做过。
三算例机制分解，共 6 次求解，一晚上能跑完：(i) pretension = 0，只有探针；(ii) probe = 0，只有预张力；(iii) 完整（现有）。然后报告 N_full 与 N_pretension + N_probe 的差（耦合项）。结论句改成："预测的卸载中 X% 来自假设预张力的松弛、Y% 来自几何缩小下的探针响应、Z% 来自二者耦合。"这条做完，整篇文章的可信度会有质的变化——你在诚实告诉读者头条数字的来源结构。
论文的输入卖点是 Reynolds 的 10 级异质染色质图。但端点是外表面的面积加权全局平均。核内异质性往外表面传递时极可能被平均掉。
如果把核换成体积平均（或 Voigt/Reuss）的均质材料，端点还是 −44.5%，那么整套数字化图像、灰度标定、10 级映射的工作对头条结论毫无贡献。审稿人一定会问。
跑均质核对照（2 次求解）。两种结果都对你有利：
不管哪种，现在不做就是留着被打。
正文写 "self-equilibrated 375 Pa nominal shear probe"、Table S4 写 probe_shear_stress_pa，但 §9 第 2 节写的是 "逐点牵引为 t₀ = P̄·N"。
t₀ = P̄N 是沿法向的牵引，即压力，不是剪切。各向同性压力加在 K/G=30 的近不可压体上基本什么也不会发生（这也许正是端点被预张力主导的原因）。
把载荷说清楚，名实一致。若确实是法向压力，全文的 "shear probe" / probe_shear_stress_pa 全部改名（如 probe_normal_traction_pa）。若本意是剪切，则载荷实现是错的，需要重做。另外要写明这是参考构型上的 dead load 还是当前构型上的 follower pressure——两组参考构型相同，dead load 意味着两组受完全相同的总外力（这其实正是你要的 "same probe"，值得明说）；follower 则外力随当前面积缩小，会系统性改变报告的差值。这个选择必须显式声明并做一次敏感性对照。
量级参照：200 mOsm 的渗透压差在 310 K 下约 515 kPa，比 375 Pa 探针大约 1 400 倍。你没有把渗透压当外压加（这是对的），但正文最好主动点明这个量级对比，否则读者会困惑。
带预张力的闭合壳需要核内有平衡压力（类似压力容器 T = pR/2）。只往壳里塞 T₀ = 0.05 mN/m 而没有配平，参考构型不处于平衡，第一次 Newton 就会把它松掉，实际实现的预张力不等于 0.05。
正文说 "residualized reference-configuration initial stress"——如果 "residualized" 指减掉 u=0 处的内力残差（常见技巧），那这个初应力就成了只进切线不进平衡的虚拟量，而它偏偏占端点的 76%。
明确写清是哪种做法，并给出初始态自检：在 u=0 和零载荷步后分别计算实现的壳张力，证明它等于（或差多少）标称的 0.05 mN/m。这个检查放 Figure S1。
从 Table S2 反推可确认最小 Je 出现在壳里：min Je(shell) × Jg(shell) = 0.69278 × 0.7255170.35 = 0.61918——与表里 minimum_total_jacobian 精确吻合，说明 coupling 按 Jg,shell = Jnat0.35 实现（这条内部一致，没问题）。
但结论是：壳被体积压缩了 31%（Je = 0.693）。真实 lamina 在面内近乎面积不可压、通过减薄响应，不会损失 31% 体积。用 ν = 0.35 的可压缩 Neo-Hookean 实体做 0.12 µm 的 lamina，会产生大量伪体积弹性应变，直接污染 resultant proxy。
壳改用膜/壳公式（面内 2D 应变能 + 由 J=1 决定的厚度减薄），或至少把壳的 ν 提到 0.49+ 并给出敏感性。同时报告 Je(shell) 的分布——0.693 这个数一旦被审稿人反推出来，很难解释。
| 表里写的 | 隐含的 | 冲突 |
|---|---|---|
| Core Poisson ratio = 0.45 | K/G = 30 → ν = 0.4835 | 同一个量两个值 |
| Core Young modulus = 800 Pa | G ∈ [170, 4000] Pa, ν=0.4835 → E ∈ [504, 11 868] Pa | 800 Pa 是什么？均质回退值？ |
Table S1 明确 core 哪些行"仅在均质对照中使用"、哪些是异质映射的实际输入；ν 只保留一个来源（建议由 K/G 派生，删掉 0.45）。
强加的是各向同性自然收缩，但贴壁细胞核在高渗下主要是变扁（高度方向优先塌陷）。各向异性 Fg 会显著改变壳 resultant 分布。而 §11 恰好提出用"核体积/轴比"做湿实验读出——那就更应该把轴比作为模型的可检验预测做出来。
补一个各向异性 Fg 敏感性算例，并把预测的三个轴比变化写进 Table S3 作为一条 supported claim + 实验读出。这是全篇少数"模型给出非平凡可测预测"的机会，不要浪费。
"面积加权的有符号 resultant 平均"，在有大面积受压区时会正负抵消（你们自己就撞上了一个穿零算例）。一个壳可以有巨大的等量拉压而平均为零。这也是为什么不得不引入"绝对差"作为对抗性目标。
主端点换成不会抵消的量：面积加权 max principal resultant、平均 |N|、或壳的总弹性能。原有的有符号平均保留为次要指标。这一改，穿零算例和 −164.35% 这种没有意义的百分数一起消失。
如果参数盒里不存在任何一点给出"加载"，那么 Flipper 的方向性检验就没有区分力：它无法把你的模型和一个平凡假设（"东西缩了就松了"）区分开。一个不可能错的预测不是检验。
构造一个真正会给出相反方向的竞争假设，并证明 FEM 在某参数集下确实能产生它。候选：
然后 Flipper 实验就变成"在两个互斥预测之间做判决"，实验价值和论文价值同时上一个台阶。
141 对（282 次求解）全在 coarse_h1p25，而已知该网格把名义点的卸载低估 6.7 pp（−37.9% vs 收敛的 −44.5%）。所以 "P5/中位/P95 = −65.5% / −27.4% / −9.1%" 整体带系统偏差。
好消息是偏差方向保守：粗网格系统性低估卸载幅度（名义点 −37.9 → −40.9 → −42.3；最不利顶点 corner_101011 从 −4.74% 到中等网格 −5.69%，幅度也变大）。
resultant N = σ·t，壳厚 t 与端点成正比，但 0.12 µm 被固定死、零不确定度，provenance 还标着 assumption。它不会在相对变化里抵消，因为两组的当前厚度不同。
同样缺席的：几何（半轴、长宽比）、核-壳界面耦合（真实是 LAD 部分耦合，比 "shell shrinkage coupling" 更重要）、K/G 比、灰度→模量映射的形状（下一条）。
把壳厚（如 0.06–0.25 µm）和界面耦合加进参数盒，或在正文显式论证为什么排除。前者更省事。
用低/中/高相比例 57%/31%/12% 做单调标定，这个想法很好，应该在正文更突出。但它只钉住 3 个分位点，箱内映射形状（S 形 vs 分段线性 vs 阶梯）完全自由，而模量跨度是 23.5 倍（170–4000 Pa）。
另外：出版图的灰度不是定量 DNA 密度（显示 gamma、亮区截断、有损压缩）；8 张 z 面板线性插值到 50 层是很强的 z 分辨率外推；遮盖 panel A 标记后那块区域填了什么值也没写。
跑 3–4 个满足同样 3 个分位约束的不同单调映射，报告端点散布；这个散布就是数字化不确定度的诚实估计。遮盖区处理方式（inpaint / 排除出积分域）写进 provenance。
两条都不贵，路线 B 尤其划算，数据已经在手上了。
离散化误差 ≥ 3.6%（同网格 P1/P2），跨离散化的组间变化散布 6.7 pp，但正文在报：
| 现在报的 | 应改为 |
|---|---|
| −42.34% | −44%（或 −44.5%） |
| 0.066064 / 0.038092 mN/m | 0.066 / 0.038 mN/m |
| 0.725517（源头 φs=0.204 只有 3 位） | 0.726 |
| 0.735031 | 0.735 |
| −4.74% / −5.69% / −164.35% | 一律降位，最后一个删掉改绝对差 |
规则：报告精度不能超过验证精度。这是审稿人最容易抓、也最伤信任的一类问题。
前者是外体积比（含收缩较少的壳），后者是核的自然体积比。拿它们说 "close" 是在比不同的东西。应该比同类：核的平衡体积比 vs 0.726。
um → µm；phi_s → φs；1e-6 → 10⁻⁶；V/V0=phi_s+(1-phi_s)C0/C 应排成公式。prospective_validation_pending、corner_101011 这类程序内部标识符不应出现在正文散文里。<w:p .../> 标签吃掉造成的，docx 文件本身格式完好，不用管。Flipper-TR 的寿命主要反映脂膜张力/脂质有序度，而你们明确声明 proxy 不是脂膜张力。那么 "lamina resultant 下降 → Flipper 寿命下降" 这个方向关系本身就是一个未经证明的假设。现在的逻辑是：用假设 X 去检验预测 Y，而 X 和 Y 的联系没人验证过。
一个方向性变化不足以判定力学假说。
模型预测的三轴轴比变化是共聚焦直接可测的，成本远低于 FLIM，而且是模型的非平凡输出（体积比是输入，轴比不是）。现在只在 §11 一句带过，既不在 claims 矩阵里，也不在任何图里。
把"预测的轴比/形态变化"提升为 Table S3 的一条 supported claim，作为第一道实验验证。配合 §2.9 的各向异性 Fg 敏感性一起做。
| # | 动作 | 求解次数 | 对应节 |
|---|---|---|---|
| 1 | 机制分解：pretension=0 / probe=0 / 完整 | 6 | §2.2 |
| 2 | 均质核对照 | 2 | §2.3 |
| 3 | production 换 P2，补 main_h0p45_P2 | 2 | §1.1 |
| 4 | 头条数字改 −44%，全文有效数字降位 | 0 | §1.1, §4.1 |
| 5 | 写清壳的单元公式 + resultant proxy 定义式 | 0 | §1.4 |
| 6 | 相对残差替代绝对残差；统一 1e-6/1e-7 | 0 | §1.5 |
| 7 | 载荷 shear/normal 命名与实现对齐 | 0–2 | §2.4 |
| 8 | 高渗组非单调收敛如实报告，重写 Table S3 对应 reason | 0 | §1.2 |
| # | 动作 | 求解次数 | 对应节 |
|---|---|---|---|
| 9 | 解析推导 dN/dλ > 0，方向性从 "UQ 发现" 改为 "解析结果"，重构 Fig 3 | 0 | §2.1 |
| 10 | 材料场网格无关投影 + 积分阶收敛检查 | ~6 | §1.4 |
| 11 | 稳定性检查（投影切线最小特征值） | 复用已有解 | §1.7 |
| 12 | 代码验证三件套（FD 切线 / 解析基准 / 独立软件交叉） | ~10 | §1.8 |
| 13 | ±10% 过渡带 6 算例 + 分层 6–10 点在 intermediate P2 重跑 | ~30 | §3.2 |
| 14 | GP/PCE 代理 + Sobol + 连续域全局最不利点搜索与验证 | ~10 | §3.5 |
| 15 | 壳厚、界面耦合加入参数盒 | ~30 | §3.3 |
| 16 | 各向异性 Fg 敏感性 + 轴比预测 | ~6 | §2.9, §5.3 |
| 17 | 壳改膜/近不可压公式 | 重跑主算例 | §2.6 |
| # | 动作 | 对应节 |
|---|---|---|
| 18 | 构造会给出相反方向的竞争假设（压实增刚 / 面积守恒皱曲），做成判决性实验 | §3.1 |
| 19 | 灰度映射形状不确定度（3–4 个满足同分位约束的映射） | §3.4 |
| 20 | Flipper 改剂量–响应 + 正交扰动 + 可逆性 + 序参量对照 | §5.2 |
| 21 | 全文统一英文、Nature 格式化、参考文献拆条 | §4.3 |