[论文解读] Domain wall of the totally asymmetric exclusion process without particle number conservation

Yunxin Zhang|arXiv (Cornell University)|Oct 7, 2008

Stochastic processes and statistical mechanics参考文献 26被引用 5

一句话总结

本文从理论上研究了具有粒子产生与湮灭的完全不对称排除过程（TASEP），利用平均场近似和Rankine-Hugoniot条件，建立了域壁存在的必要与充分条件。推导出域壁位置与高度的精确解析表达式，揭示其对边界注入率 $\alpha$、出口率 $\beta$ 以及粒子产生/湮灭率 $\Omega_a$、$\Omega_d$ 的依赖关系，结果经模拟验证。

ABSTRACT

In this research, the totally asymmetric exclusion process without particle number conservation is discussed. Based on the mean field approximation and the Rankine-Hugoniot condition, the necessary and sufficient conditions of the existence of the domain wall have been obtained. Moreover, the properties of the domain wall, including the location and height, have been studied theoretically. All the theoretical results are demonstrated by the computer simulations.

研究动机与目标

• 从理论上确定在粒子数不守恒的TASEP中域壁存在的必要与充分条件。
• 解析推导域壁位置与高度作为系统参数 $\alpha$、$\beta$、$K=\Omega_a/\Omega_d$ 和 $\Omega_d$ 的函数。
• 解决此前仅通过模拟研究的非守恒TASEP系统中域壁性质的开放问题。

提出的方法

• 应用平均场近似，推导出热力学极限下稳态密度分布的连续微分方程。
• 利用Rankine-Hugoniot条件识别密度分布中不连续性（即域壁）的位置。
• 在左右边界条件下求解平均场方程，得到密度分布 $\rho_{l\alpha}(x)$ 和 $\rho_{r\beta}(x)$。
• 推导逆函数 $\rho_{l\alpha}^{-1}(\gamma)$ 和 $\rho_{r\beta}^{-1}(1-\alpha)$ 的隐式方程，以确定域壁存在的条件。
• 解析计算域壁高度 $2\epsilon$ 和位置 $x_S$ 对 $\alpha$、$\beta$、$K$ 和 $\Omega_d$ 的偏导数，以判断其单调性。
• 通过数值模拟验证理论预测，确认解析结果的正确性。

实验结果

研究问题

• RQ1 在具有Langmuir动力学的非守恒TASEP中，域壁在何种条件下存在？
• RQ2 域壁位置 $x_S$ 如何依赖于注入率 $\alpha$、出口率 $\beta$ 以及比值 $K=\Omega_a/\Omega_d$？
• RQ3 域壁高度 $2\epsilon$ 如何随 $\alpha$、$\beta$、$K$ 和 $\Omega_d$ 变化？
• RQ4 域壁高度是否关于 $\alpha$、$\beta$、$K$ 和 $\Omega_d$ 单调？如果不是，其极值出现在何处？
• RQ5 域壁分布的精确解析形式及其在稳态下的稳定性如何？

主要发现

• 当 $0 \leq \alpha \leq 0.5$ 且 $1-\beta \leq K/(K+1)$ 时，域壁存在的充要条件为：$\rho_{l\alpha}^{-1}(\gamma) \leq 1$ 且 $\rho_{r\gamma}^{-1}(1-\alpha) \leq 0$，其中 $\gamma = \min(0.5, \beta)$。
• 当 $0.5 \leq \alpha \leq 1$ 时，无论其他参数如何，均不存在域壁。
• 域壁位置 $x_S$ 随 $\beta$ 单调增加，随 $\alpha$ 单调减小，表明随着注入率增加，域壁向右端移动。
• 域壁高度 $2\epsilon$ 在所有 $0 \leq \alpha, \beta \leq 0.5$ 范围内随 $\beta$ 单调减小，证实出口率对壁振幅具有抑制作用。
• 当 $\frac{1}{K+1} \leq \beta \leq 0.5$ 且 $K \geq 3$ 时，高度 $2\epsilon$ 随 $K$ 和 $\Omega_d$ 增加，而位置 $x_S$ 随两者减小。
• 当 $0 \leq \beta \leq \frac{1}{K+1}$ 时，高度 $2\epsilon$ 随 $\alpha$ 和 $\Omega_d$ 减小，位置 $x_S$ 对 $K$ 展现出非单调行为。