量子效应在高密度等离子体中的体现

在传统等离子体物理中，粒子往往被视为经典点电荷，统计分布采用麦克斯韦-玻尔兹曼（MB）形式。然而，当等离子体的电子密度达到或超过 (10^{23},\text{cm}^{-3}) 且温度降至几千电子伏特以下时，量子力学效应不再是微小修正，而会显著改变等离子体的宏观行为。本文系统阐述 高密度等离子体 中的主要量子效应，重点介绍简并压、量子隧穿、波函数相干以及量子流体模型（QHD），并给出典型物理情景的数值示例。

1. 量子简并效应

1.1 费米-狄拉克统计

在高密度条件下，电子的德布罗意波长 (\lambda_{!D}=h/\sqrt{2m_{e}k_{B}T}) 与粒子间距可相当，电子不再服从 MB 分布，而遵循费米-狄拉克（FD）分布：

[
f(\varepsilon)=\frac{1}{\exp!\big[(\varepsilon-\mu)/k_{B}T\big]+1},
]

其中 (\mu) 为化学势。若 (T\ll T_{F})（费米温度），电子呈强简并状态，系统的压强主要来源于 简并压：

[
P_{\text{deg}}=\frac{(3\pi^{2})^{2/3}\hbar^{2}}{5m_{e}}n_{e}^{5/3}.
]

该压强与温度无关，决定了白矮星、金属靶在惯性约束聚变（ICF）中的结构稳定性。

1.2 简并长度与临界密度

简并长度定义为

[
r_{s}=a_{0}\left(\frac{3}{4\pi n_{e}}\right)^{1/3},
]

(a_{0}) 为玻尔半径。若 (r_{s}\lesssim 2)，电子相互作用进入强耦合区，需要量子蒙特卡罗或密度泛函理论（DFT）进行精确描述。

2. 量子隧穿与波函数相干

2.1 隧穿增强的碰撞率

在经典视角下，低能电子无法克服库伦势垒进行近距离碰撞；但量子隧穿使得有效碰撞截面增大，碰撞频率可用 WKB 近似 估计：

[
\sigma_{\text{tun}}\approx\sigma_{\text{cl}}\exp!\Big[-2\int_{r_{1}}^{r_{2}}\sqrt{\frac{2m_{e}}{\hbar^{2}}\big(V(r)-E\big)},dr\Big],
]

其中 (V(r)) 为相互势能，(E) 为相对动能。该效应在 固体密度等离子体（(n_{e}\sim10^{23},\text{cm}^{-3})）的激光加热过程中尤为显著，导致能量吸收率提升数十个百分点。

2.2 相干长度与波束聚焦

波函数相干长度

[
L_{\text{coh}}=\frac{\hbar v_{e}}{k_{B}T}
]

决定了等离子体波动的相位保持范围。若激光脉冲宽度小于 (L_{\text{coh}}/c)，则 相干散射（如布里渊散射）可被观测，用于诊断高密度等离子体的电子温度与密度。

3. 量子流体模型（QHD）

经典流体方程缺少量子压强与 Bohm 势项。QHD 在连续性方程和动量方程中加入 量子势：

[
\begin{aligned}
&\frac{\partial n}{\partial t}+\nabla!\cdot!(n\mathbf{u})=0,\
&m\left(\frac{\partial \mathbf{u}}{\partial t}+\mathbf{u}!\cdot!\nabla\mathbf{u}\right)
=-\nabla!\big(P_{\text{deg}}+P_{\text{th}}\big)
+\frac{e}{m}\mathbf{E}
+\frac{\hbar^{2}}{2m}\nabla!\left(\frac{\nabla^{2}\sqrt{n}}{\sqrt{n}}\right).
\end{aligned}
]

• 第一项 为简并压（或热压）梯度。
• 第二项 为电场驱动。
• 第三项 为 Bohm 势，体现量子扩散与波动色散。

该模型能够捕捉 量子朗道阻尼、量子皮肤效应 等现象，已被广泛用于模拟固体靶在超快激光照射下的等离子体形成。

3.1 示例：线性色散关系

对均匀背景进行线性化，可得等离子体波的色散式

[
\omega^{2}= \omega_{p}^{2}+\frac{3}{5}v_{F}^{2}k^{2}+\frac{\hbar^{2}k^{4}}{4m_{e}^{2}},
]

其中 (\omega_{p}=\sqrt{4\pi n_{e}e^{2}/m_{e}}) 为等离子体频率，(v_{F}) 为费米速度。第三项即来源于 Bohm 势，导致在 短波长（(k\lambda_{D}\gg1)） 区域出现显著的 量子色散。

4. 实例与数值模拟

在实际数值实现中，常采用 分数阶有限差分 处理 Bohm 项的高阶导数，或使用 伪谱方法 提高波数空间的分辨率。下列伪代码展示了在一维 QHD 中计算 Bohm 势的核心步骤：

import numpy as np

dx = L / N
psi = np.sqrt(n)               # 取模
lap_psi = (np.roll(psi, -1) - 2*psi + np.roll(psi, 1)) / dx**2
bohm = (hbar**2 / (2*m_e)) * (lap_psi / psi)

该片段可直接嵌入时间步进循环，配合守恒的电荷连续方程即可完成全动态模拟。

5. 小结

高密度等离子体的量子效应主要体现在：

1. 简并压：支配低温高密度电子的热力学行为。
2. 量子隧穿：显著提升低能碰撞截面，影响能量沉积。
3. 波函数相干：决定散射与诊断的相位信息。
4. Bohm 势：在流体方程中引入波动色散，产生量子色散与朗道阻尼。

通过 量子流体模型 与 第一性原理数值方法 的结合，能够在实验可达的参数空间内精确预测等离子体的结构、波动以及能量传输特性，为惯性约束聚变、超快激光等前沿应用提供可靠的理论支撑。