聚合物中慢阿伦尼乌斯过程（SAP）的建模

《Soft Matter》：Modeling the slow Arrhenius process (SAP) in polymers

　　

非晶玻璃形成聚合物表现出多种弛豫过程，包括与玻璃化转变相关的结构α-弛豫以及通常遵循阿伦尼乌斯（Arrhenius）行为的较快次级弛豫。最近，在频率远低于α-过程的区域观察到了一种独特的慢阿伦尼乌斯过程（SAP）。尽管其温度依赖关系符合阿伦尼乌斯行为，但SAP涉及更长的弛豫时间，其微观起源仍不清楚。在此，研究人员扩展了双态、双时间尺度（TS2）理论，以在一个统一框架内描述α-弛豫和SAP。研究人员提出，SAP代表了动态关联簇的粗粒化流体中类α过程的温度高限。通过重整化的相互作用能和配位数参数，同一模型在没有额外可调参数的情况下定量复现了多种聚合物的α和SAP数据，并解释了观察到的迈耶-内尔德尔（Meyer-Neldel）补偿行为。该理论进一步预测，在足够低的温度下，SAP应偏离阿伦尼乌斯行为，转变为类Vogel-Fulcher-Tammann-Hesse（VFTH）动力学，从而为了解玻璃形成聚合物中的簇尺度弛豫提供了一个物理图像清晰的诠释。

论文解读：聚合物中慢阿伦尼乌斯过程（SAP）的建模

该研究论文发表于《Soft Matter》期刊，主要围绕非晶玻璃形成聚合物中近期发现的慢阿伦尼乌斯过程（Slow Arrhenius Process, SAP）展开理论建模与机理解释。目前，对于玻璃形成聚合物的弛豫行为，学界普遍认可存在与玻璃化转变相关的结构α-弛豫（α-relaxation）以及通常遵循阿伦尼乌斯行为的较快次级弛豫（如Johari-Goldstein β-弛豫）。然而，近年来的宽带介电谱（BDS）等实验技术在频率远低于α-过程的区域观察到了一种独特的SAP，其弛豫时间更长且温度依赖呈阿伦尼乌斯型，但其微观物理起源尚不明晰。为了深入理解这一普遍存在的动力学过程及其与经典α-弛豫的内在联系，研究人员开展了一项理论扩展工作，旨在将SAP纳入统一的弛豫理论框架中进行描述。

为开展此项研究，研究人员主要采用了以下关键技术方法：基于先前发展的两态、双时间尺度（Two-State, Two-Timescale, TS2）理论框架，引入动态关联簇的粗粒化（Coarse-Graining）思想，将SAP视为粗粒化流体（由簇构成的“粗粒化粒子”CGP流体）在高温下的αβ-过程极限。研究人员利用来自Thoms和Napolitano等人的13种聚合物的公开实验弛豫数据（包括P4ClS、P4MS、PBzMA、PC、PEMA、PMMA、PNPMA、PS、PTBMA、PTBS、PTBuA、PVAc、PIB），通过TS2模型的参数回归（如Excel广义简约梯度法GRG）以及阿伦尼乌斯温度（T_A）和阿伦尼乌斯时间（τ_A）的定义，在无额外可调参数的情况下对α-弛豫和SAP数据进行拟合与对比分析，并探讨了迈耶-内尔德尔（Meyer-Neldel）补偿行为的起源。

研究结果主要包括以下几个部分：

1. 1.
   引言（Introduction）
   研究人员在引言中梳理了SAP的实验发现背景，指出SAP不仅在介电谱中被观测到，还被多种互补实验技术和模拟所报道。SAP的活化能与指前因子之间通过“熵-焓补偿”即迈耶-内尔德尔规则相关联。现有理论如White、Napolitano和Lipson（WNL）的集体小位移（CSD）模型虽能描述SAP的一些特征，但α-弛豫与SAP通常由不同模型描述。研究人员提出本研究的核心假设：SAP与大簇的动力学相关，这些簇可被视为相关时间窗口内的“粗粒化粒子”（CGP），SAP对应于粗粒化材料（CGM）的αβ-过程的高温热限，且CGM的玻璃-液体热力学转变温度略低于材料本身。
