[論文レビュー] Brownian motion of free particles on curved surfaces

Ramón Castañeda-Priego, Pavel Castro-Villarreal+2 more|arXiv (Cornell University)|Nov 25, 2012

Particle Dynamics in Fluid Flows参考文献 2被引用数 5

ひとこと要約

本稿は、リーマンノーマル座標におけるラングジュアン方程式を用いて、曲がった表面における自由粒子のブラウン運動を調査し、拡散運動が多様体上の拡散方程式と整合することを示している。主な結果として、短時間では接平面での拡散が発生し、長時間では平均二乗変位が表面の幾何学的性質にのみ依存して飽和することを示しており、数値シミュレーションおよび円と球面上での解析的解により検証されている。

ABSTRACT

Brownian motion of free particles on curved surfaces is studied by means of the Langevin equation written in Riemann normal coordinates. In the diffusive regime we find the same physical behavior as the one described by the diffusion equation on curved manifolds [J. Stat. Mech. (2010) P08006]. Therefore, we use the latter in order to analytically investigate the whole diffusive dynamics in compact geometries, namely, the circle and the sphere. Our findings are corroborated by means of Brownian dynamics computer simulations based on a heuristic adaptation of the Ermak-McCammon algorithm to the Langevin equation along the curves, as well as on the standard algorithm, but for particles subjected to an external harmonic potential, deep and narrow, that possesses a "Mexican hat" shape, whose minima define the desired surface. The short-time diffusive dynamics is found to occur on the tangential plane. Besides, at long times and compact geometries, the mean-square displacement moves towards a saturation value given only by the geometrical properties of the surface.

研究の動機と目的

• 曲がった2次元多様体上における自由ブラウン運動粒子の運動を理解すること、特に細胞膜におけるリン脂質やタンパク質の拡散といったバイオフィジカルな文脈を想定する。
• 曲がった空間におけるラングジュアン方程式が、平坦空間と同様に平均二乗変位と時間の間のアインシュタイン関係を保存するかどうかを特定すること。
• 特に円や球のようなコンパクトな多様体の表面幾何学的性質が、長時間における拡散的挙動に与える影響を分析すること。
• 適応されたアルゴリズムを用いたブラウン運動シミュレーションと、拡散方程式の解析的解を併用して理論的予測を検証すること。
• 曲率と位相的制約が、長時間における平均二乗変位の飽和値をどのように決定するかを探索すること。

提案手法

• 曲がった多様体上での過減衰ブラウン運動を記述するため、リーマンノーマル座標におけるラングジュアン方程式を定式化する。
• 特にコンパクトな表面における完全な拡散的挙動を調査するため、多様体上の拡散方程式を解析的手段として用いる。
• ブラウン運動のシミュレーションを曲がった表面に沿って行うために、エルマク＝マケンモン法のヒューリスティックな適応を適用する。
• シミュレーションにおいて粒子を所望の曲がった多様体上に閉じ込めるために、『メキシカン・ハット』型の調和ポテンシャルを用いる。
• 球面上での平均二乗変位の解析的表現を、球面調和関数およびルジャンドル多項式を用いて導出する。
• 統計力学の原則に基づいて導かれる擬似補正力（pseudo-corrective force）をシミュレーションに組み込み、曲率制約を考慮する。

実験結果

リサーチクエスチョン

• RQ1 曲がった空間におけるラングジュアン方程式は、平坦空間と同様に平均二乗変位と時間の間の標準的なアインシュタイン関係を再現するか？
• RQ2 円や球のようなコンパクト表面の幾何学的性質は、ブラウン拡散の長時間挙動にどのように影響するか？
• RQ3 曲がった表面における短時間の拡散的挙動の性質は何か？ これは接平面で発生するのか？
• RQ4 多様体上の拡散方程式は、曲がった表面における自由粒子の完全な拡散的挙動を正確に記述できるか？
• RQ5 コンパクト多様体上での長時間における平均二乗変位の飽和値は、何に依存するか？

主な発見

• 短時間における拡散的挙動は接平面で発生し、これはグローバルな曲率に関係なく局所的にユークリッド的挙動が現れることを示している。
• 長時間において、平均二乗変位は表面の幾何的性質（特に全表面積や曲率）にのみ依存して飽和値に達する。
• 球面上では、長時間における平均二乗変位の飽和が解析的に導出され、シミュレーションでも確認されており、拡散方程式の解と整合している。
• 円および球面において、拡散方程式の解析的解は、有限な値に飽和する平均二乗変位を導き、これはコンパクトトポロジーと整合している。
• リーマンノーマル座標におけるラングジュアン方程式は、曲がった多様体上での拡散方程式が記述する拡散的挙動を正しく再現している。
• 適応されたエルマク＝マケンモン法と調和ポテンシャルによる閉じ込め法の両方を用いたシミュレーションが、解析的予測を検証し、長時間拡散の幾何学的依存性を確認している。