一生能有多少爱 | 量子多体中的呐喊与彷徨之三

撰文 | 孟子杨 (中国科学院物理研究所)

来源：中国物理学会期刊网

我承认，这是一个俗气的题目。

很多年前，大概是还觉得《读者》之类的鸡汤文字集合是文学的时候吧，不知在什么地方见到过以此为题目的一本书。内容自然记不清楚了，大抵都是些content-free 的大实话：人生艰难，痛苦总是多于欢欣，没有爱的扶持我们怎么能够坚持下去云云。时过境迁，已经练就了铁石心肠鸡汤不入，只有这个题目还映在脑海里。

也是很多年前，大概是还有人看《一生能有多少爱》这样的文字的时候，在凝聚态物理学量子多体研究领域中，物理学家们正在真真实实地为了非费米液体这样的人生难题而痛苦。自从先师朗夫子和他门下的AGD诸贤人创立了朗道费米液体理论框架以来，后学们格物致知的人生从未显得如此困难，从未遇到过用圈图计算不能解决的问题，从未如此需要爱的扶持。“天不生列夫，万古如长夜”，眼下我们又遇到了新的挑战，而列夫已经走远了。

那么挑战到底是什么呢？

我们接着上一篇往下说[1]。朗道费米液体理论描述了金属导体的基本性质，比如低温下，电阻满足

是一个和温度无关但是和费米面上态密度有关的常数；比热满足C~T的关系等等。这些现象的本质原因是费米面附近无穷多的电子，尽管感受到了彼此之间的量子多体相互作用，但作用的结果不过是重整化了它们，将质量变成有效质量，磁矩变成了有效磁矩。这些具有有效质量、有效磁矩的准粒子，在能量接近费米能的过程中，愈发地接近无相互作用的自由费米子，可以用量子力学单体问题成熟的计算手段来处理，故其行为可以被准确计算并解释相关实验结果。套用最近刚刚去见朗夫子的凝聚态物理学亚圣Anderson的话，此处其实不是more is different，而是more is the same，问题变得简单了。

然而， 量子多体世界一如人生，复杂才是它的本质，痛苦总是多于欢欣。当代以来许多关联电子系统的实验现象不能用费米液体理论描述。高温超导体相图中的奇异金属和赝能隙区域，重费米子量子临界区中的非费米液体等等，抛开材料的千别万化和现象的纷繁复杂，单说电阻随温度的依赖关系这一条，就不是

并且a < 2 。这怎么解释？本篇我们就从这里切入，谈谈巡游电子量子临界理论的构造。

后世的物理学家多年勤勤恳恳地整理非费米液体实验上的结果，抽象出最简单的模型，然后运用解析和数值工具，研究如是的模型中是否蕴含着非费米液体的迹象，其中一大类是费米子和临界玻色子耦合模型[2]，其微观晶格模型的拉式量可以写成如下的形式：

其中

此处的相变就是玻色子自己在(2+1)时空维度的Wilson—Fisher 普适类相变，但是耦合项的存在改变了玻色子和费米子的命运。λ 将玻色子的临界动力学涨落传递给了费米子，当m调节到相变点时，费米子在时空的路径积分框架下具有了有效相互作用，而这些相互作用的费米子还可以再反馈给玻色子，改变玻色子的临界行为， 使得整个耦合系统L[ψ,ϕ] 处在了一个新的量子临界点上。这个零温的量子临界点和其上的量子临界区域，就是孕育非费米液体的地方。如是的模型和基本的相图都在图1中展示。

