パラメータ
通常拡散のMSDは時間に比例。ログ-ログプロットの傾きで拡散型を識別。
複数ウォーカーの軌跡と平均二乗変位(MSD)をリアルタイム可視化。格子・ガウス・レヴィフライトを切り替えて拡散則を比較。
生物物理学・生態学:細胞内でのタンパク質やmRNAの拡散、微生物や動物の採餌行動の分析に使われます。レヴィフライトは効率的な探索戦略として注目されています。
材料科学:多孔質材料内でのガスや液体の浸透・拡散過程のモデル化。不規則な構造の中での粒子の動きを予測するのにランダムウォークが応用されます。
金融工学:株価や為替レートの変動をモデル化する基礎として使われます。特に、大きな価格変動(ジャンプ)を含むモデルではレヴィフライトが参考にされます。
高分子科学:ポリマー鎖の構造や動的特性を、鎖を構成するモノマーのランダムウォークの集合体として統計的に理解するアプローチがあります。
まず、「拡散係数Dは粒子の種類だけで決まる定数」と思いがちですが、それは誤解です。 例えば、同じ分子でも、水の中と油の中では拡散の速さは全く異なります。このシミュレーターで「歩行タイプ」を「連続」にし、「ステップサイズ」を固定したまま「拡散係数」のスライダーを動かしてみてください。粒子の動きの「活発さ」が変わりますよね?これは、拡散係数が「環境(粘性など)」と「温度」に強く依存することを意味しています。実務では、シミュレーションで使うDの値を実験値や文献値から正しく設定することが第一歩です。
次に、MSDのグラフを1回の実行で判断しないこと。 特にウォーカー数が少ない場合、MSDのプロットはかなりばらつきます。「格子」歩行でウォーカー数1で実行すると、MSDが綺麗な直線にならないことが確認できるはずです。これは統計的なゆらぎです。信頼性のある結果を得るには、十分な数のウォーカー(例えば100以上)でシミュレーションを走らせ、平均を取る必要があります。
最後に、レヴィフライトのパラメータ設定には注意が必要です。 べき指数αを2.0に近い値(例えば2.1)に設定すると、見た目は通常のランダムウォークとほとんど区別がつきません。超拡散の特徴を明確に見たいなら、αを1.5から2.0の間で設定しましょう。ただし、αを小さくしすぎる(例えば1.1)と、極端に長いジャンプが発生し、シミュレーション領域からすぐに飛び出してしまうことも。実現象への適用時は、観測データからαを慎重に推定するプロセスが不可欠です。
ウォーカー数N=100、ステップ数=500、格子タイプで各方向を等確率(0.25)選択する場合、t=500ステップ後のMSD期待値はおよそMSD≈500(正規ランダムウォーク:MSD∝t)。ドリフトなし、αが1.0(ガウス分布)の設定では拡散係数D≈0.5となり、δx²=2D×tの関係式を満たす。ドリフトX=0.05を加えると系全体が漂流し、平均位置がX方向に線形シフト。αを0.5に低下させるとMSD∝t^1.5の異常拡散が現れ、重い尾を持つレヴィフライトの支配領域に遷移