我爱纠缠如秋裤｜量子多体中的呐喊与彷徨之八_风闻

对弯曲空间（在数学上被称之为流形）的研究，是数学的一个大方向，数学家非常感兴趣，特别是关于这些流形所有可能的拓扑结构。其中一个研究思路就是寻找一些简单的流形，看它们是否可以作为复杂流形的基本构件：任何一个流形都可以通过这些基本构件的拼接来实现。比如二维流形的基本构件，就是“秋裤”和“西瓜帽”。任何一个二维流形都可以通过秋裤和西瓜帽的缝合来形成。这种“裁缝拓扑”是数学中的一个重要手段。

三个西瓜帽和一个秋裤可以缝制成一个球面。两个西瓜帽和两个秋裤可以缝制成一个环面。

流形的拓扑在物理上也有很重要的应用。目前物理学家非常感兴趣而且广泛研究的量子材料，其中有很强的多体量子纠缠。但多体量的纠缠是一个全新的自然现象，它们长什么样，内部有什么结构，很难说清楚，因为描写这些结构的语言（也就是数学理论）还有待（或者说正在）发明。

探索多体量子纠缠内部结构的第一个想法，就是把多体量子体系放到有各种不同拓扑性质的流形上，然后测量这个系统的一些物理性质，如基态简并度和基态所形成的纤维丛。这样具有不同拓扑性质的流形就成为了测量多体量子纠缠的探测器。有人猜想，具有所有可能拓扑结构的流形，是多体量子纠缠的完备探测器。也就是说，可以测量到所有可能的多体量子纠缠。这样，流形的拓扑结构和多体量子纠缠的内部结构，就有了一个很深刻的对偶关系。这就是为什么多体量子纠缠的内部结构被称之为拓扑序。正是因为这个深刻的联系，数学上研究拓扑结构的方法，也成为物理上研究多体量子纠缠的方法。抽象数学中的经典拓扑（就是上面说的“裁缝拓扑”）和代数拓扑在量子材料的研究中大显身手。今天这篇文章介绍了如何用“裁缝拓扑”中的“秋裤”来探测多体量子系统中的纠缠熵。

——文小刚

撰文 | 孟子杨

来源 | 本文选自《物理》2021年第3期

1 引 子

秋裤者，勤劳勇敢的中国人民御寒保腿温之神器也。从塞外北国到中原大地，从江南水乡及至岭南热土 (没暖气)，每当神州各地秋风乍起凉意渐浓的时候，妈妈总会关切地问一句“你穿秋裤了吗？”不论在北风呼啸大雪飘飘的北方，还是在阴冷潮湿室内室外同此凉热的南方，一条秋裤不仅让温度从腿部一直蔓延到心里，更唤起了你内心中一种依恋的情节，告诉你时代再变再内卷再996，国际形势再复杂再诡谲再亡我之心不死，我们漂泊的人生中总有一些不变的东西，家一样的东西，每年按时来到，让你安心，保你平安。这就是规律性和规律性给人带来的慰藉。

虽然穿不穿秋裤的争论仍在全球化的大潮中此起彼伏，但秋裤的存在，的确揭示了我们生活中的一个普遍现象，呼唤着我们心中的一种普遍心理。这样反映普遍规律的事物，愿意思考的人都喜欢琢磨一下。虽然他们中很多人并不见得会穿秋裤，比如笔者，但喜欢是真的。

科学研究当然也是揭示规律性的活动，从事这个行业的，也颇有一些愿意思考的人。有趣的是，在笔者熟悉的量子多体系统的研究中，竟也存在一个和秋裤颇为相像的事物，揭示出量子多体系统中无处不在的量子纠缠这样深刻的道理，让科研从业者中愿意思考的人安心，忘却身边的种种不顺遂，进入“此中有真意”的境界而更加起劲地探索其中的奥妙。他们发现通过秋裤的视角可以揭示量子多体系统从朗道—金兹堡对称性自发破缺，到量子相变，再到拓扑序长程纠缠和范畴对称性等等奇异的现象，秋裤之功善莫大焉。