图1 (a)费米子与铁磁量子Ising 自旋耦合模型[3]。λ = 1,2 表示正方晶格的两层自由费米子，中间一层为量子Ising 自旋，自旋之间具有铁磁相互作用J，横向的磁场h 引入Ising 自旋的量子涨落，涨落通过耦合传递给费米子，引入费米子之间的有效相互作用。(b)模型的示意相图[4]。横场h 大于临界值时，Ising 自旋量子无序，耦合不发生作用，费米子处于顺磁费米液体，费米面自旋上下简并；横场h 小于临界值时，Ising自旋铁磁有序，耦合的结果是费米面自旋简并劈裂，费米子处于Ising 铁磁金属状态，但由于此处玻色子涨落很弱，Ising 铁磁金属仍然是费米液体。只有在横场h 处在临界值时，玻色子量子涨落最强，系统在零温时处于量子临界点。费米子之间通过玻色涨落生发出强烈地相互作用，系统进入量子临界的非费米液体区，其基本的性质是费米子自能满足频率的分数幂次，对于铁磁量子临界金属这个幂次是2/3。(c)和(d)的解释见前文[1]，与(a)和(b)发生同样的物理过程，只不过是反铁磁量子临界现象。在Ising 反铁磁序形成的区域，费米面折叠，折叠出的交点是hot spots。在量子临界区中，hot spots 动量处的费米子进入非费米液体[5]，其自能的幂次为1/2

有了模型，下一步就是如何求解。此处人们又遇到了困难。从解析上可以进行一些微扰计算，仍然是运用圈图计算的方法。在这个方面，还是俄国人支配着领域的进展，足以显示朗道学派传统余威之所及。圈图计算是一种级数展开的概念，计算结果的可靠性依赖于级数是否可以在一定阶数收敛。遗憾的是，到目前为止，对于上述的耦合模型，普遍意义上的收敛形式还没有找到。但是在不同的展开阶数上，人们对于问题还是有一些基本的认识。比如在最低阶的近似下，有所谓的Hertz—Millis—Moriya 结论，这是随机相位近似(random phase approximation， RPA) 层次上的结果。在图1 中所示的2 维晶格系统量子临界点上，Hertz—Millis—Moriya对于铁磁和反铁磁量子临界点给出了高斯不动点普适类的预言。这个结论显然过于简单，基本上就是说费米子和临界玻色子耦合的结果是把量子临界问题变成了比原来纯玻色子的Wilson—Fisher 量子临界行为还要简单的问题，虽有其存在的意义， 但又是一个more is the same 精神的结论，和实验观测还是有差距。比RPA高级一些的理论是所谓的大N展开，此处N可以是费米子的内禀自由度数目 (比如自旋或者flavor数目)，或者系统中低能费米子自由度的数目 (比如反铁磁量子临界问题中，因为费米面被玻色子涨落矢量所折叠而造成的费米面上hot spots的数目)。大N展开确实可以把系统带到新的不动点和普适类，比如对于铁磁和反铁磁的量子临界点，可以得到电子的自能出现频率的分数幂次的行为，

在上一篇文章中，我们讲到相互作用下的电子格林函数，

其中

是相互作用系统中的准粒子权重，

是重整化之后的准粒子色散关系，

描述重整化之后准粒子的寿命。其实我们之前反复提到的，人们不知道非费米液体的理论描述，主要是在说非费米液体的自能应该满足什么样的表达式，尤其是做为准粒子寿命的虚部的解析表达式。而我们上面提到的大N展开的结果，就是对于临界点上费米子自能的频

time的倒数 (也就是自能虚部) 来近似也很粗糙，考虑高阶修正之后幂指数可能再发生变化，很有希望和实验观测结合起来。所以，眼下的问题就是，能否在模型的严格计算中，看到理论预言的自能虚部有偏离费米液体的行为。

说来惭愧，几十年过去了，这一步我们才刚刚看到。理论上的近似计算当然很早就得到了上面的1/N展开结果，但是后期人们主要在讨论更高阶的修正会如何改变这些低阶的结果，关于级数是否会最终收敛此事目前仍然没有定论[8]，所以不同的阶数结论不同。比如对于铁磁量子临界点，人们甚至在讨论连续相变本身是否存在，抑或是其实系统变成了1 级相变[9]；而对于反铁磁量子临界点，高阶圈图的计算似乎在指出随着费米速度向着玻色子低能波矢方向转动，甚至连临界指数本身也会发散[10]。总之，这些讨论牵扯到进阶的重整化群语言，朗学后人和业内专家们讨论起来自然是棋逢对手，口吐莲花，十分过瘾；但是在广大读者看来庶几针尖上站几个天使之类的玄学清谈，我们还是就此打住，转而介绍晶格模型数值计算的结果吧。毕竟大家都会看图和数数。

