1. 从“一团乱麻”到“有序编织”：理解相互作用玻色气体的核心挑战
在凝聚态物理和量子多体物理的研究中，相互作用玻色气体是一个经典而又充满活力的模型。它听起来可能很学术，但我们可以把它想象成一个微观世界里的“人群”。想象一下，在一个极低温的房间里，有一群完全相同的、行为非常“佛系”的粒子（玻色子），它们不仅不排斥彼此，反而倾向于聚集在同一个状态，这就是著名的玻色-爱因斯坦凝聚现象。然而，现实世界没有绝对的“佛系”，粒子之间总会有微弱的相互作用，就像人群里总有人会不小心踩到别人的脚，或者低声交谈。这些微弱的相互作用，就是“相互作用玻色气体”这个名字的由来。
那么，物理学家们关心什么呢？他们关心这个系统的“自由能”。自由能是热力学里的一个核心概念，你可以把它理解为系统在特定温度和体积下，可以用来“做有用功”的那部分能量，或者说，是系统稳定性的一个度量。计算自由能，就等于掌握了这个量子多体系统的“命脉”——我们能知道它在什么条件下会形成超流，什么时候会发生相变，它的激发谱（可以理解为系统的“声音”或振动模式）长什么样。
但计算它谈何容易。当粒子数量N极大，且它们之间两两都有相互作用时，直接求解薛定谔方程是天方夜谭。传统的微扰理论在处理强相互作用或低温区域时常常失效。这时，我们就需要一些更聪明、更强大的理论工具。标题中提到的“环路展开”与“交织变分”，正是近年来处理这类强关联量子多体问题的两柄利剑。它们不是去直接硬算那天文数字般的项，而是试图用一种更结构化的方式来“描述”或“逼近”系统的真实行为。这篇内容，我就想结合自己的一些理解，聊聊这两种方法是如何切入，并共同为我们描绘相互作用玻色气体自由能这幅复杂图景的。
2. 理论基石：为什么要用自由能，以及传统方法的瓶颈
在深入“环路”与“交织”之前，我们必须夯实基础，明白我们战斗的战场——自由能——为何如此重要，以及传统的战士为何在此处常常折戟。
2.1 自由能：量子多体系统的“全景地图”
对于处于热平衡的系统，其一切宏观可观测性质都蕴含在配分函数Z之中。而亥姆霍兹自由能F与配分函数的关系是：F = -k_B T ln Z。其中k_B是玻尔兹曼常数，T是温度。一旦我们知道了F，通过对它求各种导数，就能得到所有热力学量：内能、熵、压强、粒子数密度、超流密度、响应函数等等。因此，计算F是统计物理的核心目标。
对于相互作用的玻色气体，其哈密顿量通常写作： [ \hat{H} = \int d^3r \left[ \hat{\psi}^\dagger(\mathbf{r}) \left( -\frac{\hbar^2}{2m}\nabla^2 \right) \hat{\psi}(\mathbf{r}) + \frac{g}{2} \hat{\psi}^\dagger(\mathbf{r}) \hat{\psi}^\dagger(\mathbf{r}) \hat{\psi}(\mathbf{r}) \hat{\psi}(\mathbf{r}) \right] ] 这里，(\hat{\psi})是玻色场算符，第一项是动能项，第二项是接触相互作用项，强度由耦合常数g描述。我们的任务就是计算这个复杂量子场论模型的配分函数 (Z = \text{Tr} e^{-\beta \hat{H}})。
2.2 传统微扰论的“无力感”：当相互作用不再微弱
最直接的想法是从无相互作用的理想玻色气体出发，把相互作用项当作微扰。这就是玻戈留波夫理论所做的，它在描述弱相互作用玻色凝聚体的低能激发（即声子谱）方面取得了巨大成功。然而，这种微扰展开本质上是按相互作用强度g的幂次进行的。当系统处于强相互作用区域（例如在光学晶格中，或利用Feshbach共振将散射长度调到很大时），微扰级数收敛极慢甚至发散。更重要的是，微扰论无法很好地处理非微扰效应，例如拓扑激发（涡旋）、量子涨落导致的凝聚分数降低等。
另一种思路是平均场理论，比如Gross-Pitaevskii方程，它将场算符替换为一个经典复函数（序参量）。这相当于忽略了量子涨落，在凝聚体密度很高、涨落很弱的条件下是很好的近似。但对于低维系统（如二维玻色气体）、强关联系统或相变点附近，量子涨落至关重要，平均场理论会给出完全错误的结论，例如它预言了错误的相变类型（平均场临界指数与真实指数不符）。
这就引出了我们的核心困境：我们需要一种既能超越微扰论、又能系统性地包含涨落效应的非微扰计算方法。路径积分表述为我们提供了这样一个框架。
3. 路径积分与相干态：为“环路展开”铺平道路
要理解“环路展开”，我们必须先进入路径积分的世界。这是一种将量子力学与统计力学统一起来的强大语言，特别适合处理多体问题。
3.1 从算符到场：相干态路径积分
对于玻色系统，最方便的路径积分是基于相干态的。相干态是湮灭算符的本征态，它最接近经典的“波场”概念。通过将配分函数Z用相干态展开，我们可以将其写成一个关于复值场 (\phi(\tau, \mathbf{r})) 的泛函积分： [ Z = \int \mathcal{D}[\phi^, \phi] \exp\left( -S[\phi^, \phi] \right) ] 其中，(S) 是欧拉作用量，(\tau) 是虚时间（范围从0到 (\beta=\frac{1}{k_B T})）。这个作用量通常包含三部分： [ S[\phi^, \phi] = \int_0^\beta d\tau \int d^3r \left[ \phi^\frac{\partial \phi}{\partial \tau} + \frac{\hbar^2}{2m} |\nabla \phi|^2 - \mu |\phi|^2 + \frac{g}{2} |\phi|^4 \right] ] 这里 (\mu) 是化学势。注意，此时的场 (\phi) 不再是量子算符，而是一个经典的复值场（但满足周期性边界条件 (\phi(0)=\phi(\beta))）。然而，这个“经典”场包含了所有量子涨落的信息，因为我们对所有可能的场构型进行了积分。问题在于，这个积分是无穷维的，无法精确计算。
3.2 稳态解与涨落：分离平均场与量子修正
一个标准的策略是寻找这个作用量的稳态解（ saddle point solution）。令 (\delta S / \delta \phi^* = 0)，我们恰好得到了 Gross-Pitaevskii 方程在静态均匀情况下的解：(\phi_0 = \sqrt{\rho_0})，其中 (\rho_0) 是平均场下的凝聚体密度。这个解对应于忽略所有涨落的平均场近似。
接下来，我们将场围绕这个稳态解展开：(\phi(\tau, \mathbf{r}) = \phi_0 + \eta(\tau, \mathbf{r}))，其中 (\eta) 代表涨落场。将展开式代入作用量S，并按涨落场 (\eta) 的幂次进行整理： [ S = S_0 + S_2[\eta^, \eta] + S_{\text{int}}[\eta^, \eta] ] 这里 (S_0) 是常数（平均场自由能），(S_2) 是涨落的二次型（高斯项），而 (S_{\text{int}}) 包含了三次及更高次的相互作用项。如果忽略 (S_{\text{int}})，只对高斯项进行泛函积分，这个计算是可以精确完成的，得到的结果就是“单圈”修正。这个修正包含了零点的量子涨落（即使在绝对零度也存在）和热涨落的影响。玻戈留波夫理论实际上就对应于这个高斯近似。
那么，“多圈”修正从哪里来？答案就是来自被我们忽略的 (S_{\text{int}})。要包含这些项的影响，就需要进行“环路展开”。
4. 解密“环路展开”：用费曼图语言系统化量子修正
“环路展开”是量子场论中微扰计算的标准技术，但在凝聚态多体问题中，它提供了一种超越简单高斯近似的系统化方案。
4.1 从相互作用项到费曼规则
(S_{\text{int}}) 项（例如 (\phi^4) 理论中的 (g|\phi|^4) 项）是涨落场之间的相互作用。在路径积分中，包含这一项的贡献可以通过将其展开成级数，并利用高斯积分的威克定理来计算。这个过程被一套图形化规则——费曼图——极大地简化了。
对于我们的玻色气体，基本的费曼规则包括：
- 传播子（线）：代表自由（高斯）涨落的关联函数 (G_0(k, i\omega_n))，其中 (k) 是动量，(\omega_n) 是玻色子松原频率。
- 相互作用顶点（点）：代表相互作用强度 (g)。在 (\phi^4) 理论中，一个顶点连接四条线。
- 圈（Loop）：代表一个动量和频率的积分。图中一个闭合的环就是一个“圈”。
自由能F的对数，即 (\ln Z)，可以表示为所有连通真空费曼图的和。这里“真空”是指没有外接腿的图。
4.2 “圈”的物理意义：量子涨落的关联与散射
一个“单圈图”，比如一个简单的泡泡图，它对应着两个粒子从相互作用顶点产生，传播一段时间后又湮灭。这个图计算了涨落场自能的最低阶修正，它修正了准粒子的色散关系和寿命。玻戈留波夫理论等价于只考虑这种自能图。
“双圈图”则更为复杂。例如，它可能描述了两个粒子通过相互作用散射后，中间态粒子又发生了一次相互作用再回到原态的过程。这包含了更高阶的关联效应。每一个新增的“圈”，都意味着考虑了更多次数的粒子-粒子散射过程，从而更精确地包含了多体关联。
注意：在凝聚态语境下，“环路展开”通常不是指按耦合常数g的幂次展开（那可能不收敛），而是按“圈数”展开。即使相互作用g很强，如果系统的维度高或涨落被压制，高阶圈图的贡献也可能很小。这是一种按拓扑结构（图的连通性）进行的非微扰分类方案。
4.3 实际操作：如何执行一圈和两圈计算
假设我们已经得到了高斯近似下的传播子 (G_0)。一圈修正对自由能的贡献，来自于将所有单圈拓扑结构的费曼图（如泡泡图、蝌蚪图）求和。在均匀系统中，这通常归结为计算一个动量-频率积分： [ F_{\text{1-loop}} \propto \frac{1}{2} \sum_{\mathbf{k}, \omega_n} \left[ \ln \det G_0^{-1}(k, i\omega_n) + \text{...} \right] ] 这个积分常常是发散的，需要重整化。物理的可观测量（如密度、压缩率）是有限的。
两圈修正涉及更复杂的积分，例如两个传播子嵌套的积分。计算量急剧增加，但能捕捉到诸如阻尼效应、 beyond-mean-field 的态方程修正（如李政道-黄-杨修正）等物理。在实际计算中，我们常常需要结合数值积分和解析延拓等技术。
5. 超越微扰：“交织变分法”的核心理念
“环路展开”虽然强大，但它本质上仍是一种微扰方法（按圈数微扰）。对于强关联系统，我们需要更彻底的非微扰方法。这就是“变分法”登场的时候。而“交织变分法”是其中一种特别适用于量子多体问题的现代变分思路。
5.1 变分法精髓：在试探函数族中寻找最优解
变分法的核心思想非常直观：我们不知道真实的量子态 (|\Psi\rangle) 或自由能 (F) 是什么，但我们猜测一个由一组变分参数 ({\lambda}) 描述的试探态 (|\Psi_{\text{trial}}({\lambda})\rangle) 或试探作用量 (S_{\text{trial}})。根据量子力学或统计物理的基本原理，真实的自由能总是小于或等于任何试探方案给出的自由能估计值（对于基态是能量最低原理，对于有限温度是吉布斯-博戈留波夫不等式）。
因此，我们计算试探方案下的自由能 (F_{\text{trial}}({\lambda}))，然后通过调整参数 ({\lambda}) 来最小化它： [ F \leq F_{\text{trial}}({\lambda}), \quad F_{\text{var}} = \min_{{\lambda}} F_{\text{trial}}({\lambda}) ] 这个最小值 (F_{\text{var}}) 就是我们能找到的对真实自由能的最佳变分估计。
5.2 何为“交织”？构建包含关联的试探态
传统的平均场理论就是一种简单的变分法，其试探态是单体乘积态（或在路径积分中是高斯作用量），它完全忽略了粒子之间的纠缠（关联）。而“交织”这个词，正是指通过引入一个“交织算符”来在试探态中“编织”进多体关联。
一个经典的例子是金兹堡-朗道理论。它将自由能泛函写成序参量及其梯度的多项式展开形式。这个泛函本身就是一个变分试探泛函，其形式由对称性决定。通过最小化这个泛函，我们可以得到非均匀的解（如涡旋、畴壁），这些解本身就包含了空间关联。
更现代的“交织变分法”可能采用矩阵乘积态、投影纠缠对态或高斯态加上多体局域酉变换等形式作为试探态。例如，对于玻色系统，我们可以从平均场波函数（凝聚态）出发，对其施加一个两体（或更多体）的纠缠算符： [ |\Psi_{\text{trial}}\rangle = \hat{J} |\text{MF}\rangle, \quad \hat{J} = \exp\left( \sum_{i,j} f_{ij} \hat{a}_i^\dagger \hat{a}j^\dagger - \text{h.c.} \right) ] 这里的算符 (\hat{J}) 就是“交织算符”，它能在平均场背景上产生粒子对的关联。变分参数就包含在函数 (f{ij}) 以及平均场参数中。通过优化这些参数，我们可以得到一个比单纯平均场或微扰论包含更多强关联效应的波函数。
在路径积分框架下，“交织”可以体现在选择一个非高斯型的试探作用量 (S_{\text{trial}})，它可能包含非局部的相互作用核，从而能够描述长程关联或临界涨落。
6. 双剑合璧：环路展开与交织变分在实践中的协同
在实际研究中，“环路展开”和“交织变分”并非互斥，而是常常协同使用，互为补充。理解它们的结合方式，是掌握现代多体理论方法的关键。
6.1 变分法提供起点，环路展开进行修正
一种常见的策略是：首先使用一个精心设计的变分试探态（例如包含Jastrow因子的波函数）来捕捉最主要的非微扰关联效应。这个变分态给出了一个“优化后的平均场”背景。然后，在这个优化后的背景上，考虑剩余的、较弱的涨落，对这些涨落进行环路展开（微扰计算）。
例如，在超冷原子气体中研究强相互作用的玻色-爱因斯坦凝聚体时，可以先使用一个包含两体关联函数的变分蒙特卡洛方法获得一个近似的基态波函数。这个波函数已经包含了由于强排斥作用导致的“空穴”（即两个粒子不能靠得太近）效应。然后，在这个基础上，我们可以计算低能激发的传播子，并利用一圈展开来研究这些激发的相互作用和衰减，这比在裸的平均场背景上做微扰要准确得多。
6.2 从变分原理推导有效场论
另一种深刻的联系体现在，许多有效的场论作用量本身可以通过变分原理推导出来。例如，** Hubbard-Stratonovich变换** 是一种常用的技巧，它通过引入辅助场来将四费米子相互作用项解耦。这个引入辅助场的过程，可以看作是一种变分变换。变换后得到的新作用量，其稳态解可能对应着一个新的有序相（如超流、反铁磁），而围绕这个新稳态的涨落，就可以用环路展开来处理。这本质上是用变分的思想（引入辅助场）重新参数化了问题，为后续的微扰计算（环路展开）提供了一个更好的起点。
6.3 数值变分与解析展开的交叉验证
在解决具体模型时，数值变分方法（如变分蒙特卡洛、密度矩阵重整化群）可以提供非常精确的自由能或基态能量数据。这些数据可以作为基准，来检验各种解析的环路展开方案在多大程度上是可靠的。例如，对于二维硬核玻色子模型，变分蒙特卡洛可以给出几乎精确的基态能量。然后，我们可以尝试用包含两圈甚至三圈修正的场论计算去逼近这个数值，从而判断该场论展开的收敛性和有效性。反过来，场论计算给出的解析形式（如临界指数、标度律）可以指导我们理解数值结果背后的普适类物理。
7. 一个具体案例：均匀玻色气体的超越平均场态方程
让我们用一个相对简单但重要的例子，来具体感受一下这些方法是如何运作的。我们考虑三维均匀相互作用玻色气体在零温下的态方程，即能量密度 (E) 作为密度 (n) 的函数。
7.1 平均场与李-黄-杨修正
在平均场（Gross-Pitaevskii）近似下，能量密度很简单：(E_{\text{MF}} / V = \frac{1}{2} g n^2)。这里 (g = 4\pi\hbar^2 a_s / m)，(a_s) 是s波散射长度。
李政道、黄克孙和杨振宁在1950年代用微扰论计算了第一项量子修正。他们的计算本质上相当于在路径积分框架下，考虑了涨落场（即非凝聚部分）的单圈图贡献。这个贡献来自于凝聚态粒子与激发粒子之间的相互作用，以及激发粒子之间的相互作用。最终得到的著名结果是： [ \frac{E}{V} = \frac{2\pi\hbar^2 a_s n^2}{m} \left[ 1 + \frac{128}{15\sqrt{\pi}} \sqrt{n a_s^3} + \cdots \right] ] 括号里的第一项是1，即平均场项。第二项 (\propto \sqrt{n a_s^3}) 就是单圈修正，它被称为“李-黄-杨修正”。无量纲参数 (\sqrt{n a_s^3}) 衡量了量子涨落的强度。当气体很稀薄（(n a_s^3 \ll 1)）时，这个修正很小，平均场理论很好。但随着密度增加或散射长度增大，这一修正变得重要。
7.2 更高阶圈图与变分方法的介入
李-黄-杨公式是一个微扰结果。如果我们要研究更稠密或相互作用更强的体系（例如在Feshbach共振附近），就需要考虑更高阶的圈图（两圈、三圈……）。这些计算非常复杂，涉及多重积分和可能的重整化。一些研究通过计算两圈图，得到了 (\propto n a_s^3 \ln(n a_s^3)) 等形式的更高阶修正。
另一方面，变分方法也可以攻击这个问题。例如，可以采用包含两体Jastrow因子的试探波函数：(\Psi(\mathbf{r}_1, ..., \mathbf{r}N) = \prod{i<j} f(|\mathbf{r}_i - \mathbf{r}_j|))，其中 (f(r)) 是一个变分函数，它在短程处描述两个粒子的排斥关联（因为硬核排斥），在长程处趋于1。通过变分蒙特卡洛方法优化 (f(r)) 并计算能量，可以得到从弱耦合到强耦合整个区间内都非常精确的态方程。这个结果是非微扰的，它包含了所有阶的关联效应，可以作为检验各种场论展开的“金标准”。
在这个具体问题中，“环路展开”（微扰场论）给出了一个解析的、可以逐阶计算的级数，适用于弱耦合区。而“交织变分法”（如Jastrow变分）提供了一个全局的、非微扰的数值解。两者结合，既给了我们解析的洞察，又给了我们精确的定量结果。
8. 前沿与挑战：当前方法的应用边界与未来方向
尽管“环路展开”和“交织变分”已经是我们工具箱里的利器，但面对日益复杂的量子多体问题，它们各自仍有局限，也催生着新的发展。
8.1 环路展开的收敛性与可重整化群
对于强关联系统或低维系统，按圈数展开可能收敛得很慢，甚至不收敛。例如，在二维玻色系统中，由于涨落更强，平均场理论完全失效，即使考虑单圈修正也不足以描述其著名的Berezinskii–Kosterlitz–Thouless相变。这时，我们需要非微扰的工具，如可重整化群。RG可以看作是环路展开的一种无穷阶重求和，它通过连续地积分掉高能（或短程）自由度，来系统地追踪耦合常数随尺度的变化。现代的功能RG方法，本质上是在一个变分的截断下，对无穷多圈图进行某种求和。
8.2 交织变分法的表达能力与计算成本
变分法的威力完全取决于试探态的选择。一个简单的试探态（如平均场）计算容易但精度有限。一个高度复杂的试探态（如具有长程纠缠的张量网络态）可能能精确描述系统，但优化和计算期望值的成本极高，甚至无法进行。这就是“表达能力”与“计算可行性”之间的权衡。
对于大规模玻色系统（如数百万个原子），高精度的量子蒙特卡洛方法虽然原则上精确，但可能面临“符号问题”的困扰。而基于神经网络量子态的变分方法，近年来显示出巨大的潜力。这些神经网络可以看作是一种极其灵活、表达能力强大的“交织”试探态，能够自动学习并表征复杂的多体波函数。如何将这种基于机器学习的变分方法与传统的物理洞察（如对称性、守恒律）以及场论的解析方法（如环路展开的图式）结合起来，是一个激动人心的前沿方向。
8.3 非平衡与非均匀系统的拓展
我们以上的讨论主要针对平衡态、均匀系统。但许多有趣的物理发生在非平衡或空间非均匀的情境下，例如量子淬火后的动力学、光晶格中的玻色气体、旋转或谐振子势阱中的凝聚体。将环路展开和变分法推广到这些情况是重大的挑战。
对于非平衡，需要用到Keldysh闭路积分等更复杂的场论形式，环路展开的图技术依然适用，但传播子和顶点规则更为复杂。变分法则可能需要基于时间相关的变分原理，优化的是整个时空演化路径，计算复杂度剧增。
对于非均匀系统（如谐振子势阱），平移对称性破缺，动量空间方法不再方便，通常需要在实空间处理。这给解析的环路展开计算带来巨大困难，但为实空间的数值变分方法（如基于格点的变分蒙特卡洛、密度矩阵重整化群）提供了用武之地。两者结合的一种思路是，先用数值方法获得非均匀的背景解，再在这个背景上对涨落进行局域的场论分析。
在我自己尝试研究光晶格中强关联玻色子的相图时，就深刻体会到这种结合的必要性。单纯用微扰场论，在强耦合区域完全失效；单纯用精确对角化，又受限于系统尺寸。最终，我们采用了一种混合方案：用实空间的变分簇方法得到一个不错的初态，然后在这个态的基础上，构建低能有效哈密顿量，再用摄动方法去计算一些动态响应。这个过程充满了试错，比如如何选择变分簇的大小才能既捕捉到关联又不至于让计算不可行，以及如何判断有效哈密顿量的截断是否合理。这些经验告诉我，没有一种方法是万能的，真正的进展往往来自于根据具体问题，灵活地搭配和改造这些基础工具。
