马约拉纳零能模的非阿贝尔统计及其在拓扑量子计算的应用*

何映萍 洪健松 刘雄军3)4)†

1) (北京大学物理学院， 量子材料科学中心， 北京 100871)2) (量子物质科学协同创新中心， 北京 100871)3) (北京量子信息科学研究院， 北京 100193)4) (中国科学院大学， 中国科学院拓扑量子计算卓越创新中心， 北京 100190)(2020 年5 月31日收到; 2020 年6 月3日收到修改稿)

自1937年被预言以来， 马约拉纳费米子在粒子物理领域和暗物质领域就广受关注. 它们在凝聚态物理中的“副本”， 马约拉纳零能模(Majorana zero mode， MZM)， 被指出可以通过拓扑超导实现， 并由于满足非阿贝尔统计及可以用来实现容错的量子计算机而成为凝聚态领域最受关注的研究方向之一. 尤其在近二十年中， 马约拉纳零能模在理论和实验方面均取得了诸多重要进展， 一些综述文章对此做了较详细介绍. 本文将会重点回顾MZM的非阿贝尔统计性质以及它们在量子计算中的应用. 文章的第一部分首先简单介绍了凝聚态系统中MZM的理论发展并概述了在人工异质结体系中寻找MZM的最新理论和实验进展. 然后介绍了MZM非阿贝尔统计的基本概念， 并讨论这一性质怎样应用到量子计算中. 接下来重点讨论了利用MZM平台实现量子计算机的两个关键步骤: MZM非阿贝尔编织操作的实验实现方案和MZM量子比特的读取. 在这一部分里， 本文分别详细列举了现有的比较受关注的实现MZM编织操作和量子比特的读取实验装置. 最后， 文章介绍了在对称性保护的拓扑超导系统中实现马约拉纳的对称保护非阿贝尔统计的可能性.

1 引 言

马约拉纳费米子在20世纪30年代由意大利物理学家Ettore Majorana[1]在研究满足电子与正电子对称形式的Dirac方程解时提出. 当将Dirac矩阵用纯虚的矩阵表示， 即在Majorana表象下， Dirac方程的解为实数场， 存在的能量为零的定态解. 这个定态解满足场的共轭为自身. 用场论语言来说， 由这种场描述的粒子在实空间的产生和消失算符等于自身， 即 γ (x)=γ†(x) . 通俗来说，这相当于马约拉纳费米子的反粒子等于它自身. 在马约拉纳表象下，U (1) 规范对称性不再满足， 因此马约拉纳费米子是一种轻子数不守恒， 电中性、且自旋为 1 /2 的粒子. 这种粒子被预言以来， 在粒子领域至今未被证实. 高能物理中猜想中微子可能是马约拉纳费米子， 但是能支持这一观点的无中微子双 β 衰变尚未在实验室被观测到[2]. 由于马约拉纳费米子对进一步理解基本粒子物理有重要意义， 寻找马约拉纳费米子作为一种基本粒子的证据仍是高能物理的重要课题[3-6].

在凝聚态领域中， 马约拉纳费米子在近二十年获得大量关注和积极研究[7-12]. 凝聚态体系中的基本粒子只有电子， 马约拉纳费米子只能以类似准粒子的形式存在. 我们知道， 在电子形成的费米液体中， 基态是电子从能带底部填充到费米面而形成的费米海. 在费米面以上和费米面以下出现的激发分别是带负电的电子型激发和带正电的空穴型激发.在凝聚态物理中电子和空穴可看作彼此的准反粒子. 因此， 一个导致马约拉纳费米子的简单考虑是让电子和空穴叠加， 形成新的准粒子. 显然这类准粒子态可能在超导中存在. 超导中的伯格留波夫准粒子激发是电子和空穴的叠加. 当电子和空穴等幅值叠加时， 可以得到电中性的准粒子. 这是超导体中可能出现马约拉纳准粒子重要原因. 具体而言，在传统的s波超导体中， 库珀对是由自旋相反的两个电子配对形成的， 其配对哈密顿量可描述为

该配对哈密顿量可以等价写在Nambu空间中，即电子-空穴空间表象中， 从而变成 Hpair=， 其中如上所述， 空穴算符是电子算符的准反粒子算符 dk，↑=c†-k，↑. 这个简单形式告诉我们， 超导中形成库珀对的机制可等价看作电子和空穴的耦合. 这就使得在s波超导中出现的激发必然是电子和空穴的叠加态， 即一般满足bk=uck，↑+vdk，↓= u ck，↑+vc†-k，↓的形式. 然而这种破坏s波超导中库珀对的激发仍是自旋相反的电子和空穴的叠加， 显然不满足自共轭条件. 要出现满足自共轭的准粒子激发， 需要消除自旋自由度差异带来的影响. 因此， 马上可以想到的是由同种自旋态的电子形成库珀对的超导载体， 其对应激发态必然为自旋相同的电子和空穴的叠加， 从而可能出现马约拉纳费米子. 根据泡利不相容原理， 自旋相同的两个电子配对的空间波函数必须是奇宇称的. 其中p波超导是最简单的例子， 配对哈密顿量可 满足如下形式:

其中序参量满足奇宇称: Δp(k)=-Δp(-k) . 类似s波超导， 此时准粒子激发算符必然满足 γk=uck+vc†-k形 式. 如 写 到 实 空 间， 很 容 易 验 证 γ(x)=u(x)c(x)+v(x)c†(x) . 于是当 u =v*， 即得到马约拉纳费米子.

按照这个思路， 理论进一步发现， 在一维p波拓扑超导体的边界处和二维 p +ip 拓扑超导体的涡旋中心， 存在零能的马约拉纳零能模(Majorana zero mode， MZM)激发， 满足自共轭特征. MZM的粒子数算符满足 γ†γ=γ2=1 . MZM的自共轭性使得单个的MZM不存在粒子数空间. 两个MZM组成一个Dirac复费米子， 对应 | 0〉，|1〉 两个量子态， 因此每个MZM的量子维度是， 对应于半个费米子. 这个特殊的量子维度使得MZM满足非阿贝尔统计这一十分奇特的性质[7，13-15]. 另一方面， 组成一个复费米子态的两个MZM在空间上是可以分离的， 这样构成的复费米子态有着非局域性. 由于空间单独存在的MZM不能被局域扰动影响， 这使得非局域存在的复费米子态具有拓扑稳定特性[7，13-15]. 这两方面重要特性使得MZM有可能成为实现抗退相干拓扑量子计算的重要基本单元，因而成为近十年凝聚态物理中备受瞩目的热点研究课题.

过去三十年里， 遵循量子力学规律工作的量子计算机是物理研究的重要热点. 其中主要挑战在于量子态很容易受环境的干扰产生退相干现象， 使得难以实现大规模的量子计算机. 上述由MZM组成的非局域拓扑量子比特可望解决退相干问题， 拓扑量子计算正是因此动机而被提出. 利用非阿贝尔任意子的统计性质， 在受拓扑保护的任意子量子系统中， 量子信息能够以非局域的形式被存储起来， 这使得量子信息能抵抗局域的环境噪声影响， 从而在硬件上解决量子退相干问题. 另外， 量子信息的操控由非阿贝尔任意子的拓扑编织(braiding)来实现. 量子门操作不依赖任意子的交换轨迹， 其精度受拓扑性质保护[16]. 分数量子霍尔效应系统是最早被考虑作为实现拓扑量子计算平台的非阿贝尔任意子系统[15，17，18]. 进一步， 拓扑超导系统中的MZM被证明具有非阿贝尔统计性质后， 成为拓扑量子计算的重要研究对象[7，13，19]. 近十年多来， 在真实物理系统中实现MZM的理论工作取得突破[20-26]，寻找MZM的实验取得一系列重要进展[27-35]. 把拓扑超导作为可拓展的量子计算平台的装置也相继被提出[36-41]. 这些发展使得MZM正在成为最有潜力实现拓扑量子计算的平台[18，42].

2 凝聚态体系中的MZM

可承载MZM的最简单的系统是Kitaev[19]在2001年提出一维无自旋的p波超导链模型:

其中， 每个格点有一个无自旋的费米子态 cx，μ 是化学势， 最近邻跃迁强度为t，Δ 是最近邻p波超导配对强度，φ 是超导配对相位. 定义马约拉纳算 符 γx，A=i(c†xe-iφ/2-cxeiφ/2)，γx，B=c†xe-iφ/2+cxeiφ/2. 显然， 马约拉纳算符满足反对易关系γx，α=γx，α，{γα，x，γα′，x′}=2δαα′δx，x′， 其 中 α =A ，B . 在这组马约拉纳基矢下改写哈密顿量为:

