Domain-wall melting and entanglement in free-fermion chains with a band structure

基本信息

• 论文ID: 2507.03780
• 标题: Domain-wall melting and entanglement in free-fermion chains with a band structure
• 作者: Viktor Eisler (Graz University of Technology & University of Graz)
• 分类: cond-mat.stat-mech, cond-mat.quant-gas, quant-ph
• 发表时间: September 18, 2025 (SciPost Physics Submission)
• 论文链接: https://arxiv.org/abs/2507.03780

摘要

本文研究了具有能带结构的自由费米子链中畴壁融化现象，其中跳跃振幅的周期性变化产生能带结构。研究表明纠缠熵随时间对数增长，且前因子与初始态中填充能带数成正比。对于二聚化链，粒子密度和电流具有与均匀情况相同的表达式，仅需对速度进行重新标度。熵分布的普适贡献被加倍，而非普适部分可通过块Toeplitz矩阵数值提取。

研究背景与动机

问题背景

1. 畴壁融化问题：畴壁融化是研究一维量子链远离平衡态动力学的经典设置，通过制备阶跃状初始密度分布并观察其时间演化来理解输运性质。
2. 现有研究局限：
   • 大多数研究集中在平移不变的均匀链上
   • 对于时间演化算符本身破坏平移不变性的情况了解较少
   • 缺乏对具有能带结构系统中纠缠传播的系统理解
3. 研究动机：
   • 探索周期性跳跃振幅如何影响畴壁融化动力学
   • 理解多能带系统中纠缠熵的增长规律
   • 建立具有能带结构系统与均匀系统之间的联系

核心贡献

1. 理论框架建立：首次系统研究了具有p-格点平移不变性的自由费米子链中畴壁融化问题
2. 纠缠增长规律：证明了纠缠熵保持对数增长，但前因子与填充准粒子能带数p成正比
3. 流体力学极限：对于二聚化链(p=2)，发现密度和电流分布可通过简单的速度重标化从均匀情况获得
4. 普适性与非普适性分离：
   • 普适贡献可解析获得且相比均匀情况加倍
   • 非普适贡献需通过块Toeplitz矩阵数值计算
5. 振荡现象分析：发现并解释了熵增长中的振荡现象，频率由能隙大小决定

方法详解

任务定义

研究具有周期性跳跃振幅的自由费米子链中畴壁初态的时间演化，分析粒子密度、电流和纠缠熵的动力学行为。

模型架构

二聚化跳跃链哈密顿量

Ĥ = -∑_{m=-N/2+1}^{N/2} [(1+δ)/2 c†_{2m-1}c_{2m} + (1-δ)/2 c†_{2m}c_{2m+1} + h.c.]

对角化过程

1. 子晶格变换：引入两个子晶格的费米子算符 a_m = c_{2m-1}, b_m = c_{2m}
2. 傅里叶变换：

   (a_m, b_m)^T = (1/√N) ∑_k e^{ikm} (a_k, b_k)^T
   

3. 能带对角化：通过酉变换得到对角形式

   Ĥ = -∑_k ω_k(α†_k α_k - β†_k β_k)
   

   
   其中色散关系为：

   ω_k = √[cos²(k/2) + δ²sin²(k/2)]
   

关联函数计算

流体力学极限

在 m,n,t → ∞ 且 m/t, n/t 固定的极限下，通过驻相近似得到：

C_{m,n}(t) = ∫_{q_-}^{q_+} dq/(2π) e^{iq(2n-2m)} [(1, e^{iφ_q}); (e^{-iφ_q}, 1)]
           + ∫_{-q_+}^{-q_-} dq/(2π) e^{iq(2n-2m)} [(1, -e^{iφ_q}); (-e^{-iφ_q}, 1)]

其中费米点 q_±(m,n,t) 由 v_{q_±} = (m+n-1)/t 确定。

技术创新点

1. 块矩阵表示：利用二格点平移不变性，采用2×2块矩阵描述关联函数
2. 准粒子图像：两个积分项对应上下能带准粒子的独立贡献
3. CFT-GHD结合：将共形场论方法与广义流体力学相结合处理纠缠传播

