パラメータ
通常拡散のMSDは時間に比例。ログ-ログプロットの傾きで拡散型を識別。
複数ウォーカーの軌跡と平均二乗変位(MSD)をリアルタイム可視化。格子・ガウス・レヴィフライトを切り替えて拡散則を比較。
ランダムウォークの核心は、粒子の位置の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の傾きは小さくなる(拡散が遅くなる)でしょうか? これは、細胞内のクロマチン繊維の中でのタンパク質探索や、複合材料中のイオン伝導など、現実の複雑な系をモデル化する上で必須の考え方です。まずは紙とペンで、単純な格子に障害物をいくつか置いた場合の経路を考えてみることから始めてみましょう。