时间矩阵乘积态理论及其应用

彭 程，冉仕举

(1 中国科学院大学物理科学学院， 北京 100049； 2 西班牙光子科学研究所， 卡斯特利德费尔斯 08860)(2018年5月16日收稿； 2018年5月28日收修改稿)

在一维量子多体系统中，量子相变点附近的临界行为可以由纠缠熵关于系统尺寸或关于矩阵乘积态中截断维数的标度理论来刻画。基于(1+1)维共形场论，刻画一维量子临界性的纠缠熵标度理论已经被相对完整地建立起来[1-11]，但是在二维量子多体系统中，由于(2+1)维共形场的不稳定性，无法直接将现有的一维标度理论直接推广到二维量子系统。目前常用于刻画二维量子系统临界性的方法，例如，寻找基态长程关联或分析低温下的热力学性质，从数值计算角度来讲都对计算精度和计算复杂度提出了很高的要求。并且相对于一维而言，如何通过纠缠熵关于二维量子系统中投影纠缠对态截断维数的标度行为来刻画二维量子系统的临界性目前还不完全清楚[12]。因此，有必要发展二维量子临界系统的标度理论。

最近，针对一维无穷长链上的平移不变量子系统中的基态问题所提出的时间矩阵乘积态，可以有效地反映系统基态的物理性质，特别是刻画临界系统中关联长度发散和纠缠熵发散的行为，以及给出临界系统的中心荷。时间矩阵乘积态被定义在二维张量网络中的实(虚)时间演化方向，它是一维的连续矩阵乘积态[13]。发展关于时间矩阵乘积态的纠缠熵标度理论对刻画临界系统的性质十分有帮助。

在这个工作中，我们利用数值重整化群方法，将时间矩阵乘积态推广到二维量子系统，并用于刻画二维量子系统在量子相变点附近的临界性。无论所描述的系统是一维还是二维量子系统，时间矩阵乘积态始终是一维量子态，因此可以通过时间矩阵乘积态结合一维量子临界系统的纠缠熵标度理论描述高维量子临界系统的性质。我们还建立了时间矩阵乘积态所刻画的临界系统纠缠熵和动力学关联长度与空间方向的矩阵乘积态所描述的纠缠熵和空间关联长度之间的等价性。

我们分别在蜂窝状六角格子和kagome格子上的自旋1/2海森堡反铁磁模型中利用时间矩阵乘积态研究临界点附近的标度行为，发现二维量子临界系统中的时间矩阵乘积态反映出与一维量子临界系统中的矩阵乘积态类似的纠缠熵对数发散型标度律。

1 从数值重整化群流到时间矩阵乘积态

我们将提供一种基于数值重整化群的方法，建立起二维量子系统的基态纠缠熵与时间矩阵乘积态的纠缠熵之间的等价性。

首先，对数值重整化群做一个简单的回顾。尽管在二维或高维系统中不存在重整化流，但是在数值重整化群中通常人为地选择一条特殊路径来定义重整化流。这种做法相当于将一个二维或高维有限尺寸格点系统上的模型映射成一维有限长链上的等效模型，相应的代价是原模型中的部分最近邻相互作用在等效模型中被当作长程相互作用处理。由于该一维等效模型的存在，同样的方法可用于任意维度有限尺寸量子模型的讨论上。

考虑包含N个自旋的系统，其哈密顿量可一般性地表示为

(1)

(2)

式中β=1/T,T为温度。虚时演化的图形表示如图1(a)所示。利用数值重整化群流，系统的基态|ψ〉可表示为

(3)

(4)

(5)

方程(3)实际上是对任意多体量子纠缠态进行施密特分解，而将系统二分后左右两条半链之间的纠缠谱则被包含在Φ中。通过纠缠谱可计算出系统的纠缠熵。用方程(3)中的基态表示重写虚时演化过程，形式如下

(6)

(7)

(8)

在此定义下，图1(b)中的结构可表示为一系列有效哈密顿量的连乘。若对全部有效哈密顿量进行如方程(8)中的分解，那么将出现沿垂直于实空间方向(虚时间方向)延伸的两列新的矩阵乘积态，即时间矩阵乘积态，其表达式为