这个系统看起来很简单， 却有两个不同的拓扑相.在这里我们以两组特殊的参数说明这两个不同拓扑 相 有 完 全 不 同 的 物 理. 在 μ ／=0，Δ =t=0 时，如图1(a)所示， (4)式只有第一项不为零， 哈密顿量 简化为

此时每一个格点上的两个马约拉纳算符通过化学势 μ 耦合在一起形成能量不为零的费米子激发 cx，系统不存在空间分离的MZM， 也没有零能的激发，系统具有唯一的非简并基态. 在 μ =0，Δ =t 时， 如图1(b)所示， (4)式中只有第二项不为零， 哈密顿量简化为

边界处MZM的存在是由系统体态的拓扑性质决定的. 对Kitaev链取周期边界条件并对哈密顿量作傅里叶变换， 得到体系在动量空间可以被写为 (令 Ck=(ck，c†-k) 为Nambu空间基矢算符)

图 1 Kitaev链的两个拓扑相 (a)和(b) Majorana基矢下的哈密顿量示意图. 图(a)中 μ ／=0，Δ=t=0 ， 只有第一项 (1+iγx，Bγx，A) 不为零， 每一个格点上的两个马约拉纳算符耦合在一起， 不存在空间分离的MZM. 图(b)中μ=0，Δ=t， 只 有 第 二 项 γx，Bγx+1，A 不 为 零， 相 邻 格 点的两个马约拉纳算符耦合在一起， 超导链的两端各自剩下一个MZM. (c) Δ =0 时的色散关系. (d) Kitaev链的拓扑相图. 当化学势穿过能带时， 体系处在拓扑相， 由(d) 中的橙色区域描述; 反之， 体系处在平庸相， 由(d)中的白色区域描述Fig. 1. Two topological phases of the Kiteaev chain. (a) Schematic illustration of the Hamiltonian in Majorana basis. In(a) μ ／=0，Δ =t=0， only the first term (1+iγx，Bγx，A)survives thus Majoranas couple at the same site leaving no seperate MZMs left. In (b) μ =0，Δ=t ， only the second term γx，Bγx+1，Asurvives thus Majoranas couple at adjacent sites， leaving one MZM at each end of the chain.(c) Energy dispersion for Δ =0 . (d) Topological phase diagram of Kitaev chain. When the chemical potential crosses the nomal spectrum the system is in topological phase， as described by the orange region in (d); otherwise the system is trivial， as described by the white region in (d).

体系拓扑不变量可由一维缠绕数刻画[8，43]. 为简便，可做变换使得 φ =0， 哈密顿量变为 H(k)=hzτz+hyτy. 一维缠绕数即为当k改变 2π 时向量 h(k)=(hy，hz) 在 y -z 平面内缠绕的圈数. 对于Kitaev链，拓扑数可以进一步由高对称动量点的性质给出:N=[sgn(hz(0))-sgn(hz(π))]/2，hz=-μ-tcosk.由此可知拓扑非平凡相的区域满足 | μ|＜t . 在任何不关闭能隙的连续变换下拓扑不变量保持不变， 要发生拓扑相变必须要关闭能隙. 对角化哈密顿量H (k) 得到体态的激发能谱为

容易得到， 只有在 | μ|=|t| 时， 系统能隙关闭， 也即发生拓扑相变. Kitaev链的拓扑相图由图1(d)给出. 在 | μ|＜|t| 时， 体系处于拓扑相; 在 | μ|＞|t| 时，体系处于平庸相. 需要注意， 在有限尺寸的非平庸相下， 体系两端各存在的MZM可以存在有限尺寸耦合， 但耦合的能量随着距离的增加而指数衰减Ef∝e-L. 一方面， 可以通过控制两个MZM的耦合初始化MZM组成的非局域的费米子态f， 另一方面， 当使得两个MZM离得足够远时， 该非局域费米子f的激发能为零， 系统具有二重简并的基态.下文中会提到， MZM带来的基态简并度和非局域性使其可以成为实现拓扑量子计算的平台.

尽管内秉的Kitaev超导链在真实系统中并未被发现， 与之等价的一维拓扑超导可通过有强自旋轨道耦合效应(spin-orbit coupling， SOC)的半导体纳米线放在s波超导体上并外加与自旋轨道方向垂直的塞曼场来实现. 该条件下纳米线中可以诱导出等效的无自旋p波超导[23，24]. 此时体系由以下哈 密顿量描述:

其中 ψ =(ψ↑，ψ↓)T， m是电子的有效质量，μ 是化学势，λ 是SOC的强度， h是外加塞曼场的强度，由外加磁场强度和材料的朗德因子g决定，Δ 是s波超导体通过近邻效应在纳米线中诱导出的超导强度. 如图2所示， 由于存在SOC， 自旋和动量关系被锁定. 再加上塞曼场在 k =0 处打开能隙， 使得当费米能量处于此能隙中时， 费米面上只有一个能带， 系统的低能物理等效于“无自旋”费米子. 此时由超导近邻效应在纳米线中诱导的费米点附近±k 动量之间s配对， 再写到自旋本征态上， 等效成为一维p波超导， 因而此时体系的低能物理由Kitaev模型描述. (9)式的哈密顿量和Kitaev链的拓扑性质由同一个拓扑不变量描述. 通过对(9)式进行能谱分析， 可以得到体系的拓扑非平凡相区域为 h2＞Δ2+μ2， 由图2(c)给出[44-46]. 很快， Delft实验组[27]在2012年首先报道， 把InSb纳米线放在超导体的表面， 在一定外加磁场的条件下初步观测到了由末端MZM诱导的零偏压电导峰(Zero-Bias Conductance Peak， ZBCP).类似观测也被其他实验组看到[47-50]. 零压隧穿电导峰是MZM的间接证据之一. 特别地， 由于MZM的自共轭特性， 在零温下由MZM诱导的隧穿谱ZBCP应该是高度为 2 e2/h 的量子化电导[51-54]， 并在有限温度下会降低[55]. 此外， Yazdani研究组[28]在有SOC的铅衬底上生长的铁原子链两端也测到了ZBCP， 分析可能是MZM存在的迹象. 把图2(a)中的一维SOC纳米线换成二维拓扑绝缘体[56，57]的边缘态形成二维拓扑绝缘体边缘/s波超导体异质结， 在外加铁磁体的情况下同样可以实现等效Kitaev模型[20，58].

图 2 在一维SOC纳米线中实现马约拉纳零模激发[8] (a)装置简图； (b)一维SOC纳米线的色散关系. 在无外加塞曼场时， 系统具有时间反演对称性， 体系有偶数个费米面(红色和蓝色曲线); 在外加塞曼场时， 时间反演对称性被破坏，k =0 处被打开能隙(黑色曲线). 若化学势位于能隙中间， 该体系只有一个费米面， 其低能哈密顿量等效为Kitaev链. (c)体系的拓扑相图， 其中相边界由 h2=Δ2+μ2给出， 橙色(白色)表示拓扑(平庸)区Fig. 2. Realizing MZMs in a 1D SOC nanowire[8]. (a) Sketch of the basic setup. (b) Energy dispersion for the 1D SOC nanowire. When Zeeman field is absent， the system is timereversal symmetric(TRS) and possesses even number of Fermi surfaces (red and blue curves); when Zeeman field is introduced， TRS is broken and a gap is opened at k=0(black curves). Given that the chemical potential lies within gap， the system possesses only one Fermi surface， and the low-energy Hamilitonian is equivalent to that of the Kitaev chain. (c) Topological phase diagram of the system with the phase boundary given by h2=Δ2+μ2 . Orange(white) denotes topological (trivial) region.

此外， Pientka等[26]在2017年的工作提出一种新型的平面约瑟夫森结也产生准一维的MZM.这种约瑟夫森结是在强SOC的二维电子气体系上放置两个超导体形成约瑟夫森结， 在结区无超导处形成一维的Andreev束缚态能带. 在结区加上面内的磁场， 通过调节磁场和两端超导的相位差可以调节结区Andreev束缚态能带的拓扑与平庸的拓扑相变， 从而调控结区的垂直方向两端点的MZM出现和消失. 这个装置的优点一是通过调控结两端的超导相位差， 增大MZM出现的磁场范围， 降低实验对磁场强度的要求; 二是装置中磁场主要加在非超导区， 不影响会抑制超导序. 根据这个方案， 两个实验组分别在HgTe量子阱中的二维电子气和高迁移率的InAs电子气上面覆盖两块超导铝形成约瑟夫森结上测到了验证MZM存在的零偏压电子峰[33，34]. 基于该装置提出了MZM编织等其他性质[41，59，60].

