从冻结库仑退相干到热化：通往Balescu–Lenard和朗道动力学的经典与量子密度矩阵途径

《Annals of Physics》：From frozen Coulomb decoherence to thermalization: Classical and quantum density-matrix routes to Balescu–Lenard and Landau kinetics

　　•经典噪声与量子浴需要不同的库仑约化。
•KMS不对称性和反冲从扩散核产生摩擦。
•RPA屏蔽给出Balescu-Lenard和Landau试验粒子极限。
•冻结的库仑无序给出无热化弹性散射。
•布居动力学无法唯一确定量子相干动力学。
引言
一个带电粒子与库仑环境相互作用时，会失去相干性、交换动量，并最终趋于平衡。这些过程通常分开处理：随机介质理论采用双振幅传播子和Bethe-Salpeter核，开放系统理论使用约化密度矩阵，而等离子体动力学依赖于动态屏蔽的碰撞算子。我们考察这些描述如何相互联系，以及潜在物理假设在哪些地方不同。
规定噪声与平衡浴之间的区别至关重要。对经典随机势取平均给出随机幺正动力学，它重新分配动量，但本身并不在有限温度下选择吉布斯分布。平衡浴还具有耗散响应。在量子描述中，涨落与耗散之间的平衡通过有序浴关联子和Kubo-Martin-Schwinger（KMS）关系表达；其经典动力学对应是摩擦与扩散之间的关系 [1], [2], [3], [4], [5]。
动态结构因子为环境涨落与试验粒子跃迁率之间提供了自然联系，如van Hove所建立 [6]，并在平移协变量子线性玻尔兹曼理论中得到发展 [7], [8], [9], [10], [11]。这些理论描述碰撞退相干、动量弛豫和细致平衡；Petruccione和Vacchini还推导了微观布朗摩擦系数 [12]。这里我们聚焦库仑问题。我们将纵向介电响应与共振试验粒子Balescu-Lenard和Landau核相联系，并将所得动力学与冻结库仑场的不同弱无序极限进行比较 [13], [14], [15], [16], [17], [18]。
计算的组织方式是常见的库仑系统-环境耦合
?
{规定经典场→Born-Novikov,平衡量子浴→Born-Redfield/KMS,
?
谱动量转移核
?
{冻结无序退相干和弹性散射,FDT/Vlasov→Balescu-Lenard→Landau.
这两种约化共享一个库仑相互作用顶点和动量转移的谱描述，但不共享同一个主方程。在平衡分支中，我们从同一跃迁核导出摩擦和扩散，包括反冲和领先KMS贡献。冻结分支需要不同的长时间弱无序极限。
我们考虑弱耦合下的非相对论试验粒子。空间局域Wigner方程在领头梯度阶下使用；我们还检验了一个完全正但不唯一的布居动力学到量子相干性的扩展。
节选片段
符号、约定和假设
探针质量为m，电荷为e0，坐标为r，动量为p，能量为εp=p2/(2m)。下面考虑的环境是均匀、平稳且反演对称的；每当使用S(-q,ω)=S(q,ω)时，假定各向同性。
我们通过周期体积V的盒子归一化傅里叶振幅，然后取热力学极限：U(r)=V?1∑q Uq e^{iq·r} 和 V?1∑q→∫d3q/(2π)3。下面所有空间谱都是动量守恒δ函数的系数，而不是包含δ(0)的完整两点平均。
微观模型
共同的微观出发点是
?=?S+?B+?I,
?S=p?2/(2m),
其中库仑相互作用
?I=∫d3q/(2π)3 e^{iq·r?}??q, ?q?=?-q。
对于电势算子φ?q，有?q=e0 φ?q。平均势被吸收到?S中，因此TrB(ρ?B ?q)=0。
经典随机分支：单粒子密度矩阵
对固定实现U，单粒子密度矩阵定义为
ρU(r1,r2,t)=ψU(r1,t) ψ*U(r2,t)。
两个坐标标记同一粒子的右矢和左矢；它们不描述两个粒子。
使用式（12）及其复共轭，我们得到
i? ?t ρU = [-?2/(2m) ?12 + U(r1,t)] ρU - [-?2/(2m) ?22 + U(r2,t)] ρU。
等价地，
?t ρU = L0 ρU - (i/?)[U(r1,t)-U(r2,t)] ρU，
其中自由刘维尔算符L0 = (i?/(2m))(?12-?22)。
引入简写
Ua(t)≡U(ra,t), s1=+1, s2=-1，
使得
U(r1,t)-U(r2,t)=2∑_{a=1} s_a Ua(t)。
则式（15）变为
msub
密度矩阵方程的经典泛函变分
为了应用Novikov恒等式，我们需要泛函导数
δρU(t)/δU(r′,t′)。
这个导数衡量密度矩阵对势局部变化的响应。
经典场的Born-Novikov闭包
如果规定的经典势是零均值高斯的，则以下Novikov恒等式是精确的 [19]：
?Uc(r,t) F[U]? = ∫d3r′ ∫dt′ Kc(r-r′,t-t′) ?δF[U]/δU(r′,t′)?。
取F[U]=ρU(r1,r2,t)。然后，利用泛函导数（20）和Born近似（22）（见下文），我们得到一个闭合表达式。
经典双振幅记忆方程
将式（24）代入式（18），引入记忆时间τ=t-t′，给出
-i? ∑_a s_a ?Ua ρU? = -1/?2 ∑_{a,b} s_a s_b ∫_t0 dτ ∫d3r′ Kc(ra-r′,τ) × G0(τ) [δ(rb-r′) ρ?(t-τ)]。
因此我们得到平均密度矩阵的Bethe-Salpeter型记忆方程 [20]：
?t ρ?(r1,r2,t) = L0 ρ?(r1,r2,t) - 1/?2 ∑_{a,b=1} s_a s_b ∫_t0 dτ ∫d3r′ Kc(ra-r′,τ) × G0(τ) [δ(rb-r′) ρ?(r1,r2,t-τ)]。
方程（26）在时间上非局域，并耦合密度矩阵的两个空间自变量。即使其闭包在场耦合中是二阶的，
量子平衡浴分支：Born-Redfield方程
对于量子浴，相互作用绘景的总态满足
dρ?I_tot/dt = -(i/?)[?I_I(t), ρ?I_tot(t)]。
我们假设ρ?I_tot(0)=ρ?S(0)?ρ?B，其中[?B,ρ?B]=0且TrB(ρB Uq)=0。在HI的二阶，Born因式分解ρ?I_tot(t-τ)?ρ?I_S(t-τ)?ρ?B给出
dρ?I_S(t)/dt = -1/?2 ∫_t0 dτ TrB[?I_I(t), [?I_I(t-τ), ρ?I_S(t-τ)?ρ?B]]。
使用H_I(t)=∫d3q/(2π)3 Vq(t)?Uq(t)、均匀性和第2节定义的强度系数C>,_B(q,τ)[Vq(t), V-q(t-τ)ρ?I_S(t-τ)] + C对于p′=p+?q，式（42）和（52）给出
w_{p→p′}/w_{p′→p}=e^{-β(εp′-εp)}。
因此
feq(p)=a e^{-βp2/(2m)}
通过逐对细致平衡是平稳的。细致平衡本身并不证明收敛。如果非零跃迁率集合连接相关动量态（不可约性/遍历性）且边界项受控，则相对熵产生非正，麦克斯韦分布是唯一的归一化平稳动量
冻结弹性无序作为单独的弱无序极限
对于静态实现uq(t)=uq，其中v?1?uq u-q?=cfr(q)，一阶微扰理论给出有限盒子中一次离散转移的概率
p_{p→p+?q}(t) = cfr(q)/(v?2) |∫_t0 ds e^{iωp,q s}|2。
将无序振幅写为λu，并取弱散射（van hove）动力学极限λ→0, t→∞，且λ2t固定，得到黄金规则率。等价地，在固定二阶强度下形成vp/t后，
lim_{t→∞} (1/t) |∫_t0 ds e^{iωs}|2 = 2π δ(ω)，
给出
wfr(p,q) = (2π/?2) cfr(q) δ(ωp,q)。
因此冻结核不是
小动量转移展开和微观漂移
对动量转移δp=?q远小于f(p)明显变化特征动量尺度pf时展开增益-损失算符。保留前两个kramers-moyal项给出
?t f + (p/m)·?_r f = -?_{p_i}[a_i f] + (1/2) ?_{p_i}?_{p_j}[b_{ij} f]，
其中
a_i = ∫d3q/(2π)3 ? q_i w(p,q),
d_{ij} ≡ b_{ij}/2 = (1/2)∫d3q/(2π)3 ?2 q_i q_j w(p,q)。
将概率流写为
j_i = a_i f - ?_{p_j}(d_{ij} f) = -γ_i f - d_{ij} ?_{p_j} f
给出
γ_i = -a_i + ?_{p_j} d_{ij}。
我们现在从确定d_{ij}的同一平衡谱得到a_i。设
ω0 = (p·q)/m, r_q = q2/(2m), ω_{p,q} = ω0 + ? r_q。
使用
从fdt谱到balescu-lenard和landau方程
上面得到的跃迁率仅需要平衡谱。我们现在使用vlasov/rpa响应来评估弱耦合库仑等离子体的谱。动态屏蔽随后导致试验粒子balescu-lenard张量，并在对数精度下导致其landau极限。
总结与结论
我们遵循了两条从微观库仑耦合到试验粒子动力学的路线。经典高斯场在born闭包后导致双振幅记忆方程；平衡量子浴给出有序关联子和kms一致的跃迁率。静态无序和动力学浴因此需要不同的动力学极限，尽管两者都允许动量转移描述。
对于平衡浴，小转移展开给出扩散和漂移，来自
利益冲突声明
作者声明，他们不知道可能影响本文报告工作的已知竞争经济利益或个人关系。
致谢
导致这些结果的研究得到了高等经济大学基础研究项目（hse-br-2025-007）的资助。
yury a. budkov 对于p′=p+?q，式（42）和（52）给出 w_{p→p′} w_{p′→p}=e^{-β(εp′-εp)}。 因此 feq(p)=A e^{-βp2 (2m)} 通过逐对细致平衡是平稳的。细致平衡本身并不证明收敛。如果非零跃迁率集合连接相关动量态（不可约性 遍历性）且边界项受控，则相对熵产生非正，麦克斯韦分布是唯一的归一化平稳动量 冻结弹性无序作为单独的弱无序极限 对于静态实现uq(t)=Uq，其中V?1?Uq u-q?=Cfr(q)，一阶微扰理论给出有限盒子中一次离散转移的概率 p_{p→p+?q}(t)=Cfr(q)/(V?2) |∫_t0 ds e^{iωp,q s}|2。 将无序振幅写为λu，并取弱散射（van hove）动力学极限λ→0, t→∞，且λ2t固定，得到黄金规则率。等价地，在固定二阶强度下形成vp t后， lim_{t→∞} (1 t) |∫_t0 ds e^{iωs}|2=2π δ(ω)， 给出 wfr(p,q)=(2π/?2) cfr(q) δ(ωp,q)。 因此冻结核不是 小动量转移展开和微观漂移 对动量转移δp=?q远小于f(p)明显变化特征动量尺度Pf时展开增益-损失算符。保留前两个Kramers-Moyal项给出 ?t f + (p m)·?_r f=-?_{p_i}[A_i f] + (1 2) ?_{p_i}?_{p_j}[b_{ij} f]， 其中 a_i=∫d3q/(2π)3 ? q_i w(p,q), d_{ij} ≡ b_{ij} 2=(1/2)∫d3q/(2π)3 ?2 q_i q_j w(p,q)。 将概率流写为 j_i=A_i f - ?_{p_j}(d_{ij} f)=-Γ_i f - d_{ij} ?_{p_j} f 给出 γ_i=-A_i + ?_{p_j} d_{ij}。 我们现在从确定d_{ij}的同一平衡谱得到a_i。设 ω0=(p·q)/m, r_q=q2/(2m), ω_{p,q}=ω0 + ? r_q。 使用 从fdt谱到balescu-lenard和landau方程 上面得到的跃迁率仅需要平衡谱。我们现在使用vlasov rpa响应来评估弱耦合库仑等离子体的谱。动态屏蔽随后导致试验粒子balescu-lenard张量，并在对数精度下导致其landau极限。 总结与结论 我们遵循了两条从微观库仑耦合到试验粒子动力学的路线。经典高斯场在born闭包后导致双振幅记忆方程；平衡量子浴给出有序关联子和kms一致的跃迁率。静态无序和动力学浴因此需要不同的动力学极限，尽管两者都允许动量转移描述。 对于平衡浴，小转移展开给出扩散和漂移，来自 利益冲突声明 作者声明，他们不知道可能影响本文报告工作的已知竞争经济利益或个人关系。 致谢 导致这些结果的研究得到了高等经济大学基础研究项目（hse-br-2025-007）的资助。 yury a.>对于p′=p+?q，式（42）和（52）给出
w_{p→p′}/w_{p′→p}=e^{-β(εp′-εp)}。
因此
feq(p)=a e^{-βp2/(2m)}
通过逐对细致平衡是平稳的。细致平衡本身并不证明收敛。如果非零跃迁率集合连接相关动量态（不可约性/遍历性）且边界项受控，则相对熵产生非正，麦克斯韦分布是唯一的归一化平稳动量
冻结弹性无序作为单独的弱无序极限
对于静态实现uq(t)=uq，其中v?1?uq u-q?=cfr(q)，一阶微扰理论给出有限盒子中一次离散转移的概率
p_{p→p+?q}(t) = cfr(q)/(v?2) |∫_t0 ds e^{iωp,q s}|2。
将无序振幅写为λu，并取弱散射（van hove）动力学极限λ→0, t→∞，且λ2t固定，得到黄金规则率。等价地，在固定二阶强度下形成vp/t后，
lim_{t→∞} (1/t) |∫_t0 ds e^{iωs}|2 = 2π δ(ω)，
给出
wfr(p,q) = (2π/?2) cfr(q) δ(ωp,q)。
因此冻结核不是
小动量转移展开和微观漂移
对动量转移δp=?q远小于f(p)明显变化特征动量尺度pf时展开增益-损失算符。保留前两个kramers-moyal项给出
?t f + (p/m)·?_r f = -?_{p_i}[a_i f] + (1/2) ?_{p_i}?_{p_j}[b_{ij} f]，
其中
a_i = ∫d3q/(2π)3 ? q_i w(p,q),
d_{ij} ≡ b_{ij}/2 = (1/2)∫d3q/(2π)3 ?2 q_i q_j w(p,q)。
将概率流写为
j_i = a_i f - ?_{p_j}(d_{ij} f) = -γ_i f - d_{ij} ?_{p_j} f
给出
γ_i = -a_i + ?_{p_j} d_{ij}。
我们现在从确定d_{ij}的同一平衡谱得到a_i。设
ω0 = (p·q)/m, r_q = q2/(2m), ω_{p,q} = ω0 + ? r_q。
使用
从fdt谱到balescu-lenard和landau方程
上面得到的跃迁率仅需要平衡谱。我们现在使用vlasov/rpa响应来评估弱耦合库仑等离子体的谱。动态屏蔽随后导致试验粒子balescu-lenard张量，并在对数精度下导致其landau极限。
总结与结论
我们遵循了两条从微观库仑耦合到试验粒子动力学的路线。经典高斯场在born闭包后导致双振幅记忆方程；平衡量子浴给出有序关联子和kms一致的跃迁率。静态无序和动力学浴因此需要不同的动力学极限，尽管两者都允许动量转移描述。
对于平衡浴，小转移展开给出扩散和漂移，来自
利益冲突声明
作者声明，他们不知道可能影响本文报告工作的已知竞争经济利益或个人关系。
致谢
导致这些结果的研究得到了高等经济大学基础研究项目（hse-br-2025-007）的资助。
yury a. budkov>