聚合物溶液冷冻在低温保存、多孔材料合成和环境修复等过程中扮演着核心角色。与简单溶质不同，聚合物通过链排除、界面吸附和受限效应影响凝固过程，这些相互作用塑造了凝固前沿的结构和稳定性，导致溶质梯度、相分离或形态不稳定性。
聚合物-冰界面相互作用
当冰在聚合物溶液中形成时，前进的固-液界面会遇到具有独特物理性质的介质：长链分子、熵弹性和溶质驱动的热力学梯度。聚合物通常被排除在冰相之外，形成浓度边界层，该层影响生长前沿的形态和速度。
聚合物排除与界面富集
随着冰前沿向聚合物溶液推进，溶解的链由于其大尺寸、低扩散率和不利的焓相互作用而被排除在生长的晶格之外。这种排除导致前沿前方形成聚合物富集边界层。该层的厚度和形状由前沿运动与溶质扩散能力之间的竞争决定。这种平衡由佩克莱特数（Pe）描述，它控制着积累层的结构和浓度梯度的严重程度。当Pe ≪ 1时，扩散占主导，溶质在较宽距离内重新分布。当Pe ≫ 1时，前沿速度超过扩散速度，导致尖锐、被压缩的浓度梯度。在聚合物溶液中，扩散系数D并非常数。随着聚合物链的积累，局部粘度增加，扩散减慢，形成一个反馈回路，放大了前沿附近的拥挤效应。这种效应可以通过标度律来捕捉，其中D0是无限稀释下的扩散系数，c是重叠浓度，α是根据聚合物化学性质在1到2之间变化的参数。当局部浓度c接近或超过c时，扩散变得越来越受阻，有效佩克莱特数在冷冻过程中动态上升。
这种浓度积累不仅改变输运性质，还影响热力学和机械性质。高的局部聚合物含量会降低冰点，增加渗透压，甚至可能在溶液转变为玻璃态或相分离状态时导致玻璃化。如果前沿推进足够慢，边界层可能达到缠结阈值ce，此时链形成瞬态弹性网络，抵抗变形并进一步减慢再分布。
因此，冷冻过程中的聚合物排除是由浓度、粘度和扩散率的耦合梯度控制的。即使聚合物化学性质或冷冻速率的微小差异也可能导致显著不同的前沿行为——从平滑分层到不稳定、高度异质的凝固。这些界面积累动力学为后面讨论的形态转变奠定了基础。
界面吸附与钉扎
除了简单的排除，许多聚合物表现出对冰界面的瞬态或部分吸附。这些相互作用通常由氢键、偶极力或表面水合结构介导。即使是弱吸附聚合物也能通过局部扰动界面能和抑制结晶过程中的分子排列，为冰传播创造能量障碍。
分子动力学模拟表明，聚乙烯醇（PVA）在棱柱形冰面上表现出选择性吸附，吸附事件通过羟基与冰晶格排列而稳定。然而，这种吸附不是连续的：由于动力学限制，表面覆盖进行缓慢，导致沿前沿结合链的异质分布。这些斑块状域充当钉扎点，改变台阶流动并在纳米尺度上使界面粗糙化。
吸附的物理图景可以通过朗缪尔型模型来构建，其中表面覆盖度θ随聚合物浓度c增加。然而，在冷冻系统中，吸附传播动力学和结合位点的空间排列比平衡饱和起着更关键的作用。模拟和表面力测量表明，吸附的链可以跨越多个表面位点，产生多价锚定点，抵抗前沿推进。
与抗冻蛋白（AFPs）的比较提供了一个有用的极限情况。AFPs表现出不可逆吸附和强钉扎，通过结合到特定晶格图案完全阻止冰生长。相比之下，像PVA或PEG这样的聚合物表现出可逆的、迁移率受限的相互作用，允许部分重排但仍引起局部过冷。原子力显微镜（AFM）研究证实了这些区别：AFPs形成有序的、固定的层，而聚合物形成无序的、松散吸附的团簇，具有有限的停留时间。
吸附也改变了界面能。局部表面能γsl的变化通过吉布斯-汤姆逊关系影响曲率和冰点，突出的链或结合层的存在可以稳定凸起或钉扎台阶边缘。这些微观结构效应在定向冷冻下将前沿形态从刻面状转变为胞状，尤其是在低聚合物浓度下，吸附主导于体相积累。
总之，界面吸附提供了第二种机制，不同于溶质排除，聚合物通过它影响冰前沿动力学。其在冷冻早期阶段或低浓度区域的影响最强，即使部分覆盖也能引发形态转变并引导凝固路径。
链柔性与受限的作用
聚合物在冷冻前沿附近的行为深受其链柔性和局部排除层内受限程度的影响。这些因素控制着聚合物如何积累、取向和再分布，塑造了浓度梯度和界面动力学，为后续的形态演化奠定了基础。
柔性聚合物，如聚乙二醇（PEG）或聚乙烯醇（PVA），具有低的持续长度，通常在稀溶液中采取卷曲构型。随着冷冻进行，被排除的聚合物在前沿前方积累，其中变窄的边界层压缩了它们的构象空间。这种压缩施加了熵惩罚，表现为抵抗进一步致密化的排斥力。受限自由能可以估算。结果是界面附近出现高渗透压和空间阻力区。
当局部浓度超过重叠阈值c*时，链开始相互渗透，形成具有 dramatically 增加粘度的缠结网络。超过缠结浓度ce后，长程链约束减慢扩散并捕获异质性。这些粘弹性变化反馈到溶质再分布和界面迁移率，强化梯度并可能诱发局部不稳定性。
聚合物链的持续长度lp进一步调节这些相互作用。半刚性聚合物，包括壳聚糖和纤维素衍生物，与界面各向异性地相互作用，通常抵抗压缩并采取拉伸或排列的构型。AFM和低温成像研究表明，较硬的链可以形成有序域或跨越粗糙冰特征的桥接触点，局部改变表面能并促进界面钉扎。与变形以最小化熵成本的柔性线圈不同，刚性链约束了排除层的几何形状，放大了曲率效应并改变了局部过冷模式。
柔性也影响动力学捕获。具有短弛豫时间的聚合物更容易重排，因此被更有效地排除，而较硬或已经缠结的链可能由于其有限的构象响应时间而被部分捕获。这导致冻结相内的组成空间变化，其富集或耗尽区由链迁移率和前沿速度之间的相互作用决定。
总之，链柔性和局部尺度的受限调节了冰-溶液界面处排除层的结构、动力学和流变学。通过影响链构象、取向和弛豫行为，这些参数控制了更大尺度前沿形态和冷冻过程中形成的孔结构的最初条件。
热力学意义
聚合物富集边界层通过依数性和曲率驱动效应改变局部冷冻条件。随着溶质在前沿附近积累，由于渗透压和界面曲率，局部冰点被降低。对于稀溶液，冰点降低可以用弗洛里-哈金斯风格的表达式近似，其中ΔTf是冰点降低，R是气体常数，Tf是纯溶剂的冰点，ΔHf是熔化焓，Φ是聚合物体积分数，Mn是数均分子量。该表达式捕捉了依数效应对聚合物负载量和链尺寸的依赖性。
除了浓度驱动的冰点降低，冷冻界面的几何形状通过曲率效应起作用。当冰前沿变得粗糙或胞状时，局部曲率半径减小，通过吉布斯-汤姆逊效应放大冰点抑制，其中γsl是固-液界面能，Tm是体相熔点，ρs是固相密度，Lf是熔化潜热，r是曲率半径。
吸附聚合物主要通过改变局部表面张力和稳定表面粗糙度来影响界面能，而不是直接改变曲率。这些变化影响了沿前沿的过冷分布，并能改变控制台阶成核或小面传播的能量平衡。曲率较低或聚合物积累较少的区域优先冻结，导致前沿速度的空间变化，并在定向凝固下促进形态不稳定性。
总之，这些效应——依数性降低和曲率放大——协同作用，决定了在聚合物富集环境中冷冻在何处以及如何引发和传播。值得注意的是，这些热力学贡献的相对重要性取决于初始聚合物浓度和施加的热驱动力，这也控制了聚合物富集边界层的输运和流变状态。
溶质积累与再分布
当冰前沿穿过聚合物溶液推进时，溶质通常被排除在固相之外，并在界面前方的液体中积累。这种积累产生陡峭的浓度梯度，改变局部输运性质，并可能驱动次级现象，如微相分离、粘度增加和次级溶质的结晶。冷冻过程中溶质的分布在决定未冻结区域的组成和结构以及冻结相的形态和均匀性方面起着核心作用。
浓度边界层与非线性扩散
聚合物在冷冻前沿前方的积累形成了一个浓度边界层，其结构和演化在前沿形态中起着核心作用。如2.1节所述，这种积累是由前沿推进和聚合物扩散之间的不平衡引起的，由佩克莱特数描述。当Pe > 1时，平流主导扩散，导致前沿前方局部聚合物浓度出现陡峭梯度。
这些梯度不是静态的。在聚合物溶液中，扩散系数D随浓度急剧下降，遵循方程（2）中的标度律。随着聚合物的积累，D下降，这进一步加剧了梯度，并导致前沿遇到粘度越来越高、扩散性越来越差的物质。这种正反馈回路增强了溶质保留，并导致动态变窄的积累层。
这种反馈的形态学意义是显著的。随着D(c)下降，有效佩克莱特数在冷冻过程中上升，扩大了前沿前方的组成过冷区域，使其更容易发生形态不稳定性。在定向冷冻中，这表现为更早出现胞状或树枝状生长，以及更小的层状间距。高分子量聚合物具有更陡峭的D(c)衰减，在较低浓度下就表现出这些效应。
最近使用PVA溶液定向冷冻的实验证实了这种浓度驱动的锐化效应。这些结果强调，演化的浓度边界层及其引起的非线性输运是形态不稳定性的核心驱动因素。
溶质捕获与分配
在凝固过程中，溶质捕获指的是尽管存在热力学排除，溶质分子仍被纳入固相。在聚合物溶液中，这通常发生在冷冻前沿移动速度足够快，以至于聚合物链无法及时扩散 away 的情况下。在这些条件下，一部分聚合物被捕获在固体冰基质内，导致非平衡分配。
这种行为通常使用锐界面模型来描述，该模型定义了跨不连续界面的分配系数k = Cs/Cl，或连续捕获模型，该模型允许由局部前沿速度和溶质迁移率控制的平滑浓度梯度。在缓慢、接近平衡的冷冻极限下，高分子量聚合物经历强烈的熵和空间排除，分配系数趋近于零。
然而，在较高的前沿速度下，一些研究报告即使对于大聚合物，k值也有适度增加。这种明显的增加并不反映排除机制的失效，而是动力学捕获的开始：当界面推进速度超过聚合物链扩散或弛豫的速度时，溶质在排除完成之前就被吞没了。在这种状态下，分配变得受界面速度与聚合物扩散率之比的控制，而不是单独由平衡热力学控制。
这种解释在连续溶质捕获模型中被形式化，例如Aziz-Kaplan框架，它将分配系数表示为前沿速度的平滑函数，其中k0是平衡分配系数，VD= D/λ是特征扩散速度，D是溶质扩散率，λ是微观界面宽度。虽然最初是为金属合金开发的，但这个框架已被广泛适用于软物质和聚合物系统，以捕捉快速冷冻过程中的非平衡溶质掺入。
根据这种观点，强排除（k ≈ 0）和速度诱导捕获（k > 0）代表了同一基本过程的两种极限情况，由前沿传播和聚合物输运之间的时间尺度分离程度区分。因此，广泛共识认为高分子量聚合物在热力学上仍然被排除在冰之外，而在高冷冻速率下观察到的有限分配反映了动力学捕获，而不是真正掺入冰晶格。
在多组分聚合物溶液中，差异分配变得更加复杂。较小的添加剂（例如，冷冻保护剂或盐）可能比较长的聚合物链优先被捕获，导致局部渗透梯度或浓度不均匀性的形成。这种效应在水性双相系统中尤其相关，其中聚合物和盐之间的热力学相分离可能在冷冻过程中发生，导致具有不同溶质组成和冰点的微区。
溶质捕获也影响固相的机械和热性质。被捕获的聚合物会干扰晶体排列，诱导非晶包裹体，或降低结晶度。例如，Li等人表明，蔗糖溶液由于溶质保留和水活度降低而表现出对冰前沿的树枝状抑制。类似地，我们最近关于PVA液滴冷冻的研究表明，前沿形态——以及冷冻液滴顶部尖点最终消失——受聚合物浓度和冷冻速率的调节。我们表明，超过重叠浓度c*的溶质保留影响了热行为和固体形状的演化。
这些综合效应——溶质排斥、部分捕获和差异分配——也可能在未冻结区域诱导渗透应力，导致相分离或壳状结构。冷冻液滴的实验切片显示，当冷冻从底部向上进行时，聚合物浓度在冷冻的底部区域较低，在冷冻的顶部区域较高。这种梯度是由前沿推进过程中的动力学捕获和定向溶质位移引起的。这种不对称性在预测冷冻后的形态和输运时至关重要，特别是在冷冻保存、食品科学和冷冻铸造中的应用。
微相分离与动态异质性
当界面附近的聚合物浓度超过溶解度极限或跨越玻璃化转变阈值时，未冻结溶液可能发生微相分离。当局部聚合物浓度变得足够高，以至于未冻结相的玻璃化转变温度Tg升高到高于环境或系统零下温度时，冷冻前沿附近就会发生玻璃化。这种转变阻止了链的迁移率，并可能导致形成玻璃状的、聚合物富集的域，与生长的冰晶体相邻。这种玻璃化层限制了进一步的溶质再分布，并可以捕获浓度梯度，保留早期凝固过程中建立的异质性。Tg对聚合物化学性质、分子量和局部水合高度敏感，因此为预测系统何时从未冻结相的流体状态转变为类固态提供了一个有用的阈值。
在冷冻的PVA-壳聚糖支架中观察到的玻璃化壳聚糖富集区域证实了在冷冻过程中超过局部玻璃化转变阈值可以抑制孔洞合并并捕获浓度异质性，这与预测的Tg依赖性输运停止一致。这些次级结构通常在冷冻后持续存在，并可以模板化最终的材料形态，特别是在涉及冷冻干燥或凝固后交联的应用中。
为了说明这一进程及其空间异质性，图4展示了沿冷冻方向的聚合物浓度分布示意图。远离前沿时，溶液保持在稀区域（c < c*），但随着界面推进，溶质排除导致陡峭的浓度梯度，前沿附近的区域超过缠结阈值（ce），甚至当T < Tg(c)时跨越玻璃化线。这种演变的分布将局部热力学状态与迁移率、粘弹性以及最终印刻在冷冻材料中的形态联系起来。
冷冻前沿形态与稳定性
凝固前沿的形态在决定冷冻聚合物溶液的空间组织和异质性方面起着核心作用。随着冰前沿的推进，其形状和稳定性受热梯度、溶质排斥、界面吸附和周围溶液输运性质之间相互作用的控制。这些因素共同定义了前沿是保持平面状还是变得不稳定，产生胞状、层状或树枝状生长图案。
平面、胞状和树枝状前沿
聚合物溶液中从平面到胞状或树枝状冰前沿的转变是溶质积累驱动的形态不稳定性的一个标志。当排除层在前沿前方积聚时，局部过冷度增加，直到扰动被放大。
这种转变由Mullins-Sekerka分析很好地描述，该分析预测，当界面前方液体中的热梯度低于由溶质积累决定的临界值时，平面前沿变得不稳定，其中m是液相线斜率，C0是初始溶质浓度，k是平衡分配系数，D是溶质扩散率。这个表达式捕捉了降低的扩散率或增加的溶质排斥——这两者在聚合物溶液中都很典型——如何通过放大温度或浓度场中的扰动来 destabilize 界面。
在盐溶液中，这种经典图景通常是一阶的良好描述：离子保持高迁移率，分配主要由依数性热力学控制，胞状或糊状形态通常可以使用带有界面修正的扩散控制模型来合理化。相比之下，聚合物溶液可能违反这些假设，因为积累层可能变得高度粘稠甚至玻璃化，并且聚合物-冰相互作用（吸附、钉扎、结合水）可以直接扰动局部生长动力学。
这些差异在NaCl和PVA溶液在单向冷冻下的直接比较中很明显。虽然NaCl显示出与Mullins-Sekerka型行为一致的局部扰动，但PVA从开始就表现出全局前沿失稳。这种行为归因于聚合物链的低扩散率和强吸附，它们迅速产生陡峭的浓度梯度和微尺度钉扎点。除了前沿形状，溶