MZM也可以存在于二维无自旋 px+ipy拓扑超导体的涡旋中心[13，61]. 二维无自旋 px+ipy手性拓 扑超导体的哈密顿量可写为

二维手性拓扑超导的拓扑由第一陈数刻画 (即二维缠绕数， 为动量空间往布洛赫球面的映射)[8，43].二维的无自旋内秉 px+ipy拓扑超导体同样难以在实际材料中实现， 但可以通过s波超导和二维自旋轨道耦合电子气做成人工异质结构而通过近邻效应等效实现. 其中， 将三维拓扑绝缘体的表面和s波超导结合， 在三维拓扑绝缘体由狄拉克锥描述的表面态[62]中引入超导序而形成等效的二维px+ipy拓扑超导[20，63-66]. 类似地， 将二维半导体量子阱和s波超导形成异质结， 并在外加塞曼场于k=0 处打开能隙解除费米双重性(Fermion doubling)， 使得费米面只有一个能带时， 系统将形成等效的 px+ipy二维拓扑超导体[21，22，67，68]. 相对复杂的异质结对实验工艺提出了挑战. 2012年， 贾金锋等[69]把层状三维拓扑绝缘体 B i2Te3薄膜生长在NbSe2超导体衬底表面上， 第一次实现了超导相和拓扑相的共存. 在后续的研究中， 研究人员在体系的涡旋中心观察到了强度受自旋极化影响的零偏压峰， 实现了基于自旋探测的MZM观测[29，70-72].相比于通过异质结构实现二维拓扑超导， 更近的理论预言铁基超导体的表面态可出现天然的二维p波拓扑超导态[73，74]. 基于这个预测， 丁洪实验组[31，32，75]在2018年取得突破， 他们在 F eTe0.55Se0.45的表面观察到了超导能隙被打开， 并后续观察到了超导涡旋中心的ZBCP. 这一结果进一步被其他实验组证实[76-79]. 基于铁基超导实现二维拓扑超导有两个特点: 一是不需要生长复杂的异质结， 这降低了工艺难度并避免生长晶体过程中的杂质等由异质结生长带来的问题; 二是体态超导直接诱导表面态打开能隙， 加上铁基超导的高的转变温度能诱导出大的拓扑超导能隙.

值得一提的是， 在涡旋中实现非阿贝尔MZM并非必须实现二维拓扑超导. 近期的理论预言， 除了在拓扑超导体的涡旋中， 在体态平庸的二维超导体中的拓扑缺陷中心也可能出现MZM[80，81]. 事实上， 涡旋中的MZM并非由二维陈数刻画， 而是由定义在更高的三维空间中(二维实空间加一维涡旋相位参数空间)的Hopf拓扑数[82]或陈-西蒙斯数决定[81]. 在特定的条件下， 体态拓扑平庸的二维超导体中的涡旋中心， 也可存在满足非阿贝尔统计的MZM[81]. 这项理论结果表明MZM可能存在于更广泛的拓扑或平凡超导中. 将二维拓展到三维超导半金属中， 理论预言了第二陈数保护并满足三维圈非阿贝尔统计的手性MZM[83].

除了涡旋中的MZM， 一维手性马约拉纳边缘态存在于二维手性拓扑超导的边界[84-87]， 但实验观测存在更大的挑战. 近期报道在基于砷化镓衬底生长的单层 (Cr0.12Bi0.26Sb0.62)2Te3合金量子反常霍尔绝缘体和超导体 N b 形成的异质结结构中初步测到了由手性马约拉纳边缘态产生的半整数量子化两端隧穿电导平台[30]. 在另一基于钛酸锶( S rTiO3)衬底生长的三明治结构量子反常霍尔绝缘体和超导异质结构中没有观测到半整数量子化电导平台[88]. 应当注意， 在量子反常霍尔绝缘体/s波超导的实现中， 二维手性拓扑超导只存在于量子反常霍尔相的拓扑转变区域附近. 因为在此区域附近， 量子反常霍尔材料体内能带处于接近能隙关闭状态， 使得费米能可能穿过单一费米面， 从而通过s波超导近邻效应诱导出二维手性拓扑超导. 容易判定， 这个转变区域必须对应于材料体内磁化能量和材料上下表面耦合能可竞争比较的参数区域. 因此， 实现手性马约拉纳边缘态需要对量子反常霍尔绝缘材料体内磁化强度具有相对平缓操控的能力.并不是基于任意量子反常霍尔材料都能普遍实现手性马约拉纳边缘态. 实验观测仍需更多进步.

到目前为止， 观测MZM的绝大多数实验仍然是通过MZM诱导的隧穿电导的零偏压峰进行[27，33，35，48，89-93]. 另一方面， MZM诱导的 4 π 周期的分数约瑟夫森效应是另一个重要观测现象[7]. 在拓扑超导体形成的约瑟夫森结中， 两个MZM耦合形成Andreev束缚态， 约瑟夫森电流可以通过单电子隧穿产生. 由于单电子相位只是超导相位的一半， 拓扑超导体约瑟夫森结中的隧穿电流和超导相位形成 4π 的周期特征. 目前还没有实验能测到直流分数约瑟夫森电流， 但是已经有报道指出在InAs纳米线结构中实验观测到了交流约瑟夫森效应[94-96]. 本文接下来将会重点介绍MZM在拓扑量子计算的应用， 包括非阿贝尔统计及其在量子计算中的应用， MZM编织操作的实现方案， 以及量子比特的读取方案等.

3 MZM非阿贝尔统计及其在量子计算的应用

3.1 MZM非阿贝尔统计

全同粒子的统计性质是量子多体物理的基本性质. 熟知的量子统计包括费米和玻色统计， 分别对应在两个粒子交换后系统波函数改变或不改变符号. 更一般情形， 对于阿贝尔任意子， 交换前后的波函数满足 | Φf〉=eiφ|Φi〉，即相差一般化相位 φ .由于相应演化算符由阿贝尔的 U (1) 相 位( eiφ)描述， 这类统计性质被称为阿贝尔统计. 当多体系统存在由准粒子带来的简并基态， 比如两个分离的MZM导致 |0〉 和 |1〉 的简并基态， 准粒子间的交换操作不仅带来系统波函数的相位改变， 且会发生不同简并基态之间的演化. 一般地， 有 |Φf〉=Ufi|Φi〉.这样准粒子交换操作由在简并基态空间内的矩阵Ufi描述. 由于矩阵不可对易， 此类统计性质被称为非阿贝尔统计， 相应的全同粒子称为非阿贝尔任意子[15]. 任意子的空间-时间轨迹可以由“世界线”描述， 如图3(b)所示. 在任意子的交换的过程中， 它们的“世界线”互相缠绕如同编辫子， 每一种缠绕方式对应一种交换， 因此任意子的交换又称为编织操作. 编织操作必须满足绝热条件， 并在单次编织完成后， 系统的哈密顿量回到初始. 从绝热定理可知，要使绝热交换粒子前后系统的量子态发生变化， 系统必定具有简并的基态. 绝热交换两个任意子， 会导致在这个基态子空间中的一个幺正变换， 数学表示为[97]:

其中 e-ℏ1dtE(t)是动力学相位; B0(t)| 是贝里相位的推广， 称为贝里矩阵. 考虑2N个空间分离足够远的MZM γ1，···，2N， 把MZM组合成N个 狄 拉克费米子

图 3 MZM交换示意图 (a)含有4个空间分离足够远的MZM γ1，2，3，4的 体 系， 其 中 γ1和 γ2、 γ3和 γ4 分 别 构 成两 个 费 米 子. γ2和 γ3顺 时 针 交 换 一 次，γ2会 跨 越 γ3 所 在涡旋 的 相 位 割线获得 一 个 负 号， 而 γ3并 未 跨越 γ2 所 在 涡旋的相位割线不获得负号. 因此结果是 γ2→-γ3，γ3→γ2 . (b)描述4个MZM的时间-空间(x; t)轨迹的世界线. γ2和 γ3 被 编织 一 次， 因 此 它 们的 世 界 线 缠 绕 一 次.系统由初态 | Φi〉=|00〉 演化到 | Φf〉=|00〉+ i|11〉Fig. 3. Skecth of a MZM braiding operation. (a) A system consists of 4 MZMs far enough apart， with γ1and γ2，γ3 and γ4forming 2 fermions f1and f2. γ2and γ3are braided once clockwise. γ2 crosses the branch cut of the votex hosting γ3 and gains a minus sign， while γ3 doesn't cross the branch cut of the votex hosting γ2 and doesn't gain a minus sign. Hence the result is given by γ2 →—γ3， γ3 → γ2. (b) Worldlines in a space-time (x; t) diagram，describing four MZMs. γ2and γ3are braided once， hence their worldlines winds each other once. The initial state|Φi〉=|00〉 evolves into | Φf〉=|00〉+i|11〉 .

