用指数倾斜热力学映射在温度、压力与化学势维度上预测相变
核心概要
该工作提出指数倾斜热力学映射（expTM），在原有热力学映射的高斯先验中引入指数倾斜因子，使先验的均值与方差分别对应压力或化学势与温度，从而仅用远离相边界的少量训练数据即可在任意温度、压力或化学势条件下生成热力学正确的样本；作者在巨正则系综下用两个数据点复现了格子气模型的相变（密度差在临界化学势附近之外保持在0.05以内），并在等温等压系综下预测了CO2在压力变化下的相变，识别出Phase I与Phase III之间约在4.5–5.4 GPa出现的中间态。
illustrated in Fig. 1, with its complete theoretical formulation provided in the “Methods” section. We demonstrate its effectiveness through
· 第 2 页深度剖析
提出expTM框架，把指数倾斜因子e^{βαx}引入扩散模型的高斯先验，使先验成为均值α、方差β^{-1}的倾斜高斯，其中α对应压力（α=−P）或化学势（α=μ），β^{-1}对应温度。 原有热力学映射（TM）只能通过正则系综编码温度涨落，先验为N(0, β^{-1})，无法同时处理压力或化学势；expTM把控制参数直接映射到先验均值，从而把框架推广到等温等压（NPT）与巨正则（μVT）系综。 方法层面给出完整的随机微分方程推导（前向与反向SDE，式20、21），并在(x, η)∈R^{3D}联合空间学习联合得分函数，而非把η当作固定标签；理论推导与两个算例相互印证。
在巨正则系综的格子气模型中，仅用单一温度T=7下两个化学势（μ=−16与μ=0）的训练数据，expTM即可在μ∈[−20,4]、T∈[1,13]的625个网格点上生成粒子密度，并复现μc=−8处的临界边界。 相比原TM只用温度作为控制变量，这里首次同时控制温度与化学势两个参数；消融实验显示，去掉指数倾斜项、改用与热力学参数无关的单位高斯先验后，模型无法恢复MCMC中观察到的尖锐相边界。 以Monte Carlo Metropolis-Hastings模拟为基准，密度差Δρ=|ρ_MC−ρ_expTM|在临界化学势μc=−8附近之外保持在0.05以内；消融结果见补充材料Note 4与补充图3–5。
在等温等压系综下，expTM用P=1 GPa与P=8 GPa两个压力下的训练构型，预测了CO2在压力与温度变化下的相行为，识别出Phase I与Phase III之间的中间态出现在约P=4.5–5.4 GPa。 该算例把expTM从格点模型推广到真实分子晶体，用原子级集体变量λ1与λ3（各256维特征向量）刻画相态，展示了在MD模拟难以覆盖的高压/低压区域捕捉结构排布的能力。 训练数据取自增强采样MD（well-tempered metadynamics）在350 K、1/3/5/8 GPa下的模拟，每个压力按玻尔兹曼分布选取10^5个构型；相边界用λ1阈值0.58判定，并以0到1的“接近Phase III”指标量化，自由能差ΔG_{I-III}由三个分块平均给出误差棒。
expTM在跨多个热力学条件采样时相对MCMC具有计算优势，且其先验的均值与方差分别与压力/化学势和温度对应，保持了高斯先验的可解性与简洁性。 原TM仅把温度涨落对应到先验方差；expTM额外引入作用于先验均值的涨落变量，使先验更具表达力，同时不牺牲高斯分布的可处理性。 墙钟时间基准比较见补充材料Note 3；理论部分给出⟨x⟩=α、Var(x)=β^{-1}的解析结果。
启示与展望
该框架面向的是在远离相边界的少数热力学条件下采集到有限观测、又希望推断其他温度、压力或化学势下平衡行为的场景；作者展示的两个算例分别是巨正则系综下的格子气模型与等温等压系综下的CO2晶体，训练数据分别只有两个数据点和两个压力下的构型。对希望用少量模拟或实验数据外推相图、并需要热力学一致样本的研究者，这一方法提供了可复用的思路；其适用前提是体系存在定义良好的相边界，且所选特征（如λ1、λ3）能区分目标相态。
读者可能仍会关注：在临界化学势μc=−8附近密度差超出0.05，说明临界点附近的精度是后续可继续考察的方向；CO2算例中相边界判定依赖λ1阈值0.58这一调整，不同阈值下中间态区间的稳健性值得进一步观察；训练数据在两个压力下各自只覆盖单一相，向更多相或更极端条件的推广尚待验证；作者提到的自旋玻璃、晶体成核、固有无序蛋白与RNA等应用属于展望，本文未给出结果；此外，本文为Article in Press版本，补充材料中的消融、时间基准与MD设置细节未随正文一并加载，若需评估这些细节需查阅补充信息。