用纠缠见证探测杂质量子临界性

Probing Impurity Quantum Criticality with Entanglement Witnesses

arXiv: 2608.12317v1

论文信息

标题: Probing Impurity Quantum Criticality with Entanglement Witnesses

作者: Mateo Cárdenes Wuttig, Andrew J. Millis

发布日期: 2026-08-12

arXiv ID: 2608.12317v1

PDF 链接: 下载 PDF

3 分钟速览

  • 研究问题:这篇论文研究双杂质 Kondo 模型中量子临界性如何通过纠缠重排体现,以及可测量的自旋/电荷涨落能否作为 “纠缠见证” 来探测这种临界性。
  • 核心方法:用密度矩阵重整化群(DMRG)计算零温基态波函数、纠缠熵和量子 Fisher 信息(QFI);用数值重整化群(NRG)计算有限温动态响应函数并按公式重构 QFI。
  • 关键结果:在固定 Kondo 耦合 J=2J=2 时,DMRG 外推得到临界杂质间交换 Kc=1.620±0.005K_c=1.620\pm0.005;QFI 在临界区出现峰值,并且杂质自旋通道的 QFI 能在温度高至约 T≈0.5T\approx0.5 时仍见证纠缠。
  • 主要局限:论文自身指出,临界点对电荷转移、粒子-空穴对称破缺、不等杂质耦合和耗散等因素不稳定;QFI 若要认证纠缠需要自旋分辨测量;有限温下认证温度范围受低能标限制。
  • 适合读者:对量子杂质、非费米液体临界性、量子信息启发的凝聚态探针以及量子点/STM 实验设计感兴趣的研究者。

论文背景和研究动机

量子临界性与多体纠缠重排之间的关系是现代凝聚态物理中的一个开放问题。传统方法通常用序参量、临界指数或关联长度来刻画相变,但这些量并不直接回答 “谁与谁纠缠” 发生了什么变化。双杂质 Kondo 模型提供了一个理想的可控测试平台:两个自旋 S=1/2S=1/2 杂质分别耦合到各自的非相互作用电子库,并通过杂质间交换 KK 相互耦合。该模型实验上可通过耦合量子点、STM 表面磁性原子或金属表面有机分子等手段实现。

在反铁磁耦合 J,K>0J,K>0 且无电荷转移的情况下,模型存在一条非费米液体量子临界线 Kc(J)K_c(J),将两种基态分开:一种是由单个电子库独立 Kondo 屏蔽杂质的 “双 Kondo 单态” 相(2KS),另一种是由两个杂质相互锁定成近纯单态的 “杂质间单态” 相(IIS)。这两种相的关键差别在于每个局域磁矩的主导纠缠伙伴:在 2KS 相,每个杂质与扩展的传导电子自由度纠缠;在 IIS 相,杂质主要与另一个杂质纠缠。作者希望把这种纠缠伙伴的变化用量子信息语言表达,并通过可测量的自旋和电荷涨落来 “见证” 它。

这一研究的动机还来自更宏观的问题:重费米子物理中的 Kondo 破坏临界性、两通道 Kondo 模型中的残余熵,以及最近在 Kondo 晶格和非弹性中子散射实验中观察到的 QFI 增强。双杂质模型提供了一个微观可控的基准,可用于理解这些更复杂系统中的量子临界纠缠结构。

核心方法和技术细节

论文的哈密顿量如下:

H=He+Himp+HKondoH = H_e + H_{\rm imp} + H_{\rm Kondo}

其中 HeH_e 描述两条半填充非相互作用电子链,Himp=K S1⋅S2H_{\rm imp}=K\,\mathbf{S}_1\cdot\mathbf{S}_2 是杂质间 Heisenberg 交换,HKondo=J∑αSα⋅se,αH_{\rm Kondo}=J\sum_\alpha \mathbf{S}_\alpha\cdot\mathbf{s}_{e,\alpha} 是每个杂质与其最近邻电子位置的 Kondo 交换。电子链长度为 LL,每条链末端耦合一个杂质。

