同构惊讶 9 · 具体 9 · 可行 8

1995 年 Mark Ghiorso 与 Richard Sack 在 Contributions to Mineralogy and Petrology 119: 197-212 发表 MELTS（“Chemical Mass Transfer in Magmatic Processes IV”）。MELTS 把”输入 bulk composition（12 个氧化物比例）+ T/P/H/S” 当作 source language，把”输出矿物相组合 + 成分”当作 object code；中间的”编译过程”就是 Gibbs 自由能最小化（constrained minimization of thermodynamic potential）。从编译器视角：bulk composition 是 token stream，T/P 是 optimization flags，Gibbs 最小化是 constant-folding pass —— 因为相组合的选择在编译期就被确定，运行期（岩浆演化）只是 token-level reinterpretation。1995 vintage 让这个类比非常贴切：当年 GCC 2.7.2 是 Linux 默认编译器，MELTS 的 Fortran 77 代码也是同代产物；两者共享”精度不够靠约束补”的工程哲学 —— MELTS 用 stoichiometric conservation matrix（元素守恒）保证数值稳定，正如 GCC 2.7 用 long double + IEEE 754 模式补 32-bit 精度。

如果要做写一篇 “MELTS 1995: petrology’s GCC 2.7”——左栏画编译器三阶段（lex → parse → codegen），右栏画 MELTS 三阶段（composition input → Gibbs formulation → phase assemblage output），中间标注 “Gibbs minimization = constant-folding pass”。附录：把 MELTS 的 stoichiometric matrix 当作 type system 画出来（12 种氧化物 = 12 种 primitive types）。

同构惊讶 8 · 具体 8 · 可行 7

Mats Hillert 1988 年提出 Compound Energy Formalism（CEF），把每个晶体相描述为多个 sublattice（亚晶格）的组合，每个 sublattice 上分布着若干 species（物种/元素），sublattice 内的 mixing 用 Redlich-Kister polynomial 表达。这在结构上完全等同于 ML/Haskell 的parametric type system：phase（矿物相）= 类型构造子；sublattice = 类型参数；species 在 sublattice 上的分布 = 类型实例化时的 witness。1995 vintage 时 CEF 已经成熟，Sun-Sundman-Ågren 等人在 1980s 末完成了一般化推广。CALPHAD 用 CEF 处理有序-无序转变、间隙固溶体、离子相等复杂相的能力，等同于编译器用 type theory 处理 ad-hoc polymorphism：让”相图”（类型规则）可以表达任意 sublattice 配置而无需重新发明轮子。

如果要做写一篇 “Sublattices as types: why petrology invented type theory before Haskell”——画一张并排对照：(左) Haskell 的 `Maybe (Map String Int)` 表达式树，每个构造子带类型参数；(右) CEF 的 `(Mg,Fe²⁺,Ca)M1 (Mg,Fe²⁺,Ti,Fe³⁺,Al)M1 (Fe³⁺,Al,Si)2 TETO6` pyroxene 表达式，每个 sublattice 带 species 集。两张图视觉上完全同构 —— 都是带结构参数的复合类型构造子。附文：1988 Hillert 比 1990 Hudak-Paulson 的 Haskell 报告早 2 年。

方法论惊讶 8 · 具体 8 · 可行 8

1990 年 Robert Fourer, David Gay, Brian Kernighan 在 Bell Labs 发表 AMPL（A Mathematical Programming Language）—— 一种声明式代数建模语言，用于描述并求解大规模数学优化问题。CALPHAD 的核心是”在元素守恒约束下最小化 Gibbs 自由能”——这正是 AMPL 最擅长表达的”约束满足 + 优化”模式。1995 vintage 时 Chen 等人开始用 AMPL 重新表达 CALPHAD 问题（“Set based framework for Gibbs energy minimization”），让 petrology 从”写 Fortran 子程序”变成”声明变量 + 约束 + 目标函数”。这等同于编译器从手写 machine code 升级到使用高级宏语言：AMPL 脚本是 petrology 的 source language，编译器后端（IPOPT、MINOS）做 codegen + optimization。

