2025-11-24T06:13:17.730458

Nonequilibrium Statistics of Biased Kondo Resonance

Han
Numerical renormalization group (NRG) is formulated for nonequilibrium steady-state by converting finite-lattice many-body eigenstates into scattering states. Extension of the full-density-matrix NRG for a biased Anderson impurity model, simplified by formulating with the original orbital basis as the Hamiltonian, enables detailed studies of the sub-Kondo spectral evolution in the zero-temperature limit, confirming the double-resonance structure at bias of the Kondo energy scale $T_K$. The distribution shows distinct multi-scale spectral features at energy $ω$ below the Kondo scale ($ω\lesssim T_K$) and near the bias ($ω\gtrsim V$), leading to the nonequilibrium temperature $T_{\rm loc}$ local to the Kondo dot scaling as $k_BT_{\rm loc}\approx V$ for $V\gg T_K$. The current-voltage relation in the low-temperature limit ($T\ll T_K$) deviates from the unitary limit as the bias exceeds the Kondo scale ($V/2\gtrsim T_K$) and reaches the current saturation regime.
academic

Nonequilibrium Statistics of Biased Kondo Resonance

基本信息

  • 论文ID: 2503.14400
  • 标题: Nonequilibrium Statistics of Biased Kondo Resonance
  • 作者: Jong E. Han (State University of New York at Buffalo)
  • 分类: cond-mat.str-el (强关联电子系统)
  • 发表时间: October 14, 2025
  • 论文链接: https://arxiv.org/abs/2503.14400

摘要

本文通过将有限格点多体本征态转换为散射态，为非平衡稳态制定了数值重整化群(NRG)方法。通过使用原始轨道基作为哈密顿量的表述，扩展了全密度矩阵NRG用于偏压Anderson杂质模型，使得能够详细研究零温极限下的亚Kondo谱演化，确认了在Kondo能标TKT_K偏压下的双共振结构。分布显示了在Kondo能标以下(ω≲TK\omega \lesssim T_K)和偏压附近(ω≳V\omega \gtrsim V)的独特多尺度谱特征，导致Kondo量子点局域的非平衡温度TlocT_{\rm loc}在V≫TKV \gg T_K时按kBTloc≈Vk_BT_{\rm loc} \approx V标度。在低温极限(T≪TKT \ll T_K)下，当偏压超过Kondo能标(V/2≳TKV/2 \gtrsim T_K)时，电流-电压关系偏离幺正极限并达到电流饱和区域。

研究背景与动机

问题描述

  1. 核心挑战: 求解偏压量子系统中的非平衡统计力学问题一直是近几十年来的巨大挑战，特别是对于强关联量子杂质模型。
  2. 重要性:
    • Kondo问题作为强关联非平衡技术的试验台，对理解量子多体系统的非平衡行为至关重要
    • 量子杂质模型在通过动力学平均场理论求解凝聚态极限中发挥关键作用
  3. 现有方法局限性:
    • 非平衡杂质求解器的理论严格性和数值效率尚未达到令人满意的水平
    • 计算困难限制了对Kondo单重态在偏压下如何分解为双峰的亚Kondo区域的理解
    • 相比于模型中的谱特征，对谱的非平衡统计理解不足
  4. 研究动机:
    • 基于成熟的NRG方法构建可访问的稳态非平衡扩展
    • 提供对偏压下Kondo共振演化的清晰理解
    • 开发适用于复杂杂质相互作用的通用方法

核心贡献

  1. 方法创新: 提出了一种新的非平衡NRG方法，通过将有限格点本征态转换为散射态来处理稳态非平衡问题
  2. 技术简化: 使用原始轨道基表述哈密顿量，简化了实现过程，使其更适合扩展到广义问题
  3. 物理发现:
    • 确认了Kondo峰在V/2≈TKV/2 \approx T_K时的分裂行为
    • 发现了多尺度分布特征和居数反转现象
    • 揭示了非平衡温度的标度关系Tloc≈VT_{\rm loc} \approx V
  4. 输运性质: 验证了低温极限下的电流饱和现象，为TK≲V/2≲3TKT_K \lesssim V/2 \lesssim 3T_K区域

方法详解

理论框架

Anderson杂质模型

哈密顿量定义为： H^=∑σϵdndσ+U(nd↑−12)(nd↓−12)+∑αkσϵkcαkσ†cαkσ−tΩ∑αkσ(dσ†cαkσ+cαkσ†dσ)\hat{H} = \sum_\sigma \epsilon_d n_{d\sigma} + U\left(n_{d\uparrow} - \frac{1}{2}\right)\left(n_{d\downarrow} - \frac{1}{2}\right) + \sum_{\alpha k\sigma} \epsilon_k c^\dagger_{\alpha k\sigma} c_{\alpha k\sigma} - \frac{t}{\sqrt{\Omega}} \sum_{\alpha k\sigma} (d^\dagger_\sigma c_{\alpha k\sigma} + c^\dagger_{\alpha k\sigma} d_\sigma)