图1中的费米子与临界玻色子耦合模型，可以通过量子蒙特卡洛模拟得到确定性的结果。图2就是数值计算所得的Ising铁磁和反铁磁量子临界点模型的相图[3, 5]。相图中的数字、符号和简称大家还是要看具体的文献，但是基本的安排 (如有颜色阴影的区域为铁磁和反铁磁金属，白色区域为顺磁费米液体，零温的量子相变点之上的温度区域为量子临界区) 和图1中的示意是一致的。为了能够进行这样问题的计算，在算法的设计和优化方面也需要做相当的前期准备，比如为了克服费米子行列式蒙特卡洛的计算复杂度高和临界慢化等问题，我们发展出了自学习蒙特卡洛方法[11]；为了增大在费米面上hot spots附近的动量解析度，我们设计了鸸鹋算法[12] (鸸鹋英文名EMU，我们的算法英文名EMUS-QMC: Elective Momentum Ultra—Size Quantum Monte Carlo Method, 简称EMUS)。其实这些计算方面上的进步，很大程度上来自于从事数值计算与从事解析理论的物理学家的良性互动，这也是近年来量子多体研究的一个大趋势。对于费米面和临界玻色子耦合的量子临界问题，只要把所有的计算资源集中到低能的自由度上，增大这里的能量—动量的解析度，就可以抓住问题的本质，有效地提高对于临界行为的准确描述能力，而不用操心原本微观晶格模型中的高能自由度。其实场论模型就是这样设计的，现在的数值计算物理学家，在量子多体问题中已经开始直接模拟量子场论模型，通过模型设计和算法的进步，逐步克服以前微观模型的局限，明心见性，直指本心。

图2 (a)费米子与铁磁量子Ising 自旋耦合模型的量子蒙特卡洛计算相图[3]；(b)费米子与反铁磁Ising 自旋耦合模型的量子蒙特卡洛计算相图[4]

大概花了3 年左右的时间，正是在和解析理论物理学家的互动中，我们才意识到问题所在[4]。原来之前所讲的大N展开圈图计算，其实都是在温度严格为0的量子临界点上进行的，当年的朗学后人们也没有想到有朝一日还真有一批实心眼的数值计算物理学家会真正设计出晶格模型，并费力优化算法，就是为了能够检验场论模型的预言。所以如果现在认真看待量子蒙特卡洛晶格计算，会发现这样的计算都是在有限的晶格大小，有限的温度下进行的，虽然我们想让计算的结果尽量逼近零温的结果，但是有限温度涨落所引入的对于自能的修正是不能忽略的。这一点在当年的解析计算中没有考虑，那个铁磁量子临界点上

是纯粹的零温的结果。故而一般而言，如果考虑了有限温度的贡献，可以预期自能应该分成热涨落(thermal)和量子涨落(quantum)两部分的贡献，即

行文至此，在物理内容上可以打住了。基本信息就是现在后学们已经可以在费米子和铁磁临界玻色子耦合的模型中，通过量子蒙特卡洛计算，再扣除掉有限温度自能的贡献，明确地揭示量子临界涨落所造成的非费米液体行为，非费米液体的自能和其分数化幂律的渐进行为也可以被看到。有了这样的认识，在其他的量子临界涨落模型中，如反铁磁、电荷密度波、规范场涨落等等，当然可以进行如是的理论与数值结合的分析，逐步建立起如费米液体一般的非费米液体理论框架——当然，我们这些后学们还得如先师一样，给这个框架起一个不那么生硬拗口的名字。

话题再回到一生能有多少爱，如果说非费米液体鸡汤中的复杂现象，我们可以就着蒙特卡洛的猛火，可劲儿地搅拌翻腾理出头绪，并把先贤的道路往前推进的话。世界的变化就真是实心眼的科学家们所不能理解的了。人生的底色是痛苦、疲惫和困厄，但目之所及，还是有人忍受着生死的考验进行着体现着生命意义的创造，当然大多数人在中途放弃了。创造的内容有高下，深远者如先师朗夫子和诸多亚圣贤人为后人构造科学的框架和指出道路，后学们在这道路上前进，总会发现新的风景，得到新的认识，再为后人推广知识的边界。我们普通人的创造大多数时候，无法如先贤般穿透历史战胜时间 (不过大家也看到，朗夫子的理论也不过半个世纪，已经不足以解释新的现象了)，但就是在这样勉力为之的过程中，一个一个的普通人找到了寄托，战胜了困厄的人生，扬弃了时代和社会加诸于我们身上的枷锁。创造性的活动让我们操心劳神，忍受着生存、名利、人性中的愚蠢偏狭所造就的种种不合理，但也让我们在创造中忘记了这些额外的痛苦，和人生达到了动态的和解，找到了人生中的爱。有一个朋友，平日里写诗，虽然水平不入流却乐此不疲，最近又给我发来新作，记述他在眼下这个瘟疫横行的世界上的一次际遇。虽然仍然水平有限但其勉力为之的意思庶几近乎我这些磕磕绊绊的文字，抄在这里做为本文的结尾吧。

