[論文レビュー] Transport of Charged Particles: Entropy Production and Maximum Dissipation Principle

Chia-Yu Hsieh, YunKyong Hyon+3 more|arXiv (Cornell University)|Jul 31, 2014

Spectroscopy and Quantum Chemical Studies参考文献 18被引用数 4

ひとこと要約

本稿では、イオン種間の相対速度差に起因する追加の散逸を組み込んだ修正されたポアソン-ナーンスト・プランク（PNP）系を提案する。エネルギー的変分的アプローチを用い、元のPNPと同一の平衡状態を保ちつつ、追加のエントロピー生成に起因する異なる動的挙動を示す。ガラーキン法とシュアウトの不動点定理を用いた厳密な局所的解存在定理が証明され、有限要素法によるシミュレーションでも、同一の平衡状態を示すが一時的な挙動が発散的であることが確認された。

ABSTRACT

In order to describe the dynamics of crowded ions (charged particles), we use an energetic variation approach to derive a modified Poisson-Nernst-Planck (PNP) system which includes an extra dissipation due to the effective velocity differences between ion species. Such a system is more complicated than the original PNP system but with the same equilibrium states. Using Schauder's fixed-point theorem, we develop a local existence theorem of classical solutions for the modified PNP system. Different dynamics (but same equilibrium states) between the original and modified PNP systems can be represented by numerical simulations using finite element method techniques.

研究の動機と目的

• 標準的なPNP理論が排他体積効果や高濃度イオンによる非理想状態で失敗する、混雑した非理想環境におけるイオン輸送をモデル化すること。
• 熱力学的に整合性のある変分原理を用いて、イオン種間の相対運動に起因する追加散逸をPNPフレームワークに組み込むこと。
• 非線形性が増加した修正PNP系の厳密な局所的解存在理論を確立すること。
• 同じ境界条件と初期データのもとで、修正PNP系と元のPNP系の動的挙動を数値的に比較し、平衡状態が同一であるにもかかわらず一時的挙動に顕著な差が生じることを明らかにすること。

提案手法

• 保存力には最小作用原理（LAP）を、散逸力には最大散逸原理（MDP）を組み合わせたエネルギー的変分的アプローチを採用する。
• 散逸関数に、陽イオンと陰イオンの相対速度の二乗に比例する項を追加する：$ \frac{k_B T}{D_{n,p}} c_n c_p |\vec{u}_n^* - \vec{u}_p^*|^2 $、これは種間の抵抗を表す。
• 修正されたPNP系を、イオン濃度および静電ポテンシャルに対する非線形で退化する放物型偏微分方程式系として導出する。
• 適切な関数空間枠組みにおいて、ガラーキン法とシュアウトの不動点定理を用いて、古典的解の局所的存在を証明する。
• 時間に implicit な離散化を用いた有限要素法（FEM）スキームを実装し、安定性とエネルギー則の近似を保証する。
• 境界条件と初期データを同一にしたもとで、修正PNP系と元のPNP系の数値解を比較し、追加散逸の影響を明確に分離する。

実験結果

リサーチクエスチョン

• RQ1 相対速度に依存する散逸項の導入が、混雑したイオン系におけるイオン輸送の動的挙動にどのように影響を与えるか。
• RQ2 追加の種間ドラッグ項を含む修正PNP系が、元のPNP系と同一の平衡状態を保持できるか。
• RQ3 修正PNP系の数学的構造は何か。非線形性と潜在的な退化性が増加しても、局所的古典的解が存在するか。
• RQ4 同じ平衡状態と境界条件のもとで、修正PNP系と元のPNP系の過渡的挙動にどのような差が生じるか。
• RQ5 追加の散逸項が、イオン輸送シミュレーションにおけるエントロピー生成率およびエネルギー散逸プロファイルにどの程度の影響を与えるか。

主な発見

• 修正PNP系は、陽イオンと陰イオン間の相対速度に依存する新たな散逸項を導入し、元のPNPモデルに存在しない種間の抵抗効果を捉える。
• 複雑性と非線形性が増加しても、修正PNP系は元のPNP系と同一の平衡状態を有する。両者とも同じ自由エネルギー関数によって決定される。
• ガラーキン法とシュアウトの不動点定理を用いて、修正PNP系の古典的解の局所的存在定理が厳密に確立された。
• 有限要素法による数値シミュレーションにより、同じ境界条件下で、修正PNP系と元のPNP系が同一の平衡状態に到達することが確認された。
• 両者には顕著に異なる過渡的挙動が見られ、特に時間依存の散逸プロファイルにおいて、追加の散逸項により修正系のエントロピー生成率が高くなることが顕著に示された。
• 図2における散逸率 $\triangle^*$（修正系）と $\triangle$（元の系）の比較により、追加の散逸が平衡状態を変えることなく、緩和経路を変化させることを示している。