Fermi-Hubbard模型与量子模拟中的自旋电荷分离-开发者社区

1. Fermi-Hubbard模型基础与量子模拟价值

Fermi-Hubbard模型作为描述强关联电子系统的标准模型，其哈密顿量可表示为：

$$H = -J\sum_{\langle i,j\rangle,\sigma}(e^{i\phi_{ij}}c_{i,\sigma}^\dagger c_{j,\sigma} + h.c.) + U\sum_i n_{i,\uparrow}n_{i,\downarrow}$$

其中第一项描述电子在相邻格点间的跃迁（振幅J），第二项表示同一格点上自旋相反的电子间的库仑排斥（强度U）。这个看似简单的模型却蕴含着丰富的物理现象：

核心物理特征：

在强相互作用极限（U ≫ J）下，系统表现出Mott绝缘体行为
引入磁通量ϕ可调控能带结构，如π通量时形成狄拉克半金属态
电子分数化为自旋子和空穴子（holon）等准粒子激发

关键提示：在量子模拟中，通常将晶格格点映射为量子比特，电子占据状态编码为量子态。例如，|0⟩表示空穴，|↑⟩和|↓⟩分别对应两种自旋的电子占据，|d⟩表示双占据态。

2. 自旋电荷分离现象的物理机制

2.1 Kotliar-Ruckenstein表示与规范理论

通过Kotliar-Ruckenstein变换，可将电子算符表示为：

$$c_{i\sigma}^\dagger = f_{i\sigma}^\dagger h_i + \sigma f_{i\bar{\sigma}}d_i^\dagger$$

其中：

$f_{i\sigma}$为自旋子算符（费米型）
$h_i$为空穴子算符（玻色型）
$d_i$为双占据子算符（玻色型）

这种表示引入了U(1)规范对称性，对应的Wilson线算符：

$$V_{hd}(M) = \sum_{\substack{i,j\|i-j|=M}} \frac{h_i d_j}{N_{\text{pairs}}} \sum_{\gamma:i→j} \prod_{\langle mn\rangle\in\gamma} S_m^z S_n^z$$

可有效表征自旋与电荷的分离程度。

2.2 不同参数区域的分离特征

参数区域	自旋电荷行为	典型特征
U=0, ϕ=0	自由电子气	无分离
U>0, ϕ=0	弱约束分离	短程关联
U>0, ϕ=π	强约束分离	长程Wilson线

实验数据显示（图18-21），π通量下即使U=0时Vhd(M)也保持非零值，表明磁通增强了自旋与电荷的约束。

3. 量子计算实现方案

3.1 格点-量子比特映射

采用Jordan-Wigner变换将费米子算符映射为泡利算符：

$$c_j^\dagger = \left(\prod_{k<j} Z_k\right) \sigma_j^+$$

对于二维系统，通常使用Bravyi-Kitaev变换以减少非局域项。

3.2 Trotter分解实现

时间演化算符分解为：

$$e^{-iHt} \approx \left(\prod_{\langle ij\rangle} e^{iJ\tau K_{ij}} \prod_k e^{-iU\tau D_k}\right)^N$$

其中：

$K_{ij}$为跃迁项
$D_k$为相互作用项
τ=t/N为时间步长

关键参数选择：

单步误差~O(τ²)
总门数~O(L²N)，L为晶格尺寸
最优步长需平衡误差与噪声

3.3 误差缓解技术组合

实验采用的多层误差缓解方案：

TFLO（Training on Fermionic Linear Optics）：
- 利用高斯态可经典高效模拟的特性
- 通过对比训练校正非高斯门误差
MESR（Maximum Entropy Shot-Reweighting）：
- 构建最大熵分布： $$P(s) = \frac{1}{Z}\exp\left(\sum_i \lambda_i O_i(s)\right)$$
- 保持测量结果的对称性约束
对称性平均：
- 利用粒子数守恒等对称性
- 对破坏对称性的测量结果进行过滤

4. 关键实验结果分析

4.1 单空穴动力学（5×5晶格）

图8-11展示了中心空穴在Néel背景下的演化：

自旋传播速度：
- U=0时最快，随U增大而减慢
- π通量下出现振荡复苏现象
空穴扩散： $$ \text{holon}_{\text{RMS}} = \sqrt{\sum_i |i|^2 \langle h_i \rangle} $$ 早期呈线性增长（弹道传播），后期受相互作用抑制

4.2 空穴条纹演化

图12-13显示条纹态的特殊行为：

渗透分数（图15）： $$ P_{\text{perc}} = \frac{\text{跨越晶格的连通空穴路径数}}{\text{总样本数}} $$
- U增大时渗透率降低
- 反映强关联下空穴运动的受限性
反铁磁关联： $$ \text{AFM} = \frac{1}{N_{\text{shots}}} \sum_s \frac{\sum_{\langle ij\rangle} S_i^z(s)S_j^z(s)}{N_{\text{pairs}}(s)} $$ 随U增大衰减变缓

5. 技术挑战与解决方案

5.1 门分解优化

采用native门集实现：

# 示例：FSIM门分解 def fsim_gate(theta, phi): return [ (RZ(-phi/2), [q1]), (RZ(-phi/2), [q2]), (RX(theta/2), [q1]), (RX(-theta/2), [q2]), (CZ, [q1, q2]), (RX(theta/2), [q1]), (RX(-theta/2), [q2]), (RZ(phi/2), [q1]), (RZ(phi/2), [q2]) ]