Spin-imbalanced fermion on a dynamic lattice

技术摘要：动态晶格上的自旋不平衡费米子

问题陈述
本研究探讨了一种一维（1D）自旋-1/2 费米子动态晶格模型的基态磁性质。在传统的平衡系统（如扩展费米-哈伯德模型）中，已确立自旋密度波（SDW）序矢量的关系遵循 $k = 2k_F$；然而，在动态晶格中，这些系统在费米子自旋不平衡（$N_\uparrow \neq N_\downarrow$）情况下的行为在很大程度上仍未得到充分探索。具体而言，本文研究了通过自旋相关跳跃实现的巡游费米子与键中心定域自旋之间的相互作用，在费米子部分极化时，如何影响磁调制。

方法论
作者采用密度矩阵重整化群（DMRG）方法来模拟系统的基态。该模型哈密顿量描述了在跳跃振幅为 $t$ 的晶格位点间巡游费米子的跳跃，该跳跃通过耦合项 $\alpha$ 被位于键上的定域自旋动态调制。定域自旋受到纵向场 $\Delta$ 和横向场 $\beta$ 的作用。巡游费米子经历外部塞曼场 $h$（用于控制自旋不平衡）和原位排斥性哈伯德相互作用 $U$。

模拟聚焦于准绝热机制（$\beta \ll t$），设定 $\beta = 0.02$ 且 $t=1$。分析了尺寸范围为 $L=24$ 至 $L=60$、采用开边界条件的系统。关键观测物理量包括定域自旋结构因子 $S_z(k)$、实空间相关函数 $S_z(r)$ 以及自旋分辨的动量分布函数 $n_\sigma(k)$。利用有限尺寸缩放法将结果外推至热力学极限（$L \to \infty$）。

主要贡献与结果

1. 磁相识别：
   研究绘制了相图，识别出两个主要的物相：顺磁（PM）相和自旋密度波（SDW）相。这两相之间的转变是由自旋相关跳跃耦合（$\alpha$）与纵向场（$\Delta$）之间的竞争驱动的。当 $\alpha$ 与 $\Delta$ 量级相当时，会出现空间调制的 SDW 序。当其中一个参数占主导地位时，系统进入以均匀自旋排列和缺乏长程磁相关为特征的 PM 相。

2. 双重 SDW 序现象（SDW-I 和 SDW-II）：
   在部分自旋极化机制下，作者发现了两种截然不同的 SDW 序模式，这些模式从根本上偏离了平衡系统中标准的 $k=2k_F$ 预期：

• SDW-I： 序矢量 $k_{max}$ 表现出对多数自旋密度（$\rho_\uparrow$）的线性依赖关系，遵循 $k_{max} \approx 2\pi\rho_\uparrow$。这对应于自旋向上费米面的嵌套矢量。
• SDW-II： 相反，$k_{max}$ 表现出对少数自旋密度（$\rho_\downarrow$）的线性依赖关系，遵循 $k_{max} \approx 2\pi\rho_\downarrow$。这对应于自旋向下费米面的嵌套矢量。
  这两种现象之间的转换可以通过调节耦合参数（$\alpha, \Delta$）和填充因子来控制。

3. 微观起源：
   对自旋分辨动量分布函数的分析表明，主导性的 SDW 峰值并非由总填充密度决定，而是由自旋组分的相对占据数决定。这种双重模式的出现被解释为一种类似于皮埃尔斯（Peierls-like）转变的磁性模拟，由巡游费米子特定自旋扇区的费米面嵌套所驱动。

4. 对哈伯德相互作用的鲁棒性：
   研究发现，引入排斥性原位哈伯德相互作用（$U$）通过抑制双占据，增强了趋向于费米子自旋极化的倾向（降低了实现完全极化的临界场）。然而，排斥相互作用并未从根本上改变 SDW-I 和 SDW-II 相的性质；这些独特的序模式在其各自的参数区域内保持稳健。

意义
本文声称揭示了一种在低维关联系统中调控磁调制的全新途径。通过证明费米子自旋不平衡可以驱动系统进入具有不同嵌套矢量（$2k_{F\uparrow}$ 对比 $2k_{F\downarrow}$）的不同 SDW 相，这项工作丰富了对动态晶格磁性的微观理解。结果表明，此类系统中的磁序不仅是一个静态属性，而且是与定域自旋的构型以及巡游费米子的特定自旋组分动态耦合的。这为理解超越传统平衡极限的低维系统中自旋不平衡如何塑造量子相提供了理论框架。