[論文レビュー] Mean first-passage time to a small absorbing target in three-dimensional elongated domains

Denis S. Grebenkov, Alexei T. Skvortsov|arXiv (Cornell University)|Feb 15, 2022

Diffusion and Search Dynamics参考文献 68被引用数 12

ひとこと要約

本稿では、断面積がゆっくりと変化する3次元軸対称で細長い領域内の任意の形状をした小さな吸収的ターゲットへの平均第一通過時間（MFPT）の近似解析的公式を導出する。3次元Fokker-Planck問題を、ターゲット位置に半透膜を設けた1次元有効拡散方程式に還元することで、有効反応性を静電容量から導出した手法により、ターゲットの位置、径方向分布、形状に依存するMFPTの依存関係を捉える。モンテカルロシミュレーションによる検証が行われている。

ABSTRACT

We derive an approximate formula for the mean first-passage time (MFPT) to a small absorbing target of arbitrary shape inside an elongated domain of a slowly varying axisymmetric profile. For this purpose, the original Poisson equation in three dimensions is reduced an effective one-dimensional problem on an interval with a semi-permeable semi-absorbing membrane. The approximate formula captures correctly the dependence of the MFPT on the distance to the target, the radial profile of the domain, and the size and the shape of the target. This approximation is validated by Monte Carlo simulations.

研究の動機と目的

• 3次元細長い領域に存在する小さな吸収的ターゲットへの平均第一通過時間（MFPT）の一般化された近似を確立すること。
• ターゲットが横方向に小さく、縦方向に無視できない場合に、狭域脱出理論の限界を克服すること。
• 静電容量から導出される有効捕獲係数を用いて、ターゲットの形状と径方向位置がMFPTに与える影響を組み込むこと。
• 多様な領域形状（円柱型、円錐型、振動的プロファイル）におけるモンテカルロシミュレーションを用いて、手法の妥当性を検証すること。
• 2次元平面的手法を3次元に拡張し、細長いチューブにおける径方向の閉じ込めとエントロピー的ドリフトを考慮すること。

提案手法

• 任意の形状の小さな吸収的ターゲットを、その捕獲係数Kと同一の等価な吸収的ディスクに置き換える。Kは静電容量CからK = 4πDCの式で導出される。
• ターゲットの軸方向位置zTに、有効反応性κ = K(rT)/S(zT)を有する半透膜（半吸収膜）を導入する。ここでS(z)は断面積を表す。
• MFPTを1次元問題としてモデル化し、後向きFokker-Planck方程式：d/dz [S(z) dT/dz] = -S(z)/D を用いる。境界条件はz=0およびz=ℓでノイマン境界条件。
• ターゲットサイズν = a/r(zT)と径方向オフセットη = rT/r(zT)の関数として、軸方向および壁隣接極限の間を補間する関数Ψ(ν,η)を用いて、位置依存捕獲係数K(rT)の解析的近似を導出する。
• zTにおける整合条件を満たすように、区間(0,zT)と(zT,ℓ)で1次元MFPT方程式を別々に解き、全MFPTプロファイルを求める。
• 多様なプロファイル（円柱型、円錐型、振動的プロファイル）を有する領域において、広範なモンテカルロシミュレーションと比較して、解析結果の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 狭域脱出理論の限界を超えて、3次元細長い領域内の小さなターゲットへのMFPTが、径方向位置と形状にどのように依存するか？
• RQ2 任意のターゲット幾何形状を有する、非等方的で断面積がゆっくりと変化する3次元領域において、有効1次元膜モデルがMFPTを正確に捉えられるか？
• RQ3 任意の形状の小さなターゲットを置き換える膜の正しい有効反応性κは何か？また、そのターゲットの静電容量および位置とどう関係するか？
• RQ4 径方向の閉じこめとエントロピー的ドリフトが、細長い領域におけるMFPTに与える影響は何か？また、これを低次元化された1次元モデルで捉えることができるか？
• RQ5 非球形ターゲットや複雑な領域プロファイルに対しても、提案された近似はどの程度有効性を保っているか？

主な発見

• MFPTは、領域の曲率の影響を反映する位置依存捕獲係数K(rT)に強く依存するため、ターゲットの径方向位置rTに依存する。
• 有効反応性κ = K(rT)/S(zT)は、1次元モデルにおけるターゲットの吸収効率を的確に捉え、K(rT)は軸方向および壁隣接極限を補間するキャリブレーション関数Ψ(ν,η)から導出される。
• Fick-Jacobs型方程式 d/dz [S(z) dT/dz] = -S(z)/D を用いた1次元還元は、断面積がゆっくりと変化する領域でもMFPTを正確に再現する。
• モンテカルロシミュレーションとの比較により、円柱型、円錐型、振動的プロファイルを有する多様な領域形状において、本手法が高い精度を示すことが確認された。
• ターゲットが横方向に小さく、縦方向に無視できない場合でも、古典的狭域脱出理論の限界を克服でき、近似が有効であることが示された。
• 本手法により、チューブ状または糸状構造を有する生物学的・物理的系における拡散的探索のMFPTを高速かつ解析的に推定できるようになった。