[論文レビュー] Diffusion on a lattice: transition rates, interactions and memory effects

M. A. Di Muro, M. Hoyuelos|arXiv (Cornell University)|Apr 3, 2022

Advanced Thermodynamics and Statistical Mechanics参考文献 35被引用数 1

ひとこと要約

本稿では、過剰化学ポテンシャルから導かれる遷移率を用いて記憶効果を考慮することで、格子拡散における相関係数 f の一般化された解析的表現を提案する。ソフトコアおよび拡張されたハードコア制約による粒子相互作用のモデル化を通じて、濃度に依存する f を導出し、平均場トレーサー拡散係数の予測を改善した。この方法は、粒子の排除サイズや密度が異なる状況において、モンテカルロシミュレーションと強い一致を示した。

ABSTRACT

We analyze diffusion of particles on a two dimensional square lattice. Each lattice site contains an arbitrary number of particles. Interactions affect particles only in the same site, and are macroscopically represented by the excess chemical potential. In a recent work, a general expression for transition rates between neighboring cells as functions of the excess chemical potential was derived. With transition rates, the mean field tracer diffusivity, $D^ ext{MF}$, is immediately obtained. The tracer diffusivity, $D = D^ ext{MF} f$, contains the correlation factor $f$, representing memory effects. An analysis of the joint probability of having given numbers of particles at different sites when a force is applied to a tagged particle allows an approximate expression for $f$ to be derived. The expression is applied to soft core interaction (different values for the maximum number of particles in a site are considered) and extended hard core.

研究の動機と目的

• 2次元格子上におけるトレーサー拡散の相関係数 f の一般化された解析的フレームワークを構築すること。
• 粒子相互作用に起因する空間相関を組み込むことで、平均場拡散係数理論を拡張すること。
• 過剰化学ポテンシャルと熱力学的要因を用いて、ソフトコアおよび拡張ハードコア相互作用を巨視的にモデル化すること。
• 異なる粒子排除領域および密度において、理論的相関係数をモンテカルロシミュレーションと照合すること。

提案手法

• 詳細なバランスとウィドム挿入公式を用いて、格子サイト間の遷移率を導出し、過剰化学ポテンシャルでパrameter化する。
• 平均場トレーサー拡散係数 DMF = D0 / Γ を導入し、Γ は粒子密度および化学ポテンシャルから導かれる熱力学的要因である。
• 外部力が小さい条件下での遷移率非対称性を用いて、粒子数の同時確率分布を用いて相関係数 f を関数として定式化する。
• ソフトコア（1サイトあたり Ω 粒子）および拡張ハードコア（b を用いた排除領域）相互作用にこの手法を適用し、有効な充填率 ξ を定義する。
• 2次近似 f = 1 − c₁ξ + c₂ξ² を用い、異なる b 値に対してモンテカルロデータにフィットするようにパラメータ c₁, c₂ を決定する。
• 全トレーサー拡散係数を D/D₀ = f / Γ = (1 − c₁ξ + c₂ξ²)(1 − ξ/ξ_max) として導出する。ここで ξ_max = ρ_max b / 4 である。

実験結果

リサーチクエスチョン

• RQ1 格子拡散における記憶効果を、相関係数 f の理論的表現に体系的に組み込む方法は何か？
• RQ2 2次元正方形格子上での粒子密度および相互作用タイプ（ソフトコア対拡張ハードコア）に応じて、相関係数 f はどのように変化するか？
• RQ3 巨視的熱力学的量のみを用いて、異なる相互作用モデルに一般化可能な f の式を導出できるか？
• RQ4 調整可能なパラメータを有する理論的 f は、異なる粒子排除領域（b）および密度において、シミュレーション結果とどの程度一致するか？

主な発見

• 理論的式 D/D₀ = (1 − c₁ξ + c₂ξ²)(1 − ξ/ξ_max) は、b（排除領域サイズ）の異なる値において、シミュレーションデータに合わせて c₁ と c₂ を調整することで、トレーサー拡散係数を正確に予測している。
• ソフトコア相互作用では、Ω が増加するにつれて相関係数 f は 1 に近づき、粒子の移動に伴うサイトの擾乱が最小限であるため、記憶効果が小さくなる。
• 拡張ハードコア相互作用では、最大充填率 ξ_max は b の増加に伴い減少し、b=1 のとき 1 から b=37 のとき 1/12 にまで低下する。これは、より大きな排除領域を反映している。
• 相関係数 f は充填率 ξ の増加に伴い減少し、特に b が小さい場合（例：b=1,5）に高密度で強い記憶効果が顕著に現れる。
• モンテカルロシミュレーションは理論モデルと優れた一致を示し、すべてのテストされた相互作用タイプおよび密度において、このアプローチの妥当性が検証された。
• この手法は、f の非単調的挙動、特に負の c₂ 値（例：b=21, c₂=−0.41）を正しく捉えており、高密度領域での f(ξ) の曲率反転を示している。