Zhang Rui, Schweizer Kenneth S
Department of Materials Science and Frederick Seitz Materials Research Laboratory, University of Illinois, 1304 West Green Street, Urbana, Illinois 61801, USA.
J Chem Phys. 2015 Oct 14;143(14):144906. doi: 10.1063/1.4932679.
We heuristically formulate a microscopic, force level, self-consistent nonlinear Langevin equation theory for activated barrier hopping and non-hydrodynamic diffusion of a hard sphere penetrant in very dense hard sphere fluid matrices. Penetrant dynamics is controlled by a rich competition between force relaxation due to penetrant self-motion and collective matrix structural (alpha) relaxation. In the absence of penetrant-matrix attraction, three activated dynamical regimes are predicted as a function of penetrant-matrix size ratio which are physically distinguished by penetrant jump distance and the nature of matrix motion required to facilitate its hopping. The penetrant diffusion constant decreases the fastest with size ratio for relatively small penetrants where the matrix effectively acts as a vibrating amorphous solid. Increasing penetrant-matrix attraction strength reduces penetrant diffusivity due to physical bonding. For size ratios approaching unity, a distinct dynamical regime emerges associated with strong slaving of penetrant hopping to matrix structural relaxation. A crossover regime at intermediate penetrant-matrix size ratio connects the two limiting behaviors for hard penetrants, but essentially disappears if there are strong attractions with the matrix. Activated penetrant diffusivity decreases strongly with matrix volume fraction in a manner that intensifies as the size ratio increases. We propose and implement a quasi-universal approach for activated diffusion of a rigid atomic/molecular penetrant in a supercooled liquid based on a mapping between the hard sphere system and thermal liquids. Calculations for specific systems agree reasonably well with experiments over a wide range of temperature, covering more than 10 orders of magnitude of variation of the penetrant diffusion constant.
我们试探性地建立了一种微观的、力水平的、自洽的非线性朗之万方程理论,用于描述硬球渗透剂在非常致密的硬球流体基质中的活化势垒跳跃和非流体动力学扩散。渗透剂动力学受渗透剂自身运动引起的力弛豫与集体基质结构(α)弛豫之间丰富竞争的控制。在没有渗透剂 - 基质吸引力的情况下,根据渗透剂 - 基质尺寸比预测了三种活化动力学 regime,它们在物理上的区别在于渗透剂跳跃距离以及促进其跳跃所需的基质运动性质。对于相对较小的渗透剂,当基质有效地充当振动无定形固体时,渗透剂扩散常数随尺寸比下降最快。由于物理键合,增加渗透剂 - 基质吸引力强度会降低渗透剂扩散率。对于尺寸比接近 1 的情况,出现了一种独特的动力学 regime,与渗透剂跳跃对基质结构弛豫的强烈从属相关。在中间渗透剂 - 基质尺寸比处的交叉 regime 将硬渗透剂的两种极限行为连接起来,但如果与基质有强吸引力,基本上会消失。活化渗透剂扩散率随基质体积分数强烈下降,且这种下降随着尺寸比增加而加剧。我们基于硬球系统与热液体之间的映射,提出并实施了一种用于刚性原子/分子渗透剂在过冷液体中活化扩散的准通用方法。特定系统的计算在很宽的温度范围内与实验结果相当吻合,涵盖了渗透剂扩散常数变化超过 10 个数量级的范围。