液相色谱的多尺度模拟：流动相流速如何调节随机与有序多孔颗粒床层中的板高—保留因子关系

《Journal of Chromatography A》：Multiscale simulation of liquid chromatography: How the mobile-phase flow velocity modulates the plate height–retention factor relationship in random and ordered beds of porous particles

研究人员开展了液相色谱（LC）的多尺度模拟研究，探讨了流动相平均流速（uav）如何调节总板高（H）与保留因子（k）之间的关系。研究发现，H对k的依赖性随流速变化显著：在低流速下，H随k的增加而增加；而在高流速下，H则随k的增加而减小。在中等保留范围内（k约为1至10），板高可降低40%（从5降至3微米），这表明分离效率得到了实质性提升。此外，研究人员发现，在高流速下，H-k曲线的下降趋势反映了床层异质性的程度。研究还表明，分析物分子被多孔颗粒摄取和释放的过程会中断涡流扩散，增加保留效应从而增强这种效果，进而减少谱带展宽。这些发现揭示了多孔颗粒床层中复杂的质量传递机制，并为优化液相色谱分离条件提供了理论依据。

研究背景与意义

高效液相色谱（HPLC）柱的分离效率通常由板高（H）与平均流动相流速（u_av）的关系来表征。然而，分析物保留因子（k）对板高的精确影响，特别是对随机填充床中涡流扩散（eddy dispersion）贡献的影响，仍未得到充分理解。传统的H-u_av曲线分析忽略了流速对H-k关系的系统性调制作用。本研究旨在通过多尺度模拟，揭示不同流速下随机与有序多孔颗粒床层中H-k关系的演变规律，阐明保留因子在色谱动力学中的微观机制，为色谱柱填充质量评估和分离条件优化提供理论支持。

关键技术方法

研究人员采用自下而上的多尺度模拟方法。首先，通过分子动力学（MD）模拟建立反相液相色谱（RPLC）狭缝孔模型，获取分析物密度和表面平行扩散系数分布。其次，将这些分布嵌入到介孔空间重构模型中，利用随机游走模拟计算有效介孔扩散系数（D_meso）。随后，在无流动条件下进行随机游走模拟以获得有效床层扩散系数（D_bed）。最后，利用格子玻尔兹曼方法（LBM）计算流场，并结合随机游走粒子追踪法模拟有流动条件下的纵向弥散过程，提取长时间渐近纵向弥散系数（D_L），进而计算总板高。研究对比了物理重构的随机填充床与计算机生成的面心立方（fcc）有序床两种形态。

研究结果

3.1 无流动条件下的多尺度扩散模拟

研究结果表明，D_meso是乙腈（ACN）体积分数（即保留因子）的函数，且随着k的增加而单调下降。当按各自本体扩散系数（D_bulk）归一化后，D_meso/D_bulk的值始终小于1，这推翻了RPLC中D_meso大于D_bulk的传统假设。在床层尺度上，随机填充床和fcc结构的D_bed差异很小（最大约5-10%），主要受比表面积影响，且D_bed随k增加而降低。

3.2 H-k关系的流速依赖性

研究发现，流速决定了H-k关系的基本形态。在低流速下，传质受纵向扩散项（B-term）主导，H随k增加而单调增加。在中等流速下，fcc结构的H随k线性增加，而随机填充床的H在k约为0.5至10范围内基本保持不变。在高流速下，fcc结构的H趋于平稳，而随机填充床的H随k增加而显著下降，并在k约为0.5至1之间出现明显的最大值。这表明在高流速下，增加保留因子能显著提高随机填充床的分离效率。

3.3 涡流扩散板高贡献的分析

通过扣除B项和固定相区传质阻力项（H_Cs），研究人员分离出涡流扩散项（H_A）。结果显示，H_A是随机填充床中H随k下降的主要原因。在高流速下，随机填充床的H_A在k约为0.5至1之间出现峰值，随后在k约为1至10范围内下降45%至46%，而fcc结构的H_A则保持较低且平坦的曲线。这归因于随机填充床特有的短程通道间（short-range interchannel）微结构异质性。

3.4 多孔颗粒床层H-k关系背后的物理学机制

研究揭示了高流速下H随k下降的物理机制。一方面，随着保留增加，更多的分析物分子从发生涡流扩散的流动区域（mobile zone）转移到纯粹的扩散传输的固定相区域（stationary zone），减少了涡流扩散的程度。另一方面，多孔颗粒对分析物分子的反复摄取和释放充当了“黑盒混合器”，中断了分析物在流动路径中的速度记忆积累，加速了纵向速度自协方差的衰减，从而降低了涡流扩散。

讨论与结论翻译

研究人员得出结论，低流速下的谱带展宽由B项主导，而高流速下的展宽则由随机填充床中的涡流扩散主导。短程通道间涡流扩散与颗粒内传输的耦合导致了弱保留条件下的H-k关系最大值。此外，H-k曲线在高流速下的下降趋势反映了实际的涡流扩散程度和填充床层异质性程度。这一机制不仅适用于RPLC，也是HPLC中普遍存在的物理现象。

本研究的重要结论包括：传统RPLC中关于D_meso> D_bulk的假设是错误的；在k约为1至10的实际相关范围内，高流速下H随k增加而降低，这源于分析物分子向固定相区的重分布及其对涡流扩散记忆效应的中断；k约为1至10区域内H-k曲线的斜率可作为评估色谱柱填充质量的互补指标，斜率越陡表明填充均一性越差。