Extended dynamical density functional theory for nonisothermal binary systems including momentum density

Mori-Zwanzig-Forster 射影演算子法を用いて、運動量密度とエネルギー密度を新たな変数として取り込んだ非等温二成分系の拡張動的密度汎関数理論（EDDFT）を導出することで、拡散と対流の両方のダイナミクスを記述し、硬球系における正確なエントロピー・自由エネルギー汎関数を導き、音速の正しい値を得ることを示した。

要約

この論文は、**「複雑な流体（液体や気体）の動きを、より正確に、より詳しく予測するための新しい計算ルール」**を作ったという研究です。

専門用語を避け、日常の例えを使ってわかりやすく解説します。

1. 従来のルール（DDFT）の限界：「歩行者」しか見えていない

これまで科学者たちは、**「動的密度汎関数理論（DDFT）」**というツールを使って、液体やコロイド（泥水のようなもの）の動きを予測していました。
しかし、このツールには大きな弱点がありました。

• 従来のルール： 粒子が「どこにいるか（濃度）」と「どう拡散するか（ゆっくり広がる動き）」しか計算できませんでした。
• 現実の問題： 実際には、風が吹いたり、水流があったりして、粒子が**「流されて動く（対流）」ことがあります。また、「温度の違い」**で動きが変わることもあります。
• 例え： 従来のルールは、**「歩行者の足跡だけを見て、その人がどこへ向かうか予測する」ようなものです。しかし、実際にはその人が「自転車に乗って風に乗って走っている」**場合、足跡だけでは予測できません。

2. 新しいルール（EDDFT）の登場：「自転車とエンジン」を追加

この論文の著者たちは、この弱点を克服するために、**「拡張された動的密度汎関数理論（EDDFT）」**という新しいルールを開発しました。

この新しいルールでは、以下の 4 つの要素を同時に計算に組み込みました。

1. 全体の重さ（質量密度）： どれくらい重いか。
2. 成分の割合（濃度）： コロイドと溶媒（水など）の混ざり具合。
3. 勢いと方向（運動量密度）： **「ここは自転車に乗って走っている！」**という動き。
4. 熱エネルギー（エネルギー密度）： **「ここは暑くて、粒子がバタバタしている！」**という熱の動き。

例え話：

• 従来のルール： 「歩行者（粒子）」が「ゆっくり歩く（拡散）」ことしか考えない。
• 新しいルール： 「歩行者」が「自転車に乗って風に乗って走る（対流）」ことと、「夏の日差しで汗だくになっている（熱）」ことも同時に計算する。

これにより、「金属合金の製造」や「セメントの流れ」、さらには**「ウイルスが空気中でどう広がるか（COVID-19 など）」**といった、複雑で激しく動く現象を、よりリアルにシミュレーションできるようになります。

3. すごいところ：「音の速さ」を正しく計算できる

この研究の最大の功績の一つは、**「音の速さ」**を正しく計算できる点です。

• なぜ重要か？ 音は、空気や液体の中を「圧縮と膨張」の波として伝わります。このとき、温度の変化が重要な役割を果たします（断熱過程）。
• 従来の失敗： 以前のツールは「温度は一定」と勝手に決めていたため、音の速さを計算すると、実際の値とズレていました。まるで「氷点下の雪の上を走る車の速さ」を「真夏のアスファルトのデータ」で計算してしまうようなものです。
• 新しい成功： この新しいルールは「温度が場所によって変わる」ことを考慮しているため、「音の速さ」を物理的に正しい値で導き出すことができました。

4. 具体的な仕組み：「モリ・ツワンジグ・フォースター」という魔法のフィルター

この新しいルールを作るために、著者たちは**「モリ・ツワンジグ・フォースターの射影演算子法」**という高度な数学のテクニックを使いました。

• 例え： 想像してください。数千億個の粒子が暴れ回っている部屋（ミクロの世界）があるとします。それをすべて追いかけるのは不可能です。
• 魔法のフィルター： このテクニックは、「重要な情報（重さ、動き、熱）」だけを抽出して、無関係なノイズを捨て去るフィルターのようなものです。
• これによって、複雑すぎるミクロな世界を、私たちが理解できるマクロな「流体の方程式（ナビエ・ストークス方程式の進化版）」に変換することに成功しました。

5. まとめ：何ができるようになるの？

この新しい理論（EDDFT）を使うと、以下のようなことが可能になります。

• よりリアルなシミュレーション： 温度差がある場所での流体の動き（熱対流など）を正確に描ける。
• 産業への応用： 金属の鋳造、セメントの混合、薬の製剤など、工業的なプロセスの最適化。
• 感染症対策： 換気や空調によるウイルスの飛沫の動きを、温度や気流を考慮して予測する。
• 音響の理解： 液体や気体の中を音がどう伝わるかを、微視的な視点から説明する。

一言で言うと：
「これまでのルールは『静かな池の水面』しか見られなかったが、この新しいルールを使えば、『暴風雨の中を走る川』や『熱い蒸気の中で動く粒子』の動きまで、驚くほど正確に予測できるようになった」という画期的な研究です。

技術概要

この論文「Extended dynamical density functional theory for nonisothermal binary systems including momentum density（運動量密度を含む非等温二成分系の拡張動的密度汎関数理論）」は、コロイド懸濁液や二相流などの非等温ダイナミクスを記述するための新しい理論的枠組みを提案したものです。以下に、問題提起、手法、主要な貢献、結果、および意義について詳細な技術的サマリーを日本語で記述します。

1. 問題提起 (Problem)

標準的な動的密度汎関数理論（DDFT）は、コロイド粒子の局所濃度の拡散ダイナミクスを記述する上で非常に成功していますが、以下の点に限界がありました。

