第七章 分野（上）
与 Kamenev 教授的一席谈，如同在闷热的夏日午后注入了一股清冽的泉水，不仅涤荡了李叶心头的些许迷茫，更在理论构建的荒原上，为他标出了几条可能通向绿洲的路径。回到机房，他做的第一件事，便是将交流中的关键要点、思路启发、以及 Kamenev 教授提到的参考文献，事无巨细地记录下来。那些关于“辅助场平均场”、“涨落修正”、“阻挫可能的拓扑性”的提示，如同散落的珍珠，被他小心翼翼地拾起、串联，重新审视自己之前的工作。
他首先按照 Kamenev 教授的建议，深入研读了关于用 Hubbard-Stratonovich 变换处理自旋系统、特别是存在外场和阻挫情况的相关文献。这是一个他之前虽有涉猎但未曾深入的系统性技巧。现在，带着明确的问题和目标，他重新学习，理解如何将交错磁场项和反铁磁交换项，分别用辅助的“磁场”场和“配对”场来表示，写出配分函数，然后在平均场近似下，求解这些辅助场的自洽方程。
这个过程充满了复杂的代数运算和场论技巧。他需要在草稿纸上推导大段的公式，处理复杂的矩阵对角化，理解不同近似（如静态近似、忽略涨落等）的物理含义和适用范围。这比他之前尝试的、更为直接的强耦合展开要系统得多，但也复杂得多，常常为了一个符号、一个极限的处理，要反复查阅资料、推导验证。
“辅助场的引入，本质上是对相互作用项的解耦，”李叶在笔记上写道，试图理清思路，“交错磁场是外场，可以直接处理。反铁磁相互作用是四算符项，通过HS变换引入辅助配对场，变成二次型，就可以在平均场下处理。但关键是，这两个场是耦合的，交错磁场会影响配对场的空间分布，阻挫项又会引入额外的、可能是非局域的耦合……”
推导到具体形式时，他遇到了第一个难关。对于一维链，交错磁场是周期性的（比如，格点i上的磁场为 h_i = h_0 * (-1)^i），这使得辅助的“磁场”场也具有了相同的空间调制。而配对场，如果要描述最近邻自旋单态的形成，也自然地定义在键上。如何写出一个统一的、包含这两种辅助场，并满足自旋算符约束（每个格点自旋算符的模固定）的场论表述？这需要仔细选择辅助场的表示，并正确处理约束条件。
他翻出 Kamenev 教授提到的一篇关于“自旋液体平均场理论”的经典综述，仔细研读其中关于用复数希格斯场（plex Higgs field）表示自旋单态配对，用规范场（这里可能是U(1)或Z2规范场）来处理约束的部分。这让他意识到，问题可能比他最初想的更复杂。在阻挫存在的情况下，简单的平均场处理可能不够，需要考虑规范涨落的影响，甚至可能出现 deconfined 的分数化激发——这正是刘逸正在研究的领域。
“看来，即使是最简单的平均场起点，也需要小心处理规范结构。”李叶在草稿纸边缘标注。他决定先采用一个相对简单的处理：忽略可能的内禀规范场，专注于实数配对场和周期性磁场场的耦合，在平均场近似下求解，看看能否得到非零的配对场期望值，作为“二聚体化”的序参量。这当然是一个粗糙的近似，但可以作为起点。
他花了将近一周的时间，才初步搭建起这个简化的理论框架，写出了平均场自由能的表达式。这是一个关于周期性磁场场 φ_i 和最近邻配对场 Δ_{＜i,j＞} 的泛函。自洽方程要求自由能对这两个场的变分为零。由于空间的周期性，可以借助傅里叶变换到动量空间处理，但交错磁场破坏了平移不变性，使得问题必须在实空间或包含两个子格的原胞中处理，方程变得耦合且非线性。
解析求解希望渺茫。李叶转向数值。他编写了MATLAB程序，尝试用自洽迭代的方法求解这个方程组。他需要猜测初始场位形，然后迭代更新直到收敛。然而，由于非线性，迭代常常不收敛，或者收敛到平庸解（Δ=0，即无配对）。他需要调整参数（交换耦合J、阻挫系数α、交错磁场强度h0），尝试不同的初始猜测，小心翼翼地在参数空间中“巡逻”，寻找可能的有序相。
这又是一项繁琐而耗时的任务。屏幕上，迭代误差的曲线时而振荡发散，时而缓慢趋于一个值，但往往是零。李叶不断调整算法，加入阻尼，尝试不同的迭代格式。有时，在某个参数区间，当交错磁场h0足够大，阻挫α在一个特定范围时，迭代似乎开始收敛到一个非零的 Δ 值，而且 Δ 的分布呈现出空间调制——在某些键上更强，暗示着可能的“二聚体”形成模式。这让他精神一振。
然而，这只是第一步。得到非零的平均场解，只是说明在这个近似下，系统可能存在某种有序倾向。但平均场解是否稳定？量子涨落会如何修正它？更重要的是，这个平均场解对应的低能激发是什么？是否与数值模拟中看到的、暗示分数化自旋子连续谱的特征相符？要回答这些问题，他需要回到完整的路径积分表述，将辅助场的涨落积分回去，得到有效作用量，然后分析其低能激发谱。这涉及到计算高斯涨落行列式，分析其本征模式……又是一个庞大的工程。
李叶感到既兴奋又疲惫。兴奋的是，Kamenev 教授指出的这条路径似乎是可行的，他至少找到了一个切入点，并且得到了初步的、定性的迹象。疲惫的是，这条路径如此漫长，充满了技术细节，每一个步骤都可能出现新的问题。但他知道，这是从猜想到理论必须跨越的桥梁。他决定双线推进：一方面继续优化和改进这个平均场计算，探索更广泛的参数空间，并尝试将阻挫项更细致地纳入辅助场框架（Kamenev教授提到阻挫可能引入拓扑项，这需要更精巧的场论表述）；另一方面，他必须开始着手处理涨落修正，这需要他进一步学习泛函积分、高斯积分、行列式计算等场论技术。
就在李叶沉浸在新思路的探索和旧工作的深化中，几乎进入一种“闭关”状态时，317宿舍里，其他人的轨道也在悄然发生着变化，并且开始呈现出越来越明显的分野。这种分野，不仅是研究方向的差异，更体现在工作状态、思维习惯乃至日常生活的节奏上。
最显着的变化发生在刘逸身上。自从得到方文教授关于数值合作的可能线索后，刘逸仿佛被打通了任督二脉。他不再像过去几个月那样，被一种无形的焦虑和茫然所笼罩，而是进入了一种高度聚焦、目标明确的工作状态。他的书桌一角，整齐地码放着方教授推荐的和自己查找到的关于Z2规范场、自旋子平均场、随机相位近似（RPA）计算的文献，每一本都贴满了彩色标签，空白处写满了密密麻麻的笔记和推导。
他首先集中精力攻克方教授指出的第一个难点：平均场解的稳定性分析。这需要计算高斯涨落下的有效势，判断平均场解是局域极小值还是鞍点。这涉及到计算复杂的海森矩阵（Hessian matrix），并分析其本征值。刘逸花了整整三天，推导了描述规范场和自旋子场耦合系统涨落的二次型作用量，然后尝试对角化。计算中，他遇到了如何处理规范固定（gauge fixing）的棘手问题。Z2规范理论虽然比连续规范理论简单，但规范冗余仍然存在，需要在路径积分中小心处理，否则会导致发散或无意义的结果。他不得不重新翻阅关于离散规范场量子化的资料，学习如何引入适当的规范固定项（如类似于连续理论的Faddeev-Popov鬼场，但更简单），并确保物理可观测量与规范选择无关。
这个过程极其烧脑，刘逸常常在宿舍里一坐就是到深夜，草稿纸用掉一张又一张，上面布满了诸如“∑＜ij＞ σ{ij} (ψ_i^? ψ_j + h.c.)”、“积分掉ψ场得到有效作用量 Seff[σ]”、“在平均场解 σ_0 附近展开 σ = σ_0 + δσ”之类的式子。他的眉头紧锁，嘴唇无意识地抿着，完全沉浸在抽象符号的世界里。只有偶尔推导出关键的一步，或者发现之前的错误时，他才会抬起头，眼神放空片刻，然后迅速低头继续。
“怎么样，逸哥，规范固定搞定了没？”张海峰有时晚上回来，看到刘逸还在苦战，会凑过来看一眼那如同天书般的草稿，咋舌问道。
刘逸摇摇头，揉了揉发酸的眼睛：“还在弄。离散规范场的鬼场有点不一样，要引入Grassmann数但又有Z2约束……比较绕。不过快有眉目了。”
“牛逼。”张海峰佩服地拍了拍他肩膀，“这玩意我看着就头疼。还是我的蒙特卡洛实在，虽然也头疼，但至少是跟随机数打交道。”
刘逸苦笑一下，没接话。他知道张海峰有自己的苦处。他重新低下头，继续与那些顽抗的符号搏斗。终于，在又一天的深夜，他成功推导出了规范固定后的有效作用量形式，并写出了计算海森矩阵所需的主要表达式。虽然距离完整的稳定性分析还有几步（如具体计算矩阵元，对角化），但最艰难的规范处理部分，总算啃下来了。他长舒一口气，靠在椅背上，感到一种久违的、源自纯粹智力挑战被攻克后的疲惫与满足。
紧接着，他开始着手撰写给那位做张量网络（Tensor Network）的博士后的邮件。这封邮件他反复修改了多遍。他深知，对于一个潜在的合作者，特别是忙碌的博士后来说，邮件的清晰、扼要和吸引力至关重要。他必须用最简洁的语言，说明白三件事：1. 研究的问题是什么（阻挫晶格上Z2规范场耦合自旋子模型，探索可能的量子自旋液体相）；2. 他已经做了什么（平均场解，初步的一圈图修正，发现可能的新不稳定性迹象）；3. 合作能给对方带来什么（一个有趣且具有挑战性的具体模型，理论预言可供数值验证，潜在的共同发表成果的可能）。
他还精心准备了一个简短的摘要幻灯片，用图示化方式展示了模型、平均场相图、以及一圈图修正对自旋子能谱的可能影响。在邮件中，他诚恳地表达了合作的意愿，并提出了初步的合作设想：由对方利用张量网络方法（如无限 projected entangled-pair states, iPEPS 或无限 matrix product states, iMPS）计算该二维模型的基态相图和低能激发，与自己的场论预言进行对比；如果结果吻合，可以进一步深入探究特定相的性质（如拓扑序、分数化激发谱等）；如果存在差异，则可以共同分析原因，深化理论理解。
邮件发出后，便是忐忑的等待。那几天，刘逸变得有些心神不宁，每隔一会儿就要刷新一下邮箱。他既期待回复，又害怕被拒绝。这种期待与焦虑交织的心情，直到几天后收到对方的肯定回复时才得以平复。虽然只是初步的、试探性的合作意向，但这对刘逸而言，意义重大。这意味着他的工作，不再仅仅是理论上的“纸上谈兵”，开始有了接受独立数值检验的可能，也意味着他正式踏入了“理论-数值”结合的研究范式，这是现代凝聚态理论物理的主流路径之一。
他立刻与那位名叫陆云峰的博士后约了线上讨论。讨论持续了一个多小时，刘逸详细解释了模型的哈密顿量、对称性、以及他所关注的物理问题。陆云峰问了很多具体的技术问题，比如模型的具体参数范围、自旋子表示是复费米子还是马约拉纳费米子、规范场是紧致的还是非紧致的等等。刘逸一一作答，展现出了对模型细节的充分把握。陆云峰表示，这个模型确实有研究的价值，特别是其中可能涌现的Z2拓扑序和分数化自旋子，是张量网络方法可以很好处理的问题。他同意先做一些小规模的计算试试水，看看基态能量、近邻自旋关联等基本性质，能否与刘逸的平均场结果大致吻合。
合作的大门，就这样推开了一条缝。刘逸的生活节奏再次变得异常紧张。他需要不断与陆云峰沟通，提供更详细的模型参数和理论预期，同时还要继续完善自己的理论计算，特别是按照方教授的建议，开始着手用随机相位近似（RPA）分析规范涨落可能诱导的集体不稳定性。他几乎将所有时间都投入到了文献、推导、代码和邮件往来中，但这次，紧张中充满了期待和动力。他不再像以前那样，独自在理论的迷宫中绝望徘徊，而是看到了前方有同伴，有可以互相印证、共同前行的可能。这种转变，不仅体现在工作内容上，更体现在他的精神状态上——眼神中重新有了光，偶尔在宿舍与李叶讨论问题时，语气也变得更加自信和主动。
