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