第10章 裂 隙（1/3）

第十章 裂 隙

盛夏的蝉鸣在梧桐叶间做最后的挣扎，声嘶力竭，却难掩季节更迭的脚步。科大校园里，暑假的气息早已被紧张的研究氛围驱散。对317宿舍的年轻人们而言，时间不再是线性的流淌，而是被切割成一个个项目节点、代码运行周期、实验测量轮次、以及导师约定的汇报日期。每个人都像上了发条的齿轮，在自己选定的轨道上高速旋转，试图在既定的节奏中，捕捉那一丝可能带来突破的灵光，或是至少，不被甩出轨道。

李叶的生活进入了理论推导与数值验证深度交织的“双螺旋”模式。白天的大部分时间，他沉浸在辅助场平均场和低能有效理论的数学世界中。草稿纸和笔记本上，布满了复杂的公式、积分符号、以及不断被划掉重写的推导步骤。他成功地在忽略阻挫的简化模型中，得到了平均场方程的非平庸数值解，证实了交错磁场确实可以“诱导”出非零的配对场期望值，并且这种期望值呈现出空间调制，在强磁场区域配对更强，直观上支持“二聚体”倾向。这给了他初步的信心。

然而，将涨落积分回去，得到低能有效作用量的过程，则充满了陷阱。他需要围绕平均场解展开，计算高斯涨落行列式，然后追踪出对低能物理有贡献的软模。这个过程涉及复杂的泛函积分和矩阵对角化，他不得不反复查阅场论教材和文献，学习如何处理有空间调制背景下的涨落。在“静默连接”的帮助下，他能够快速掌握技术细节，但计算本身的繁琐和容易出错，仍然消耗着他大量的精力。常常为了一个符号的正负、一个积分的收敛性，要反复检查推导数遍。

与此同时，他也在加紧分析新的DMRG动力谱数据。他将从简化模型平均场计算中初步猜测的 Luttinger 参数 K 的范围，代入标准 Luttinger 液体理论预言的关联函数幂律衰减形式，与DMRG计算得到的实际关联函数进行拟合比较。结果喜忧参半：在某些参数区间，拟合得相当不错，幂律指数与猜测的 K 值大致吻合；但在另一些参数区间，特别是阻挫较强时，拟合效果很差，关联函数呈现出更复杂的衰减行为，甚至带有振荡。

“涨落修正，或者阻挫诱导的拓扑项，可能改变了简单的 Luttinger 液体行为。”李叶在实验记录中分析道，“平均场近似过于粗糙，可能漏掉了关键的低能模式，或者 Luttinger 参数本身是依赖于波矢 k 的？” 这迫使他回头重新审视理论推导，考虑是否需要在有效作用量中引入更高阶的项，或者用更非微扰的方法来处理阻挫的影响。理论预言的不足，反过来为数值分析提供了更具体的目标：他需要更精确地测量低能激发谱的精细结构，看看是否能分辨出不同的激发分支，或者是否存在能隙。

这种理论猜想与数值证据之间的不断对话、相互修正，构成了李叶研究的主旋律。进展是缓慢的、螺旋式上升的，每一次小的领悟，都建立在对无数次失败尝试的反思之上。他开始撰写论文的初稿，试图将现有的数值结果、初步的平均场分析、以及基于 Luttinger 液体猜想的物理图像整合起来。写作的过程异常痛苦，因为他必须清晰地陈述哪些是坚实的数值结果，哪些是基于这些结果的合理猜想，哪些是尚待证实或可能被证伪的理论尝试。他必须诚实面对工作中的不完美和未解之处，这种诚实，有时比得到一个漂亮的结果更让人煎熬。

刘逸那边的进展则呈现出不同的节奏。与张量网络博士后陆云峰的初步合作，在谨慎的乐观中推进。陆云峰利用无限 projected entangled-pair states (iPEPS) 方法，对刘逸的二维阻挫Z2规范场耦合自旋子模型进行了小规模的计算。初步结果表明，在某些参数区域，系统的基态能量、近邻自旋关联函数与刘逸平均场计算的结果大致趋势相符，这给了刘逸不小的鼓舞。然而，在另一些参数区域，特别是刘逸理论预言可能出现新不稳定性的区域，数值结果显得比较模糊，或者呈现出复杂的竞争，难以明确判定。

“iPEPS 方法本身在捕捉某些高度纠缠的量子自旋液体相时，可能需要非常大的 bond dimension 才能收敛，计算量很大。”陆云峰在邮件中解释道，“而且，你的模型中包含阻挫和规范场，基态可能接近简并，或者具有拓扑简并，这对张量网络优化也是挑战。我需要更多时间优化算法，尝试更大的 bond dimension。”