其中α=L,R\alpha = L,R表示源极和漏极库，UU为库仑排斥参数，ϵd\epsilon_d为杂质能级。

非平衡密度矩阵

稳态非平衡密度矩阵构造为： ρ^neq=1Zneqe−β(H^−VY^)\hat{\rho}_{\rm neq} = \frac{1}{Z_{\rm neq}} e^{-\beta(\hat{H} - V\hat{Y})}

其中Hershfield的Y算符定义为： Y^=12∑kσ(ψLkσ†ψLkσ−ψRkσ†ψRkσ)\hat{Y} = \frac{1}{2}\sum_{k\sigma}(\psi^\dagger_{Lk\sigma}\psi_{Lk\sigma} - \psi^\dagger_{Rk\sigma}\psi_{Rk\sigma})

散射态构造

关键创新在于通过以下关系构造Y算符： Y^=iη−L+iηY^0\hat{Y} = \frac{i\eta}{-\mathcal{L} + i\eta}\hat{Y}_0

其中L\mathcal{L}是Liouvillian算符，Y^0\hat{Y}_0是断开杂质时的初始Y算符。在能量本征态表象中： ⟨α∣Y^∣β⟩=iη−Eα+Eβ+iη⟨α∣Y^0∣β⟩\langle\alpha|\hat{Y}|\beta\rangle = \frac{i\eta}{-E_\alpha + E_\beta + i\eta}\langle\alpha|\hat{Y}_0|\beta\rangle

数值实现

Wilson链表示

  • 使用Wilson链(WC)模拟库，链长NWC=100N_{\rm WC} = 100，重整化参数Λ=1.2\Lambda = 1.2
  • 环境轨道处理：考虑偏压下环境轨道的自由能变化
  • 密度矩阵权重在Wilson链长度n=Mn = M处达到峰值，满足DM≈V/2D_M \approx V/2

全密度矩阵NRG (FDM-NRG)

  • 保持态数Nkept=2000N_{\rm kept} = 2000
  • z平均化使用8个z值减少Wilson链的离散性
  • 每次NRG迭代的Hilbert空间维度为16×Nkept16 \times N_{\rm kept}

实验设置

模型参数

  • 强Kondo区域: U/Γ=10U/\Gamma = 10, Γ=0.1\Gamma = 0.1, TK=0.0044T_K = 0.0044
  • 粒子-空穴对称: ϵd=0\epsilon_d = 0
  • 粒子-空穴非对称: ϵd=U/4\epsilon_d = U/4
  • 温度: T=0.1TKT = 0.1T_K (低温极限)

计算参数

  • 重整化参数: Λ=4\Lambda = 4
  • 谱展宽参数: δ=0.2×max⁡(TK,T,V/2)\delta = 0.2 \times \max(T_K, T, V/2)
  • 无穷小参数: η=0.1×max⁡(TK,T,V/2)\eta = 0.1 \times \max(T_K, T, V/2)

评价指标

  • 谱函数: ρΣ(ω)\rho_\Sigma(\omega) (使用自能技巧计算)
  • 分布函数: f(ω)=Gd<(ω)/[2πiρ(ω)]f(\omega) = G^<_d(\omega)/[2\pi i\rho(\omega)]
  • 非平衡温度: Tloc2=6π2∫−∞∞ω[f(ω)−Θ(−ω)]dωT_{\rm loc}^2 = \frac{6}{\pi^2}\int_{-\infty}^{\infty} \omega[f(\omega) - \Theta(-\omega)]d\omega
  • 电流: 通过图解技术计算

实验结果

谱演化分析

非相互作用验证

  • 在非相互作用共振能级模型中验证了方法的正确性
  • 小偏压下电流遵循I=G0VI = G_0 V，其中G0=1/(2π)G_0 = 1/(2\pi)为幺正量子电导
  • 大偏压下符合精确关系I=(Γ/π)tan⁡−1(V/2Γ)I = (\Gamma/\pi)\tan^{-1}(V/2\Gamma)

Kondo共振分裂

  1. 零偏压: 谱函数在ρΣ(0)≈1/(πΓ)\rho_\Sigma(0) \approx 1/(\pi\Gamma)处满足幺正极限
  2. 临界偏压: 当V/2≈TKV/2 \approx T_K时，Kondo共振开始显示峰分裂
  3. 高偏压: 峰位置对齐在±V/2\pm V/2，但线宽远大于TKT_K，应视为弱关联