其中每个费米子 fj的能量为零 ϵj=0 . 因此这N个费 米 子 构 成 2N重 简 并 的 基 态 空 间 |Φi〉=|n1，···，nN〉， 其 中 nj=0，1 是 占 据 数 n ˆj=fj†fj的本征值. 由于超导满足费米宇称(Fermion parity)守恒， 编织操作只能使得量子态在 2N-1-维度的简并空间演化. 由MZM导致系统的基态简并， 这满足非阿贝尔统计的必要条件.

MZM的非阿贝尔统计性质最早在二维系统中被提出来. Read和Green[13]在研究 5 /2 -分数量子霍尔效应的Moore-Read Pfaffian态时发现它们adf与二维 p +ip 超导系统有着对应关系， 前者中的非阿贝尔统计在二维 p +ip 超导的涡旋中心的零能激发将同样成立， 这一点由Ivanov[14]给出直接证明. 此外， Teo和Fu[98-100]提出三维系统中满足特定条件的拓扑缺陷形成的MZM激发也可以满足非阿贝尔统计. 2011年， Alicea等[101]通过构造一维网络证明了准一维拓扑超导系统中的MZM同样满足非阿贝尔统计性质， 这个工作开启了以一维MZM系统作为实现主要拓扑量子计算研究平台的广泛研究.

为方便描述， 我们以二维 p +ip 拓扑超导中的MZM介绍非阿贝尔统计性质. 如图3(a)， 考虑二维拓扑超导中的4个相互离得足够远的涡旋， 每一个涡旋中心有一个MZM， 其中 γ1和 γ2按(12)式规则构成费米子 f1，γ3和 γ4构成 f2. 超导相位沿着每个涡旋绕一周变化 2π . 一种方便的选择是让超导相位为关于割线的阶跃函数， 在割线以外的区域， 超导相位恒定， 越过割线时， 超导相位跃变 2π .这条割线可以是任意一条从涡旋中心延伸到无穷远处的线(或延伸到一个反涡旋的中心). 由于超导序参量是2个电子的乘积， 而MZM 是电子空穴的线性叠加， 所以MZM 的相位是超导相位 φ 的一半， 即 γ ∝(feiφ/2+ f†e-iφ/2)[14]. 现 在 让 γ2和 γ3所在的涡旋顺时针交换一次， 如图3(a)所示，γ2所在涡旋会跨越 γ3所在涡旋的相位割线， 超导相位改变2π， 则 γ2相应获得 π 相位， 即变负号， 而 γ3并未跨越 γ2所在涡旋的相位割线不获得负号. 因此 γ2和γ3顺时针交换一次结果是

系统态在编织过程的演化用世界线的语言来描述，如图3(b). 编织导致不同简并态之间的转化， 表明编织操作属于矩阵操作， 因而对应为非阿贝尔编织. 需要注意的是， 编织结果与MZM具体走过的路径无关， 只由路径的拓扑性质决定， 是受拓扑保护的. 如果初始态是， 则编织后的结果为. 注意这里得到的是单次编织操作的结果. 一次完整的编织包含两次上述的单次编织， 其结果导致和复费米子空间的变化为.类似地可以得出， 顺时针交换 γ1和 γ2操作算符为. 容易知道， 编织操作 B12不改变 f1和 f2的宇称， 因而粒子数态表象中为对角化 矩 阵， 并 可 得 到 B12|00〉=|00〉 和 B12|11〉=-i|11〉 . 因此， 在偶费米宇称子空间 (|00〉，|11〉) 中编织操作 B23和 B12可以写成如下矩阵:

为方便， 矩阵 B12乘了一个整体相位因子. 这不改变非阿贝尔统计物理. 显然， 贝里矩阵满足

这重复验证了MZM之间的编织操作满足非阿贝尔特征. 事实上， 2N个MZM的编织操作在基态子空间的贝里矩阵对应着编织群 B2N的一个投影表示[14]. 通过编织实现对量子态的操控， 可用于量子计算.

需要注意， 关于MZM的非阿贝尔统计常存在误解. MZM的非阿贝尔统计核心在于: 我们在复费米子态空间看MZM的交换性质. 事实上， 不同MZM的算符满足费米子对易关系 {γi，γj}=0( i ／=j ). 但是， 单个MZM没有完好定义的复费米子态空间， 而有量子维度. 必须由MZM两两组合出复费米子态空间. 这导致即使不同MZM算符满足费米子反对易关系， 在复费米子态空间来看， 其编织相当于对非局域的量子态进行“分数化”操控， 从而导致量子态在简并空间的演化. 由此可知， MZM的非阿贝尔统计是其在复费米子态空间的特殊量子维度所致. 这导致我们在复费米子态空间看MZM的编织出现非阿贝尔特征. 由此得知，零能量并非非阿贝尔统计的本质因素， 即使非零能(但简并的)马约拉纳费米子， 只要能够以某种方式实现编织， 也必然满足非阿贝尔统计行为. 但是MZM在实现单个MZM之间的编织上， 显然要自然和方便， 并受拓扑能隙的保护.

3.2 量子计算基础

1982年， 物理学家费曼[102]提出了按照量子力学规律工作的计算机的概念， 这是最早的量子计算机的思想. 在讲量子计算机以前， 首先回答: 什么是计算机? 计算机是一种处理信息的工具， 所有的计算系统都要有两大基本功能: 信息的存储和操控. 经典二进制计算机是把信息存储在二进制数组里， 通过对输入二进制数组实现特定的布尔运算得到新的二进制数组， 从而实现对信息的操控. 而量子计算机则是把信息存储在量子态里， 通过对量子态在特定条件下按照量子力学演化得到新的量子态， 从而实现对量子信息的操控. 显然， 量子物理具有很多经典物理所不具有的特点， 如量子态的叠加原理， 量子纠缠和量子相干性等性质. 通过充分利用这些奇特的量子特性制造出的量子计算机将有可能在一些问题的计算上比经典计算机有很大的优势. 目前已被证明的大大优于经典计算机的量子算法主要有三大类[16]: 用于做指数分解的Shor算法[103]， Grover的搜索算法[104]， 和模拟真实的物理系统[102，105].

经典计算的最小信息单位是比特， 每个比特可以是0或1两个状态. 量子计算的最小信息单位则包含两个量子态， 记为 | 0〉，|1〉 ， 称为量子比特. 和经典比特不同的是， 量子态是概率性的、满足叠加原理的， 也就是说， 每一个比特 | ψ〉 都可以是 |0〉 和 |1〉的叠加态， 即

因此一个寄存器能同时存储 2N个信息状态. 从一个包含N个量子比特的量子态函数出发， 做一次计算实际上可相当于做了 2N次计算， 这被称为量子并行计算[106，107]. 量子并行计算是人们认为量子计算有可能大大超越经典计算的主要原因之一. 除了态的叠加原理外， 量子力学中的量子态纠缠和干涉都是量子算法区别于经典算法的重要性质.

经典计算实现信息操控的是逻辑门[107]， 每一个经典逻辑门对应实现一个布尔函数

每一个经典门由一个 2m×2n的真值表描述. 而量子计算的操控是实现量子门， 每一个量子门由一个2N×2N的幺正矩阵描述， 每一个量子门对应一个演化算符，

与经典量子门对应的量子门是指按照经典门的真值表实现量子态矢量转化的量子门. 比如， 经典的非门是一个单比特门， 由一个2 × 2的真值表表示: 0→1，1→0， 对应的量子门使得 |0〉→|1〉，|1〉→， 矩阵形式为

实现通用的量子计算要求， 对于一个N比特的存储器， 必须要实现2N维希尔伯特空间的任意幺正变换 U (2N) . 正如在经典计算中实现通用计算只需实现一个完备的逻辑门集， 其他的逻辑门都可以由这个基本的逻辑门集里的门组合给出， 实现通用量子计算同样只需要实现一个通用量子门集， 其他的量子门都可以由这个集合里的量子门通过有限次组合得到. 下一节讨论基于MZM的拓扑量子计算中将给出具体形式.