实验设置

数值计算参数

• 链长：N = 300 (开边界条件，避免前沿区域加倍)
• 时间演化：t = 200, 300
• 二聚化参数：δ = 0.02, 0.05, 0.1, 0.2, 0.3, 0.5

初始态设置

畴壁初态：左半链完全填充，右半链完全空：

C(0) = [1 0; 0 0]

评价方法

1. 密度和电流：与解析流体力学结果对比
2. 纠缠熵：通过关联矩阵本征值计算
3. 标度行为：验证对数增长规律

实验结果

主要结果

1. 密度和电流分布

• 密度：ρ(n,t) = (1/π)arccos(n/(v_max t))，其中 v_max = 1-δ
• 电流：J(n,t) = (v_max/π)√1-(n/(v_max t))²
• 验证：数值结果与解析表达式完美吻合

2. 纠缠熵增长

• 半链熵：S ~ (1/3)ln t + const（相比均匀情况的1/6 ln t加倍）
• 普适贡献：S_u = (1/3)ln2v_max t(1-ζ²)
• 非普适贡献：需通过块Toeplitz矩阵数值计算

3. 振荡现象

发现熵增长中存在衰减振荡：

S - (1/3)ln t ≈ α + β cos(2δt + φ)/(δt)

振荡频率恰好等于能隙大小2δ。

消融实验

δ → 1 极限

在完全二聚化极限下，解析得到：

lim_{δ→1} S_0(ζ) = (1/3)ln√(1-ζ²) + s_0

三能带模型

对于三格点周期性，验证了：

• 半链熵：S = (1/2)ln t + const
• 熵分布出现多个标度区域，反映不同能带的速度差异

案例分析

时间标度行为

• t ≪ δ^{-1}：表现与均匀链相同
• t ≫ δ^{-1}：过渡到加倍的对数增长
• 交叉时间尺度由能隙倒数决定

相关工作

畴壁融化研究

1. 均匀链：Antal等人首次研究XX链中畴壁融化10
2. 相互作用系统：通过GHD方法扩展到相互作用情况11
3. 纠缠研究：CFT方法理解纠缠传播21,22

非均匀系统

1. 局域缺陷：单个杂质导致线性熵增长38-44
2. 周期性调制：本文首次系统研究周期性跳跃振幅的影响

结论与讨论

主要结论

1. 普适标度：纠缠熵保持对数增长，前因子与能带数成正比
2. 速度重标化：流体力学量可通过简单速度重标化获得
3. 能隙效应：产生特征振荡，频率由能隙决定

局限性

1. 非普适部分：需要数值计算，缺乏解析表达式
2. 边缘结构：超越驻相近似的精细结构分析困难
3. 相互作用扩展：二聚化破坏可积性，GHD方法不适用

未来方向

1. 任意填充：研究非满填充初态的情况
2. 边缘物理：分析t^{1/3}标度区域的普适行为
3. 相互作用系统：探索是否存在有效流体力学描述

深度评价

优点

1. 理论完整性：从哈密顿量对角化到纠缠计算的完整理论框架
2. 方法创新性：巧妙结合CFT、GHD和矩阵分析方法
3. 结果普适性：建立了多能带系统的一般性规律
4. 数值验证：解析结果得到精确数值验证

不足

1. 技术局限：块Toeplitz矩阵常数项缺乏解析结果
2. 物理解释：振荡现象的深层物理机制需进一步阐释
3. 实验可及性：理论结果的实验验证方案未讨论

影响力

1. 理论贡献：为非均匀量子系统动力学提供新视角
2. 方法价值：块Toeplitz方法可推广到其他系统
3. 应用前景：为量子模拟实验提供理论指导

适用场景

1. 冷原子系统：光晶格中可实现周期性跳跃调制
2. 固体物理：二聚化材料中的激子动力学
3. 量子信息：多能带系统中纠缠传播的理解

参考文献

本文引用了76篇相关文献，涵盖：

• 广义流体力学理论3-8
• 畴壁融化经典工作9-17
• 纠缠熵理论18-27
• 共形场论方法47-49
• Toeplitz矩阵理论53-55,72-76