GitLab

背景

氢原子最轻，在晶体里像量子小球一样大幅"抖动"（零点运动）。压到极高压、堆出富氢结构时，这种抖动加上势能面的强非谐性，可以让常规（简谐）声子计算给出的"稳定结构"和"超导转变温度 T_c"发生定性的错——这正是 SSCHA（Errea / Calandra / Mauri 等）反复证明的事：LaH_10、H_3S、PdH 的相界与 T_c 在计入量子非谐后才对得上实验。如今人们更进一步，用生成式 AI 直接"造"出候选晶体去找常压（而非百万大气压）的高温超导体——但这些候选的"高 T_c"几乎都建立在"氢老老实实坐在平衡位"这一简谐 / 准简谐假设之上。问题来了：当氢的量子抖动大到让它离开平衡位、扩散甚至两两成键时，这套"高 T_c"还成立吗？

核心问题

对一个由生成模型预测、且在简谐与标准 SSCHA下都被判为常压高 T_c（λ≳ 2）的氢化物，计入完整量子核运动（超越 SSCHA 高斯 ansatz 的路径积分 / 神经网络变分）后，其氢亚晶格是保持定域、T_c 稳健（发现为真），还是扩散 / 二聚、λ 与 T_c 跌落数倍（"高 T_c"是钳住原子核的假象）——二者必居其一；且存在一个可从静态结构直接读出的判据预言落入哪一类，并可由常压合成与 H→D 同位素效应符号证伪。

为什么值得做（PRL 级意义）

• 这是整条"生成式常压氢化物超导发现"路线的地基是否结实的问题，而非又一次大海捞针。 生成模型（InvDesFlow、CrystalFormer 等）+ Eliashberg/McMillan 已经在批量产出"常压、超 McMillan 极限、甚至破液氮温区"的候选；但这些判据默认核是经典钳住的点。本提案要回答的是机制：在哪些情形下量子核运动只是"温和重整"T_c，在哪些情形下它会"灾难性摧毁"超导——这决定了该路线产出的候选有几成是真的。
• 已有一个"同一结构、相反结论"的尖锐活靶子，分歧纯粹在方法对量子核运动的处理。 Li_2AuH_6：生成式发现 + 各向异性 Eliashberg 给出常压 T_c≈140 K、λ=2.84，且标准 SSCHA 判其动力学稳定（无虚频）；但 2026 年的路径积分计算（不依赖 SSCHA 的"存在平衡位"高斯 ansatz）发现它动理学失稳——氢亚晶格塌缩、部分二聚成分子，λ 掉到 0.84、T_c 掉到 22 K。输入是同一个结构，结论差一个数量级，差别只在"氢被允许跑多远"。 把这件事从"两篇打架"升级为"可预先判定、可被合成证伪的普适机制"，就是 PRL 级的贡献。
• 这个组是少数同时握有两个半边的团队。 一边是生成式常压超导发现（曹振东 CrystalFormer / CrystalFormer-RL、欧仕刚"常压超导体 Pareto 前沿"+ DFT/DFPT 流水线）；另一边是氢的量子核效应的世界级方法储备（李子航稠密氢神经网络变分、柴神都 NN-VMC 研究氢的核量子效应、乔贝 PIMC）。Errea/Zurek 阵营有 SSCHA 却没有生成 + NN-VMC；生成发现阵营有候选却缺超越高斯 ansatz 的核量子机器。"生成候选 × 超越 SSCHA 的量子核运动判定"恰好落在该组两套独门工具的交集上。
• 判据 + 同位素 = 干净可证伪。 常压合成远比 187 GPa 可及；常规 BCS 同位素效应 T_c∝M^-0.5（α=+0.5），而"量子脆弱"体系因氢的量子弥散驱动失稳，应表现出反常 / 反同位素效应（如 PdH 的 α<0）——这给了一个即使不精确测到 T_c 也能判生死的实验旋钮。

现状与争议

• 生成式发现 + Eliashberg（"常压 140 K"）：High temperature superconductivity in Li_2AuH_6 mediated by strong electron-phonon coupling under ambient pressure, arXiv:2501.12222 (2025)。InvDesFlow 生成 + 各向异性 Eliashberg：常压 λ=2.84、T_c≈140 K，Au–H 八面体呼吸模 + 含 Li 低频模主导耦合。这是"生成模型造出常压高温超导"的旗舰式声明。
• 路径积分反驳（"其实只有 22 K"，最尖锐靶子）：Kinetic instability and superconductivity in Li_2AuH_6 and Li_2AgH_6 at ambient pressure, arXiv:2604.12367 (2026)。明确指出 SSCHA 依赖"存在平衡位 + 有效声子"的高斯 ansatz，对强扩散体系会失效；改用随机路径积分（量子 + 热涨落、无高斯假设）发现氢亚晶格扩散 / 二聚，λ(0)=0.838、T_c≈22 K，远低于先前 80–140 K。与上一条用近似不同、结论差一个数量级——这正是机制问题的核心。
• "a priori 判非谐符号"的 SOTA 判据（但不覆盖塌缩）：F. Belli, E. Zurek, I. Errea, A chemical bonding based descriptor for predicting the role of anharmonicity induced by quantum nuclear effects in hydride superconductors, npj Comput. Mater. (2026), DOI 10.1038/s41524-026-01973-7。用 iCOBI / 键价构造判据，把氢化物分为"对称键合（量子效应压低 T_c）"与"非对称键合（量子效应抬高 T_c）"两类。但其框架默认结构有良定义平衡位，预言的是修正的"符号"，并不刻画 Li_2AuH_6 那种"氢亚晶格直接塌缩"的灾难性情形——这正是本提案要补的缺口。
• "简谐选错多形体、唯非谐才对"的高压存在性证明：P. Hou, Y. Ma, H. Xie, M. Li, Y. Cai, Y. Shen, X. Wang, M. Pang, Significant Impact of Quantum and Anharmonic Effects on the Structural Stability and Superconductivity of NbH_3 at High Pressures, arXiv:2410.08606 (2024)。简谐预测 I4̄3d 为基态（33–400 GPa），实验在 187 GPa 见到的却是 Fm3̄m；SSCHA 把 Fm3̄m 的动力学稳定边界下推约 145 GPa、判 I4̄3d 在 187 GPa 失稳，并给出 Fm3̄m 在 187 GPa 的 T_c≈44 K（McMillan，µ^*=0.15）与实验吻合。确立了"量子非谐定性改判结构与 T_c"在富氢体系真实存在。
• 奠基性机制与同位素判据：I. Errea, M. Calandra, F. Mauri, First-Principles Theory of Anharmonicity and the Inverse Isotope Effect in Superconducting Palladium-Hydride Compounds, Phys. Rev. Lett. 111, 177002 (2013), DOI 10.1103/PhysRevLett.111.177002（PdH 反同位素效应 α<0）；I. Errea et al., Quantum crystal structure in the 250-kelvin superconducting lanthanum hydride, Nature 578, 66–69 (2020), DOI 10.1038/s41586-020-1955-z（量子效应稳定 LaH_10，否则简谐失稳）。