纠缠度量 使用 von Neumann 熵。论文定义了两个关键分区:S(AL)S(A_L) 表示左右两条单杂质子链之间的纠缠,S(Aimp)S(A_{\rm imp}) 表示两个杂质作为一个子系统相对全部电子的纠缠。在 SU(2) 不变的单态扇区中,S(Aimp)S(A_{\rm imp}) 可以直接由杂质-杂质关联函数 CII=⟨S1⋅S2⟩C_{\rm II}=\langle \mathbf{S}_1\cdot\mathbf{S}_2\rangle 表示(见附录 A.3.1)。这建立了传统关联函数与纠缠熵之间的解析联系。

量子 Fisher 信息(QFI) 是核心工具。对于热态 ρT\rho_T 和算符 O\mathcal{O},QFI 定义为

FQ[ρT,O]=4π∫0∞dω tanh⁡(ω2T)χOO′′(ω,T)F_Q[\rho_T,\mathcal{O}]=\frac{4}{\pi}\int_0^\infty d\omega\,\tanh\left(\frac{\omega}{2T}\right)\chi''_{\mathcal{O}\mathcal{O}}(\omega,T)

其中 χOO′′\chi''_{\mathcal{O}\mathcal{O}} 是动态磁化率的耗散部分。纯态情形简化为 FQ=4 Var(O)F_Q = 4\,\mathrm{Var}(\mathcal{O})。对左右对称/反对称组合的集体算符,可分态的 QFI 上界为 22,超过 22 即见证左右纠缠。作者重点使用杂质自旋差 SIIz,−=S1z−S2zS^{z,-}_{\rm II}=S_1^z-S_2^z、杂质自旋和 SIIz,+S^{z,+}_{\rm II}、最近邻电子自旋差 seez,−s^{z,-}_{ee} 以及最近邻电子密度差 nee−n^{-}_{ee} 等生成元。

数值方法 包括 DMRG 和 NRG。DMRG 用 ITensor 以 MPS 形式获得零温基态,最大键维数 χmax⁡=4000\chi_{\max}=4000,链长 LL 从 50 到 500。NRG 使用 TRIQS/NRG Ljubljana 实现,保留最多 Nkeep=1500N_{\rm keep}=1500 个多体态,对数离散化 Λ=2\Lambda=2。有限温 QFI 通过 NRG 计算动态响应并按上述积分重构。

创新点和贡献

论文的第一项创新是把纠缠见证从抽象的量子信息概念变为实验中可测量的动态响应量。此前的工作常把 QFI 与谱函数联系,但本文明确用于双杂质 Kondo 模型的临界区,并给出了杂质自旋和邻近电子密度两个独立通道的可操作性方案。对纯 SU(2) 单态,杂质自旋 QFI 与 CIIC_{\rm II} 之间有简单关系 FQ(SIIz,±)=2±83CIIF_Q(S^{z,\pm}_{\rm II})=2\pm\frac{8}{3}C_{\rm II},因此 “QFI 之差” 即可直接读出杂质间自旋关联。

第二项创新在于发现邻近电子密度差 nee−n^{-}_{ee} 的 QFI 也能追踪临界点,且不需要直接测量杂质间自旋关联。图 3(d) 显示 ∂KFQ(nee−)\partial_K F_Q(n^{-}_{ee}) 在有限尺寸临界耦合 Kc(L)K_c(L) 附近出现明显峰。这说明杂质临界性的纠缠重排会传播到周围传导电子的局域密度涨落中,为实验提供了更灵活的探测路径。自旋分辨的 FQ(n↑,ee−)F_Q(n^{-}_{\uparrow,\mathrm{ee}}) 还能同时见证纠缠,因为它超过了可分性上界 22(图 3(d) 插图)。

第三项贡献是把零温 QFI 分析扩展到有限温度。NRG 结果(图 4)表明,临界响应在低温下表现为尖锐峰,随温度升高变宽且被抑制;杂质自旋 QFI 在 T≈0.5T\approx0.5 时仍能超过可分性上界,而相邻电子自旋分辨 QFI 的认证温度约在 T≈0.1T\approx0.1。论文指出这些温度量级与耦合常数相当,对实验探测是有利的。

最后,论文补充了一个可精确求解的四自旋强耦合极限模型,它能在局域层面复现纠缠结构,并解释有限温下 CIIC_{\rm II} 和 CIEC_{\rm IE} 的行为。这为理解更复杂的数值结果提供了直观基础。

实验结果分析

