2025-11-24T22:16:17.465368

Molecular Dynamics Simulation of Hydrodynamic Transport Coefficients in Plasmas

Damman, LeVan, Baalrud

Molecular dynamics (MD) simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green-Kubo formalism over a broad range of Coulomb coupling strength, $0.01 \leq Î\leq 140$. Emphasis is placed on testing standard results of the Chapman-Enskog solution in the weakly coupled regime ($Î\ll 1$) using these first-principles simulations. As expected, the results show good agreement for $Î\lesssim 0.1$. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green-Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal and thermoelectric coefficients computed in the Chapman-Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime ($Î\gg 1$). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

academic

Molecular Dynamics Simulation of Hydrodynamic Transport Coefficients in Plasmas

基本信息

论文ID: 2511.09787
标题: Molecular Dynamics Simulation of Hydrodynamic Transport Coefficients in Plasmas
作者: Briggs Damman, Jarett LeVan, Scott D. Baalrud (University of Michigan)
分类: physics.plasm-ph
发表时间: 2025年11月14日
论文链接: https://arxiv.org/abs/2511.09787

摘要

本研究使用分子动力学(MD)模拟计算双组分等离子体在排斥性库仑势下的输运系数。通过Green-Kubo形式在广泛的库仑耦合强度范围(0.01 ≤ Γ ≤ 140)内计算了热导率、电导率、电热系数、热电系数和剪切粘度。研究重点测试了Chapman-Enskog解在弱耦合区域(Γ≪1)的标准结果。结果显示在Γ≲0.1时与理论有很好的一致性，但这种一致性只有在仔细定义线性本构关系时才能实现。研究表明，在强耦合区域(Γ≫1)，只有输运系数的动力学分量是有意义的。

研究背景与动机

1. 核心问题

等离子体流体力学方程的闭合需要准确的输运系数模型。传统上这由Chapman-Enskog(CE)求解Boltzmann方程提供，但这些表达式缺乏充分的实验验证，特别是在广泛的等离子体条件范围内。

2. 问题的重要性

实验困难: 在等离子体中进行精确的输运系数测量极其困难，现有测量精度有限、覆盖系数少、条件范围窄
理论局限: 许多实验测试与碰撞输运理论预测不一致，可能因为难以在热稀等离子体中实现真正的碰撞区域
应用需求: 准确的输运系数对磁流体动力学(MHD)方程的闭合至关重要

3. 现有方法的局限性

一组分等离子体(OCP)研究: 已有大量OCP和Yukawa OCP的MD模拟，但因无扩散贡献，与电子-离子混合物的输运性质关联性差
定义不一致: 动力学理论和非平衡热力学对输运系数的定义不同，导致比较困难
量子效应: 物理致密系统的MD模拟需考虑量子力学效应(如Pauli阻塞)，密度泛函理论分子动力学(DFT-MD)计算代价昂贵

4. 研究动机

使用第一性原理MD模拟提供从牛顿运动方程的形式精确计算
测试等离子体动力学和输运理论的假设
为强耦合条件提供基准数据
量化标准碰撞等离子体理论的有效条件范围

核心贡献

系统性基准测试: 首次系统地使用MD模拟验证双组分等离子体CE解在弱耦合区域(Γe≲0.1)的准确性
定义一致性框架: 明确建立了动力学理论(CE方法)和非平衡热力学(Green-Kubo方法)中输运系数定义之间的对应关系，解决了长期被忽视的定义不一致问题
广泛参数空间数据: 提供了跨越0.01≤Γe≤140范围的五个输运系数(σ, λ, ϕ, φ, η)的完整数据集
强耦合区域洞察: 揭示了在强耦合区域，排斥性库仑系统中热流的势能和维里分量发散，只有动力学分量有物理意义
拟合公式: 提供了便于使用的经验拟合公式，适用于弱耦合和强耦合区域
OCP与TCP关系: 证明了剪切粘度可以从一组分系统准确获得，而其他输运系数需要双组分模拟

方法详解

任务定义

输入:

库仑耦合参数Γe (表征等离子体耦合强度)
粒子数N、数密度ne、温度T

输出:

五个输运系数: 电导率σ、热导率λ、电热系数ϕ、热电系数φ、剪切粘度η

约束条件:

使用排斥性库仑势避免非物理束缚态
经典MD框架(不考虑量子效应)

理论框架

1. Green-Kubo关系

输运系数通过平衡态流量涨落的自相关函数计算:

$\sigma = \frac{V}{3k_BT}\int_0^\infty dt \langle j(t)\cdot j(0)\rangle$

$\lambda = \frac{V}{3k_BT^2}\int_0^\infty dt \langle q(t)\cdot q(0)\rangle$

$\eta = \frac{V}{6k_BT}\sum_{i=1}^3\sum_{j\neq i}^3\int_0^\infty dt \langle\hat{P}_{ij}(t)\hat{P}_{ij}(0)\rangle$

其中流量通过Irving-Kirkwood程序表达为粒子轨迹:

$j = \frac{1}{V}\sum_{i=1}^N q_iv_i$

$q = \frac{1}{V}\sum_{i=1}^N\left[v_i\frac{1}{2}m_i|v_i|^2 + \frac{1}{2}v_i\sum_{j\neq i}\phi_{ij} + \frac{1}{2}\sum_{j\neq i}(r_i\cdot v_i)\frac{\partial\phi_{ij}}{\partial r_i}\right]$

热流包含动力学、势能和维里三个分量。

2. Chapman-Enskog理论

提供输运系数的显式表达式，基于电子库仑碰撞时间τe:

$\tau_e = \frac{3}{2\sqrt{2\pi}}\frac{(4\pi\epsilon_0)^2\sqrt{m_e}(k_BT)^{3/2}}{n_eq_e^4\ln\Lambda}$

输运系数为:

$\sigma_{CE} = 1.93\frac{n_eq_e^2\tau_e}{2m_e}$
$\lambda_{CE} = 1.02\frac{n_ek_B^2T\tau_e}{m_e}$
$\eta_{CE} = 0.96n_ik_BT\sqrt{\frac{m_i}{2m_e}}\tau_e$

3. 定义对应关系

关键发现是两种形式主义的输运系数定义不同。通过重新组织线性本构关系，建立了对应关系:

$\sigma = \sigma_{CE}$

$\phi = T\phi_{CE} + \frac{5k_BT}{2q_e}\sigma_{CE} = 3.19\frac{k_BTn_eq_e\tau_e}{m_e}$

$\lambda = \lambda_{CE} + T\frac{\phi_{CE}^2}{\sigma_{CE}} + \frac{5k_BT}{q_e}\phi_{CE} + \frac{25k_B^2T}{4q_e^2}\sigma_{CE} = 11.58\frac{n_ek_B^2T\tau_e}{m_e}$