DiVincenzo[108]在2000年总结出了实现量子计算的物理系统必须满足的五大要求: 1) 具有良好特征的量子寄存器的可扩展物理系统; 2) 在简单的基准状态下初始化量子寄存器的能力; 3) 足够长的相关的相干时间， 远长于门控操作时间;4) 可以实现一组“通用”的量子门; 5) 具有系统的测量能力. 然而， 量子系统都是“脆弱”的， 易于受环境干扰而且精确操控难度极高， 所以可扩展的量子计算机的研究碰到了很大的障碍. 目前而言， 对量子计算机的探求中遇到的两个最大的“绊脚石”是: 环境噪声导致的量子退相干问题和量子操作精度有限导致的量子门误差. 这两个困难使DiVincenzo给出的五大要求中的第三和第四个要求不能满足.拓扑量子计算正是为了从根本上解决这两个问题而提出. 在受拓扑保护的非阿贝尔任意子量子系统中以非局域的形式存储量子信息， 并通过任意子的拓扑交换来实现量子信息的操控， 从硬件上可以实现抗退相干效应的量子计算[16]. 为此， 拓扑超导中的MZM系统是近十年来被研究最多的可用来作为拓扑量子计算平台的非阿贝尔任意子量子系统.

3.3 基于MZM实现量子计算

由于量子计算的优越性， 量子计算机被广泛研究， 但是由于量子态非常易于受环境干扰产生退相干效应和量子操作精度有限导致的量子门误差， 人们至今仍没有实现通用量子计算机. 而MZM形成的Dirac复费米子态具有非局域特性. 这使得由MZM组成的量子态受到拓扑保护， 不会被局域扰动所改变; 另一方面， MZM的非阿贝尔统计性质又保证了通过交换MZM的位置可以在保持量子态非局域性的前提下高精度地实现量子态操控， 这些性质使得MZM成为实现量子计算的重要物理平台[109].

前文中提到， 含有 2 N+2 个MZM的系统具有 2N+1重简并的基态， 但如前所述， 由于超导系统中的费米子宇称守恒(即费米子奇偶数守恒)， 这一系统实际上只能形成N个独立的量子比特. 图4中我们以密集编码方式的2比特的量子计算为例，说明用MZM怎样实现信息的存储和操控[110].2比 特 系 统 由6个MZM构 成 γ1，2，3，4，5，6， 按 照构成3 个费米子， 对应的占据数可以标记基态 | n1n2n3〉 . 现选奇数基态子空间作为计 算空间， 选取我们的计算基矢为

图 4 通过MZM编织操作实现量子门[110] (a)-(d) 分别为实现单比特量子门H-门、Z-门、双比特门CNOT门和双比特门CZ门的编织操作示意图Fig. 4. Quantum gates realized by MZM braiding operations[110]. (a)-(d) The elementary braids corresponding to the single-qubit gates H-， Z-gates on the first qubit as well as the 2-qubit CNOT and CZ gates.

类似可以给出其他两个MZM之间的编织算符. 通过组合不同的编织操作可以实现不同功能的量子门. 其中常见的一些量子门的实现可以由图4中的“世界线”给出.

理论证明， 三个单比特量子门， 即H-门， Z-门和T-门， 加上一个纠缠的双比特门， 如CNOT-门，可以组成完备的通用量子门集[42，111，112]. 基于这4个基本的量子门， 其他的量子门都可以由这4个门的有限次组合近似得到， 进而实现通用量子计算. 这几个量子门表示如下

其中， H-门， Z-门和CNOT-门都可以在图4 中的世界线给出. 然而， T-门并不能由受拓扑保护MZM的编织组合给出. 因此， 单纯通过MZM的非阿贝尔编织还不能实现通用量子计算. 一种补充方案是通过非拓扑保护的方式实现T-门 ， 从而结合拓扑和非拓扑的手段相结合来实现通用量子计算[113-117].

4 MZM编织操作的物理实现

4.1 一维T-型结

在严格的一维体系中， MZM不可能绕开彼此完成编织过程. 因此， Alicea等[101]提出用T型几何结构(图5)的拓扑超导体 (T-junction) 实现MZM的编织. 通过调节门(gate)电压移动拓扑相和平庸相之间的畴壁， MZM可以在T-型结上移动， 借助竖直部分避免了MZM在编织过程中相遇. 利用 f =(γ1+iγ2)/2 定义简并的多粒子基态波函数

原则上可以通过计算贝里相位

验证编织结果的非阿贝尔统计性质. 详细计算可以证明两个态之间的贝里相位差为 π /2 . 具体可参考文献[101]. 此处通过如下直观定性分析快速得到编织结果.

从图5的编织过程可以看到， 基于T-型结实现的一次编织实际上相当于 γ1和 γ2相对于结点逆时针转圈. 与二维系统不同， 一维链中MZM并不伴随涡旋， 因此没有涡旋割线描述. 但是在T-型结的结点处， 由三端结点往纳米线的三个方向(左、右、下)， 以及虚拟的往上方向， 超导配对参数 Δp必须存在等效的相位涡旋特征， 方式如下. 考虑横向纳米线中为px-波超导， 则从结点往右(图5(a)箭头方向)和往左方向配对参量分别正负符号(奇宇称px-波配对往左往右反号)， 因此对应为 0 和 π相位. 另一方面， 竖直方向的纳米线中， 超导配对参数必须有ipy成分， 即必须和横向纳米线中配对参数有相位差. 否则， 如果竖直方向配对参数完全为实函数， 比如为py-波超导， 则这段纳米线和横向的左端， 或者右端部分必将形成 π -约瑟夫森结， 从而在编织过程中在结点出现额外的MZM， 导致破坏编织. 因此， 最简单且满足条件的构型是在竖直方向纳米线中(朝上方向， 即图5(b)竖直箭头方向)形成ipy配对. 这样一来， 从结点往上和往下竖直方向的配对相位分别为 π /2 和 3 π/2 . 由此， 绕结点逆时针一圈等效改变 2π 相位. 于是可以假想从结点出引出一割线. MZM在编织中则会跨过割线，获得 π 相位， 从而得到二维中涡旋编织一样的结果.其中图5 描述的单次编织后， MZM的变换满足γ1→γ2，γ2→-γ1， 对应非阿贝尔统计.

图 5 T-型结进行编织操作以及键盘门的操作方式[101]. T-型结由两条水平链和一条竖直链组成， 深蓝色部分为拓扑相， 浅蓝色部分为平凡相， 超导体上向上和向右的箭头分别代表 φ =π/2 和 φ =0， 用向左和向下的箭头代表 φ=π和 φ =3π/2 . T-型结外的箭头代表MZM运动的方向. 黑色和灰色方块分别代表局域门的关闭和打开状态， 对应平凡相和拓扑相. (a)—(d) 给出 γ1 先运动到竖直链上， 然后γ2从水平链右端运动到左端， 最后 γ1 运动到右端的过程，过程结束后箭头反向. (e) 中局域门的存在可以保证在不关闭能隙的前提下逐渐移动MZMFig. 5. A T-junction allows for braiding process and the keyboard gates[101]. The T-junction consists of two horizontal segments and one vertical segment. Dark blue segments are in topological phase， and light blue lines trivial phase.φ=0 or π /2 is represented with rightward or upward pointing arrows， while φ =π or 3 π/2 represents the leftward or downward pointing arrows. MZMs are transported according to the arrows around the T-junction. Black and gray blocks denote different states of tunable gates in accordance with trivial and topological phases. (a)—(d) sketch the process which γ1is transported to vertical line firstly，then γ2travels from the right end to the left end and at lasTγ1is transported to the right end. After this process，the arrow points to the opposite direction. Local gates in(e) ensure that the MZMs can be manipulated gradually without closing the gap.

另一种等价描述方式即将T-型结看作三段Kitaev链并如下选择参数. 如图5所示， 用向上和向右的箭头分别代表 φ =π/2 和 φ =0 的部分， 用向左和向下的箭头代表 φ =π 和 φ =2π/2 的部分. 在编织过程中， 为了保证结点处只有一个零模， 拓扑非平凡区域箭头的指向必须保持一致， 并且在结点处， 箭头指向必须一进一出. 这样在编织前后， 横向整条链的箭头必将反向. 此时对横向链所有的费米子算符做规范变换