零温相图与纠缠重排。图 2 展示了 CIIC_{\rm II}、CIEC_{\rm IE}、两种纠缠熵和两种 QFI 在 JJ-KK 平面上的分布。在 J,K>0J,K>0 区域,CIIC_{\rm II} 从 2KS 相约 0 变到 IIS 相约 −3/4-3/4;CIEC_{\rm IE} 从约 −3/4-3/4 的 AFM 屏蔽变到约 0。对应的 S(Aimp)S(A_{\rm imp}) 在 2KS 区接近最大值 2ln⁡22\ln2,在 IIS 区趋近 0;S(AL)S(A_L) 则呈现相反趋势。QFI 图 2(c) 显示 FQ(SIIz,−)F_Q(S^{z,-}_{\rm II}) 在单态关联区域超过可分性边界 22,直接认证左右纠缠。这些结果说明,当经过临界线时,杂质的主要纠缠伙伴从 “各自电子库” 切换为 “另一个杂质”,并且这一切换可由 QFI 捕获。

临界点附近的奇异响应。图 3 固定 J=2J=2,改变 KK 并比较不同链长 LL。∂KCIE\partial_K C_{\rm IE}、−∂KS(Aimp)-\partial_K S(A_{\rm imp}) 和 ∂KFQ(SIIz,−)\partial_K F_Q(S^{z,-}_{\rm II}) 都在尺寸依赖的 Kc(L)K_c(L) 附近形成峰,且峰高随 LL 增大而增强。论文用有限尺寸外推 KcK_c 到热力学极限,得到 Kc=1.620±0.005K_c=1.620\pm0.005(DMRG),与 NRG 在 T=0.001T=0.001 时的 Kc=1.59±0.03K_c=1.59\pm0.03 一致(考虑不同红外截断后)。这验证了 QFI 的导数可以作为临界点的实用刻度。

有限温实验窗口。图 4 的 NRG 数据显示,静态磁化率 χ′(SIIz,−)\chi'(S^{z,-}_{\rm II}) 在低温 T=0.001T=0.001 时出现尖锐峰,随着温度升高而变宽。FQ[ρT,SIIz,−]F_Q[\rho_T,S^{z,-}_{\rm II}] 在 KK 的一定范围内仍高于可分界 22,直到 T≈0.5T\approx0.5。图 4(d) 表明邻近电子密度 QFI 的导数在同一温度窗口内仍可见临界结构。作者指出,T∗≪T≪TKT^*\ll T\ll T_K 是非费米液体临界的生存窗口,超出此窗口后 QFI 结构被热涨落抹平。

局限与待解决问题

本文的理论架构在真实实验面前仍有一些需要处理的不稳定性。作者明确指出,双杂质 Kondo 临界点对粒子-空穴对称破缺、两电子库之间的电荷转移、杂质耦合不对称以及耗散都敏感;一旦这些扰动存在,尖锐临界点会被圆滑成交叉行为。因此,实验上必须仔细评估耦合量子点或表面磁性原子系统中这些扰动的幅度。

有限温 QFI 的认证窗口也受限。虽然杂质自旋通道可以见证到 T≈0.5T\approx0.5,但相邻电子自旋分辨 QFI 的认证温度只有约 T≈0.1T\approx0.1,这意味着如果实验分辨率只能看电子密度通道,认证纠缠的温度范围会显著缩小。论文没有给出这些温度在具体器件中的物理温度标度,只以 tt 为单位展示,因此需要结合具体材料的 tt 才能判断实际可行性。

另一个开放问题是空间分辨的纠缠结构。本文使用 S(AL)S(A_L) 和 S(Aimp)S(A_{\rm imp}) 这类全局分区熵,无法回答 “Kondo 屏蔽云内部不同距离的电子如何与杂质分享纠缠”。作者推测,结合总自旋算符和动量分辨的自旋/电荷 QFI,可能把这种杂质尺度上的重排连接到重费米子准粒子的形成与消失,但论文没有在双杂质模型中给出相应的动量分辨数据。此外,与两通道 Kondo 模型的残余熵 12ln⁡2\frac12\ln2 的关系也值得进一步研究。

最后,DMRG 和 NRG 本身存在有限尺寸/离散化误差。论文通过改变键维数、链长、Λ\Lambda 和 NkeepN_{\rm keep} 做了收敛检查,但导数峰的绝对高度仍有几个百分点的系统变化。因此,文中的临界峰值行为应理解为在该数值设置下的收敛趋势,而不是解析证明的标度形式。