2. 2.
   模型（The model）
   2.1 通用TS2建模与α-转变的参数化（The general TS2 modeling and parameterization of the α-transition）
   研究人员简要回顾了TS2理论的基本假设：非晶材料中存在“域”（可理解为协同重排区CRR），每个域可处于“液态”或“固态”。快速Johari-Goldstein（JG）β-弛豫时间由阿伦尼乌斯函数给出，慢α-弛豫时间则为超阿伦尼乌斯函数（如Ginzburg在TS2模型中提出的表达式）。TS2模型原本有五个参数，但其中三个在适当无量纲化后是通用的，仅两个（如热力学玻璃-液体转变温度T_X和基本材料时间τ_el）具有材料依赖性。为方便，研究人员引入阿伦尼乌斯温度T_A（α-弛豫时间偏离高温阿伦尼乌斯渐近线5%处的温度）和阿伦尼乌斯时间τ_A（T=T_A时的α-弛豫时间）来重新表述方程，且该比值T_A/T_X等材料无关。

2.2 TS2模型在慢阿伦尼乌斯过程中的应用（Application of the TS2 model to slow Arrhenius process）

研究人员将SAP理解为粗粒化流体（CGP流体）中的弛豫过程。在分子水平上考虑单体和结合能ε_α及配位数Z_α；经过粗粒化，一个域视为一个粗粒化粒子（CGP），CGP流体的配位数为Z_SAP。研究人员假设在CGP流体中，SAP有效地替代了高温α-弛豫，且在T > T_A,SAP时，SAP活化能对应E_1,SAP，转变温度T_0,SAP对应相关参数，阿伦尼乌斯截矩对应log(τ_el,SAP)。研究人员预测若在低温下继续测量，SAP最终可能变得非阿伦尼乌斯。

慢阿伦尼乌斯过程的迈耶-内尔德尔补偿定律（The Meyer–Neldel compensation law for the slow Arrhenius process）

研究人员指出SAP数据满足迈耶-内尔德尔补偿律，即ln(τ_A,SAP) = A - kE_a,SAP。在TS2结合Sanchez-Lacombe（SL）状态方程框架下，活化能与相互作用能ε_SS和配位数Z有关：E_a,SAP∝ aZ_SAPε_SAP。由此重新写为ln(τ_A,SAP) = A - bN_SAP，其中N_SAP与CGP流体中阻碍空穴扩散的键数目有关。研究人员解释：未实现的潜在键（N_max- N_SAP）迫使扩散空穴花费更多时间，较低阿伦尼乌斯时间下限约为10^-10s。

1. 3.
   结果与讨论（Results and discussion）
   3.1 结果（Results）
   研究人员汇编了基于Thoms和Napolitano数据的α-过程和SAP的TS2回归参数（表1，涉及13种聚合物）。除PIB和PNPMA外，均有T_A,SAP< T_A,α且log(τ_A,SAP) > log(τ_A,α)。以聚苯乙烯（PS）为例，SAP曲线在高温下是阿伦尼乌斯的，近乎平行于TS2 JG线；模型预测SAP线在约250 K（1000/T = 4.0，低于T_g约100 K）附近会从阿伦尼乌斯转变为VFTH行为，“SAP玻璃化转变温度”估计约为337 K。对于较强的聚异丁烯（PIB，脆性m≈35），α-转变本身仅在近T_g时偏离阿伦尼乌斯，α和SAP曲线近乎平行，两者差异较其α-过程曲线差异小。其余11种聚合物的数据也表明SAP模型曲线在无可调参数下生成。图5显示了α-和SAP的ln(τ_A) vs T_A，SAP满足迈耶-内尔德尔补偿的直线关系，α-过程则无。图6显示T_A,SAP与T_A,α及T_g,SAP与T_g,α的关系，由于τ_A,SAP更高推高T_g,SAP而T_A,SAP更低推低它，部分抵消致T_g,SAP≈T_g,α。研究人员估计了动态减慢参数λ≈1.52，并计算了α-过程和SAP的有效配位数Z和相互作用能ε。图7显示α-过程对应链段/分子弛豫，ε_α具材料特异性（1.8至3.4 kJ mol^-1），Z_α约3.6至5.5；SAP对应域的弛豫，ε_SAP几乎普适（≈3.7 kJ mol^-1），Z_SAP较小（约1.8至3.2）。

3.2 讨论（Discussion）

