动力学淬火过程中的不动点及衍生拓扑现象*

邓天舒 易为

(中国科学技术大学,中国科学院量子信息重点实验室,合肥 230026）

本文对近两年来有关淬火动力学过程中拓扑现象的研究做简要综述.这些动力学拓扑现象被动力学过程中的衍生拓扑不变量保护,与淬火前后体系的拓扑性质有密切关系.基于人工量子模拟平台的高度可控性,已在诸如超冷原子、超导量子比特、核磁共振、线性光学等众多物理体系中,通过对人工拓扑体系动力学过程的调控,观测到如动力学涡旋、动量-时间域的Hopf映射及环绕数、拓扑保护的自旋环结构、动力学量子相变、动量-时间斯格明子等诸多动力学拓扑现象.其中某些拓扑结构还可以在非幺正动力学淬火过程中稳定存在.这些研究将人们对拓扑物相的认识和研究从平衡态推广到非平衡动力学领域,具有重要的科学价值.

1 引 言

对拓扑相与相变的研究是现代物理的重要方向之一.在传统的金斯堡-朗道相变中,量子相与体系的对称性有深刻的关系,量子相变由局域序参量刻画,同时伴随着对称性的破缺.拓扑相和相变与此截然不同:拓扑相变一般不伴随对称性的破缺,拓扑相则与体系基态波函数在希尔伯特空间的几何性质相关并由拓扑不变量刻画[1,2].从人们发现量子霍尔效应至今已有三十余年了,此间人们对无相互作用体系中拓扑相的理解和分类日趋完善,对有相互作用体系中拓扑相的认识也有了长足进展.近年来拓扑物态研究领域的一个新兴方向是在凝聚态物质之外的体系中寻找、刻画拓扑相及与拓扑相关的现象.随着超冷原子、线性光学等量子模拟平台的蓬勃发展,具有拓扑能带及拓扑特性的人工拓扑体系相继在这些平台上得以实现[3−11].此外,人工拓扑体系的高度可控性使对动力学过程及非厄米体系中拓扑现象的实验研究成为可能[7−11],而这也驱动了相关领域的理论研究[12−16].值得强调的是,动力学过程及非厄米体系中的拓扑现象在传统凝聚态体系中一般很难被实现.因此各量子模拟平台上的人工拓扑体系具有重要的价值.

在这一背景下,本文简要综述近年来人们对淬火动力学过程中拓扑现象的理论和实验研究.所谓淬火,即体系由初始哈密顿量的本征态出发,在末态哈密顿量作用下演化的动力学过程.对于拓扑体系,一个重要的问题是该体系的淬火动力学过程中是否也会出现受拓扑保护的拓扑现象?更进一步,如果初态和末态哈密顿量具有不同拓扑性质的能带,体系淬火过程中的拓扑现象是由初态还是末态哈密顿量的拓扑性质决定?对这些问题,一个最直接的回答来自于对体系瞬时波函数几何性质的刻画.文献[17,18]研究了几种简单拓扑体系在淬火过程中瞬时波函数的拓扑不变量,结果表明瞬时波函数的拓扑不变量不随时间变化,即淬火过程的拓扑性质完全由初始哈密顿量决定.由于这些工作中考虑的体系都是无相互作用(或者平均场意义下无相互作用)的厄米体系,淬火动力学时间演化相当于对初始态做局域幺正变换.此时,由于局域幺正变换不破坏保护拓扑性质的对称性,瞬时态的拓扑不变量必然不变.然而,人们发现在淬火过程中会出现衍生的动力学拓扑现象,这些拓扑现象或与初态和末态哈密顿量的拓扑性质相关,或由衍生的动力学拓扑数保护,并在动力学过程中表现出可观测效应[19−31].有趣的是,即使基于同一淬火过程中观测到的拓扑现象,也可以构造不同的衍生拓扑结构加以解释.这一方面说明了淬火动力学过程中丰富的物理内涵,另一方面也会引发关于动力学拓扑现象普适描述的思考.

冷原子体系的高度可控性使其成为研究淬火动力学过程中拓扑性质的理想平台之一.2017年初，德国汉堡大学的研究组通过周期驱动光晶格,在冷原子中实现了基于Floquet Haldane模型的淬火过程,并利用动力学量子相变理论对淬火过程中的动力学涡旋进行了刻画[24].此后,清华大学的研究组从理论上研究了淬火动力学过程中的衍生拓扑结构,利用Hopf映射和环绕数重新解释了汉堡实验,并把环绕数与淬火前后哈密顿量的静态拓扑性质联系了起来[19].这一理论随后被汉堡大学研究组的另一实验所证实[25].与此同时,来自中国科学技术大学和北京大学的研究组在对拉曼晶格上的冷原子气体进行淬火时发现了动力学过程中原子自旋分布中的环形结构,并通过能带反转面的理论把环的出现与静态哈密顿量的陈数联系了起来[22,23,26].这些实验充分反映了淬火动力学中衍生拓扑现象的普遍性和丰富性.除冷原子气体之外,淬火过程中的动力学拓扑现象也在超导量子比特、单光子量子行走等量子模拟平台上被实验观测.目前这些模拟平台都还是通过操控单量子比特模拟费米体系拓扑模型的淬火过程,其优势是更易于进行实验操控和态重构,以及可以方便的实现耗散以研究非幺正动力学淬火过程.

我们将从几个不同的方面介绍这一领域的一些最新进展,着重描述与拓扑保护的动力学不动点相关的衍生拓扑结构.所谓动力学不动点,即在这些参数点上体系不随时间演化.而拓扑保护的动力学不动点的存在性与淬火前后静态哈密顿量的拓扑性质直接相关.由此,拓扑保护的动力学不动点为这些动力学拓扑现象提供了一个普适的基础,并成为连接动力学拓扑现象和体系静态拓扑性质的桥梁.

2 动力学量子相变及动力学拓扑序参量

我们首先介绍动力学量子相变及其与拓扑保护不动点的关系.先回顾动力学过程中Loschmidt回波的定义.Loschmidt回波是量子体系动力学演化中的一个重要可观测量,它描述体系初态在演化过程中的重现行为.如体系的时间演化由决定,则Loschmidt回波可表示为

从上述定义可知,动力学相变发生在瞬时态和初态正交的时刻,此时动力学自由能非解析.与热力学相变不同,动力学量子相变的临界点可以在体系演化时周期性地出现,而不是通过调节体系哈密顿量参数达到.动力学量子相变在多大程度上可以定义为相变仍然是有争议的话题.比如动力学量子相变是否继承了传统热力学相变中重要的性质,如在临界点附近的标度行为,相变普适类等.

与此同时,人们发现动力学量子相变有深刻的几何溯源,它的出现可以由动力学拓扑序参量联系起来[35].在这个意义上,动力学量子相变可以看作是动力学过程中的衍生拓扑现象.动力学序参量的定义与所谓的Pancharatnam几何相位有密切联系.Pancharatnam 相位是 Loschmidt振幅的辐角中与动力学相位无关的部分,源于体系动力学演化过程中的几何性质.以无相互作用的半满晶格模型为例,由于晶格平移对称性,不同晶格动量子空间的动力学演化解耦.此时,Pancharatnam 几何而为的 辐 角 .体 系 的 动 力 学 相 位 为的Bloch哈密顿量(假设不含时).

在拓扑不动点的基础上,可以进一步定义动力学拓扑序参量.对于一维体系,动力学拓扑序参量可以写作

从上述讨论可以看出,动力学量子相变及动力学拓扑序参量与拓扑保护的不动点有深刻联系.基于这一认识,下面我们将进一步揭示动力学量子相变与其他衍生动力学拓扑不变量的关系.

3 动力学拓扑不变量

我们在第2节讨论了拓扑体系淬火过程中的拓扑不动点.在拓扑不动点上,Pancharatnam相位在时间演化过程中恒为零.相位在拓扑不动点消失的本质原因在于这些不动点同时也是动力学不动点,即在这些参数点上体系的态不随时间演化.基于这一认识,人们发现可以基于不动点在动量-时间的二维参数空间里定义衍生动力学拓扑不变量[20,21,27,28,31].这些动力学拓扑不变量跟同样基于不动点的动力学量子相变、动力学拓扑序参量以及淬火前后哈密顿量的拓扑性质等紧密关联.目前文献里这方面的工作主要集中在一维两能带拓扑体系的淬火过程,我们也将具体介绍这类体系中的不动点结构和动力学拓扑不变量.高维多带体系的动力学淬火过程中的拓扑结构会更加丰富,但仍有待进一步研究.

对于无相互作用且具有晶格平移对称性的体系,不同动量部分的动力学解耦.因此,一维两能带拓扑体系的普适Bloch哈密顿量可表示为

如图1所示,淬火动力学过程可以形象的表示为瞬时态矢量从初态出发,在Bloch球面上绕轴的转动.这里有两种值得注意的特殊情况.

图1 Bloch 球上的动力学演化 (a)态矢量在 Bloch 球上绕 运动;(b)动力学不动点对应于 ;(c)临界点对应于.实线代表 (绿色) 与 (红色),虚线代表态矢量;假设初态处于 基态上,即 时态矢量与 方向相反Fig.1.Visualizing dynamics on the Bloch sphere:(a)State vector revolving around the axis;(b)illustration of fixed points when ;(c)illustration of critical points with .

由此可见,拓扑不动点的存在与动力学量子相变临界点的存在有直接关系.同时,基于这些动力学不动点,我们可以在衍生的动量-时间空间内定义动力学拓扑不变量.由于每个动量上的态演化都是周期性的,态矢量的演化在动量-时间参数空间内的轨迹是垂直于动量轴的圆.而在动力学不动点,态的轨迹是动量轴上的一个点.这样,任意两个动力学不动点间的态演化都可以映射到一个Bloch球上.我们由此可以定义衍生的动力学拓扑陈数

图2 淬火前后哈密顿量具有不同拓扑数时的典型斯格明子结构.黑色箭头为自旋在平面内方向,背景颜色对应自旋在与平面垂直方向上的分量大小,蓝色对应向内,黄色对应向外.竖直虚线为不动点位置,红色实线表示不同动量 点的周期Fig.2.Momentum-time skyrmions when pre- and postquench Hamiltonians possess different winding numbers.

4 非厄米体系中的非幺正淬火过程

前述的动力学拓扑结构也存在于非厄米体系的非幺正淬火过程中,前提是体系的动力学演化仍然是相干的[41].这一前提满足的条件体系具有宇称-时间对称性.一般而言,非厄米哈密顿量的本征值不是实数.但如果该哈密顿量有宇称-时间对称性,即,且本征态也为算符的本征态时,该本征态对应的哈密顿量本征值为实数[42−44].这里为宇称算符,为时间反演算符.如果宇称-时间守恒的哈密顿量所有本征值均为实数,则体系处于宇称-时间对称守恒相;反之,如哈密顿量某些本征值非实数,体系处于宇称-时间对称自发破缺相.

这里我们考虑最简单的宇称-时间守恒的拓扑体系,一维宇称-时间守恒的Su-Schieffer-Heeger(SSH)模型[45,46].

如图3(a)所示,该模型与标准SSH模型的区别在每个格点上均存在增益或损耗项.体系的Bloch哈密顿量仍然可记为,其中扑数可由全局贝里相对应的绕数刻画

图3 非厄米 SSH 模型及其拓扑相图 (a)非厄米 SSH 模型.在厄米SSH模型的基础上,每个格点上均有增益或损耗;(b)体系拓扑相图.实线为拓扑边界,虚线为宇称-时间对称与对称破缺区域的边界. , 为SSH模型的跃迁系数, 为增益损耗系数, 为绕数Fig.3.Non-Hermitian SSH model and its topological phase diagram:(a) Non-Hermitian SSH model;(b) topological phase diagram.

当末态哈密顿量处在宇称-时间守恒区域时,瞬时态矢量在Bloch球上绕着南北极旋转.此时,Bloch球南北极对应的动量就是动力学不动点.当末态哈密顿量处在宇称-时间对称自发破缺区域时,瞬时态矢量在Bloch球上会沿着大圆渐近趋近北极.此时不存在动力学不动点.基于这一图像,我们可以在动力学不动点存在的时候构造动力学拓扑不变量.事实上,可以证明当体系在宇称-时间守恒的不同拓扑相间淬火时,不同种类的不动点(即分别对应于南北极的不动点)会成对出现.此时,所有幺正淬火演化过程中的动力学拓扑构造和相关可观测量也会存在.而只要淬火前后的哈密顿量之一处于宇称-时间对称自发破缺区域,不动点以及动力学拓扑构造就不一定存在了.需要强调的是,在非幺正动力学淬火过程中,不动点和动力学量子相变的临界点不再出现在布里渊区的高对称点.这会导致多个临界时间尺度的出现.如图4所示,此时动力学量子相变的发生会存在多个周期,同时也会存在多个对应的动力学拓扑序参量.

图4 非厄密SSH模型淬火中的典型动力学自由能与动力学拓扑序参量 (a)动力学自由能 ;(b)动力学拓扑序参量 .在非厄米淬火过程中存在两个临界时间尺度及两个动力学拓扑序参量Fig.4.Dynamic free energy and dynamic topological order parameter in the quench dynamics of non-Hermitian SSH model: (a) Dynamic free energy ;(b)dynamic topological order parameter .

实验上,具有损耗-增益的宇称-时间守恒非厄米体系在光波导,微波腔阵列,线性光学等多中物理体系中实现.在一些难以实现增益的物理体系如冷原子气体、单光子量子行走中,可以实现基于纯损耗的被动宇称-时间守恒非厄米体系[49,50].被动宇称-时间守恒非厄米体系虽然不会有纯实数的本征谱,但体系的能谱及动力学性质可以方便地映射到基于损耗-增益的主动宇称-时间守恒非厄米体系上.最近,动力学量子相变、动力学拓扑序参量、动力学拓扑不变量等衍生动力学拓扑现象已经在具有被动宇称-时间守恒性质的单光子非幺正量子行走实验里被探测到[31].

5 总结与展望

本文集中讨论了几种基于拓扑保护动力学不动点的动力学拓扑结构的构造和关联.由于相关动力学拓扑不变量的定义都基于不同类拓扑不动点的存在,而后者又与静态哈密顿量的拓扑性质密切相关,拓扑保护的动力学不动点成为有机联系各种动力学衍生拓扑现象及平衡态拓扑相的核心要素.目前对于动力学拓扑现象的讨论多数集中在低维(一维或者二维)两能带五相互作用体系的突变淬火过程.高维多能带拓扑系统淬火过程中会出现更丰富的拓扑结构,如何系统的刻画它们是今后研究的一个方向.除本文讨论的衍生拓扑现象以外,其他已知的动力学拓扑结构如Hopf映射、环绕数、基于能带反转面的拓扑环等与动力学不动点及动力学量子相变临界点的关系则需要进一步研究讨论.因此,关于动力学拓扑结构的更普适的描述是个亟待解决的问题.