溶剂调控的自由能景观主导沥青质聚集路径

《Energy & Fuels》：Solvent-Regulated Free-Energy Landscapes Govern Asphaltene Aggregation Pathways

　　

沥青质聚集控制着石油体系中许多流动保障挑战，然而溶剂环境在调控聚集路径中的微观作用仍未得到完全理解。研究人员将全原子分子动力学模拟与最小统计聚集框架相结合，研究芳烃溶剂和脂肪烃溶剂如何重塑沥青质从二聚到介观尺度团簇形成的组装过程。尽管甲苯和庚烷均表现出强烈的局部芳环堆叠，模拟揭示了由溶剂介导的界面效应驱动的根本不同的聚集路径。在甲苯中，芳烃溶剂化稳定了暴露的沥青质核心并阻碍持续生长，导致形成小型、开放、弱连接纳米聚集体。相比之下，庚烷通过更大的溶剂配位损失、更宽的构象可达性以及向近乎空间填充团簇的渐进致密化来促进聚集。通过在一个统一框架内定量关联溶剂配位、溶剂接触自由能、聚集体形态和团簇统计，研究人员表明溶剂质量不仅通过溶解度控制聚集，还通过重塑底层聚集自由能景观来调控聚集。所得框架提供了溶剂依赖性聚集的分子解释，并为通过溶剂环境或添加剂设计控制沥青质缔合提供了途径。

**溶剂调控的自由能景观主导沥青质聚集路径——论文解读**

**一、研究背景与科学问题**

沥青质作为原油中最重且极性最强的组分，其自缔合与沉淀倾向对石油开采、运输和炼制过程中的结垢、沉积、黏度增加及流动保障问题具有决定性影响。Yen-Mullins框架将沥青质聚集描述为从单体到纳米聚集体再到更大团簇的分级自组装过程，该模型虽提供了有价值的介观图像，但未能完全解析连接分子结构、溶剂环境与聚集体形态的微观路径。大量实验与模拟研究已确立溶剂质量是沥青质聚集的一阶变量，如芳烃溶剂（甲苯）稳定分散沥青质，而脂肪烃溶剂（庚烷）促进缔合与沉淀，但溶剂依赖性相互作用如何控制从二聚体到更大团簇的进展尚未被系统阐明。现有研究多聚焦于特定相互作用模式或选择性热力学条件，缺乏从二聚体尺度到团簇尺度的统一、溶剂分辨的聚集图像。因此，溶剂质量如何重塑聚集自由能景观并导致从二聚体形成到介观团簇形态的不同路径，仍是未解的关键问题。

**二、研究内容与主要结论**

研究人员采用全原子分子动力学模拟，系统比较了沥青质模型分子Asphaltene 29（C₆₉H₅₃NS₂；摩尔质量960.32 g mol^–1）在甲苯和庚烷中的聚集行为。每个体系包含40个沥青质分子，浓度为7 wt%，分别在两种溶剂中进行长达1000 ns（甲苯）和4000 ns（庚烷）的模拟。研究发现，两种溶剂均表现出强烈的局部芳环堆叠，但聚集路径存在根本性差异：甲苯中的芳烃溶剂化稳定了暴露的沥青质核心并阻碍持续生长，导致形成小型、开放、弱连接的纳米聚集体；庚烷则通过更大的溶剂配位损失、更宽的构象可达性和渐进致密化促进聚集，最终形成近乎空间填充的致密团簇。通过定量关联溶剂配位、溶剂接触自由能、聚集体形态和团簇统计，研究人员证明溶剂质量不仅通过溶解度控制聚集，更通过重塑底层聚集自由能景观来调控聚集路径。相关成果发表在期刊《Energy》。

**三、主要技术方法**

