We develop a convergent reaction-drift-diffusion master equation (CRDDME) to facilitate the study of reaction processes in which spatial transport is influenced by drift due to one-body potential fields within general domain geometries. The generalized CRDDME is obtained through two steps. We first derive an unstructured grid jump process approximation for reversible diffusions, enabling the simulation of drift-diffusion processes where the drift arises due to a conservative field that biases particle motion. Leveraging the Edge-Averaged Finite Element method, our approach preserves detailed balance of drift-diffusion fluxes at equilibrium, and preserves an equilibrium Gibbs-Boltzmann distribution for particles undergoing drift-diffusion on the unstructured mesh. We next formulate a spatially-continuous volume reactivity particle-based reaction-drift-diffusion model for reversible reactions of the form $\textrm{A} + \textrm{B} \leftrightarrow \textrm{C}$. A finite volume discretization is used to generate jump process approximations to reaction terms in this model. The discretization is developed to ensure the combined reaction-drift-diffusion jump process approximation is consistent with detailed balance of reaction fluxes holding at equilibrium, along with supporting a discrete version of the continuous equilibrium state. The new CRDDME model represents a continuous-time discrete-space jump process approximation to the underlying volume reactivity model. We demonstrate the convergence and accuracy of the new CRDDME through a number of numerical examples, and illustrate its use on an idealized model for membrane protein receptor dynamics in T cell signaling.
翻译:我们提出了一种收敛的反应-漂移-扩散主方程(CRDDME),以促进研究在一般区域几何形状中空间输运受到单体质势场引起的漂移影响的反应过程。广义CRDDME通过两个步骤获得。我们首先推导了可逆扩散过程的非结构化网格跳跃过程近似,从而能够模拟漂移由偏置粒子运动的保守场引起的漂移-扩散过程。利用边平均有限元方法,我们的方法在平衡时保持了漂移-扩散通量的细致平衡,并为在非结构化网格上进行漂移-扩散的粒子保留了平衡吉布斯-玻尔兹曼分布。接下来,我们为形式为$\textrm{A} + \textrm{B} \leftrightarrow \textrm{C}$的可逆反应建立了一个基于粒子的空间连续体积反应性反应-漂移-扩散模型。采用有限体积离散化来生成该模型中反应项的跳跃过程近似。该离散化的构建确保了组合的反应-漂移-扩散跳跃过程近似与平衡时反应通量的细致平衡一致,并支持连续平衡态的离散版本。新的CRDDME模型代表了底层体积反应性模型的连续时间离散空间跳跃过程近似。我们通过多个数值算例展示了新CRDDME的收敛性和准确性,并说明了其在T细胞信号传导中膜蛋白受体动力学的理想化模型上的应用。