2Dランダムウォーク 戻る
確率・拡散

2Dランダムウォーク

複数ウォーカーの軌跡と平均二乗変位(MSD)をリアルタイム可視化。格子・ガウス・レヴィフライトを切り替えて拡散則を比較。

パラメータ

ウォーカー数10
ステップサイズ5
Xドリフト0.0
歩行タイプ
境界条件
速度5
0
ステップ
0.0
現在MSD
0.0
理論MSD=2Dt
0.0
最大変位
$$\langle r^2(t)\rangle = 2Dt$$
通常拡散のMSDは時間に比例。ログ-ログプロットの傾きで拡散型を識別。

2Dランダムウォークとは

🧑‍🎓
ランダムウォークって、ただ酔っぱらいがフラフラ歩くだけのモデルですよね?これが何の役に立つんですか?
🎓
ざっくり言うと、不規則な動きを理解するための「ものさし」なんだ。例えば、細胞内で栄養分子がどう拡散するか、株価の変動、昆虫の餌探しの動きまで、幅広く応用されているよ。このシミュレーターで「歩行タイプ」を「格子」から「連続」に変えてみると、動きの滑らかさが全然違うのがわかる。
🧑‍🎓
なるほど。でも、グラフにある「平均二乗変位(MSD)」って何を見てるんですか?
🎓
MSDは、粒子が平均的にどれだけ遠くまで広がったかを数値化したものだ。実務では、実験で計測した粒子の動きから拡散の速さ(拡散係数)を求めるのに使うことが多いね。シミュレーター右側のログ-ログプロットを見て。通常の拡散(格子・連続)なら、MSDの直線の傾きはだいたい1になる。でも「レヴィフライト」に切り替えて「べき指数α」を小さく(例えば1.5)すると、傾きが1より大きくなるのが観測できるよ。これが「超拡散」だ。
🧑‍🎓
超拡散って面白い!でも、そんなに大きくジャンプする動きって、現実にありますか?
🎓
あるんだ。例えば、海のサメやアルバトロスが餌を探す動きは、小刻みに動いたあと、まれに長距離を移動する。これがまさにレヴィフライトで近似できる。シミュレーターで「ウォーカー数」を増やして「レヴィフライト」で実行してみて。ほとんどの粒子は狭い範囲にいるけど、たまに遠くまで飛び出す粒子がいるのが視覚的にわかるはずだ。

物理モデルと主要な数式

ランダムウォークの核心は、粒子の位置の2乗平均(平均二乗変位:MSD)が時間とともにどう増加するかです。通常の拡散(ブラウン運動)では、以下の関係が成り立ちます。

$$\langle r^2(t) \rangle = 2d D t$$

ここで、$\langle r^2(t) \rangle$は時刻$t$における平均二乗変位、$d$は次元数(このシミュレーターでは2)、$D$は拡散係数です。この関係から、MSDの時間に対する傾きが拡散の「型」を示すことがわかります。

レヴィフライトでは、ステップ長$l$の確率分布がべき乗則 $P(l) \sim l^{-\alpha}$ に従います。これにより、MSDは時間のべき乗に比例する超拡散を示します。

$$\langle r^2(t) \rangle \sim t^{\beta}, \quad \beta > 1$$

ここで、$\alpha$はべき指数(シミュレーターのパラメータ)、$\beta$はMSDの時間依存性の指数です。$\alpha \leq 3$のときに超拡散($\beta > 1$)が現れ、$\alpha$が小さいほど大きなジャンプが起きやすくなります。

実世界での応用

生物物理学・生態学:細胞内でのタンパク質やmRNAの拡散、微生物や動物の採餌行動の分析に使われます。レヴィフライトは効率的な探索戦略として注目されています。

材料科学:多孔質材料内でのガスや液体の浸透・拡散過程のモデル化。不規則な構造の中での粒子の動きを予測するのにランダムウォークが応用されます。

金融工学:株価や為替レートの変動をモデル化する基礎として使われます。特に、大きな価格変動(ジャンプ)を含むモデルではレヴィフライトが参考にされます。

高分子科学:ポリマー鎖の構造や動的特性を、鎖を構成するモノマーのランダムウォークの集合体として統計的に理解するアプローチがあります。