可将编织后的哈密顿量变回到编织前的形式， 由于f†=(γ1-iγ2)/2， 该规范变换等价于 γ1→γ2，γ2→-γ1， 同 样 得 到 非 阿 贝 尔 统 计. 很 容 易 将 两 个MZM的情形拓展到四个以及更多MZM的情况.基于T-型结原则上可以进一步构造纳米线网络，实现MZM的编织和拓扑量子计算.

实现MZM在一维拓扑超导体上的运动可采用如图5(e)所示的键盘门(keyboard gate)的方式，即局域地调节纳米线的化学势使得拓扑相和平庸相之间的畴壁逐渐移动. 但是该方案在实验上难以操作而且不能保证体系绝热演化[118，119]， 后续基于T-型结以及改进的编织方案更多考虑利用调节MZM之间的相互作用， 如隧穿效应或库仑相互作用， 实现MZM位置的交换[36，37，120-124].

4.2 基于测量的编织

通过T-型结或类似装置实现MZM编织的最大困难在于保证移动MZM过程中的绝热条件. 尤其在T型结点处， 实际构型复杂， 容易出现非MZM的低能模式. 在调节MZM通过结点时，MZM和其他复杂模式的耦合容易破坏绝热条件的成立， 从而影响MZM的编织操作. 避免这类困难的一种解决方案是基于测量实现MZM的编织.这类方案的主要好处是不需要真实通过结点移动MZM， 而是通过对MZM做一定顺序的测量，导致对MZM的teleportation从而实现编织. 但另一方面， 这种方法需要使用辅助MZM(ancilla Majorana). 这会带来对测量的一系列额外要求.基于测量的拓扑量子计算一般理论参见文献[125].具体方案实现介绍如下.

首先考虑两个MZM情形， 即对 γ1和 γ2进行编织. 为此， 需要引入两个辅助MZM， 计为 χ1和χ2. 注意到引入两个辅助MZM会将需要操作的量子态空间增大一倍， 因此编织过程必须保证在固定的子空间进行. 编织过程分为如下步骤. 第一步，测 量 两 个 辅 助MZM χ1和 χ2， 并 设 想 测 量 值 为iχ1χ2=+1 . 则下面的测量中， 为确保编织过程系统在固定的子希尔伯特空间进行， 每次测量均选择相同的测量值. 第二步， 测量 γ1和 χ1， 选择测量值iγ1χ1=+1 . 第 三 步， 测 量 i γ2和 χ1， 同 样 选 择iγ2χ1=+1 测量值. 第四步， 再次测量 χ1和 χ2并选择测量值 i χ1χ2=+1 . 完成这些步骤， 即实现了MZM γ1和 γ2的一次编织操作. 这些操作过程没有移动MZM， 但测量操作会产生量子态隐形传输.在第二步结束(完成测量 γ1和 χ1)后， 把初始 γ1和γ2的拓扑量子比特转移到 γ2和 χ2. 第三步操作则导致进一步转移到 χ2和 γ1. 第四步转移回 γ1和 γ2的拓扑量子比特， 完成编织. 如果把这个过程和图4中的T-型结对比， 其第二步相当于把左边的γ1移动到垂直纳米线上， 第三步把右边的 γ2移动到左端， 第四步把 γ1从垂直纳米线移到右端.

上述两个MZM编织的方案很容易应用到四个MZM γ1，2，3，4的编织中. 这时同样需要两个辅助MZM χ1，2. 通过和上面相同的操作过程， 对 γ2，3和χ1，2进 行 四 步 测 量， 从 而 完 成 对 γ2，3的 编 织.MZM的测量通过调节MZM之间耦合强度实现[126].基于测量编织方案进一步构造量子线路在近期也有相应研究[39].

可以看到， 基于测量的编织方案具有不必移动MZM的优势， 但同时必须引入辅助MZM. 这导致拓扑量子态希尔伯特空间增大. 因此要保证准确编织， 必须对测量值做选择， 确保在固定的子空间进行. 这带来两类潜在挑战. 首先， 量子态可出现向期望的子空间以外的态演化， 造成信息损失(information leakage); 另外， 测量中总可能出现不是待选择值的其他结果， 比如在第二步中测量得到iγ1χ1=-1 . 这时需要重复第一步使系统回到初态，对应 i χ1χ2=+1 ， 然后再重做第二步测量， 检测是否得到期望结果 i γ1χ1=+1 . 如仍为非待选择值，则再次重复前两步直到得到期望结果， 然后做进一步测量. 显然这样完成一次编织， 实际可能需要进行 大量的测量， 从而带来更多的复杂性.

4.3 基于自旋操作实现编织

以上介绍的两类MZM编织操作分别要求移动MZM位置和操控不同MZM间的耦合来实现.前者要求使MZM跨过T型结点， 后者要求引入额外MZM、扩大简并空间而导致操作的复杂性.因此这些在实验工艺上均存在大的难度， 使得MZM编织的实验实现仍是巨大挑战. 而拓扑量子计算理论指出， 对于有内禀空间的任意子， 描述它们演化的世界线应被扩展为“世界带”. 在世界带的语言中， 交换两个任意子的位置一次在拓扑上等价于两个任意子的世界带各自扭转半圈[127]. 根据这个思想， 最近提出， 通过转动MZM自旋自由度来扭转MZM世界带可以实现MZM的等效编织， 从而实现既不移动MZM位置也不引入额外辅助MZM以及它们的耦合的编织方案[128]. 该方案应用到实验上常见的超导/二维拓扑绝缘体/铁磁绝缘体(superconductor/2D topological insulator/ferromagnetic insulator， SC/2 DTI/FI)的异质结系统中， 提出了一个可以操控MZM自旋进而实现MZM编织的实验装置.

如图6(b)所示， SC区形成一维拓扑超导，FI区为平庸绝缘体， 两区的交界出现孤立的MZM， 此时MZM的自旋方向可由平庸区的铁磁体的磁化方向操控. 在这个实验装置中， 编织两个由铁磁区连接的MZM一次(单次编织操作)需要绝热地转动相应铁磁体的铁磁方向半圈， 即转过 π角， 编织两个MZM两次(一个完整的编织操作)需要绝热地转动铁磁体的铁磁方向一整圈， 即转过 2π 角. 编织结果不受系统中的局域杂质和磁化不规则转动路径等因素的影响， 具有拓扑稳定性. 这个稳定性的物理原因在于编织操作引起的奇偶性变化和转动磁场引起的拓扑电荷泵之间的联系. 这个方案的优点一是不需要辅助MZM， 这使得编织操作只涉及最少的MZM; 二是可以在保持基态简并度(无需测量操作)同时不移动MZM位置就等效地实现MZM编织操作. 但是这个装置不能实现超导两端的MZM的编织， 在做到编织任意两个MZM上存在限制.

图 6 通过转动铁磁体的磁化方向实现铁磁体两端的MZM[128] (a)交换世界带两次在拓扑上等价于世界带各自扭转 2π . 箭头表示MZM的自旋; 世界带的红蓝边界标记MZM的内部自由度的时间演化. (b) SC/QSH/FI异质结装置中的MZM. 黄色(红色)箭头表示MZM的局域自旋方向(铁 磁 磁 化 方 向). 下 方 的 红 色 箭 头 转 π 角，γ1和 γ2 编 织一次; 转 2 π 角，γ1和 γ2 编织两次. 逆方向转动代 表逆交换.Fig. 6. Braiding operation via winding FI magnetization[128].(a) The monodromy operator can be realized by either braiding two MZMs or twisting each worldribbons by 2 π .The arrows indicate the MZM spin. The blue and red edges of the ribbon denote the evolution of internal degree of freedom. (b) MZMs in the SC/QSH/FI hybrid system. The yellow (red) arrows represent the directions of local spin polarizations for MZMs (FI magnetization). Winding the red arrow at the bottom by π，γ1and γ2are braided once; by 2π they are braided twice. A reverse rotation leads to an inverse braiding operation.

5 MZM量子比特的读取

零压电导峰(ZBCP)是MZM的一个间接证据， 但是只是单个MZM的信息， 并不能读出量子比特. 量子比特是两个MZM组合成的Dirac复费米子的占据状态， 要读出量子比特， 必须把两个MZM耦合到一起， 考虑两个耦合的MZM作为一个整体的性质. 前比较广泛接受的读出MZM量子比特的方法主要有三种: 分数约瑟夫森效应， 库仑阻塞(Coulomb blockade)方法， 及干涉方法.

5.1 分数约瑟夫森效应

