[論文レビュー] A note on phase transitions for the Smoluchowski equation with dipolar potential

Pierre Degond, Amic Frouvelle+1 more|arXiv (Cornell University)|Dec 17, 2012

Nonlinear Dynamics and Pattern Formation参考文献 14被引用数 7

ひとこと要約

本稿では、球面上の修正されたスモラッチョフスキー方程式における相転移を、双極子ポテンシャルと、秩序パラメータに依存する整列および拡散強度を用いて研究する。第二種および第一種相転移の条件を同定し、解析的および数値的手法を用いてヒステリシスを示し、安定平衡状態への指数的収束率を確立する。

ABSTRACT

In this note, we study the phase transitions arising in a modified Smoluchowski equation on the sphere with dipolar potential. This equation models the competition between alignment and diffusion, and the modification consists in taking the strength of alignment and the intensity of the diffusion as functions of the order parameter. We characterize the stable and unstable equilibrium states. For stable equilibria, we provide the exponential rate of convergence. We detail special cases, giving rise to second order and first order phase transitions, respectively. We study the hysteresis diagram, and provide numerical illustrations of this phenomena.

研究の動機と目的

• 秩序パラメータに依存する整列および拡散係数を有する修正されたスモラッチョフスキー方程式における相転移を分析すること。
• 整列および拡散関数に一般な条件を課した場合の安定および不安定な平衡状態を特徴づけること。
• 初期密度をゆっくり変化させた際に、系の長期的ダイナミクスにヒステリシスが現れる仕組みを調査すること。
• 安定平衡状態への明示的な指数的収束率を提示し、相転移の臨界閾値を同定すること。
• 理論的予測であるヒステリシスと第一種相転移を数値的に検証すること。

提案手法

• 秩序パラメータの大きさ $|J_f|$ に依存する整列および拡散係数を有する球面上の修正されたスモラッチョフスキー方程式を定式化する。
• 自由エネルギー汎関数と勾配流構造を用いて平衡状態の安定性を分析し、$\Phi(|J|)$ を比 $h(|J|) = \nu(|J|)/\tau(|J|)$ の原始関数として用いる。
• 初期密度 $\rho$ をゆっくり変化するパラメータとして扱うことで、時間スケール分離の議論を用い、ヒステリシスサイクルを模擬する。
• 空間方向に有限差分法、時間方向に陰解法を用いた数値スキームを採用し、不安定性閾値を超える収束を促進するための摂動戦略を導入する。
• 臨界値に近い場合の均一状態付近での数値的停滞を回避するため、閾値に基づく摂動を導入する。
• 数値的シミュレーションと理論的ヒステリシスループを比較し、第一種相転移の存在を検証する。

実験結果

リサーチクエスチョン

• RQ1 修正されたスモラッチョフスキー方程式が第一種相転移と第二種相転移を示す条件は何か？
• RQ2 秩序パラメータに依存する整列および拡散の依存性が、平衡状態の安定性に与える影響は何か？
• RQ3 初期密度 $\rho$ が系の長期的ダイナミクスにおけるヒステリシスの発生に果たす役割は何か？
• RQ4 初期密度 $\rho$ をゆっくり変化させた場合、系はヒステリシスループを示すことができるか？また、標準的なスモラッチョフスキー方程式と比較するとどうなるか？
• RQ5 安定平衡状態への指数的収束の明示的レートは何か？また、系のパラメータにどのように依存するか？

主な発見

• 比 $\nu(|J|)/\tau(|J|)$ が線形である場合、系は連続的（第二種）相転移を示し、古典的スモラッチョフスキー方程式と整合する。
• 非線形な比の場合は、臨界値 $\rho_c$ および $\rho_*$ において秩序パラメータに不連続なジャンプが生じ、ヒステリシスループを形成する第一種相転移を示す。
• 解は安定平衡状態へ指数的かつ高速に収束し、線形化作用素の固有値特性から明示的な収束レートが導出される。
• 数値的シミュレーションにより、$\rho$ をゆっくり変化させた際、濃度パラメータ $c = |J_{\widetilde{f}}|$ にヒステリシスループが存在することが確認され、理論的予測と一致する。
• 臨界値 $\rho_c$ は均一分布の不安定化の開始を、$\rho_*$ は非等方的平衡状態の再安定化をそれぞれ示す。
• 自由エネルギーが均一状態と非等方的状態で等しくなる値 $\rho_1$ は、均一状態の強い準安定性のため、ダイナミクスに顕著な影響を及ぼさない。