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

本文利用 Mori-Zwanzig-Forster 投影算符技术，通过引入总质量密度、组分浓度、总动量密度和能量密度，推导了包含对流与扩散动力学的非等温二元系统扩展动态密度泛函理论（EDDFT），并给出了硬球系统的精确熵与自由能泛函，同时探讨了其流体动力学极限、与玻璃化转变模式耦合理论的关系以及声速的正确预测。

要点

这篇文章介绍了一种名为**“扩展动态密度泛函理论”（EDDFT）的新数学模型。为了让你轻松理解，我们可以把这篇论文的研究对象想象成“一锅正在被搅拌的、冷热不均的浓汤”**。

1. 以前的模型：只关注“哪里人多”，忽略了“怎么动”和“温度”

想象一下，以前科学家研究这锅汤（比如由胶体粒子和溶剂组成的混合物）时，用的是一种叫“标准动态密度泛函理论”（DDFT）的工具。

• 它的局限： 这个工具就像是一个只数人头的统计员。它只能告诉你：“在这个位置，有多少个粒子？”（浓度）。
• 它的问题： 它假设粒子只是像蚂蚁一样漫无目的地慢慢爬行（扩散）。它完全忽略了两个关键因素：
  1. 流动（对流）： 如果有人在搅拌汤，粒子会跟着汤流一起跑，而不是慢慢爬。
  2. 温度（热量）： 如果汤的一边热一边冷，粒子会怎么跑？热量怎么传递？
  3. 惯性： 粒子是有质量的，它们动起来有冲劲，停下来也有惯性，以前的模型把这些都简化了。

这就好比你在看一场足球赛，以前的模型只能告诉你“球在哪里”，却完全不知道球是被人踢飞的（有速度、有惯性），也不知道风（温度梯度）是怎么吹动球的。

2. 新模型（EDDFT）：给模型装上了“大脑”和“感官”

这篇论文的作者们（来自德国多所大学的物理学家）开发了一个升级版的模型。他们给这个“统计员”装上了更高级的感官，让它能同时关注四个核心变量：

1. 总质量密度： 汤有多稠？
2. 浓度： 某种特定粒子（比如胶体）有多少？
3. 动量密度（速度）： 汤流得有多快？方向在哪？（这是最关键的新增项，让模型能处理“流动”和“惯性”）。
4. 能量密度（温度）： 哪里热？哪里冷？

核心比喻：
以前的模型像是在看一张静态的照片，只能数数人。
现在的模型像是在看4K 高清实时直播，不仅能数人，还能看到人跑得多快、往哪跑，甚至能感觉到空气的冷热变化。

3. 他们是怎么做到的？（Mori-Zwanzig-Forster 投影算子技术）

听起来很吓人，其实原理很简单。
想象你有一万个粒子在疯狂运动，你不可能追踪每一个粒子的每一个动作（那太累了，计算机也跑不动）。

• 以前的做法： 强行简化，忽略大部分细节。
• 新做法（投影技术）： 就像是用一个智能过滤器。它把那些无关紧要的、瞬间变化的微观细节“过滤”掉，只保留那些慢速的、宏观的、重要的信息（比如整体的流动和温度）。
• 这就好比你在嘈杂的酒吧里，虽然周围有几千人在说话（微观粒子），但你只关心“乐队在唱什么”和“舞池里的人在怎么跳舞”（宏观变量）。这个技术就是帮你从噪音中提取出旋律的“降噪耳机”。

4. 这个新模型有什么用？（实际应用场景）

作者们不仅推导出了公式，还证明了它在几个关键方面非常厉害：

• 算得准“声速”：

  • 比喻： 想象你在汤里扔一块石头，产生的波纹（声波）传播速度是多少？
  • 以前的错误： 以前的模型因为忽略了温度变化（声波传播其实是绝热过程，温度会变），算出来的声速是错的。
  • 现在的正确： 新模型考虑了热量和温度的耦合，算出来的声速和真实物理世界完全一致。这就像以前算声音传播像算“在真空里”，现在算的是“在真实大气中”。