这听起来好像有点离奇，待我为君细细道来。

2 量子纠缠

量子纠缠是一个深刻的概念，其全面的外延与内涵，当然不是这篇小文可以承担的。但是仅就在凝聚态物理学量子多体问题的研究中，人们逐渐认识到，量子纠缠的重要性体现在对其的测量能够反映量子多体系统的规律性——尤其是这样的规律性无法从常规测量中得到时——以至于在许多新奇量子物质形态的探索中，纠缠成了一锤定音的判据。

那么纠缠作为一个概念，怎样在量子多体系统的计算中进行量化呢？这就是纠缠熵。为了计算纠缠熵，需要定义约化密度矩阵 (reduced density matrix)。如图1所示的一个量子多体系统，我们将其分为A 和B 两个子系统 (也可以记为A 与其补集

)，那么就子系统A 而言的约化密度矩阵就是将系统波函数中属于B 的自由度积分掉，写成公式就是

然后随之而来的纠缠熵 (此处主要介绍雷尼(Renyi)熵) 就是

其中 q 是一个整数。q → 1 的极限就给出了冯·诺依曼熵，而在量子蒙特卡洛的晶格计算中，人们常常计算 q = 2 的Renyi熵。对于空间维度 d = 1 的量子多体系统，纠缠熵在量子场论的理论框架之下，以及密度矩阵重正化群 (density- matrix renormalization group，DMRG) 的严格数值计算之中，目前人们已经获得了十分深入近乎完善的理解。比如对于 s =1/2 海森伯模型的自旋链，假设链长为L，而其中的子系统A长为lA，其 q = 2 的Renyi熵就是

其中的常数c就是系统的中心荷 (central charge) ，是一个普适量(universal constant)(对于 s = 1/2 海森伯自旋链 c = 1)，而b是一个非普适量。这方面的文献已然汗牛充栋，比如文献[1]，笔者也不赘述了。总之纠缠熵的计算在1维的量子多体系统中是如此多见，以至于眼下做DMRG计算的人，在进行量子多体计算时，第一个要看的就是纠缠熵如何，“纠缠熵怎么样？”，“纠缠熵大不大？”，“纠缠熵和D怎么标定？”这样的话，在DMRG从业人员的日常交流中，已经平常到“你吃了吗？”，“你家孩子几岁了？”，“今天天气，哈哈哈……”这样的见面寒暄用语的程度了。

不过在这篇文章里，出于笔者个人的行业积习，我们谈谈纠缠熵在空间维度 d = 2 的量子多体系统中，如何通过蒙特卡洛计算得到，以及得到的结果反映了什么普适规律。当然，最重要的是，纠缠熵与秋裤到底是什么关系。

图1 空间维度为d 的量子多体系统，可以分为A，B 两个子系统，计算其纠缠熵。A 和B的交界，记为lA，就是下文中所说的面积律里面的面积。注意此处的面积是指A 和B 交界的面积，即d-1维度面积，对于d 维来说，lA是A 和B 交界的周长

3 秋裤的算法

其实纠缠熵之所以能够在DMRG这类计算中生根发芽，主要是因为这些计算可以得到量子多体系统的波函数 (起码是波函数很高程度的近似)，有了波函数，自然可以套用上面所说的公式，直接计算系统的约化密度矩阵和纠缠熵。但是对于空间维度 d = 2 和更加高维的晶格系统，如正方、三角、Kagome等晶格的Hubbard模型和海森伯模型等等，面对系统中自由度指数增加的问题 (常常也称为“指数墙”问题)，此时还想要普遍性地写出系统的波函数，就勉为其难了。取而代之的，是通过统计物理的方法，在 d+z 维度的相空间中 (d 为系统的空间维度，z 为系统的“时间”维度，其实z 有一个学名，叫做动力学临界指数，在我们讨论的问题里，很多时候 z = 1，也就是大家常常听到的那句其实是有错误的话“量子系统就是d+1维的经典系统”的由来)，进行满足系统统计规律的蒙特卡洛抽样，才有可能以代数增长的计算复杂度，克服指数增长的物理问题自由度，获得对于量子多体系统的严格结论。