よくある誤解と注意点

まず、「拡散係数Dは粒子の種類だけで決まる定数」と思いがちですが、それは誤解です。 例えば、同じ分子でも、水の中と油の中では拡散の速さは全く異なります。このシミュレーターで「歩行タイプ」を「連続」にし、「ステップサイズ」を固定したまま「拡散係数」のスライダーを動かしてみてください。粒子の動きの「活発さ」が変わりますよね?これは、拡散係数が「環境(粘性など)」と「温度」に強く依存することを意味しています。実務では、シミュレーションで使うDの値を実験値や文献値から正しく設定することが第一歩です。

次に、MSDのグラフを1回の実行で判断しないこと。 特にウォーカー数が少ない場合、MSDのプロットはかなりばらつきます。「格子」歩行でウォーカー数1で実行すると、MSDが綺麗な直線にならないことが確認できるはずです。これは統計的なゆらぎです。信頼性のある結果を得るには、十分な数のウォーカー(例えば100以上)でシミュレーションを走らせ、平均を取る必要があります。

最後に、レヴィフライトのパラメータ設定には注意が必要です。 べき指数αを2.0に近い値(例えば2.1)に設定すると、見た目は通常のランダムウォークとほとんど区別がつきません。超拡散の特徴を明確に見たいなら、αを1.5から2.0の間で設定しましょう。ただし、αを小さくしすぎる(例えば1.1)と、極端に長いジャンプが発生し、シミュレーション領域からすぐに飛び出してしまうことも。実現象への適用時は、観測データからαを慎重に推定するプロセスが不可欠です。

関連する工学分野

このツールで扱う計算は、「輸送現象」を扱うほぼ全ての工学分野の基礎として登場します。例えば、化学工学では、反応器内での試薬の混合や、カラムクロマトグラフィーでの分離プロセスをモデル化する際、物質の拡散が核心となります。多孔質触媒内部での反応物質の拡散をランダムウォークでシミュレートし、反効率を予測するのです。

半導体工学では、イオン注入や拡散プロセスで不純物原子がシリコン中をどう移動するかを予測するのに、ランダムウォークの考え方が使われます。ここでは「格子」ランダムウォークが結晶格子の中の原子のホッピングを抽象化したモデルとして直接対応します。

また、土木・環境工学における地下水や土壌中の汚染物質の拡散予測も重要な応用先です。地層の不均質性を反映させるため、ステップサイズや方向の確率分布を調整したランダムウォーク(連続タイプ)が用いられます。さらに、バイオエンジニアリングでは、創薬において薬剤が生体組織内をどう拡散して標的に到達するか、その効率を評価する際の基礎モデルとなっています。

発展的な学習のために

まず次の一歩は、「フォッカー-プランク方程式」を調べてみることです。 これは、ランダムウォークのような確率過程を、粒子の確率密度分布の時間発展を記述する「決定論的な」偏微分方程式で表現したものです。拡散方程式 $$\frac{\partial P(\mathbf{r}, t)}{\partial t} = D \nabla^2 P(\mathbf{r}, t)$$ はその最も単純な形です。シミュレーターで多数の粒子のランダムな歩行を追跡する(モンテカルロ法)ことと、この方程式を解くことは、表裏一体の関係にあると理解できるでしょう。

数学的背景を深めたいなら、「確率過程論」や「ストキャスティック微分方程式(SDE)」が鍵になります。通常のランダムウォークは「ウィーナー過程」へ、レヴィフライトはより一般的な「レヴィ過程」へと発展します。ここでは、シミュレーターの「ステップ」が数学的には「確率変数の実現値」であることを厳密に学べます。

実践的な次のトピックとしては、「障害物がある場合のランダムウォーク」を推奨します。このシミュレーターに、動けない障害物を配置したら粒子の動きはどう変わるか? MSDの傾きは小さくなる(拡散が遅くなる)でしょうか? これは、細胞内のクロマチン繊維の中でのタンパク質探索や、複合材料中のイオン伝導など、現実の複雑な系をモデル化する上で必須の考え方です。まずは紙とペンで、単純な格子に障害物をいくつか置いた場合の経路を考えてみることから始めてみましょう。