• 模拟复杂流动：

  • 它可以用来模拟水泥的流动、金属合金的凝固，甚至是病毒在空气中的传播（比如新冠气溶胶）。
  • 以前这些情况太复杂，因为涉及到了“冷热不均”和“快速流动”，新模型能同时处理这些，就像给工程师提供了一把万能瑞士军刀。

• 连接微观与宏观：

  • 它像一座桥梁，一头连着微观的原子碰撞，另一头连着宏观的流体力学方程（纳维 - 斯托克斯方程）。这让科学家既能看到“森林”（整体流动），也能看到“树木”（粒子相互作用）。

5. 总结

这篇论文的核心贡献是：我们终于有了一个数学工具，既能描述粒子“怎么扩散”，又能描述它们“怎么流动”，还能描述“温度如何影响这一切”。

• 以前： 只能看静态的、慢吞吞的、温度恒定的世界。
• 现在： 可以模拟动态的、有惯性的、冷热不均的真实世界。

这就好比从黑白默片升级到了3D 彩色有声电影，让我们能更真实地预测和理解自然界中复杂的流体和混合物的行为。这对于工业制造（如更好的合金、水泥）和科学研究（如理解玻璃化转变、疾病传播）都有着巨大的潜力。

技术摘要

这是一篇关于**包含动量密度的非等温二元系统扩展动力学密度泛函理论（EDDFT）**的学术论文详细技术总结。

1. 研究背景与问题 (Problem)

• 现有理论的局限性： 标准的动力学密度泛函理论（DDFT）主要描述胶体悬浮液中胶体粒子的局部浓度演化，成功应用于扩散动力学。然而，标准 DDFT 存在两个主要缺陷：
  1. 它忽略了对流动力学（convective dynamics），即无法描述溶剂流动场对粒子的输运作用，这在存在显著流动场（如工业应用、多相流）时至关重要。
  2. 它通常假设系统是等温的，无法处理温度梯度引起的热输运现象（如 Ludwig-Soret 效应和 Dufour 效应）。
• 工业与科学需求： 在金属合金、水泥流动、气溶胶传播（如 COVID-19 病毒传播）以及多相流等实际应用中，必须同时考虑惯性效应（动量）、温度梯度以及不同组分（如胶体与溶剂）之间的相互作用。
• 核心挑战： 如何从微观层面出发，构建一个统一的理论框架，能够同时描述质量密度、浓度、动量密度（速度场）和能量密度的非等温演化，并正确包含热力学涨落。

2. 方法论 (Methodology)

本文采用Mori-Zwanzig-Forster 投影算子技术 (MZFT) 来推导新的理论框架。

• 微观基础： 考虑由 $N_c$ 个胶体粒子和 $N_s$ 个溶剂原子/分子组成的二元混合物。系统由哈密顿量描述，包含动能、外部势能和粒子间的对相互作用势（硬球模型）。
• 相关变量选择： 选取四个慢变量作为宏观描述的基础：
  1. 总质量密度 ($\rho_m$)
  2. 胶体粒子的局部浓度 ($c$)
  3. 总动量密度 ($\vec{g}$)
  4. 能量密度 ($\varepsilon$)
• 投影过程： 利用 MZFT 将微观相空间动力学投影到上述相关变量空间。
  • 定义相关概率密度 $\rho(t)$，使其仅依赖于宏观变量的平均值。
  • 引入共轭变量（基于熵泛函 $S$ 的变分导数），包括温度倒数 $\beta$、化学势等。
  • 推导精确的输运方程，包含可逆部分（对流项）和不可逆部分（耗散项，由记忆核描述）。
• 近似处理：
  • 马尔可夫近似 (Markovian Approximation)： 假设相关变量演化较慢，将记忆核近似为局域时间的耗散项，从而得到广义的 EDDFT 方程。
  • 绝热近似 (Adiabatic Approximation)： 利用自由能泛函的变分导数来近似压力张量，从而封闭动量方程。
  • 硬球模型： 针对硬球系统，利用基本测度理论 (FMT) 构建精确的熵泛函和自由能泛函。