MZM的分数约瑟夫森效应最早是Kitaev[7]在讨论Kitaev超导链模型的时候提出来的， 这是存在MZM的奇宇称拓扑超导系统的一个重要特性. 当两个MZM通过约瑟夫森结耦合在一起时，通过测量穿过约瑟夫森结的零偏压的电流响应， 也就是直流(DC)约瑟夫森电流， 我们可以推出这两个MZM组成的费米子的占据状态， 也就是读出MZM量子比特.

图7(a)所示为利用约瑟夫森效应的量子比特读出装置的简示图. 约瑟夫森结两端的超导体由绝缘区微弱地耦合起来， 当结两端各自存在一个MZM时， 流过结区的无阻超导电流包含两部分I=I2e+Ie. 其中第一项是库珀对隧穿贡献的传统约瑟夫森电流 I2e， 第二项是两个MZM耦合在一起贡献的单电子隧穿电流 Ie. 正是通过测量 Ie[101]，可以给出MZM耦合的量子比特的信息. 下面讨论Ie的形成机制和特性.

假设形成约瑟夫森结的两端拓扑超导相位分别为 φL和 φR. 结区存在的两个零能的MZM会通过量子隧穿耦合在一起， 形成能量非零的费米子态. 此时约瑟夫森结的低能有效哈密顿量为

其中，φ =φL-φR是结两端超导相位差，λ 是两个MZM耦合的强度，γ1，2是马约拉纳算符， f是 γ1和γ2耦合成的费米子 f =(γ1-iγ2)/2 . 保持系统的粒子数奇偶守恒时， f的占据数不变， 流过结区的电 流

其中 n0为系统初态的费米子f的占据数. 由电流的表达式可以看出， 和传统的库珀对隧穿形成的超导电流 I2e不同的是，I2e是关于关于结两端超导相位差 φ 以 2π 为周期变化的， 而 Ie是关于结两端超导相位差 φ 以 4π 为周期变化的. 因此这个现象被称为分数约瑟夫森效应. I2e是由库珀对隧穿形成的， 是一个二阶过程， 载流子电荷为 2e ; Ie是由于MZM透过结区耦合引起的单电子隧穿形成的电流， 是一个一阶过程， 载流子电荷为e. 进一步，Ie中的 (-1)n0因子说明， 费米子的占据状态改变时，Ie的方向也会改变. 因此通过测量 Ie， 可以推断出费米子占据状态， 即读出MZM量子比特[20].

图 7 约瑟夫森结的量子比特读出装置[129] (a)为MZM约瑟夫森结的示意图. 结两端各自存在一个MZM.蓝色区域为一维拓扑超导区， 绿色区域为绝缘体区. 拓扑区的长度足够长， 使得两个MZM通过拓扑区的耦合可以忽略. γ1和 γ2 通过足够短的绝缘区耦合起来， 形成一个能量非零的费米子. 通过改变通过线圈的磁场通量 Φ 可以改变结两端的超导相位差 φ . (b) 直流分数约瑟夫森电流随 φ的变化关系. 其中红色虚线代表空占据态 | 0〉 ， 蓝色实线代表占据态 | 1〉 . 和传统的以 2 π 周期的约瑟夫森电流不同的是， 分数约瑟夫森电流关于 φ 的变化周期是 4 π . 量子比特|0〉 和 | 1〉 相 对 应的电 流方向 相 反， 因 此通过 测 量直流 约瑟夫森电流可以读出MZM量子比特Fig. 7. Basic set-up for qubit readout using a Josephson junction(JJ)[129]. (a) The schematic of a JJ with 2 MZMs residing at the junction. The blue region denotes 1d TSC and the green denoted trivial insulator. The TSC region should be long enough so that the coupling of the MZMs through TSC is negligible. γ1and γ2couple weakly at the junction， forming a non-zero energy fermion. The phase difference φ can be varied by changing the magnetic flux Φ .(b) The d.c. fractional Josephson current flowing across the junction versus φ . Instead of conventional 2 π -periodic JJ current induced by Cooper pair tunneling， the fractional JJ current induced by MZMs exhibits 4 π periodicity. The red dashed line denotes | 0〉 and the blue solid line denotes | 1〉 .The direction of the current is inverse for the qubiT| 0〉 and|1〉 ， which enables the readout of the qubit by measuring the direct Josephson current.

虽然理论上测量出 Ie就可以读出MZM量子比特， 但是实验上要测到 Ie还有很多难点要克服.一是怎么把 Ie从 I2e和其他由于工艺不完善引入的直流信号中分离出来. 比如， 当TSC 区不够长的时候，γ1和 γ2通过TSC区的耦合不能被忽略， 这时两个MZM形成的 Ie会既包含 4π 周期成分， 也包含 2π 周期成分，4π 周期的电流信号会被压制[129];当结区的绝缘区有超导存在的时候， 绝缘区的库珀对分裂成两个电子， 分别进去结两边的超导区， 引入更复杂的电流成分[130]. 二是调制 φ 的时间长度要适合. 若 φ 变化太慢， 系统将保持在基态上， 约瑟夫森结内部的局域费米子宇称将不守恒， 粒子数奇偶发生变化. 这样电流的 4π 周期会还原为 2π 周期;若 φ 变化太快， 会在系统中激发很多准粒子， 影响测量. 目前还没有实验能提取出直流分数约瑟夫森电流的符号以读出量子比特， 但是交流约瑟夫森效应已经有报道[94，95，131]. 直流分数约瑟夫森效应的测量需要在已有装置的基础上进一步改良工艺. 最近提出基于 4π 约瑟夫森效应的编织方案[132].

5.2 库仑阻塞

在考虑拓扑超导由于马约拉纳带来的基态简并时， 已经假设超导处于接地状态， 或者与一个宏观系统处于强耦合状态. 这样超导的电容无穷大，改变超导的费米子宇称， 即改变系统电子数并不需要能量. 当考虑有限尺寸的非接地超导， 即使存在MZM形成的 |0〉 和 |1〉 不同费米宇称态， 改变超导的费米子数(电子数)会存在电容充电能(可由电平控制). 这使得奇偶费米子宇称态不再简并， 从而可以从通过参数操控的能谱判定奇偶费米子宇称态.事实上， 此时超导不同的费米子数态可以变得不再简并， 具体由外加电平控制[39，96，133]. 一个研究得比较多的装置是通过电平开关把MZM和量子点耦合起来[37，39，96，124，133，134]， 自旋量子比特的测量也适用这个方法[135-137]. 装置如图8 所示， MZM通过量子隧穿效应和量子点弱耦合起来， 调节电平改变势垒高度可以改变MZM波函数的展宽， 从而调节MZM和量子点的有效耦合强度t， 控制电平开关的状态. 改变电平开关状态， 系统的基态能量将发生变化， 并受到MZM量子比特奇偶性的调控.如图8(a)所示的两MZM系统中，ti=1，2由关闭状态 (ti=0) 到 打开状态 (ti／=0) 会引起哈密顿量的基态能量的变化， 根据二阶微扰理论， 系统基态能 量为

其中 E0(j=0，1)分别对应为 T=0 时整体奇宇称和偶宇称基态的能量， 与量子点和超导充电能有关; EC为超导充电能; Ng为受电平控制的超导感应电荷.由方程(31)可见， 当 (t*1t2-t1t*2)／=0 时， 系统在调节量子点电平开关打开t的前后的基态能量变化ΔEg.s.受 量子比特的奇偶性 i γ1γ2调控. 如前所述，这个能量的变化可以通过测量系统的能谱、量子点电荷或者微分电容等方式测出， 从而读取出量子比特. 这个装置的优点一是MZM与量子点的耦合强度可以通过量子点的宏观电平参数来调节; 二是可拓展性好， 可通过拓展U型网络构造大规模的量子线路[39，124].

图 8 MZM和量子点耦合的装置简图[39](a) ， (b)分别是2个和4个MZM系统的装置; (c) 在调节量子点与MZM之间跃迁矩阵元的局域电平开关Fig. 8. Sketch of MZMs coupled to quantum dots[39]， with 2-MZM and 4-MZM system shown in (a) and (b)， and in (c) the local gates controlling the coupling between MZMs and quantum dots are shown.

5.3 马约拉纳干涉仪

Fu[138]提出利用MZM可以实现电子的隐态传输， 且传输相移与拓扑量子比特的宇称相关. 因此， 通过干涉装置探测传输后电子的相位变化可以得到MZM的宇称[38，133，139].