如果要做写一篇 “AMPL: petrology’s answer to Lisp macros”——给出 AMPL 写的一个 toy Gibbs minimization model（5 个元素 + 3 个相），左栏对比 Fortran 77 等价代码（300 行），右栏是 AMPL 版本（30 行）。两个版本输出完全相同的相图，但 AMPL 版本把”什么是物理约束”和”怎么求解”解耦 —— 这正是宏语言的本质。下一步：把 OpenCalphad 包装成 AMPL callable library。

文化隐喻惊讶 8 · 具体 8 · 可行 7

1995 年原版 MELTS 只能模拟压力 0-2 GPa（地壳浅部），因为 1995 vintage GCC 2.7.x 的 -O2 -ffast-math 在 64-bit 双精度浮点上仍有 edge case bug —— 高压 (>2 GPa) 相平衡的 Gibbs 自由能差异极小（mJ/mol 量级），需要严格 IEEE 754 舍入模式。Ghiorso 等 2002 年的 pMELTS（Geochemistry Geophysics Geosystems 3: 2001GC000217）把压力上限推到 3 GPa，跨越上地幔 —— 这次扩展的关键不在物模型，而在编译器基础设施升级：g77 / GCC 3.0 / Intel Fortran 7 让严格 IEEE 754 浮点变得可靠（不再需要 -ffloat-store 这类 GCC 2.x hack）。这意味着 petrology 的”mantle 全段模拟”被编译器版本间接解锁 —— 一个”低层 bug 修复 → 高层科学发现”的清晰链路。Side note：GCC 3.0 2001 年发布，正是 pMELTS 同年发表；GCC 2.95 到 GCC 3.0 的浮点重写（“Voltaire” 重构）让 mantissa 表达完全符合 IEEE 754。

如果要做写一篇 “When the compiler version unlocked the mantle: pMELTS 2002 and GCC 3.0”——做一个小实验：用 GCC 2.95 和 GCC 3.0 分别编译同一份 Gibbs 自由能计算子程序，对比 >2 GPa 区间的输出差异；预期 GCC 2.95 在某些 input ranges 有 mJ/mol 量级偏差（pMELTS 论文 supplementary 给了参数）。再写一段”为什么地球物理学家不该假设编译器’就是对的’”。

同构惊讶 9 · 具体 9 · 可行 8

Hurrian hymn h.6（约公元前 1400 年 Ugarit 出土的 cuneiform 板）用楔形符号同时记录歌词（Hurrian 语言）+ 9 弦琴的弦名 + 音程间距 —— 它是最早的”音乐源文件”。West 1994 在 Music and Letters 75: 161-179 详细分析了其音律学含义。公元前 3 世纪，希腊化亚历山大的 Ctesibius 发明 hydraulis（水风琴）—— 这是已知最早用物理装置”读取”音乐指令的机器：水泵产生稳定气压，键盘控制气流进入不同长度的笛管。Hydraulis 的键盘本质上就是把”音程编号”翻译成”管长”的物理 decoder —— 与 h.6 楔形符号把”音程编号”翻译成”琴弦张力”的功能完全对应。两者间隔 1100 年，恰好是 cuneiform 记号传统被希腊化世界吸收、并被工程师重新实现为机械的过渡期。

如果要做写一篇 “Hydraulis = Hurrian hymn’s physical decoder, 1100 years late”——画一张三阶段对照：(左) h.6 cuneiform 板（楔形符号 + 9 弦记号，平面 clay substrate），(中) Ctesibius hydraulis 键盘 + 笛管（凸面物理按键 + 长度映射到音高），(右) 现代 MIDI keyboard（电压信号 + 采样器）。三者在”符号 → 物理动作 → 声音”的链路上同构。附文：1948 年希腊考古出土的 Aquincum hydraulis 副本是已知最早的物理重建实例。

同构惊讶 8 · 具体 8 · 可行 7

al-Jazari（1136-1206 CE）的浮动乐师船是已知最早的”可编程自动乐队”：4 个机械乐师在池塘中央的小船上演奏，旋转的钉柱圆筒（cam cylinder）触发不同杠杆、对应不同乐器的动作。Noel Sharkey 论证其本质是”可编程自动机”—— 每个 cam 上的钉 = 一条指令，cam 转动一周 = 一段乐谱。Sumerian cuneiform 的革命性在于它是第一个 logographic script（一个符号对应一个完整词/概念，不是 phoneme）。al-Jazari 浮动乐师船的 cam cylinder 完美复现了这个性质：每个钉 = 一个完整动作（不是基础 trigger），整周 cam = 一首完整乐曲。这是把 cuneiform “sign-to-meaning binding” 哲学物理化的最佳案例 —— 楔形在 clay 上、钉在 cam 上，两者都是”压一个印记 = 启动一个完整语义单元”。

