[论文解读] The blockage problem

Ovidiu Costin, Joel L. Lebowitz+2 more|arXiv (Cornell University)|Jul 27, 2012

Stochastic processes and statistical mechanics参考文献 15被引用 8

一句话总结

本文研究了在整数集 Z 上的 totally asymmetric exclusion process (TASEP) 中，当位置 0 处的跃迁速率降低为 r 时的阻塞问题。通过电流 j(r) 的 16 阶幂级数展开，作者提供了强有力的启发式证据，表明临界值 rc = 1，即最大电流仅在 r ≥ 1 时达到 1/4。他们进一步表明，当 r ≲ 1 时，j(r) ≈ 1/4 − γexp[−a/(1−r)]，其中 a ≈ 2，并证明了 j(r) = J_max(r)，从而确立了电流函数的连续性。

ABSTRACT

We investigate the totally asymmetric exclusion process on Z, with the jump rate at site i given by r_i=1 for i nonzero, r_0=r. It is easy to see that the maximal stationary current j(r) is nondecreasing in r and that j(r)=1/4 for r>=1; it is a long outstanding problem to determine whether or not the critical value r_c of r such that j(r)=1/4 for r>r_c is strictly less than 1. Here we present a heuristic argument, based on the analysis of the first sixteen terms in a formal power series expansion of j(r) obtained from finite volume systems, that r_c=1 and that for r less than 1 and near 1, j(r) behaves as 1/4-γ\exp[-{a/(1-r)}] with a approximately equal to 2. We also give some new exact results about this system; in particular we prove that j(r)=J_max(r), with J_max(r) the hydrodynamic maximal current defined by Seppalainen, and thus establish continuity of j(r). Finally we describe a related exactly solvable model, a semi-infinite system in which the site i=0 is always occupied. For that system, the critical r is 1/2 and the analogue j_s(r) of j(r) satisfies j_s(r)=r(1-r) for r<=1/2; j_s(r) is the limit of finite volume currents inside the curve |r(1-r)|=1/4 in the complex r plane and we suggest that analogous behavior may hold for the original system.

研究动机与目标

• 为长期悬而未决的问题提供解答：Z 上 TASEP 的临界阻塞速率 rc 是否严格小于 1。
• 理解跃迁速率的局部变化对非平衡系统中稳态电流的全局影响。
• 确定 r 接近 1 时电流 j(r) 的函数形式，并判断其是否表现出本性奇点。
• 建立最大电流 j(r) 作为 r 的函数的连续性。

提出的方法

• 通过有限体积系统推导出电流 j(r) 的 16 阶幂级数展开。
• 利用 Padé 逼近和二阶差分最小化方法分析系数，以估计奇点的位置。
• 应用分式线性变换将 r=1 处的奇点映射至无穷远，改善收敛性，并估计极点位于 r₀ ≈ -1.5437。
• 使用鞍点法推导 exp[a/(1−r)] 的泰勒系数在大 k 时的渐近行为。
• 证明 j(r) = J_max(r)，即 Seppäläinen 定义的流体动力学最大电流，从而确保 j(r) 的连续性。
• 研究了一个相关且可精确求解的半无限系统，其中位置 0 始终被占据，得到 j^{s-i}(r) = r(1−r) 对于 r ≤ 1/2。

实验结果

研究问题

• RQ1 临界阻塞速率 rc 是否严格小于 1，还是如平均场理论所暗示的 rc = 1？
• RQ2 当 r 接近 1 时，电流 j(r) 的函数形式是什么，特别是其是否以本性奇点方式趋近于 1/4？
• RQ3 电流 j(r) 在 r=1 处是否保持连续，且能否严格证明这一点？
• RQ4 j(r) 的幂级数展开是否可用于推断复平面上奇点的位置与性质？
• RQ5 相关且可精确求解的半无限模型的行为是否暗示了无限系统具有类似行为？

主要发现

• 基于 j(r) 前 16 项幂级数的分析，临界值 rc 被确定为 1，而非小于 1。
• 当 r ≲ 1 时，电流表现为 j(r) ≈ 1/4 − γexp[−a/(1−r)]，其中 a ≈ 2，表明在 r=1 处存在本性奇点。
• 证明了最大电流 j(r) 等于 J_max(r)，即流体动力学最大电流，从而确保 j(r) 作为 r 的函数是连续的。
• 通过两种独立的数值方法，高精度确定了 j(r) 的幂级数在 r₀ ≈ -1.5437 处存在极点。
• 在相关半无限模型中，位置 0 始终被占据，得到 j^{s-i}(r) = r(1−r) 对于 r ≤ 1/2，且 j^{s-i}(r) 是复平面上 |r(1−r)| = 1/4 曲线内部有限体积电流的极限。
• 发现 exp[a/(1−r)] 的泰勒系数 b_k 的渐近行为为 b_k ∼ const·k^{-3/4}e^{2√(ak)}，当 k→∞ 时，与 j(r) 的指数形式一致。