(9)

(10)

图1 虚时演化和时间矩阵乘积态的示意图Fig.1 Diagrams of the imaginary time evolution and the time matrix product state

在时间矩阵乘积态中，指标{cn}的作用相当于空间矩阵乘积态中的物理指标，而{an}和{bn}的作用相当于辅助指标。因此可以发现，方程(7)中所定义的有效哈密顿量实则为时间矩阵乘积态内积的转移矩阵。在时间矩阵乘积态的框架下，可以将方程(6)中的虚时演化过程重新表述为

(11)

该式等号成立的条件是Φ为式(7)中有效哈密顿量的最大本征态。

2 时间矩阵乘积态的纠缠熵与关联长度

由于时间矩阵乘积态是一维的平移不变矩阵乘积态，因此可以很方便地求出纠缠熵和关联长度，并且建立在一维矩阵乘积态基础上的纠缠熵标度理论对时间矩阵乘积态同样有效。由于时间矩阵乘积态的纠缠熵和关联长度分别等价于量子系统的基态纠缠熵和动力学关联长度，这为研究二维量子临界系统提供了便利。

需要说明的是，一方面上述等价性建立在数值重整化群流能够在多大程度上反映真实的二维量子多体系统的基态；另一方面，方程(8)中分解的对称性也制约着二者的等价性。我们通过数值模拟发现，即使这两方面要求并不能够严格地满足，时间矩阵乘积态依然能够反映出基态关联长度和纠缠熵的特性。

3 结果和讨论

首先，在蜂窝状六角格子上的自旋1/2各向异性海森堡反铁磁模型中，应用时间矩阵乘积态研究基态临界行为。系统的哈密顿量为

(12)