如果要做写一篇 “al-Jazari’s floating band: cuneiform’s logographic principle, mechanically realized”——画两张并排：(左) Sumerian 楔形符号表（每个符号对应一个完整词，如 “lú” = man, “šà” = heart），(右) al-Jazari cam cylinder 展开图（每个钉 = 一个完整动作，如拉弓/敲鼓/吹笛），两图下方并列：“1 symbol = 1 word” vs “1 pin = 1 musical phrase”。附文：两者都是”atomic semantic units”，区别仅在 substrate（clay vs metal）。

方法论惊讶 8 · 具体 8 · 可行 8

Banū Mūsā 兄弟（约 850 CE，Baghdad House of Wisdom）发明自动笛手——使用钉柱滚筒（peg cylinder）触发不同笛孔。这是有史以来第一个 known 可编程机器（Koetsier 2001 “On the prehistory of programmable machines”）。从符号学角度：滚筒上钉的分布 = 一段乐谱的物理表示；钉旋转到触发位置 = 读取一个 music instruction；不同钉的高度 = 不同 pitch。这与 cuneiform 楔形在 clay 上的分布完全同构：cuneiform wedge = pin；clay tablet surface = cylinder surface；Hymn to Nikkal h.6 的楔形序列 = 滚筒上的钉序列。“可编程性”的关键不在电路，而在”机械载体能承载任意的、有结构的符号序列”—— Banū Mūsā 物理化了 cuneiform 的”刻在 substrate 上的有序符号序列”这一核心思想。Koetsier 明确论证 Banū Mūsā 的自动笛手是 Jacquard 提花机（1804）的直接祖先 —— 从 clay tablet 到 perforated card 的物质技术演化链条上，Banū Mūsā 是中间环节。

如果要做写一篇 “From clay wedges to brass pins: Banū Mūsā’s peg cylinder as the physical transpiler of cuneiform”——画一张演化时间轴：(公元前 3200) Sumerian 楔形 on clay → (公元前 1400) Hurrian 楔形音乐 notation on clay → (公元 850) Banū Mūsā peg cylinder on brass → (1804) Jacquard perforated card → (1940s) Hollerith punch card → (1970s) paper tape。每个阶段都是”ordered symbols on substrate” 的物理实现，载体从 clay 演到 metal 到 paper 到 mylar。每个符号承载的语义从”logogram” 演到”machine instruction”。

方法论惊讶 8 · 具体 9 · 可行 8

1970 年代以来，Anne Kilmer、Richard Crocker、Martin West 等学者尝试重建 h.6 的实际音乐 —— 至少有 5 种相互竞争的重构版本（Duchesne-Guillemin 公开演奏版最知名）。每个版本必须同时做出两个决定：(1) cuneiform 楔形符号到底对应哪些具体音高（存在多义性，例如 “idu” 既可指一种音程也可指一种调式）；(2) 9 弦琴的具体调音（West 推断为 diatonic 七声音阶，但具体音高有 ±50 cent 的不确定性）。实体装置约束强迫我们把这两个不确定性同时解决：物理装置一旦造成，就锁定了两个决定 —— 琴的调音决定哪些音高可用，cuneiform 翻译规则决定如何分配音高给符号。这是 instrument design 的”反问题”：传统乐器先有乐谱再调音；这里必须先造装置再确认乐谱能演奏。结果是：物理装置的选择反向验证了 cuneiform 翻译规则 —— 哪些规则能让装置演奏出合理的旋律？这是一个全新的”考古乐理学方法论”。

如果要做写一篇 “Reverse engineering h.6: when the physical instrument constrains the cuneiform decipherment”——做一个 prototype：用 3D 打印造一个 9 弦现代版 hurrian sammûm（lyre），让它能演奏 West 1994 提出的所有 5 种音高分配方案；在每个方案下，让 musicologist 给旋律打分（流畅度、调式稳定性）；统计哪个音高分配方案在物理装置上”听感最好”。这是把 instrument design 变成 cuneiform 翻译规则的判别准则。