本研究采用的主要技术方法包括：全原子分子动力学模拟（力场参数基于GAFF，静电势经量子化学计算并采用RESP电荷拟合法，长程静电相互作用采用PME方法处理）；二维势均力（PMF）计算（以单体间距d、相对取向cos θ和溶剂配位数Nsolv为集合变量构建自由能景观）；径向分布函数（RDF）分析（用于揭示溶剂特异性短程有序和优先相互作用模式）；团簇尺寸分布分析（利用GROMACS的gmx cluster模块，以4 ?截断距离判定分子间核心距离）；全局形状参数分析（包括回旋半径Rg、非球度Δ和形状参数S）；标度系数分析（R_g ∝ N^ν）；以及统计聚集模型构建（从模拟可观测量参数化得到临界团簇尺寸N*）。

**四、研究结果**

**二聚体形成的自由能景观与溶剂配位**：二维PMF投影表明，两种溶剂中二聚化均通过部分溶剂化构型而非完全去溶剂化接触进行。甲苯中低自由能区域在溶剂配位数上较为局域化，且结合区集中于|cos θ| ≈ 1，表明强偏好对齐堆叠二聚体；庚烷中低自由能区域在Nsolv上延伸更宽，cos θ范围更大，表明取向相关性较弱且倾斜结合构型更多。RDF分解分析显示，两种溶剂中沥青质-沥青质相互作用均由核-核接触主导（芳环堆叠），而甲苯的芳环优先接触沥青质芳核（表面润湿），庚烷则涉及脂肪族片段的核-尾和尾-尾接触模式。

**聚集动力学与标度行为**：甲苯体系中最大团簇尺寸N_max和平均团簇尺寸?N?在约1 μs内趋于稳定，生长在纳米聚集体区间饱和；庚烷体系中两者在整个4 μs轨迹中持续增长并达到更大值。标度分析得到ν_tol ≈ 0.74 ± 0.04和ν_hep ≈ 0.37 ± 0.02，对应有效分形维数D_f≈2.6（庚烷，接近致密三维物体）和D_f≈1.35（甲苯，更开放延伸生长）。

**聚集体形态的时间演化**：甲苯体系中聚集体在整个模拟期间保持分散，Rg分布局域在较小值，Δ和S分布居中于小值，形态主要为弱各向异性；庚烷体系则显示持续生长和显著聚集，Rg分布显著拓宽并移向更大值，Δ和S分布随时间演化为更小且更窄的值，表明结构逐渐收敛为致密近各向同性团簇。

**溶剂依赖性聚集的统计模型**：研究人员构建了最小统计自由能模型F_s(N) = ??_sN + γ_sN^αs，其中内聚参数?_s从PMF图提取，界面参数γ_s≈z_surfΔw_s由RDF导出的溶剂接触自由能尺度与溶剂配位损失估算，指数α_s由标度关系确定。模型预测的临界团簇尺寸N*_theory为甲苯约2–3、庚烷约9–18，与模拟值（甲苯约3、庚烷约20）定量吻合。分析表明溶剂依赖性不能仅由局部核-核吸引力解释，溶剂介导的界面效应（γ_s）和聚集体形态差异（α_s）是主导因素。

**五、结论与展望**

研究结果表明，溶剂质量通过重塑整个聚集自由能景观控制沥青质聚集，而非仅通过改变溶解度。甲苯中的聚集受挫源于芳烃溶剂化对暴露芳环表面的更强稳定化，庚烷中的聚集则归因于更大的溶剂排除自由能增益和生长过程中的致密化。该框架定量重现了溶剂依赖性平均团簇尺寸，捕捉了从甲苯中有限纳米聚集到庚烷中大规模致密团簇的转变，为溶剂依赖性聚集提供了物理动机明确的分子解释，并为通过溶剂环境或添加剂设计控制沥青质缔合提供了机制路线。研究局限在于仅考虑单一代表型岛式沥青质在纯溶剂中的行为，未来工作应扩展至多分散沥青质库和多组分SARA环境，并借助机器学习粗粒化模型探索更大体系和更长时标。