刘逸理解其中的困难。数值模拟，尤其是处理强关联二维系统，从来不是一蹴而就的。他一方面继续与陆云峰保持沟通，提供更详细的理论预期，帮助对方理解物理图像，以便更有针对性地进行数值搜索；另一方面，他按照方文教授的指导，开始系统地用随机相位近似（RPA）分析规范涨落可能诱导的集体不稳定性。

本小章还未完，请点击下一页继续阅读后面精彩内容！

RPA计算比之前的一圈图修正更为复杂，它涉及到计算规范场的极化率 bubble 图，然后将其代入Dyson方程，得到重整化的规范场传播子。不稳定性出现在重整化传播子出现极点（即发散）时。这需要计算复杂的费曼图积分，并自洽地求解方程。刘逸再次被淹没在繁复的代数中，但他已经逐渐习惯了这种强度。推导、编程、调试、分析结果……周而复始。他偶尔会和李叶交流场论技术上的心得，比如如何处理发散积分，如何数值求解自洽方程。两人虽然模型不同，但在场论技巧上有很多共通之处，这种交流往往能互相启发。

方文教授对刘逸的进展保持着密切关注。每周的讨论，刘逸都需要汇报新的推导结果和遇到的困难。方教授依然严厉，会尖锐地指出推导中的漏洞、近似的合理性、以及物理解释的模糊之处。但刘逸能感觉到，导师的态度在发生微妙的变化。批评依然存在，但不再仅仅是斥责，更多是建设性的、引导式的提问。有时，当他提出一个不错的想法，或者解决了一个技术难点时，方教授甚至会微微点头，简短地评价一句“有进步”。这对刘逸而言，已是莫大的鼓励。他逐渐明白，方教授的严格，并非针对他个人，而是源于对研究工作本身严谨性的极高要求。当你展现出与之匹配的努力和潜力时，你就能赢得他的尊重和指导。

张海峰依然在他那片名为“负符号问题”的沼泽中艰难跋涉，但似乎，在无尽的黑暗摸索中，前方终于出现了一丝极其微弱的光亮。在经历了无数次复 Langevin 方法的不稳定和错误结果后，他将目光转向了另一个更为复杂、但也理论上更严格的方法： Lefschetz 硫柱 (thimble) 方法。这个方法的核心思想，是将原本在实轴上进行、因快速振荡而导致符号问题的路径积分，通过解析延拓，变形到复空间中的稳定流形（即 Lefschetz 硫柱）上进行。沿着这些硫柱，被积函数的相位变化平缓，从而缓解甚至消除符号问题。

然而，实现这一想法异常困难。首先，需要找到被积函数（即作用量）在复空间中的临界点（即梯度为零的点）。对于多变量复杂作用量，这本身就是一个高维非线性方程的求解问题。其次，需要从这些临界点出发，构造出对应的稳定流形（硫柱），这涉及到解一组常微分方程（梯度下降流方程）。最后，还需要计算沿这些流形积分的雅可比行列式，并处理多个硫柱贡献的叠加（如果存在多个相关临界点的话）。

张海峰几乎是硬着头皮开始学习这个方法的数学基础。他阅读了提出该方法的原始论文，以及后续一些应用于场论和统计物理的文献。复杂的复分析、微分几何概念让他头晕目眩，但别无选择。他选择了一个尽可能简单的模型——一维 Hubbard 模型在有限化学势下——作为“试验田”，尝试实现硫柱方法。即使对于这个简化模型，计算也异常繁重。他需要编写程序寻找复空间中的临界点，这通常需要用到诸如牛顿-拉夫森法等数值优化算法，并且对初始猜测非常敏感。然后，他需要数值求解梯度流方程，以构造硫柱。

过程充满了挫折。程序常常无法收敛到正确的临界点，或者梯度流方程数值求解不稳定，导致构造出的硫柱不准确。他需要不断调整算法参数，尝试不同的数值积分方法。有时，程序能跑出看似合理的结果，但与已知的精确解或没有符号问题的其他方法结果一对比，却发现存在明显偏差。这意味着他的实践中还有隐藏的错误。

但这一次，张海峰没有像以前那样轻易放弃。或许是多次失败磨砺了他的耐心，或许是他内心深处的不甘被彻底激发。他将每一次失败都详细记录在日志中，分析可能的原因，然后逐一尝试改进。他开始在相关的学术论坛和邮件列表上小心翼翼地向国际同行请教，尽管大多石沉大海，但偶尔一两条有价值的回复，都能让他兴奋半天，并给他新的调试思路。

经过不知多少次的尝试和修改，终于，在一个凌晨，当他再次运行优化后的程序，将得到的一维 Hubbard 模型某些物理量（如双占据数）的结果，与基于贝特拟设的精确解进行比对时，两组曲线在误差范围内基本重合了！