式中：J和J′分别是沿着如图2(a)所示的不同方向的耦合系数。在反铁磁相互作用情况下，哈密顿量中的耦合系数大于零。在热力学极限下，模型的临界点位于Jc=0.5(4)[17]。在临界点左侧(J/J′

图2 有限尺寸格点系统示意图Fig.2 Diagrams of the finite size system

图3 时间矩阵乘积态的纠缠熵ST和动力学关联长度ξT随观测位置的变化趋势图Fig.3 The position dependences of entanglement entropy STand dynamic correlation length ξT

利用密度矩阵重整化群方法研究该模型的基态性质，取有限尺寸的蜂窝状六角格子，如图2(a)所示，每个自旋占据格子中的一个交点，自旋与自旋之间的相互作用由黑色实线表示。图中阴影部分为选取的有限尺寸格子，相邻阴影部分之间的黑色实线代表边界相互作用。在后面的计算中，考虑两个方向的周期边界条件。为便于程序设计，将图1(a)中虚线圆圈中的两个自旋用一个等效的格点代替，因此蜂窝状六角格子等效为一个方格子，如图2(b)所示，每个等效格点在矩阵乘积态中用一个张量表示，矩阵乘积态以图中序号所表示的路径，遍历方格子中的所有格点。因此很容易发现，某些最近邻相互作用在矩阵乘积态中以长程相互作用的形式出现。由于矩阵乘积态中不包含格点的平移不变性，因此应在每个格点处观测时间矩阵乘积态的纠缠熵ST与关联长度ξT。ST和ξT与空间位置的关系如图3所示。为尽量减小边界的影响，选择远离边界的格点计算ST与ξT，并对被选中的格点取平均值。

图4 ST和ξT关于Trotter-Suzuki步长的收敛趋势示意图Fig.4 Trotter-Suzuki step dependencs of STand ξT

为排除Trotter-Suzuki分解的误差对结果的影响，令步长τ逐渐减小并趋于零，以观察纠缠熵和关联长度关于步长τ的收敛情况。尽管Trotter-Suzuki误差不能被完全避免，但是如图4所示，我们发现纠缠熵ST和动力学关联长度ξT在τ<10-4以后几乎达到饱和值，因此在后面的计算中固定取值为τ=10-6，以平衡Trotter-Suziki误差与机器精度之间的关系。

通过调节耦合系数J的大小，蜂窝状六角格子上的自旋1/2各向异性海森堡反铁磁模型会经历两个不同的相区，分别是从基态到第一激发态之间存在激发能隙的二聚体相以及无激发能隙奈尔相[17]。如图5(a)所示，在两个不同的相区内，纠缠熵ST和关联长度ξT关于时间矩阵乘积态的截断维数χ存在不同类型的标度关系。当J=0.1时，系统处于二聚体相，纠缠熵和关联长度随着截断维数的增加都快速收敛到有限大小的饱和值。由于时间矩阵乘积态实质上是一维矩阵乘积态，并且已知对于固定截断维数的时间矩阵乘积态，有限大小的截断维数将限制系统的关联长度[7]。对于基态有激发能隙的系统，基态关联函数呈指数衰减，也就是说基态只存在短程关联，关联长度接近晶格常数。因此随着截断维数的增加，系统关联长度会出现短暂的上升，当关联长度上升至接近晶格常数(饱和值)时，即使继续增加截断维数也不会再观察到关联长度的上升。而受到关联长度的限制，纠缠熵同样随着截断维数的增加而短暂上升，当关联长度达到饱和值的同时，纠缠熵也收敛到一个有限大小的饱和值。

图5 蜂窝状六角格子上自旋1/2海森堡反铁磁体的纠缠熵ST和关联长度ξT标度律示意图Fig.5 The scaling laws of the entanglement entropy STand the correlation length ξT

在各向同性点J=1处，系统处于奈尔相中，纠缠熵在出现有限尺寸效应以前，随着截断维数的增加呈对数型发散，关联长度与截断维数之间呈指数型发散，具体表示为

ST=αlnχ+const,ξT∝χκ.

(13)

上述行为与一维临界系统中的纠缠熵和关联长度关于截断维数的标度行为一致。我们在不同的尺寸下观察临界系统纠缠熵和关联长度的标度关系，二者满足

ST=cTlnξT+const.

(14)

其中，因子cT=α/κ类似于一维量子临界系统的中心荷，随着系统尺寸的增加，有cT～1。

我们还在kagome格子上的自旋1/2各向同性海森堡反铁磁模型中，利用时间矩阵乘积态的纠缠熵标度行为研究系统的基态。如图6所示，时间矩阵乘积态的纠缠熵和关联长度清晰地反映出系统基态的临界行为。我们的结果与支持kagome格子上的自旋1/2各向同性海森堡模型中的基态是无激发能隙的量子自旋液体这一结论的实验和理论预言[18-30]相一致。

图6 kagome格子上自旋1/2海森堡反铁磁体的基态纠缠熵ST和关联长度ξT关于标度律示意图Fig.6 The scaling laws of the entanglement entropy STand the correlation length ξTfor the spin-1/2 HAFM on the kagome lattice

4 总结与展望

本文提出一种用于刻画二维量子临界系统的一维量子态——时间矩阵乘积态。时间矩阵乘积态的纠缠熵和关联长度与二维量子系统中通过数值重整化群流所定义的纠缠熵和动力学关联长度是等价的。在蜂窝状六角格子上的自旋1/2各向异性海森堡反铁磁模型中展示了时间矩阵乘积态所给出的纠缠熵和关联长度在2个相区内的不同标度行为，在有能隙的二聚体相中纠缠熵和关联长度关于增加的截断维数先短暂上升后快速收敛，而在无能隙的奈尔相中纠缠熵关于增加的截断维数呈对数型发散，关联长度呈指数型发散。还通过kagome格子上自旋1/2各向同性海森堡模型对应的时间矩阵乘积态的纠缠熵标度行为，发现纠缠熵关于截断维数的变化满足临界系统的纠缠熵标度律。本文提出的方法很好地给出了有能隙的二聚体相与无能隙的奈尔相内的纠缠熵标度行为，但是在靠近相边界的区域，由于计算误差的影响(如裁剪误差、有限尺寸效应等)，标度行为的特点(收敛或对数发散)在数据上变得不明显，这需要在未来使用更精确的计算进一步研究，或发展更好的物理性质提取方法。同时，我们的结果暗示了通过(1+1)维共形场论可以对时间矩阵乘积态的临界标度律进行进一步研究。