化学平衡与物理平衡塑造有机晶体中双重冰成核(ice nucleation)路径——以间苯三酚(phloroglucinol, PGL)为例

《Communications Chemistry》：Chemical and physical equilibria shape dual ice-nucleation pathways in an organic crystal

　　

有机晶体在自然环境的冰形成过程中起关键作用，但其冰成核活性的分子机制仍不清楚。本研究表明，间苯三酚（PGL，一种简单的多羟基芳香化合物）表现出的优异冻结效率，源于化学平衡与物理平衡的协同作用，二者共同产生两条并行的成核路径。一条路径源自PGL晶体表面，其作为强效冰模板(ice template)；另一条发生在溶液中，溶解的PGL分子组装为纳米尺度聚集体(aggregate)，即使在低于溶解度极限的浓度下也能诱导冰成核。研究人员发现，碱性pH通过诱导互变异构(tautomeric shift)改变氢键模式并抑制分子组装，从而消除上述两条路径。溶液中的老化(aging)可能通过氧化或聚合过程降低冰成核活性，而冻融循环(freeze-thaw cycling)通过生成新鲜PGL晶体可部分恢复活性。这些结果表明分子结构、溶解度、pH依赖的互变异构及相行为是共同控制冰成核的关键决定因素，为理解和预测酚类及其他有机冰成核剂(ice nucleator, IN)在大气及低温生物学复杂环境中的活性提供了可推广的框架。

本研究发表于《Communications Chemistry》，围绕简单多羟基芳香化合物——间苯三酚（phloroglucinol, PGL）——探究有机冰成核剂(ice nucleator, IN)的冰成核分子机制。已有研究表明有机INs广泛存在于大气气溶胶、土壤和生物流体中，但传统无机INs的晶格匹配(lattice matching)与表面极性概念无法解释有机分子柔性结构、质子化平衡、多相态下的冰成核效率差异，且PGL文献报道的冻结温度高度分散，无法用单一静态成核模型解释，因此有必要从化学平衡（互变异构、质子化）与物理平衡（溶解度、相变、聚集）耦合的角度揭示其成核机理。

研究人员采用的主要关键技术方法为：（1）高通量液滴冻结实验——采用双板冰成核检测(Twin-plate Ice Nucleation Assay, TINA)，以连续降温速率1 °C min^?1测定不同浓度及pH下PGL溶液的结冰分数f_ice及半数冻结温度T₅₀，并以Vali公式计算单位质量冰成核位点数N_m(T)；（2）拉曼光谱(Raman spectroscopy)——表征不同pH下PGL的酮-烯醇互变异构及芳香性变化，关联~990 cm^?1（烯醇芳香环面内振动）与~1640 cm^?1（酮式C=O伸缩）信号强度与T₅₀；（3）动态光散射(dynamic light scattering, DLS)——测定溶液中PGL纳米聚集体的水合半径R_h及其随pH的变化；（4）粉末X射线衍射(powder X-ray diffraction, XRD)与光学显微术——确认PGL二水合物晶体物相及晶体粒径分布；（5）异质基底分析法(Heterogeneous Underlying-Based, HUB method)——从f_ice谱中提取冰成核温度微分分布，解析IN亚群；(6) 经典成核理论(classical nucleation theory, CNT)的异质冰成核温度(Heterogeneous Ice Nucleation Temperature, HINT)模拟——关联成核温度与冰结合域尺寸；(7) 基于van't Hoff方程的PGL溶解度-温度关系拟合。

浓度依赖性冻结揭示PGL溶液中存在不止一种冰成核物种

液滴冻结实验显示PGL溶液冻结温度随稀释向更低温偏移；浓度≥5 mg mL^?1时出现双峰冻结模式（~?7 °C与~?12 °C），高浓度时较暖峰占主导并接近晶体PGL活性，稀释后暖峰消失仅剩冷峰。饱和溶液（50 mg mL^?1）含未溶微晶（平均直径约100 μm），分子模拟表明PGL二水合物(010)晶面是高效冰结合位点，与暖峰对应。Nm(T)谱随浓度不重叠且各浓度具独特冻结曲线，说明稀释改变了IN结构种群，表明PGL在溶液中形成浓度依赖的动态物种组合，而非固定IN群体，可能涉及聚集、化学平衡、溶解度或相变调控。

互变异构与质子化平衡调控冰成核活性