研究人员提出SAP的新诠释：通过TS2形式化，将其视为“簇”的高温αβ-过程；粗粒化流体（CGP流体）具有比材料本身更低的热力学转变温度T_x和更低脆性。因此实验中SAP呈阿伦尼乌斯过程；预测若在更低温度和频率测量同一峰，可能看到类VFTH的对阿伦尼乌斯偏离。研究人员探讨CGP流体热力学与材料本身的差异：通用SL-TS2的两个参数可 reinterpret 为ε和Z；对于粗粒化流体（SAP），范德华参数相当恒定（与Lipson等发现一致），其普适性及是否限于聚合物/小分子或见于其他材料尚待研究。研究人员指出阿伦尼乌斯时间下限约10^-10s（高脆性材料），而“绝对”下限约10^-13–10^-14s；因此脆弱材料的阿伦尼乌斯区域较窄（3-4数量级），强材料较宽（6-15数量级），SAP属后者。低脆性可能与低“活跃”或“强”接触数有关。研究人员也讨论了方法潜在局限：从阿伦尼乌斯方程“映射”到通用TS2方程虽总是可行，但若预测的阿伦尼乌斯温度过高则会失败（如PBVCl的SAP、水的α-过程）。SAP弛豫时间的潜在超阿伦尼乌斯行为在低温下可能不遵从通用TS2或甚至VFTH，若其他相互作用变得重要；或者DMA、BDS或吸附测量的SAP弛豫时间在极低温继续服从阿伦尼乌斯甚至变成亚阿伦尼乌斯（由于微观与宏观动力学解耦）。总体框架为综合建模SAP及由其决定的现象（如聚合物吸附）提供新途径，未来工作应包括进一步阐明这种“自然粗粒化”的分子机制及产生SAP的“簇”形状和大小，以及SAP低温延伸更像JG β-转变（继续阿伦尼乌斯）还是 bifurcate 出CGM的α-和β-转变。

3.3 簇形成和时间尺度分离的物理起源（Physical origin of cluster formation and timescale separation）

研究人员强调此处粗粒化簇不应解释为永久结构聚集体或热力学不同相，而是玻璃形成系统中弛豫时间尺度分离产生的动态涌现实体。当温度接近玻璃化转变，动力学变得空间异质且关联在有限长度尺度上，这些关联区常指协同重排区（CRR）。在TS2图像中，α-弛豫对应此类域的集体重排。但在τ_β? t ? τ_SAP时间窗口，新的有效动力学层级浮现：域内运动与弛豫已平衡，而域间重排（关联SAP）尚未发生，这些关联域组表现为有效刚性单元。时间尺度分离自然导致层级粗粒化：短时间相关自由度是分子链段；中间时间是域；更长是域的簇。SAP可解释为这种重整化粗粒化流体的结构弛豫，即SAP不是独立弛豫机制，而是其基元单元被底层能景异质性动态重定义的系统的类α弛豫。动态关联簇的形成可理解为集体约束和连接性的结果；温度降低，域间有效强接触数增加，导致瞬态连接的域网络，其集体重排需许多单元协同位移，重排活化能垒随有效接触数缩放，这自然解释了SAP的大阿伦尼乌斯活化能。由于粗粒化水平的有效相互作用能在所检聚合物系统内近似材料无关，遂使阿伦尼乌斯参数表现出迈耶-内尔德尔补偿行为。簇由动力学准则操作性地定义，其大小和寿命由域内弛豫和域间重排竞争决定，故取决于温度和观测时间尺度；该观点将SAP置于无序材料层级弛豫更宽框架中，并支持玻璃态动力学源于连续动态粗粒化层级的看法。

1. 4.
   结论（Conclusions）
   研究人员利用“通用”SL-TS2理论框架分析了慢阿伦尼乌斯过程（SAP），提出SAP代表了“簇”粗粒化流体的α-过程高温区；该粗粒化流体相较材料本身具有更低脆性和更低转变温度。基于此分析，研究人员预测“慢阿伦尼乌斯过程”在低于玻璃化转变温度的温度下会显示出向上的、类VFTH的偏离阿伦尼乌斯行为。研究人员还提出了为SAP（而非α-过程）观察到的熵-焓补偿的机制——粗粒化范德华能对材料特性的独立性。理解这种独立性的分子机制及潜在违反情形将是未来研究的课题。