非平衡统计特征

分布函数演化

  • 低偏压: 分布函数接近Fermi-Dirac函数
  • 中等偏压: 在∣ω∣≳V|ω| \gtrsim V处发展出大而宽的峰
  • 居数反转: 在分裂Kondo峰之间(−V/2≲ω≲V/2-V/2 \lesssim \omega \lesssim V/2)出现f′(0)<0f'(0) < 0

非平衡温度标度

关键发现：在高偏压极限下 Tloc≈VforV/2≫TKT_{\rm loc} \approx V \quad \text{for} \quad V/2 \gg T_K

这表明局域电子温度与偏压成正比。

输运性质

电流-电压关系

  1. 低温极限 (T=0.1TKT = 0.1T_K):
    • 小偏压下满足幺正极限G0=(2π)−1G_0 = (2\pi)^{-1}
    • V/2≳TKV/2 \gtrsim T_K时开始偏离幺正极限
    • TK≲V/2≲3TKT_K \lesssim V/2 \lesssim 3T_K区域出现电流饱和
  2. 高温 (T=TKT = T_K):
    • 不出现电流饱和现象
    • 证明电流饱和源于强关联效应

相关工作

非平衡量子杂质方法

论文系统回顾了处理Kondo问题的各种非平衡技术：

  • 解析理论: Bethe ansatz方法
  • 图解方法: Landauer公式、微扰理论
  • 数值方法: 实时图解模拟、虚时模拟、tensor网络方法
  • 重整化群: 实时RG、微扰RG、函数RG

与现有NRG方法比较

  • SNRG方法: Anders提出的稳态NRG使用非相互作用散射态构造Y^0\hat{Y}_0
  • 本文方法: 使用原始哈密顿量基，实现更直接，适合扩展到复杂杂质相互作用

结论与讨论

主要结论

  1. 方法验证: 成功将Hershfield散射态理论与NRG方法结合，为非平衡强关联系统提供了可靠的数值工具
  2. 物理图像: 确认了Kondo峰在V/2≈TKV/2 \approx T_K时的分裂，揭示了亚Kondo区域的复杂多尺度结构
  3. 非平衡特征: 发现了居数反转和非平衡温度的线性标度关系，为理解强关联非平衡态提供了新洞察
  4. 输运现象: 验证了电流饱和现象，证实了理论预测

局限性

  1. 计算成本: 由于需要处理两个Wilson链，Hilbert空间增长迅速，限制了保持态数
  2. 参数依赖: 方法对无穷小参数η\eta的选择敏感，需要仔细调节
  3. 温度限制: 主要适用于低温区域T≪TKT \ll T_K

未来方向

  1. 扩展应用: 适用于复杂杂质相互作用模型
  2. 格点理论: 通过动力学平均场理论应用于非平衡格点模型
  3. 实验对比: 与量子点输运实验进行详细比较

深度评价

优点

  1. 理论严谨性: 基于Hershfield的严格理论框架，数学推导清晰
  2. 方法创新: 巧妙地将有限格点本征态转换为散射态，避免了显式构造散射态的复杂性
  3. 数值可靠: 通过非相互作用极限验证，确保了方法的正确性
  4. 物理洞察: 提供了对Kondo共振非平衡演化的深入理解

不足

  1. 计算复杂度: 相比平衡态NRG，计算成本显著增加
  2. 参数调节: η\eta参数的选择需要经验，可能影响结果的普适性
  3. 实验验证: 缺乏与实际量子点实验的直接比较

影响力

  1. 方法学贡献: 为非平衡强关联系统提供了新的数值工具
  2. 理论价值: 深化了对Kondo效应非平衡行为的理解
  3. 应用前景: 可扩展到更复杂的杂质模型和格点系统

适用场景

  • 量子点输运实验的理论分析
  • 强关联电子系统的非平衡研究
  • 动力学平均场理论的杂质求解器
  • 分子结电子学的理论建模

参考文献

论文引用了该领域的重要文献，包括：

  • Kondo原始工作和Hewson的经典教科书
  • Hershfield的非平衡稳态理论
  • Wilson的数值重整化群方法
  • Anders的稳态NRG方法
  • 各种非平衡量子杂质理论和数值方法

这篇论文在非平衡量子多体理论领域做出了重要贡献，提供了一个既严格又实用的数值方法来研究强关联非平衡系统，特别是对理解Kondo效应的非平衡行为具有重要价值。