pH扫描实验表明酸性条件（pH 1.4）与中性（pH 5.5）冻结行为相近，碱性（pH≥12）完全丧失IN活性；T₅₀在pH≈pK_a1(~8)和pK_a2(~9.2)逐步下降。拉曼光谱显示中性至弱碱条件PGL以芳香烯醇式(I)为主（~990 cm^?1强峰），一次去质子化得单阴离子仍芳香(II)，IN活性略降；二次去质子化(pK_a2)破坏芳香性生成非芳香双负离子酮式(III)（~1640 cm^?1出现C=O峰），此时IN活性完全丧失；更高pH(pK_a3>12)恢复芳香性但无—OH，亦无IN活性。将pH回调可恢复IN活性。研究表明PGL的IN能力严格依赖烯醇式芳香结构与羟基质子化状态，去质子化诱导酮-烯醇互变异构破坏匹配冰的氢键几何，pH波动可显著调控酚类有机IN的环境冰成核行为。

化学老化与相变塑造冰成核剂种群

新鲜与长期储存（老化）PGL溶液对比显示老化使两冻结峰均向低温移动——暖峰中心由?7.8 °C移至?8.9 °C（效率降约56%），冷峰由?12 °C移至?15.2 °C，溶液变黄提示酚类氧化/聚合使PGL分子失活。冻融循环可使老化溶液中约19%暖峰活性恢复，并伴有新生微晶析出，说明冷却驱使PGL过饱和析出具IN活性的新鲜微晶，补偿部分因不可逆化学改性损失的活性，此过程可逆且与生物IN经冻融失活不同，证实PGL IN活性受可逆物理化学过程（溶解度、再结晶）与不可逆化学老化共同塑造。

溶解度平衡决定PGL的冰成核路径

PGL水中溶解度约12 mg mL^?1（20 °C）。浓度≥溶解度限（≥5 mg mL^?1）降温时超出溶解度驱动PGL二水合物晶体析出，对应暖温成核模式（晶体(010)面提供有利晶格匹配与氢键模板）；浓度低于溶解度限（≤1 mg mL^?1）无晶体检出（无沉降、无冻融引起的晶体记忆效应），DLS检测到Rh~86 nm的纳米聚集体，其散射强度与尺寸随pH升高（去质子/失去芳香性）减小至消失，表明聚集体靠π-π相互作用与氢键自组装，介导较低温成核模式。van't Hoff拟合显示≥5 mg mL^?1在测试温区跨越溶解度线，≤1 mg mL^?1始终未饱和。冻融循环中近饱和溶液出现"记忆效应"——暖温IN减少、冷温IN增强，反映新晶体的动态形成与再溶解。综上存在双路径：①溶解度临界以上——既有或新析PGL晶体表面异相成核（高效）；②溶解度以下——溶液中PGL分子自组装纳米聚集体诱导成核（较弱）。

经典成核理论建模量化尺寸-活性关系

HINT-CNT计算表明：微米级晶体具大于~18×18 nm²的有效冰结合面，成核温度趋近平台对应高T_het；溶液中纳米聚集体等效冰结合域约4×4至7×7 nm²时可解释观测到的较低T_het。两家族IN分别与晶体大尺度平整面及纳米聚集体受限结合域相符，从热力学上支撑化学-物理平衡共同调控的双路径成核模型。

总结与讨论（结论部分翻译与归纳）

研究人员得出结论：PGL并非单一静态冰成核剂，而是受化学平衡（酮-烯醇互变异构、电离状态）与物理平衡（溶解度、相行为、聚集）动态调控的系统。当浓度超过溶解度限时，PGL二水合物晶体(010)面提供大面积冰模板，主导高温异相成核；低于溶解度限时，溶解PGL分子自组装为纳米聚集体，以较低效能在更低温诱导冰成核。碱性pH及老化通过破坏芳香烯醇结构与促进氧化/聚合抑制两种路径；冻融循环通过再结晶部分恢复活性。经典成核理论定量吻合两家族IN所对应的冰结合域尺寸。该结果挑战了将IN简单分为"晶体型"与"可溶型"的传统分类，证明相态与浓度可逆切换成核模式，对多羟基芳香/酚类化合物具普适参考意义。研究强调在大气与低温生物学等变动环境中，须同时考量化学与物理平衡方能准确预测有机冰成核剂的冻结行为，为理解复杂体系中有机IN的分子构效关系提供了可推广框架。