3. 关键贡献 (Key Contributions)

1. 构建了通用的 EDDFT 方程组： 推导了一组包含总质量密度、浓度、动量密度和能量密度的耦合偏微分方程（方程 49-52）。这组方程是标准 DDFT 和纳维 - 斯托克斯（Navier-Stokes）方程的推广。
2. 推导了非等温硬球系统的精确熵泛函： 针对硬球混合物，推导出了包含动量密度和能量密度的非等温熵泛函 $S[\rho_m, c, \vec{g}, \varepsilon]$ 的解析形式（方程 81）。这使得理论可以直接利用成熟的 DFT 方法处理硬球相互作用。
3. 建立了速度与动量密度的热力学联系： 证明了动量密度的热力学共轭变量 $\vec{g}^\natural$ 直接对应于流体的速度场 $\vec{v}$，即 $\vec{v} = \vec{g}/\rho_m$，从而在微观理论中自然导出了流体速度。
4. 导出了正确的声速： 证明了该理论能够正确描述声波传播（绝热过程），并给出了与热力学一致的正确声速公式，解决了早期惯性 DDFT 因假设等温而导致的声速计算错误问题。
5. 建立了 H 定理： 证明了在该框架下，系统的熵产生率非负，满足热力学第二定律，确立了理论的热力学一致性。

4. 主要结果 (Results)

• 广义输运方程：
  • 质量守恒： $\dot{\rho}_m = -\nabla \cdot \vec{g}$
  • 浓度演化： 包含扩散项（由浓度梯度驱动）和热扩散项（由温度梯度驱动，即 Soret 效应）。
  • 动量演化： 包含对流项、压力梯度项、外力项以及粘性耗散项。压力张量 $\Pi$ 通过自由能泛函的变分导数给出。
  • 能量演化： 包含对流输运、热传导以及由浓度梯度引起的热扩散耦合项（Dufour 效应）。
• 流体动力学极限 (Hydrodynamic Limit)： 在小波矢和小频率极限下，方程退化为经典的悬浮液流体动力学方程（方程 124-127），其中扩散系数、粘度和热导率由微观记忆核积分给出。
• 声速计算： 通过对线性化方程的分析，推导出声速 $c_s$ 的表达式：
  $$c_s^2 = \frac{\gamma}{\rho_{m,0} \chi_T}$$
  其中 $\gamma$ 是比热容比，$\chi_T$ 是等温压缩率。这证明了 EDDFT 能够正确捕捉声波作为绝热过程的物理本质。
• 与模式耦合理论 (MCT) 的关系： 讨论了 EDDFT 与玻璃化转变模式耦合理论的联系。指出 EDDFT 通常采用马尔可夫近似（忽略记忆效应），而 MCT 保留了记忆核的非马尔可夫特性，两者在描述玻璃化转变时是互补的。

5. 意义与展望 (Significance)

• 理论突破： 该工作填补了微观统计力学与宏观流体力学之间的空白，提供了一个统一的框架，能够同时处理扩散、对流、热输运和相变。
• 工业应用潜力： 该理论适用于描述复杂的多相流系统，如金属合金凝固、水泥浆体流动、以及气溶胶在空气中的传播（特别是涉及温度梯度和惯性效应的场景）。
• 物理准确性： 通过包含能量密度和热涨落，修正了以往惯性 DDFT 在声速计算和热力学一致性上的缺陷。
• 未来方向： 论文指出未来的扩展方向包括引入取向自由度（用于液晶和活性物质）、结合化学反应动力学，以及进一步研究非马尔可夫效应在玻璃化转变中的作用。

总结： 这篇文章通过严谨的投影算子技术，成功构建了一个包含动量和能量变量的非等温二元系统扩展动力学密度泛函理论。它不仅统一了扩散和对流动力学，还通过精确的熵泛函推导和声速验证，确立了该理论在描述复杂软物质和流体系统非平衡态演化中的准确性和广泛适用性。