底特律机场与爵士乐

瘟疫正在世界上横行，动荡的时代更让人时时追问生命本来的目的。阴差阳错中穿行于大洋彼岸，在底特律机场碰到现场爵士乐钢琴表演，黑人乐手自信忘我的神情和清脆激越的琴声片刻间安抚了我焦灼的心，但这表演却被大多行色匆匆的旅客忽略。遂想起生命旅途中的迷惑亦如此，美好常见于疲惫之中，匆匆的行旅又有多少是不自觉的盲从。站在琴声中既悲伤又欣慰，记录下来示不忘也。

在底特律机场里听爵士乐

冬日的阳光从云层中洒出

抬头见大厅钢窗外的天

轻快如少年

琴声珍珠般灵动

清脆激越

黑人乐手迷人的投入

被来往的旅人所忽略

美好总是飘渺的

又有谁会驻足？

瘟疫肆虐

只求不要再放大人们心中的恶

争夺、欺诈

宁可毁了自己只为让别人难受的撕裂

这几年来，可曾少见？

生命的意义真如这琴声

只有想听的人才能获得慰藉

机场里的街头音乐家

自遣自娱自徘徊

只有一个异国的旅客

片刻间被他熨平了满心的焦渴

停下来，纵然只有一小会儿

激赏和忘情

真的只有一小会儿

又要开始下一段疲惫的旅程

少年的心事毕竟会远去

只愿人性的恶不要胜于瘟疫

真的只有一小会儿

在琴声中看到了善良

想起来有谁说过

不知该原谅些什么

诚觉世事尽可原谅

参考文献

[1] 孟子杨. 白马非马，非费米液体—非—费米液体. 物理，2020，49(3)：186

[2] Abanov A，Chubukov A V，Schmalian J. Quantum- critical theory of the spin- fermion model and its application to cuprates: Normal state analysis. Advances in Physics, 2003，52：119

[3] Xu X Y，Sun K，Schattner Y et al. Non-Fermi-liquid at (2+1)d ferromagnetic quantum critical point. Phys. Rev. X，2017，7：031058

[4] Xu X Y，Klein A，Sun K et al. Extracting non-Fermi liquid fermionic self-energy at T=0 from quantum Monte Carlo data. 2003，arXiv：2003.11573

[5] Liu Z H，Pan G P，Xu X Y et al. Itinerant quantum critical point with fermion pockets and hot spots. PNAS，2019，116 (34): 16760

[6] Altshuler B L，Ioffe L B，Millis A J. Lowenergy properties of fermions with singular interactions. Phys. Rev. B，1994，50：14048[7] Chubukov A V. Self- generated locality near a ferromagnetic quantum critical point. Phys. Rev. B，2005，71:245123

[8] Lee S S. Low-energy effective theory of Fermi surface coupled with U(1) gauge field in 2 + 1 dimensions. Phys. Rev. B，2009，80：165102

[9] Metlitski M A，Sachdev S. Quantum phase transitions of metals in two spatial dimensions. I. Ising-nematic order. Phys. Rev. B，2010，82：075127

[10] Metlitski M A，Sachdev S. Quantum phase transitions of metals in two spatial dimensions. II. Spin density waveorder. Phys. Rev. B，2010，82：075128

[11] 孟子杨. 从德尔斐箴言到自学习蒙特卡罗. 物理，2017，46(4): 248[12] 孟子杨. 我们的壮游. 物理，2019，48(10)：683