• 対流ダイナミクスの欠如: 標準 DDFT は慣性や対流（流れ場）を無視しており、溶媒の顕著な流れ場がある場合や、金属合金、セメント流、感染症の空気中伝播（COVID-19 など）のように慣性と温度勾配が重要な工業的・物理的現象を記述できません。
• 非等温性の不足: 従来の慣性を含む拡張 DDFT の多くは等温系を仮定しており、熱輸送や熱対流、ルードヴィッヒ・ソレー効果（温度勾配による濃度勾配）などを正しく扱えていませんでした。
• 音速の誤差: 既存の慣性を持つ DDFT 手法では、音波が断熱過程であることを正しく反映できず、音速の値が不正確になる傾向がありました。

これらの課題を解決し、質量密度、濃度、運動量密度、エネルギー密度をすべて変数として含む、より包括的な理論が必要とされていました。

2. 手法 (Methodology)

著者らは、モーリ＝ツワンジグ＝フォスター射影演算子法（Mori-Zwanzig-Forster projection operator technique, MZFT） を用いて、古典的なハミルトニアン力学系から新しい理論を導出しました。

• 変数の選択: 以下の 4 つの「遅い」関連変数（relevant variables）を射影演算子法に適用しました。
  1. 全質量密度 ($\rho_m$)
  2. 一方の種の局所濃度 ($c$)
  3. 全運動量密度 ($\vec{g}$)
  4. エネルギー密度 ($\varepsilon$)
• 理論的枠組み: 微視的な多粒子系のダイナミクスを、これらのマクロ変数に射影することで、閉じた運動方程式を導出しました。この過程で、散逸項（拡散テンソル）と可逆項（対流項）が自然に現れます。
• 硬球モデルへの適用: 具体的な応用可能性を高めるため、硬球系（Hard spheres）に対して、エントロピー汎関数と自由エネルギー汎関数の明示的な形を導出しました。これには、基礎的測定理論（FMT）の結果が活用されました。
• 断熱近似: 圧力項を自由エネルギー汎関数の汎関数微分を用いて近似する「断熱近似」を導入し、実用的な計算を可能にしました。

3. 主要な貢献 (Key Contributions)

この論文の主な貢献は以下の通りです。

1. 新しい EDDFT の導出: 質量、濃度、運動量、エネルギーの 4 つの変数を含む、非等温な二成分系（コロイドと溶媒など）のための拡張動的密度汎関数理論（EDDFT）を初めて体系的に導出しました。
2. 厳密なエントロピー・自由エネルギー汎関数の導出: 硬球系に対して、運動量密度と温度依存性を考慮した厳密なエントロピー汎関数と自由エネルギー汎関数を導出しました。これにより、既存の DFT 手法（FMT など）を非等温・慣性系へ拡張することが可能になりました。
3. H 定理の証明: 導出した EDDFT 方程式系に対して、エントロピー増大則（H 定理）が成り立つことを示しました。これは、理論が熱力学第二法則と整合的であることを保証します。
4. 音速の正確な再現: 運動量密度とエネルギー密度（温度変動）を同時に扱うことで、音波が断熱過程であることを正しく記述し、文献値と一致する正しい音速の値を得られることを示しました。
5. 一般流体力学方程式との接続: 長波長・低周波の極限（流体力学極限）において、この EDDFT が二成分系の一般化されたナビエ・ストークス方程式および熱伝導方程式に帰着することを示しました。

4. 結果 (Results)

• 運動方程式: 導出された EDDFT 方程式（式 49-52）は、連続の式、運動量保存則（ナビエ・ストークス型）、エネルギー保存則、および濃度拡散方程式を拡張した形をしています。これらは、拡散項だけでなく、対流項（$\vec{v} \cdot \nabla$）や圧力テンソル、粘性項、熱拡散項を含みます。
• 硬球系での具体化: 硬球系に対して、圧力テンソルや拡散テンソルを具体的な関数形として記述可能であることを示しました。
• 音速の計算: 一成分系における線形化解析を行い、導出された式から音速 $c_s$ を計算した結果、熱力学関係式（比熱比 $\gamma$ を含む）と完全に一致する値が得られました。これは、温度一定を仮定した従来の慣性 DDFT との決定的な違いです。
• モード結合理論（MCT）との関係: この理論がガラス転移のモード結合理論（MCT）と自然に接続されることを議論しました。MCT は記憶効果を重視するのに対し、EDDFT はマルコフ近似（記憶項の無視）を行う点で補完的ですが、両者は MZFT という共通の枠組みに基づいています。

5. 意義 (Significance)

この研究は、軟物質物理学および流体力学の分野において重要な進展です。

• 工業応用への道筋: 金属合金の凝固、セメントの流動、二相流（気泡の形成など）など、慣性と温度勾配が支配的な工業プロセスのシミュレーションを、微視的な粒子レベルからマクロな流体力学まで一貫して記述できる基盤を提供します。
• 感染症伝播モデルへの応用: 空気中のウイルス飛沫（COVID-19 など）の伝播は、流体の慣性と温度勾配（呼吸による熱）が重要であり、この EDDFT はそのような現象のモデル化に有効です。
• 理論的整合性の向上: 従来の DDFT が抱えていた「音速の誤り」や「熱力学非整合性」を解決し、熱力学第二法則（H 定理）を満たす堅牢な理論体系を構築しました。
• 将来の拡張: この枠組みは、液晶（配向自由度の導入）やアクティブマター（自走粒子）、化学反応を伴う系への拡張も可能であり、複雑な非平衡系の研究における強力なツールとなります。

要約すると、この論文は、微視的な粒子相互作用から出発し、マクロな流体力学と熱力学を自然に包含する、運動量とエネルギーを考慮した次世代の密度汎関数理論を確立した画期的な研究です。