如图9 所示考虑拓扑超导体与普通金属的连接， 线路中穿过可调节的磁通 Ψ ， 线路中的电流将随磁通变化而改变. 调节超导体的电容能使之小于超导能隙但是大于电子的隧穿效应， 从而抑制电子直接从金属跃迁到超导体中. 借助于空间上分隔很远的MZM， 可以存在一端电子进入超导体， 另一端电子离开超导体进入金属的过程， 且跃迁过程中保持长距离的相位相干. 跃迁的哈密顿量可以写成如 下一般形式

其中a， b表示两个MZM所处的位置，Tab表示跃迁振幅， 对于费米子宇称 Pab=±1 的量子比特而言， 电子的传输相移分别为 - iTab和 i Tab. 当电流分别经过拓扑量子比特和具有足够长相位相干长度的普通金属后发生干涉时， 由于不同宇称的量子比特会给电子的传输相移带来 π 的相位差， 因此电导率表示为

图 9 马约拉纳干涉仪装置[38] (a)电流从上方经过两个MZM也即一个拓扑量子比特， 从下方经过相位相干长度足够长的金属， 电路中间穿过可调节的磁通. 两条路径的相位差由电子的传输相移以及穿过线路的磁通决定.(b)通过测量干涉后的电导可以得到拓扑量子比特的宇称信息. 实线和虚线分别对应宇称为1和—1的量子比特的电导信号Fig. 9. Majorana inferometry[38]. (a) One path goes through two MZMs i.e. a topological qubit while the other path goes through a normal metal with a sufficiently long phase-coherence length. Ψ is the applied magnetic flux enclosed by the two paths. The phase difference of two paths is determined by the phase transition shift and the magnetic flux，which can be measured by the conductance. (b) Majorana interferometer provides a projective measurement of the fermion parity. Solid line and dotted line represent the conductance signals corresponding to qubits with parity 1 and—1 respectively.

从电导的测量数据中可以读取费米子宇称的信息，进而探测编织的结果. 这个方案中费米子宇称的读取不需要将两个MZM相互融合， MZM之间始终保持较远的距离， 波函数的交叠很少， 因此会较少受到热涨落的干扰[140].

6 对称性保护的非阿贝尔统计

对拓扑超导体加上额外的对称性要求如时间反演对称[84，141-149]或空间点群对称性[150-154]等，可以丰富其体态的拓扑分类. 这对应着边界可以存在多个MZM， 且这些零模受到对称性的保护而不会互相耦合打开能隙. 对这样的体系进行编织操作将意味着对多个MZM进行交换. 一般而言这可能出现平庸的结果. 但是在对称性的保护下， 如果编织过程可以分解成多对独立的部分， 每对MZM的编织将满足对称保护的非阿贝尔统计[155，156]. 其中单次编织满足

其中 γiL代表编织前处于左边的N个MZM，γiR为之前处于右边的N个MZM. 此式描述的核心含义是， 在编织过程中， 一边的任意一个MZM只能看到另一边的其中一个MZM， 而看不到其他MZM.显然这种情况对于二维手性拓扑超导中的涡旋是无法成立的. 比如考虑左右分布且空间分离的多对涡旋， 对左右进行编织操作， 显然处于一边的MZM将受到另一边的所有MZM的影响.

MZM的对称保护非阿贝尔统计最早提出于时间反演对称的拓扑超导中[155]. 如图10(a)所示， 考虑一维时间反演对称的拓扑超导体， 超导体的每一端都存在两个MZM， 它们在时间反演算符T 作用下组成马约拉纳克拉默斯对 (Majorana Kramers pair， MKP)对称保护统计将使得一次编织后MZM满足. 然而， 对于时间反演对称保护的拓扑超导， 其MKP的非阿贝尔统计是否同样只需要哈密顿量满足时间反演对称， 这是值得探讨的问题. 事实上， 在编织过程中即使哈密顿量每个时刻满足时间反演对称， 边界的MZM与体态模的耦合仍可能诱导出MKP内部的有效耦合， 造成一端MKP的局域旋转操作，这对上述非阿贝尔统计带来破坏[157，158].

完善MKP对称保护非阿贝尔统计理论涉及一个深刻的问题， 即时间反演对称的动力学破缺[156，159]. 马约拉纳编织实际可看作一个时间从t=-T/2 到 T=T/2 的含时动力学演化过程， 其中T是单次编织时间. 相应的编织矩阵是一个幺正算符 Bij=U(T)，等于编织过程中哈密顿量 H (t) 的编时(time-ordered)积分. 由于时间反演对称是反幺正算符， 即使哈密顿量 H (t) 每个时刻满足时间反演对称， 其时间积分得到的幺正演化算符可能和时间反演算符不对易， 从而动力学破缺时间反演对称. 这是出现MKP局域操作的根本原因， 这一性质首先在文献[156]中被揭示. 具体而言， 这样即使对于每一时刻的哈密顿量满足时间反演对称[ H(t)，T]=0 ， 可发现U(T)=e-iHET

图 10 时间反演对称性保护的编织过程及非阿贝尔统计的结果[156] (a) 时间反演不变的拓扑超导体两wlxb 端各有一对MZM， 以T-型结的方案完成编织. (b)， (c) 不同的无条 序件强 下度 均 W 有0 下经过完整编织后MKP的演化， 结满果 足，非在 阿不 贝同尔统计. 经过一个周期后 给出绝热演化的条件Fig. 10. Time-reversal symmetry protected braiding process and the results of non-Abelian statistics[156]. The TSC in (a)hosts a pair of MZMs at each end， and the braiding is fulfilled by the T-junction scheme. (b)， (c) The evoluation of MKPs after the full braiding in the presence of different disorder strengh W0 . The non-Abelian statistics is confirmed by . The adiabatic conditi on is satisfied in that.

如果将编织演化算符写成， 则对于编 织过程对应的等效哈密顿量 HE， 由(35)式有

这便是时间反演对称在编织过程中发生动力学破缺. 因此， 在时间反演对称拓扑超导体中， 对称保护非阿贝尔统计需要: 1)拓扑超导哈密顿量在每个时刻满足时间反演对称; 2)整个编织产生的动力学演化满足类似的时间反演对称. 为了使MKP满足非阿贝尔统计， 需要在编织的动力学过程中重建时间反演对称性. 为此， 文献[156]定义了一个交换对称性(swapping symmetry)， 使得

即分解为两独立的编织[156]. 进一步研究证明， 该交换对称性S在一维的纳米线中总是成立. 即使存在静态无序， 由于左右两端的MZM等效走过互逆的路径， 系统的无序对两端MZM造成的局域旋转是相反的， 此时交换对称性的定义总可以成立. 同时由于对称性保护， 噪声带来的影响为二阶项， 且能被噪声的随机分布所抑制(如图10(b)(c)所示)[156].MKP的编织带来系统量子态的变化， 而由时间反演对称拓扑超导中马约拉纳组成的不同量子比特同样可以通过约瑟夫森电流读取[155，160]. 对称性保护的非阿贝尔统计预期能很好地在实际体系中实现并为拓扑量子计算提供更丰富的操作[40].

7 总结与展望

本文对拓扑超导中MZM的非阿贝尔统计， 以及在拓扑量子计算中的应用做了介绍. 在对承载MZM的拓扑超导做简单回顾后， 对MZM的非阿贝尔统计特性， 构造拓扑量子比特和量子门进行介绍. 然后， 重点讨论MZM的几种主要编织方案，以及由MZM构成的拓扑量子比特测量方案， 并对这些方案做了比较. 最后对新发展的马约拉纳克拉默斯对满足的对称保护非阿贝尔统计进行了介绍，这可能对MZM的研究带来新的方向.

从MZM被提出可以被用来做量子计算到今天， 拓扑量子在这约二十年间取得了大量包括理论和实验上的重要进展， 已经有大量的实验证据可以支持MZM的存在， 但是这些实验证据都是对局域的单个MZM的性质直接测量操作. 要论证MZM的存在， 最直接也是最确凿无疑的证据就是验证MZM的非阿贝尔统计性质， 也就是实验验证: 编织两个MZM前后系统的基态发生了变化.这包括两方面的难题: 一是读取两个MZM融合后的量子态; 二是实现MZM的非阿贝尔编织操作、从而实现对量子态的操控. 目前还没有任何实验可以证明两个MZM组成的复费米子态的非局域性，并直接读出两个MZM的量子态. 因此， 当前研究的最大热点还是为两个MZM融合后的量子态读取寻找实验可操作的平台， 进而实现MZM的编织操作并读取出编织操作前后系统量子态的变化. 一旦实现了MZM的编织操作、并验证了相应的非阿贝尔统计， 进一步应用于拓扑量子计算则将成为可能.