第十章 裂 隙
盛夏的蝉鸣在梧桐叶间做最后的挣扎，声嘶力竭，却难掩季节更迭的脚步。科大校园里，暑假的气息早已被紧张的研究氛围驱散。对317宿舍的年轻人们而言，时间不再是线性的流淌，而是被切割成一个个项目节点、代码运行周期、实验测量轮次、以及导师约定的汇报日期。每个人都像上了发条的齿轮，在自己选定的轨道上高速旋转，试图在既定的节奏中，捕捉那一丝可能带来突破的灵光，或是至少，不被甩出轨道。
李叶的生活进入了理论推导与数值验证深度交织的“双螺旋”模式。白天的大部分时间，他沉浸在辅助场平均场和低能有效理论的数学世界中。草稿纸和笔记本上，布满了复杂的公式、积分符号、以及不断被划掉重写的推导步骤。他成功地在忽略阻挫的简化模型中，得到了平均场方程的非平庸数值解，证实了交错磁场确实可以“诱导”出非零的配对场期望值，并且这种期望值呈现出空间调制，在强磁场区域配对更强，直观上支持“二聚体”倾向。这给了他初步的信心。
然而，将涨落积分回去，得到低能有效作用量的过程，则充满了陷阱。他需要围绕平均场解展开，计算高斯涨落行列式，然后追踪出对低能物理有贡献的软模。这个过程涉及复杂的泛函积分和矩阵对角化，他不得不反复查阅场论教材和文献，学习如何处理有空间调制背景下的涨落。在“静默连接”的帮助下，他能够快速掌握技术细节，但计算本身的繁琐和容易出错，仍然消耗着他大量的精力。常常为了一个符号的正负、一个积分的收敛性，要反复检查推导数遍。
与此同时，他也在加紧分析新的DMRG动力谱数据。他将从简化模型平均场计算中初步猜测的 Luttinger 参数 K 的范围，代入标准 Luttinger 液体理论预言的关联函数幂律衰减形式，与DMRG计算得到的实际关联函数进行拟合比较。结果喜忧参半：在某些参数区间，拟合得相当不错，幂律指数与猜测的 K 值大致吻合；但在另一些参数区间，特别是阻挫较强时，拟合效果很差，关联函数呈现出更复杂的衰减行为，甚至带有振荡。
“涨落修正，或者阻挫诱导的拓扑项，可能改变了简单的 Luttinger 液体行为。”李叶在实验记录中分析道，“平均场近似过于粗糙，可能漏掉了关键的低能模式，或者 Luttinger 参数本身是依赖于波矢 k 的？” 这迫使他回头重新审视理论推导，考虑是否需要在有效作用量中引入更高阶的项，或者用更非微扰的方法来处理阻挫的影响。理论预言的不足，反过来为数值分析提供了更具体的目标：他需要更精确地测量低能激发谱的精细结构，看看是否能分辨出不同的激发分支，或者是否存在能隙。
这种理论猜想与数值证据之间的不断对话、相互修正，构成了李叶研究的主旋律。进展是缓慢的、螺旋式上升的，每一次小的领悟，都建立在对无数次失败尝试的反思之上。他开始撰写论文的初稿，试图将现有的数值结果、初步的平均场分析、以及基于 Luttinger 液体猜想的物理图像整合起来。写作的过程异常痛苦，因为他必须清晰地陈述哪些是坚实的数值结果，哪些是基于这些结果的合理猜想，哪些是尚待证实或可能被证伪的理论尝试。他必须诚实面对工作中的不完美和未解之处，这种诚实，有时比得到一个漂亮的结果更让人煎熬。
刘逸那边的进展则呈现出不同的节奏。与张量网络博士后陆云峰的初步合作，在谨慎的乐观中推进。陆云峰利用无限 projected entangled-pair states (iPEPS) 方法，对刘逸的二维阻挫Z2规范场耦合自旋子模型进行了小规模的计算。初步结果表明，在某些参数区域，系统的基态能量、近邻自旋关联函数与刘逸平均场计算的结果大致趋势相符，这给了刘逸不小的鼓舞。然而，在另一些参数区域，特别是刘逸理论预言可能出现新不稳定性的区域，数值结果显得比较模糊，或者呈现出复杂的竞争，难以明确判定。
“iPEPS 方法本身在捕捉某些高度纠缠的量子自旋液体相时，可能需要非常大的 bond dimension 才能收敛，计算量很大。”陆云峰在邮件中解释道，“而且，你的模型中包含阻挫和规范场，基态可能接近简并，或者具有拓扑简并，这对张量网络优化也是挑战。我需要更多时间优化算法，尝试更大的 bond dimension。”
刘逸理解其中的困难。数值模拟，尤其是处理强关联二维系统，从来不是一蹴而就的。他一方面继续与陆云峰保持沟通，提供更详细的理论预期，帮助对方理解物理图像，以便更有针对性地进行数值搜索；另一方面，他按照方文教授的指导，开始系统地用随机相位近似（RPA）分析规范涨落可能诱导的集体不稳定性。
这章没有结束，请点击下一页继续阅读！