[论文解读] Critical temperature for fermion pairing using lattice field theory

Matthew Wingate|arXiv (Cornell University)|Feb 15, 2005

Cold Atom Physics and Bose-Einstein Condensates参考文献 1被引用 5

一句话总结

本论文首次采用蒙特卡罗模拟对强相互作用两组分费米气体在有限温度下的超流相变进行了格点场论研究。该研究在普遍性区域确定了临界温度 $T_c$，在单位率下得到 $T_c/T_F \approx 0.035 \pm 0.004$，证明了在无符号问题和无基态假设条件下，第一性原理计算在强耦合极限下的可行性。

ABSTRACT

Dilute gases of 2-component fermions are of great interest in atomic and nuclear physics. When interactions are strong enough so that a bound state is at threshold, universal behavior is expected. Lattice field theory provides a first principles approach to the study of strongly interacting systems such as this through Monte Carlo simulation. Results of exploratory simulations are presented here. In particular, the finite temperature phase transition between superfluid and normal states is studied. We present first results for the critical temperature T_c and describe the future work necessary to determine T_c as a function of interaction strength.

研究动机与目标

• 开发并应用格点场论结合蒙特卡罗模拟，研究强相互作用两组分费米气体在有限温度下的超流相变。
• 在 $1/a \to 0$ 的普遍性区域计算临界温度 $T_c$，避免对基态或波函数的假设。
• 证明在大正则系综中使用辅助场方法进行有限温度研究时，可避免符号问题，具有可行性。
• 在探索性模拟中，系统地外推至 $J=0$，并校正有限体积和格点间距的影响。
• 绘制临界温度 $T_c$ 随相互作用强度的变化曲线，尤其关注BCS与BEC区域之间的交叉区域。

提出的方法

• 为具有接触相互作用的两组分费米子构建哈密顿量，并通过 Hubbard-Stratonovich 变换将其映射为路径积分形式，利用辅助标量场 $\phi$ 分解四费米子相互作用。
• 采用大正则系综，费米子在虚时方向采用反对称边界条件，从而可直接访问有限温度热力学性质。
• 通过混合蒙特卡罗（HMC）方法进行随机模拟，采用多个 $J$ 值（0.07、0.1、0.14）外推至 $J=0$，确保无费米子符号问题。
• 将费米子凝聚能 $\Sigma$ 作为超流性的序参量，监测其在相变点附近的急剧下降。
• 利用 $T_c/4$ 温度下的有限温度密度，通过关系式 $T_c = \xi \cdot T_F / (n_0 b_s^3)$ 计算 $T_c/T_F$，其中 $n_0$ 由 $N_t=48$ 的模拟结果确定。
• 采用数据分箱和误差估计方法以考虑自相关性，$N_{\text{bin}}=5$ 足够控制统计误差。

实验结果

研究问题

• RQ1 格点场论结合HMC模拟能否准确确定强相互作用费米气体在有限温度下的超流-正常相变？
• RQ2 在普遍性区域（$1/a \to 0$）下，采用无基态假设的第一性原理方法，$T_c/T_F$ 的值是多少？
• RQ3 有限体积和格点间距效应对提取的 $T_c$ 有何影响？能否系统地控制这些效应？
• RQ4 在相变路径上，散射长度 $a$ 如何变化？需要何种参数路径才能在改变化学势 $\mu$ 的同时保持 $a$ 恒定？
• RQ5 在大正则系综中使用辅助场方法能否在强耦合极限下可靠地计算 $T_c$，即使在无符号问题的情况下？

主要发现

• 费米子凝聚能 $\Sigma$ 在相变点附近表现出急剧下降，表明在参数空间的四条路径上均存在清晰的超流到正常相的转变。
• 临界温度估计为 $T_c/T_F = 0.035 \pm 0.004$，对应 $m^2 = 0.155$ 和 $0.175$，即接近单位率的普遍性区域。
• 在相变区域，散射长度 $a$ 显著变化，$m^2=0.155$ 时 $b_s/a$ 从 $-0.10$ 变化至 $0.10$，$m^2=0.175$ 时从 $-0.11$ 变化至 $0.08$，因此需谨慎选择路径以保持 $a$ 恒定。
• 零温密度 $n_0 b_s^3 = 0.2612(6)$ 和 $0.2893(7)$ 在 $N_t=48$ 条件下计算得出，$T_c/4$ 温度下的有限温度密度对 $N_t$ 或 $J$ 的依赖性极小。
• 有限体积和格点间距误差远大于统计误差，表明在交叉区域实现高精度需更大规模模拟。
• 该方法成功避免了费米子符号问题，且无需对基态做任何假设，验证了其在第一性原理研究强耦合超流性方面的适用性。