把分子"看得见"：一文读懂 HS-PPM 算法
从硬球模型到拟颗粒并行模拟，连接微观碰撞与工程尺度传递过程
如果想理解气体在催化剂颗粒等多孔介质材料中的扩散、吸附以及反应过程，单靠传统宏观模型往往不够。基于连续介质假设的CFD方法可以描述反应器与颗粒尺度上的流场、浓度场和温度场，但只能把分子碰撞、孔道受限扩散、吸附位点作用等微观细节平均化；传统分子动力学虽然能够追踪原子/分子的运动轨迹，但计算成本高，模拟尺度和时间跨度都受到明显限制。HS-PPM算法 正是为连接这两类尺度而生：它用硬球碰撞描述基础微观运动，同时引入拟颗粒方法提升并行效率，让"看得见的分子运动"和"算得动的大体系"之间有了一个实用桥梁。
一、HS-PPM是什么？
HS-PPM，Hard-Sphere/Pseudo-Particle Modeling的缩写，即硬球-拟颗粒模型。所谓"硬球"，是指把真实分子间复杂的势能作用，近似为粒子之间的自由飞行和遵循动量守恒、能量守恒的瞬时弹性碰撞；所谓"拟颗粒"，则强调用更适合并行计算的时间推进和区域交换机制来组织这些粒子的运动。这样做不会试图复现所有电子结构细节，而是抓住稀薄气体、孔道扩散和传递过程中的关键几何约束与碰撞统计。
图1 HS-PPM耦合算法框架图：PP子区像河流，HS子区像岛屿
从算法框架看，HS-PPM并不是简单把整个模拟区域"一刀切"。它会按照计算模型和并行通信需要，把体系划分为若干子区：在大部分内部区域，采用硬球事件驱动模型以保持较高的碰撞精度；在进程边界和耦合区域，采用PPM进行同步更新和信息交换，这是理解HS-PPM可扩展性的关键。
二、为什么不用传统MD直接算？
传统分子动力学能够给出原子级轨迹，但每一步都要计算相互作用势，体系一大、时间一长，计算成本就会迅速上升。对于稀薄气体扩散、多孔介质传质、扩散-反应耦合等问题，研究者常常更关心的是粒子如何在复杂几何中迁移、碰撞、被吸附或参与反应，而不是每一个短程相互作用的精细势能曲线。HS-PPM用硬球有效直径替代具有连续排斥特征的软球势函数相互作用，把核心计算集中在碰撞事件（包括吸附、催化反应、脱附）的处理上，因此更适合大规模、长时间的介尺度模拟。
图2 硬球碰撞事件示意：粒子自由飞行与瞬时碰撞更新
硬球思想的优势在于清晰：粒子在两次碰撞之间做直线运动，只有当距离达到碰撞条件时才更新速度。这种处理方式天然适合描述稀薄气体中的碰撞控制过程。与此同时，PPM部分让算法避免完全依赖全局事件列表，在边界区域按固定时间步推进，从而更容易与并行计算、复杂边界和多区域耦合结合起来。
三、算法核心：精度与并行的折中
HS-PPM最重要的设计不是"简化"，而是"把复杂性放在合适的位置"。在碰撞频繁、需要准确处理粒子相互作用的区域，HS模型负责事件判断和速度更新；在跨区域通信、边界缓冲和同步推进的位置，PPM负责把局部计算组织成可并行的大系统。这样既保留了硬球分子动力学对碰撞过程的物理表达，又避免了纯事件驱动方法在大规模并行时的通信瓶颈。
图3 HS-PPM在边界区域耦合HS与PPM的计算流程
四、它能算什么？从扩散到吸附反应耦合
HS-PPM最早被用于气体流动与传递性质模拟，随后逐渐进入多孔材料、催化剂孔道、吸附剂颗粒和多组分气体扩散研究。
图4 HS-PPM计算气体自扩散系数并与经典理论/实验对照
在多孔催化剂中，HS-PPM可以把"扩散"和"反应"放进同一个孔道模型中统一考察。相关反应-扩散耦合研究引入了扩散因子D、反应因子R和单位时间产率Y，用来刻画孔道形状、支路数量与反应活性之间的相互制约关系：支路增多会带来更多反应表面，但也可能让分子迁移路径变得更曲折、削弱扩散通畅性。因此，最优孔结构并不等同于"表面积越大越好"，而往往取决于扩散能力与反应能力之间的平衡。
图5 基于HS-PPM算法的催化剂孔道内的反应-扩散耦合模拟：获取实验难以观测的微观动力学演化过程
对于真实催化剂失活问题，HS-PPM还可用于观察焦炭沉积如何改变分子通行路径。甲烷在ZSM-5中的扩散模拟显示，均匀分布的少量焦炭影响相对有限；但当焦炭在部分T位点附近集中并逐步堵塞孔道时，自扩散系数会明显下降。这说明HS-PPM不仅能给出扩散快慢，还能把"哪里被堵住、哪条方向更难走"转化为可统计的结构效应。
图6 基于HS-PPM算法观测催化剂积碳失活
图7 基于HS-PPM用于扩散-吸附耦合的吸附剂模型
五、HS-PPM的价值在哪里？
对化工、材料和能源系统而言，HS-PPM的价值在于它提供了一种"中间尺度"的计算语言。它比宏观经验模型更微观，能看到粒子如何穿过孔道、如何在活性位点附近形成局部浓度差；又比全原子MD更轻量，能够把体系做到更大，把时间拉得更长，并服务于工程参数的筛选。
因此，HS-PPM不是为了替代所有分子模拟方法，而是补上一个重要位置：当问题的关键在于气体输运、孔道几何、碰撞统计和扩散-反应/吸附耦合时，它能以较低成本给出具有物理解释力的结果。对于多级孔催化剂、吸附剂优化、复杂气体混合物传质和反应器多尺度模拟，HS-PPM都可以作为连接微观机制与宏观性能的重要工具。
结语
HS-PPM把真实分子世界适度"粗粒化"，再用高效并行算法把它重新组织起来。它让我们不必在"看得清"和"算得动"之间二选一，而是用硬球碰撞、拟颗粒耦合和区域并行，搭建起理解气体扩散与多孔材料传质的新窗口。
参考资料
(1) Li, Y. P. et al. Simulation of the effect of coke deposition on the diffusion of methane in zeolite ZSM-5. Chem. Eng. J. 2017, 320, 458-467.
(2) Zhang, C. L. et al. Hard-sphere/pseudo-particle modelling (HS-PPM) for efficient and scalable molecular simulation of dilute gaseous flow and transport. Mol. Simul. 2016, 42(14), 1171-1182.
(3) Qiu, T. H. et al. HS-PPM simulation for diffusion coefficients of binary and multi-component gas mixtures. AIChE J. 2025, 71(9).
(4) Li, Y. P. et al. Simulation Study on the Reaction-Diffusion Coupling in Simple Pore Structures. Langmuir 2017, 33(42), 11804-11816.
(5) Zhu, Y. F. et al. Simulation of Gas Diffusion-Adsorption Coupling in Adsorbents for Optimizing Pore Structures. Ind. Eng. Chem. Res. 2024, 63(40), 17402-17412.
