[論文レビュー] Open reaction-diffusion systems: bridging probabilistic theory and simulations across scales

Mauricio J. del Razo, Kostré, Margarita|arXiv (Cornell University)|Apr 10, 2024

Gene Regulatory Network Analysis被引用数 5

ひとこと要約

この論文は粒子レベルの開放反応-拡散系を記述する化学拡散マスター方程式（CDME）を開発し、次にこの確率的枠組みを巨視的な濃度ベースのモデルと貯蔵体（リザーバー）との相互作用へ結びつけ、マルチスケールシミュレーションを実現する方法を示す。

ABSTRACT

Reaction-diffusion processes are the foundational model for a diverse range of complex systems, ranging from biochemical reactions to social agent-based phenomena. The underlying dynamics of these systems occur at the individual particle/agent level, and in realistic applications, they often display interaction with their environment through energy or material exchange with a reservoir. This requires intricate mathematical considerations, especially in the case of material exchange since the varying number of particles/agents results in ``on-the-fly'' modification of the system dimension. In this work, we first overview the probabilistic description of reaction-diffusion processes at the particle level, which readily handles varying number of particles. We then extend this model to consistently incorporate interactions with macroscopic material reservoirs. Based on the resulting expressions, we bridge the probabilistic description with macroscopic concentration-based descriptions for linear and nonlinear reaction-diffusion systems, as well as for an archetypal open reaction-diffusion system. Using these mathematical bridges across scales, we finally develop numerical schemes for open reaction-diffusion systems, which we implement in two illustrative examples. This work establishes a methodological workflow to bridge particle-based probabilistic descriptions with macroscopic concentration-based descriptions of reaction-diffusion in open settings, laying the foundations for a multiscale theoretical framework upon which to construct theory and simulation schemes that are consistent across scales.

研究の動機と目的

• 生物学や社会現象に存在する開放的な反応拡散系の動機づけと、一貫した開放系モデリングの必要性。
• 反応による系の次元変化を扱う確率論的な粒子レベルの記述を定式化する。
• 粒子レベルのダイナミクスを巨視的貯蔵体と結びつけ、巨視的方程式を導出する枠組みを開発する。
• 線形および非線形反応に対してCDMEから濃度ベースのRD方程式を回収するためのワークフローを確立する。

提案手法

• CDMEを、位置 x^(n) にある n 個の粒子の確率を表す密度 rho_n(x^(n))の階層として提示する。
• 拡散は標準的なブラウン運動項と、喪失および獲得成分を含む反応演算子を用いて記述する。
• 単純な分解/創生のCDMEの明示形と、二分子反応（多種反応を含む一般化）を提供する。
• 離散的または連続的なレートを用いた境界層境界条件を介して貯蔵体交換をモデル化し、開放境界を取り入れる。
• 線形・非線形の場合について、CDMEの平均場極限として巨視的RD-PDEが出現することを示す。
• 拡散影響下・開放反応における粒子レベルのパラメータを限界手順で巨視的速度定数へ結びつける。

Figure 1: Structure of the phase space of the CDME. a. Phase space for a system with one chemical species $A$ . The structure of the phase space consists of discrete sets of continuous diffusion domains $\mathbb{X}$ depending on the number of particles. The transitions between sets depend on the spe

実験結果

リサーチクエスチョン

• RQ1 粒子数が変化する反応拡散ダイナミクスを粒子レベルで確率論的にどう定式化できるか？
• RQ2 巨視的貯蔵体との相互作用を粒子レベルのCDMEに一貫して組み込むにはどうするか？
• RQ3 線形・非線形反応について、CDMEから濃度ベースの開放RD方程式をどう回収するか？
• RQ4 拡散影響下/開放系における微視的パラメータと巨視的速度定数との関係は何か？

主な発見

• CDMEは粒子数が変化する拡散と反応を結ぶ無限次元相空間での線形発展枠組みを提供する。
• 開放境界はリザーバー連結境界層としてモデル化され、ディラックデルタ表現と整合する出入り速度の明示形を得る。
• CDMEの平均場極限は、分解・創生の古典的RD-PDEを回収し、開放系における拡散影響下の反応項を生じさせる。
• 非線形の二分子様反応の場合、巨大コピー数仮定の下で微視的から巨視的への対応付けを確立し、反応速度定数を巨視的動力学へ結びつける。

Figure 2: Diagram representing loss from and gain into the $n$ particle state due to the reaction $kA\rightarrow lA$ assuming $k>l$ . The loss $\mathcal{L}_{n}\rho_{n}$ must depend on the reactant’s positions, so it is a function of $x^{(n)}$ , while the gain $\mathcal{G}_{n}\rho_{n+k-l}$ depends on