24小时热门版块排行榜

返回列表

新手已上路

木虫 (职业作家)

应助: 82 (初中生)
金币: 4153.3
红花: 422
帖子: 4933
在线: 616.2小时
虫号: 3474904
注册: 2014-10-14
性别: GG
专业: 元素化学

[交流] 相互作用拓扑绝缘体中拓扑不变量的计算取得进展已有2人参与

拓扑物质形态，如具有拓扑序的量子霍尔态以及拓扑绝缘体，是近年来凝聚态物理领域最大的进展之一。在拓扑绝缘体中加入电子之间的强关联相互作用，研究可能涌现的新的拓扑物态，以及拓扑物态之间由于相互作用驱动的拓扑相变，是人们关心的问题。前期拓扑绝缘体的理论和实验研究，主要集中于无相互作用或弱相互作用系统，而对于强关联电子系统中的拓扑相变和拓扑不变量的精确计算却付之阙如，并且往往存在把无相互作用系统中的拓扑不变量计算方法简单推广到相互作用系统中，造成概念和实践上的混乱等现象。因此，如何正确计算相互作用电子系统中的拓扑不变量，探测相互作用驱动的拓扑相变和拓扑物态，是十分基础的问题。
中国科学院物理研究所/北京凝聚态物理国家实验室(筹)副研究员孟子杨，中国人民大学博士研究生何院耀、邬汉青，教授卢仲毅组成的研究团队，在相互作用拓扑绝缘体中拓扑不变量的数值计算方面，开展了系统的工作。该团队首先发展了计算相互作用电子系统的拓扑不变量的数值方法，并运用大规模量子蒙特卡罗计算在常见的拓扑绝缘体模型中进行了数值实现；随后，他们找到并严格计算了几种相互作用驱动的、没有对称性自发破缺的拓扑量子相变，发现目前被大量使用的、通过单粒子格林函数计算拓扑不变量的方法对于这类拓扑相变完全失效，指出了现有的拓扑不变量理论在向强关联电子系统拓展时，不能依靠单粒子格林函数来计算，而需要在系统的集体激发中寻找拓扑信息，澄清了理论上的混乱和误区。
如图1 和图2 所示，在上述研究工作的前半部分中，该团队运用相互作用系统中的零频率单粒子格林函数，设计了准确计算Z2 不变量和自旋陈数的计算方法。由于量子蒙特卡罗方法都是计算有限尺寸的系统，自旋陈数因为有限尺寸效应得到的都是非整数的结果，该团队借鉴了团簇动力学平均场理论中的周期化方法，提出一种针对周期函数的差值方案，有效地克服了自旋陈数结果的有限尺寸效应，使得即使对于有限尺寸系统也能计算得到整数化的自旋陈数。
相互作用拓扑绝缘体中拓扑不变量的计算取得进展

图1 (a)相互作用拓扑绝缘体中的Kane—Mele—Hubbard模型，系统具有最近邻跃迁t，自旋轨道耦合λ，三近邻跃迁t3，以及可以驱动系统进入能带绝缘体的最近邻跃迁td，和库仑相互作用U；(b)系统的布里渊区
相互作用拓扑绝缘体中拓扑不变量的计算取得进展

图2 由三近邻跃迁t3/t 导致的系统拓扑相变(a)Z2不变量的变化，(b)和(c)显示系统自旋陈数的变化。在自旋陈数的变化中，L=12 和L=18 的数据为有限尺度蒙特卡罗计算的结果，可以看到有限尺度系统的自旋陈数并没有整数化；而(b)和(c)中的插值结果为我们提出的周期函数差值所得的拓扑不变量，可以看到有限尺度效应被克服，系统自旋陈数整数化
如图3 和图4 所示，在上述研究工作的后半部分中，该团队使用自己发展的计算自旋陈数的方法探测了一些相互作用驱动的、没有对称性自发破缺的拓扑量子相变，发现在某些直积态(拓扑平庸态，自旋霍尔电导率为0)系统中，数值计算得到了非零的自旋陈数。显然，这个非零自旋陈数并不对应物理上可观测的自旋霍尔电导率，因为后者在直积态为0。这说明用单粒子格林函数构造的自旋陈数在这样的直积态中失效，而该团队发现这样的直积态的共同特点是，没有无相互作用的能带绝缘体能与之不经历相变地、没有对称性破缺地连接。
相互作用拓扑绝缘体中拓扑不变量的计算取得进展

图3 相互作用拓扑绝缘体的双层Kane—Mele—Hubbard—Heisenberg 模型，模型中层间电子的海森伯J 相互作用将会导致新型的相互作用驱动的拓扑相变，相变中没有对称性自发破缺，系统进入拓扑平庸的二聚绝缘体直积态。这样的拓扑相变是用单粒子格林函数构造的拓扑不变量所无法描述的，系统的拓扑信息要在集体激发中寻找
相互作用拓扑绝缘体中拓扑不变量的计算取得进展

图4 相互作用拓扑绝缘体双层Kane—Mele—Hubbard—Heisenberg 模型中，层间电子的海森伯J 相互作用导致的拓扑相变，单粒子格林函数所构造的自旋陈数在相变的左右两边都是整数化的，而在大J 的二聚绝缘体相中，系统的自旋霍尔电导率应该为零，故其物理的自旋陈数应该为零。这就说明了，在这个相互作用驱动的拓扑相变中，单粒子格林函数不能正确反应系统的自旋霍尔电导率
这样现象背后的原因在于，在无相互作用电子系统中，拓扑不变量有确定的物理意义，即对应系统可观测的(自旋)霍尔电导率。但是在相互作用电子系统中，用单粒子格林函数构造的拓扑不变量并不对应于系统中真实的(自旋)霍尔电导率。长期以来人们一直有这样的怀疑，但是并没有找到真正的例子。该团队的研究工作首次通过严格的数值计算给出了用单粒子格林函数构造的拓扑不变量在相互作用拓扑绝缘体中失效的证据，给出了用单粒子格林函数构造的拓扑不变量的适用条件，指出了在一大类相互作用的拓扑物态中需要在系统的集体激发中寻找拓扑信息。
从现有的研究工作来看，相互作用电子系统中拓扑不变量的构造和分类还远不完备，需要使用什么样的拓扑不变量来描述相互作用量子相的拓扑性质是正在快速发展的领域。该团队这一系列工作推动了相互作用拓扑物态中拓扑不变量的数值计算和拓扑相变的研究。该系列工作分为两篇长文，背靠背地连载在最近的Physical Review B上。
致谢量子蒙特卡罗模拟所需的大规模的并行计算在中国科学院物理研究所量子科学模拟中心、中国人民大学物理系高性能并行计算物理实验室、天津国家超算中心天河1 号、广州国家超算中心天河2 号等计算平台上完成，计算过程中得到了天津国家超算中心的有力配合，在此感谢。
本文选自《物理》2016年第9期

e2700076a91e5862e30.jpg