那么问题就来了，既然量子蒙特卡洛是解决 d = 2 的量子多体系统的不二法门，怎么用蒙特卡洛的方法计算纠缠熵呢？办法自然是有的，这就是秋裤登场的地方。不过我们还是以 d = 1 ，z = 1 的量子多体系统为例 (比如上文中的 s = 1/2 海森伯自旋链) 来说明秋裤的算法吧，如图2所示 (此处当然也可以画一个空间维度 d = 2 的晶格，不过那样卷成的秋裤就太后现代了，还是一维比较清楚)，长度为L，温度为T (“时间”长度 β = 1/T ) 的量子多体系统的配分函数，可以用图2分母中的一个长为L 的圆筒表示，圆筒的周长就是系统的时间维度β，之所以卷起来是周期性边界条件的意思。如此系统的配分函数和种种物理可观测量，其蒙特卡洛计算自然是不在话下[2]。问题是纠缠熵需要把系统两分 (bipartition，当然实际计算中可能不止两分而是多分) 为子系统A和其补集，就是

然后刻画系统内子系统之间的纠缠。

图2 d =1，z =1的量子多体系统，2阶Renyi纠缠熵S2(A )的秋裤表示。系统本身的配分函数Z就是一个长为L 的圆筒，筒的周长是β。现在为了计算纠缠熵，需要把Z 平方一下，就是ln里面分母中的

，其中的∅是空集的意思，就是系统不两分，子系统A 为空集。而ln里面分子上的

，就是我们所说的秋裤构型或者秋裤配分函数。这个配分函数的裤裆就是子系统A，具有时间周期2β；裤腿就是A 的补集

，

，时间周期为β 。2阶Renyi纠缠熵

其实通过简单的推导 (比如在文献[1]中)，如是的以约化密度矩阵写出的纠缠熵：

可以转化为用配分函数写成

其中 q = 2 就是图2所示的，我们关心的2阶Renyi纠缠熵，

而ln里面分母上的

就是系统不两分时配分函数的平方 (此时子系统A为空集∅)。ln里面分子上的

就是我们喜欢的秋裤构型的配分函数。其中裤裆部分的空间维度就是子系统A，其时间上的周期为2β；裤腿部分的A的补集，其时间周期为β。因为我们此处讨论的是2阶Renyi纠缠熵，所以秋裤有两条腿。如果是3阶、4阶的话，那就会有3条、4条腿甚至更加复杂联通的秋裤了，是给外星生物穿的了。

图2中描述的，就是这样的一个2阶Renyi纠缠熵的秋裤表示。一旦可以写出配分函数，不管其形状多么奇怪，后面的蒙特卡洛抽样等等就有章可循了。重要的是有秋裤这样一个清晰的意象。当然算法的发展，也经历了颇为曲折的过程，而且还在不断优化之中，从最早的、略显尴尬的在蒙特卡洛中硬要DMRG上身的Swap算符[3]，到后来的运用非平衡抽样的手段，优化配分函数比例的统计质量[4]等等，更加简便和鲁棒的计算窍门，正在不断向前推进。但是正是拜秋裤的意象所赐，蒙特卡洛从业人员可以从配分函数入手计算纠缠熵 (其实也不止纠缠熵了，而是以其为代表的广泛的非局域测量方法，如威尔逊环 (Wilson loop)， 无 序 算 符 (disorder operator)等等，后文会有提及)，这样的视角在技术上对蒙特卡洛从业人员来说，比纠结于波函数的方法更加可亲，也更容易在不同的模型中程序实现。可以想见在不久的将来，随着计算方法的普及，量子蒙特卡洛的从业人员讨论起纠缠熵的时候，也可以像我们的 DMRG同行那样稀松平常，“今天天气，哈哈哈……”了。

4 秋裤的应用

有了通过秋裤配分函数计算Renyi纠缠熵的方法，我们就可以谈谈纠缠熵作为物理可观测量，在量子多体问题中揭示了什么普遍性的规律。其实很多读者都听说过纠缠熵要满足面积律 (area law)，就是说