第九章 余波
与Kamenev教授的交流，如同在平静的湖面投下一颗石子，激起的涟漪在李叶的研究中持续扩散，其影响远超一次简单的学术讨论。在接下来的几周里，李叶的生活被彻底重组，仿佛一台精密的仪器，被注入了新的指令和动力，开始以更高的效率和更明确的目标全速运转。
他首先面对的，是Kamenev教授提到的那些参考文献。这并非仅仅是几篇论文，而是通往一个更系统、更强大的理论工具库的钥匙。他花了一周多的时间，几乎不眠不休，沉浸在这些文献的海洋中。从经典的Hubbard-Stratonovich变换在自旋系统中的应用，到更现代的关于调制磁场诱导的量子相变、阻挫系统中的量子无序与分数化激发的综述文章，再到具体处理一维系统中自旋子Luttinger液体理论的场论技巧。他像一块干燥的海绵，贪婪地吸收着其中的思想和方法。
“辅助场平均场”，这个之前只是模糊概念的名词，现在在他的脑海中逐渐清晰、立体起来。他理解了如何将复杂的四算符相互作用，通过引入辅助玻色场（如配对场、磁化场）来“解耦”，转化为更容易处理的二次型。他明白了在存在外场（如他的交错磁场）的情况下，如何将外场与辅助场耦合，写出统一的场论作用量。他也开始领会，平均场近似只是第一步，更重要的是随后的“涨落修正”——将积分掉的高能模式（或辅助场的涨落）重新纳入考虑，得到低能有效理论。这个过程，就像是先用简单的线条勾勒出物理图像的轮廓（平均场），再用细腻的笔触添上阴影和细节（涨落），使其变得丰满、真实。
然而，理解原理与具体实现之间，横亘着巨大的鸿沟。当李叶真正开始动手，将他的一维阻挫自旋链模型“翻译”成场论语言时，复杂性立刻显现。交错磁场的周期性调制，使得辅助场也必须具有空间依赖性，平均场方程从简单的代数方程变成了耦合的非线性差分方程组。阻挫项的引入更是棘手，它不仅带来了次近邻相互作用，其特殊的符号结构（有些键是反铁磁，有些是铁磁？取决于具体模型）可能对辅助场的相位配置产生微妙影响，甚至如Kamenev教授暗示的，可能诱导出某种“拓扑性”。
李叶决定从简化版本开始。他先忽略阻挫，只考虑交错磁场下的反铁磁海森堡链。即使这样，问题也不简单。他需要写出包含交错磁场项和反铁磁交换项的配分函数，然后引入辅助配对场 Δ_{i,i+1} 来解决交换项。由于交错磁场的存在，系统不再具有平移不变性（平移一个晶格常数，磁场方向反转），因此必须处理两个子格（A和B子格）。这导致辅助场、自旋算符都需要在扩大了的原胞中标记。他花了好几天，与傅里叶变换、矩阵对角化、以及繁琐的代数运算搏斗，才得到了平均场自由能的表达式——一个关于子格依赖的辅助场 Δ_A, Δ_B 的复杂函数。
下一步是求解自洽方程：让自由能对 Δ_A, Δ_B 的变分为零。这通常没有解析解，需要数值迭代。李叶编写了MATLAB脚本，尝试在不同交错磁场强度 h0 和交换耦合 J 下，寻找非零的 Δ 解。这又是一个需要耐心和技巧的过程。迭代可能不收敛，可能收敛到平庸解（Δ=0），也可能依赖于初始猜测。他需要小心翼翼地在参数空间中扫描，寻找可能的有序相边界。
就在他埋头于这些繁复的数值计算和公式推导时，动力DMRG的后续计算也陆续返回了结果。他之前提交的、针对更大系统尺寸（L=80, 100）和更精细能量、动量网格的计算终于完成了。处理这些海量数据本身就是一个挑战。他编写了新的脚本，对得到的动力自旋结构因子 S(k,ω) 进行更细致的分析，包括尝试用各种函数形式（如幂律、洛伦兹峰等）去拟合低能部分的谱权重，进行有限尺寸标度分析，试图提取出更清晰的能谱特征。
结果喜忧参半。在某些参数下，低能区域的连续谱特征似乎更加明显，弥散的范围更宽，尖锐的峰更少。有限尺寸标度分析也显示，随着系统尺寸增大，连续谱的特征并未消失，反而有趋于某个极限的趋势，这似乎支持无能隙激发（或能隙极小）的图像。然而，拟合得到的“自旋子”连续谱边界（spinon continuum）依然模糊，难以精确确定。而且，在另一些参数区域，低能激发似乎又呈现出某些“软化”但尚未完全弥散的特征，可能暗示着接近量子相变点，但序参量尚未完全消失。
“看来，系统可能处于一个‘临界’或‘近临界’的区域，”李叶在分析笔记中写道，“并非严格的标准 Luttinger 液体，可能受到某些微弱相关的影响，或者平均场图像需要修正。” 这与他正在尝试的平均场计算中观察到的现象有些呼应——在某些参数下，平均场解对初始猜测非常敏感，似乎处于有序和无序的边界附近。
这章没有结束，请点击下一页继续阅读！