切入路径

工具（手段，非卖点）：

• 造 / 选候选：CrystalFormer（Z. Cao, X. Luo, J. Lv, L. Wang, Sci. Bull. 70, 3522 (2025); arXiv:2403.15734）+ CrystalFormer-RL（Z. Cao, L. Wang, arXiv:2504.02367 (2025)，原文已点名 T_c 为下一步奖励）/ CrystalFormer-CSP（Z. Cao, S. Ou, L. Wang, arXiv:2512.18251 (2025)），生成 / 整理常压、目标 λ 的候选与其竞争多形体；通用 ML 力场做廉价弛豫。
• 简谐基线：DFPT 声子 + Eliashberg/McMillan（该组 DFT/DFPT 流水线）。
• 高斯非谐：SSCHA（L. Monacelli et al., J. Phys.: Condens. Matter 33, 363001 (2021)），可由 ML 势加速（F. Belli, E. Zurek, npj Comput. Mater. 11, 87 (2025)）。
• 超越高斯（该组独门优势）：把组内 PIMC / NN-VMC / 神经正则变换用于氢亚晶格的完整量子核运动，判定平衡位是否成立——这是 SSCHA 在 Li_2AuH_6 上失手的那一步。

判据（卖点）：寻找一个可从静态结构直接读出的"量子核脆弱度"序参量，事先预言"温和重整"还是"灾难性塌缩"。候选物理量（择优、可证伪）：(i) 最近邻 H–H 软模 / Born–Oppenheimer 面上的双阱特征；(ii) 氢零点弥散幅度与最近邻 H–H 间距之比 σ_ZP/d_HH（类 Lindemann）；(iii) 键不对称度（iCOBI，承接 npj 2026 但延伸到"塌缩"区）。机制主张：当该序参量越过某阈值，量子涨落把单阱压成双阱 / 抹平势垒，氢离位 → N(E_F) 与 λ 崩塌 → "高 T_c"消失。

最小可交付第一步：以 Li_2AuH_6（必要时加 Li_2AgH_6）为基准，复现"分裂"——简谐 / SSCHA 端（λ≈2.84、T_c≈140 K、无虚频）对完整量子核运动端（氢二聚、λ≈0.84、T_c≈22 K）；并在同一体系上算上述脆弱度序参量。再取一个已知稳健的对照（如某团簇笼合物 CaH_6 / LaH_10 在其稳定压强，或一个"对称键合"SB 类体系）做标定。

成功判据：(a) 在误差范围内同时复现 Li_2AuH_6 的两端极限（这验证 pipeline）；(b) 脆弱度序参量在"脆弱"（Li_2AuH_6/Li_2AgH_6）与"稳健"（笼合物 / SB 类）之间干净分开；(c) 把判据用到一个该组 CrystalFormer-RL 常压候选上，输出二选一（稳健 / 脆弱）+ 实验判别量：可合成性 + 预测的 H→D 同位素效应符号（脆弱 ⇒ 反常 / 反同位素效应，类 PdH）。

风险与证伪

• 超越高斯的核量子动力学很贵：对整块氢化物晶体做 PIMD / NN-VMC 成本高。缓解：先在 Li_2AuH_6 小基准上验证再放大（与组内"先小规模验证"约定一致）；用 ML 力场承担采样、DFT 抽检鞍点。
• 脆弱度序参量可能不普适 / 非单值：若它在标定集上无法把已知"稳健 vs 脆弱"分开，则判据本身被证伪——这是诚实的负结果，仍对场内有约束价值（至少给出"哪些静态特征不可靠"）。
• "动理学失稳"与"有限温热涨落"需区分：必须固定在合成相关温度报告；零温量子塌缩与热致扩散要分别给出。
• ML 力场在二聚鞍点附近易偏置：靠近塌缩路径处需 DFT 校核，避免势面偏差伪造 / 抹平塌缩。
• 命题级证伪（最关键）：若完整量子核运动显示 Li_2AuH_6（及所选候选）其实保持定域氢亚晶格、λ≈2（即 arXiv:2501.12222 对、arXiv:2604.12367 错），则常压高 T_c 对这一类是稳健的，"量子脆弱"主张对该类不成立——这同样是一个干净、决定性的结论，且直接为生成式常压超导发现"松绑"。

👍 可做 · 👎 不行 · 🤔 不确定 · 🚀 优先 · ❓ 需补充