-
早在上世纪30年代, Majorana求解了Dirac相对论协变的电子运动方程, 发现了一种不带电荷的费米子, 它的反粒子是其自身. 人们为了寻找它的踪迹一直在不懈地努力, 然而最终Majorana零模在凝聚态物理中被发现, 并成为重要的研究课题[1-5]. 超导体系中
$ U(1) $ 规范对称性的破缺为Majorana费米子的产生提供了可能性, 人们已经在具有强自旋-轨道耦合的半导体纳米线[6-10], 磁性原子链[11-13], 平面约瑟夫森结[14-16]以及常规超导体和拓扑绝缘体[17-19]的界面等体系中发现了它的存在. 另一方面由于Majorana费米子具有局域性且满足非阿贝尔统计[20-22]等特性, 使得它成为实现容错拓扑量子计算[5,23]最有力的竞争者. 由于拓扑量子计算的巨大应用前景, 使得Majorana费米子相关性质的研究越来越被人们重视. 特别是近年来, 随着冷原子技术的发展, 人们发现通过周期驱动光格子可以实现物质拓扑态[24-26], 通过周期驱动具有p波配对的超导量子线, 有可能会产生额外的π模[27]. 通过多个时间周期驱动的Kitaev链产生了可以支持Majorana零模的新区域, 对Majorana费米子的寻找提供了理论基础[28]. 拓扑相最初是在厄密系统中发现的, 但人们对非厄密系统中拓扑相的研究也存在很大的兴趣[29-33]. 由于Majorana零模可以在非厄密体系中出现且可以持续存在, 其对环境具有很强的鲁棒性, 为更好地研究Majorana费米子提供方法. 最近, Wu等[34]阐述了实现非阿贝尔编织的一种新途径, 利用Jackiw-Rebbi零模也可以实现非阿贝尔编织, Jackiw-Rebbi零模不具有Majorana零模的自共轭特性, 其可以出现在非超导体系中. Jackiw-Rebbi零模的研究为拓扑量子计算提供了新的思路. 有趣的是, Majorana零模可以被认为是Jackiw- Rebbi零模在具有粒子-空穴对称性时的特例[35,36].Kitaev链是研究Majorana费米子的重要模型, 在此基础上人们意识到通过对Kitaev链的调制可以极大地改变系统的拓扑相变过程. 如Lang和Chen[37]研究了周期性调制对Majorana费米子产生的影响, 他们发现随着调制强度的增大, 拓扑非平庸超导相可能会被破坏. 由于Majorana零模的稳定性受到超导能隙的保护, 因此在加入周期调制化学势的情况下Majorana费米子可能是不稳定的, 会随着调制化学势强度的增大而消失. 然而在某些特殊参数下, 调制强度无法改变Majorana费米子的存在性. 与此同时, Cai等[38]讨论了非公度调制对拓扑相变的影响, 发现随着非公度调制强度的增加系统将经历从拓扑非平庸相向平庸的安德森局域相的转变. 随后相当多的工作对调制的Kitaev链进行了深入的研究[39-41]. 本文将讨论(准)周期调制的p波超导量子线系统中的拓扑量子相变.
-
考虑一维具有(准)周期调制的p波超导量子线, 其哈密顿量可以写为
$\hat{H} = \sum\limits_{i}\!\left[\!\left( -t\hat{c}_{i}^{\dagger }\hat{c}_{i+1}+\varDelta \hat{c}_{i}\hat{c}_{i+1}+ h.c.\right) +V_{i}\hat{n}_{i}\!\right], $ 其中
$ \hat{c}_{i}^{\dagger} $ ($ \hat{c}_{i} $ )是费米子的产生(湮灭)算符,$ \hat{n}_{i} = \hat{c}_{i}^{\dagger }\hat{c}_{i} $ 是粒子数算符, t是最近邻格点间的跃迁强度, 被设为能量单位($ t = 1 $ ),$ \varDelta $ 为超导配对项中的对产生或对湮灭的强度. 化学势项可以写为$ V_{i} = V\frac{\cos \left( 2{\text{π}} i\alpha + \delta \right) }{1-b\cos \left( 2{\text{π}} i\alpha+\delta \right) }, $ 其中V是化学势的强度,
$ \delta $ 是任意的相移.$ \alpha $ 控制系统的调制周期, 若$ \alpha = p/q $ 是有理数(p和q是互质的整数), 则$ V_{i} $ 是公度势; 若$ \alpha $ 是无理数, 系统则具有非公度调制. 化学势是参数b的连续函数, 其中$ b\in [0, 1) $ . 在没有超导配对的情况下, 即$ \varDelta = 0 $ , 当$ b = 0 $ 时, 若$ \alpha $ 为有理数, 系统处于拓扑非平庸相, 由非零整数的陈数所标记[42]; 若$ \alpha $ 为无理数, 模型退化为著名的AA模型[43], 此时如果$ V < 2 t $ , 系统中所有的单粒子本征态为扩展态并且具有非平庸的拓扑性质, 而当$ V > 2 t $ 时, 所有的本征态都为局域态,$ V = 2 t $ 是扩展到局域相的转变点, 此时所有的本征态展现多分形的特性, 而这一系统中并不存在迁移率边[44]; 对于$ b\ne 0 $ 且$ \alpha $ 为无理数的情况[45], 系统具有能量依赖的自对偶特性, 其迁移率边可以解析地表示为$ E_{\rm c} = (2 t-V)/b $ . 对于存在超导配对的情况, 即$ \varDelta\ne 0 $ , 若$ \alpha $ 为有理数, 模型哈密顿量为周期调制的p波超导量子线, 已经被广泛地研究[37], 文献[37]中指出此系统的拓扑相变依赖于相移$ \delta $ , 而在某些特殊$ \delta $ 点系统一直处于拓扑非平庸相不会受周期调制强度V所控制; 对于非公度调制, Cai等[38]指出随着非公调制强度的增大, 系统经历一个由拓扑非平庸相到安德森局域相的转变, 转变点在$ V_{\rm c} = 2 t+2\varDelta $ 处. 由此可见, 在$ b = 0 $ 的情况, 模型具有丰富的拓扑及局域化特性, 已经引起了广泛的兴趣. 在这篇文章中我们关注$ \varDelta\ne 0 $ ,$ b\ne 0 $ ,$ \alpha $ 分别为有理数和无理数情况下系统的拓扑相变, 以及在$ \alpha $ 为无理数时系统的局域化特性.通过 Bogoliubov-de Gennes (BdG)变换[46-48]把系统的哈密顿量(1)对角化, 定义一组准粒子算符:
$ \hat{\eta}_{n}^{\dagger} = \sum\limits_{i = 1}^{L}\left[ u_{n,i}\hat{c}_{i}^{\dagger }+\nu_{n,i}\hat{c}_{i}\right], $ 其中L是系统的格点数, n是能级指标且
$ n = $ $ 1, \cdots, L $ . 由于在哈密顿量(1)中所有的参数都选为实数,$ u_{n, i} $ 和$ \nu_{n, i} $ 也均为实数. 哈密顿量可以用准粒子算符表示为$ \hat{H} = \displaystyle\sum\nolimits_{n = 1}^{L}E_{n}\left( \hat{\eta}_{n}^{\dagger }\hat{\eta}_{n}-{1}/{2}\right) $ , 其中$ E_{n} $ 是准粒子的本征能量. 由对角化关系$ \left[ \hat{\eta}_{n}^{\dagger }, H\right] = -E_{n}\hat{\eta}_{n}^{\dagger} $ , 得到下面的BdG耦合方程:$ \left( {\begin{array}{*{20}{c}} {\hat h}&{\hat \varDelta }\\ { - \hat \varDelta }&{ - \hat h} \end{array}} \right)\left( {\begin{aligned} {{u_n}}\\ {{\nu _n}} \end{aligned}} \right) = {E_n}\left( {\begin{aligned} {{u_n}}\\ {{\nu _n}} \end{aligned}} \right), $ 其中
$\begin{aligned} & \hat{h}_{ij} = -t\left( \delta_{j, i+1}+\delta_{j, i-1}\right) +V_{i}\delta_{ji},\\ & \hat{\varDelta}_{ij} = -\varDelta\left( \delta_{j, i+1}-\delta_{j, i-1}\right), \\ & u_{n}^{\rm T} = \left( u_{n, 1}, ..., u_{n, L}\right) ,\\ & \nu_{n}^{\rm T} = \left( \nu_{n, 1}, ..., \nu_{n, L}\right). \end{aligned}$ 通过求解BdG方程, 可以得到准粒子的本征能量及其相应的本征波函数. 由于BdG方程满足电子-空穴对称性, 即
$ \hat{\eta}_{n}\left(E_{n}\right) = \hat{\eta}_{n}^{\dagger}\left(-E_{n}\right) $ , 系统的能谱关于零点对称. 系统的基态对应于所有负的准粒子的能级被填满的情况. 在下面的分析中取$ b = 0.5 $ . -
这一小节讨论
$ \alpha $ 为有理数情况下, 系统的拓扑相变. 在开边界条件下, 我们通过数值求解BdG方程(4)得到准粒子的本征能量$ E_{n} $ , 若系统处于拓扑非平庸相, 能谱中会出现零能的Majorana边缘态, 而当系统处于拓扑平庸相, Majorana零模将消失. 图1计算了在$ b = 0.5 $ ,$ \varDelta = 0.2 $ ,$ V = 1.5 $ 和$ \delta = 0 $ 时, 能谱随参数$ \alpha $ 变化的情况, 即Hofstadter蝴蝶谱[49,50], 其中红色点表示非平庸的零模. 随着$ \alpha $ 的增加, 系统表现出复杂的拓扑相变过程. 作为具体的例子, 我们将分别讨论$ \alpha = 0 $ ,$ 1/2 $ ,$ 1/3 $ 的情况.图 1 Hofstadter蝴蝶谱: 随
$\alpha$ 变化的能谱, 红色点是零能Majorana费米子b = 0.5, L = 120,$\varDelta=0.2$ , V = 1.5, δ = 0Figure 1. Hofstadter butterfly: the energy spectrum varying with
$\alpha$ . The red dotted point denotes the Majorana Fermion.$b=0.5,\; L=120,\; \varDelta=0.2, \;V=1.5$ and$\delta=0$ 在
$ \alpha = 0 $ 时, 哈密顿量退化为标准的Kitaev模型[5], 系统在$ \left|\dfrac{V\cos{\delta}}{1-b\cos{\delta}}\right| = 2 t $ 处经历一个拓扑相变, 在$ \left|\dfrac{V\cos{\delta}}{1-b\cos{\delta}}\right| < 2 t $ 区域处于由Majorana零模所标记的拓扑非平庸相. 可以看出当$ \delta $ 取$ {\text{π}}/2 $ 奇数倍时, 系统将一直处于拓扑非平庸相, 并不依赖于V的取值.我们知道, 非平庸的Majorana零模可以由
$ Z_2 $ 拓扑不变量来表征[5,37]. 对于$ \alpha = 1/2 $ 和$ 1/3 $ 的情况, 可以通过计算$ Z_2 $ 拓扑不变量, 解析地得到系统的相变点. 考虑具有周期性边界的系统并对其进行傅里叶变换,$ \hat c_{i} = \hat c_{s, l} = \sqrt{{q}/{L}}\sum_{k}\hat c_{s, k}{\rm e}^{{\rm i}kql} $ . 其中,$ i = s+\left( l-1\right) q, s = 1, \cdots, q $ 表示一个超导元胞内的格点数,$ l = 1, \cdots, L/q $ 是第l个超导元胞的位置, k表示动量, 其取值范围为$ [0, 2{\text{π}}/q] $ . 哈密顿量(1)进行傅里叶变换之后可以写为$\begin{split} \hat{H}_{k}= \;&\sum\limits_{k}\bigg[ \sum\limits_{s = 1}^{q-1}\left( -t\hat{c}_{s,k}^{\dagger } \hat{c}_{s+1,k}+\varDelta \hat{c}_{s,k}\hat{c}_{s+1,-k}\right)\bigg. \\ & \bigg.-t\hat{c} _{q,k}^{\dagger }\hat{c}_{1,k}{\rm e}^{{\rm i}kq}+\varDelta \hat{c}_{q,k}\hat{c} _{1,-k}{\rm e}^{-{\rm i}kq}+{\rm {h.c.}}\bigg] \\ & -\sum\limits_{k}\sum\limits_{s = 1}^{q}V_{s}\left( \hat{c}_{s,k}^{\dagger }\hat{c}_{s,k}- \frac{1}{2}\right). \end{split} $ 在动量空间下, 我们定义一组准粒子算符为[51]:
$ \hat{\gamma}_{2 s-1}\left( k\right) = \hat{c}_{s, k}+\hat{c}_{s, -k}^{\dagger },\; \hat{\gamma}_{2 s}\left( k\right) = \left( \hat{c}_{s, k}-\hat{c} _{s, -k}^{\dagger }\right) /i $ , 它满足反对易关系:$ \left\{ \hat \gamma_{m}^{\dagger}\left( k\right) , \hat \gamma_{n} ( k^{^{\prime }}) \right\} = 2\delta_{mn}\delta_{kk^{\prime}} $ 以及$ \hat \gamma_{m}^{\dagger}\left( k\right) = \hat \gamma_{m}\left( -k\right) $ . 可以看出只有$ \hat \gamma_{m}\left( 0\right) $ 和$ \hat \gamma_{m}\left( {\text{π}}/q\right) $ 满足Majorana费米子算符的定义, 即$ \hat\gamma_{m}^{\dagger}\left( 0\right) = \hat\gamma_{m}\left( 0\right) $ ,$ \hat\gamma_{m}^{\dagger}\left( {\text{π}}/q\right) = \hat\gamma_{m}\left( -{\text{π}}/q\right) = \hat\gamma_{m}\left( {\text{π}}/q\right) $ . 在新的算符基矢下, 可以把哈密顿量重新写成如下形式:$ \hat{H}_k = \frac{i}{4}\sum\limits_{k}\sum\limits_{m,n}B_{m,n}\left( k\right) \hat\gamma_{m}\left( -k\right) \hat\gamma_{n}\left( k\right). $ 对于
$ s = 1, \cdots, q $ 时,$ B_{2 s-1, 2 s}\left( k\right) = -B_{2 s, 2 s-1}\left( k\right) = -V_{s};$ 对于
$ s = 1, \cdots, q-1 $ 时,$\begin{aligned} & B_{2 s-1, 2 s+2}\left( k\right) = -B_{2 s+2, 2 s-1} = \varDelta-t,\\ & B_{2 s, 2 s+1}\left( k\right) = -B_{2 s+1, 2 s}\left( k\right) = \varDelta+t,\end{aligned}$ 对于
$ s = q $ , 有$\begin{aligned}& B_{2 s-1, 2 s}\left( k\right) = -B_{2 s, 2 s-1}\left( k\right) = -V_{s},\\ & B_{2, 2 q-1}\left( k\right) = -B_{2 q-1, 2}^{\ast}\left( k\right) = -\left( \varDelta-t\right) {\rm e}^{-{\rm i}kq}.\end{aligned}$ $ { B}\left( k\right) $ 是一个$ 2 q\times2 q $ 的矩阵, 并且只有$ { B}(0) $ 和$ { B}({\text{π}}/q) $ 是反对称矩阵. 系统的$ Z_2 $ 拓扑不变量可以定义为[5,51]:$ M \!=\! {\rm{sgn}}[{\rm{Pf}}( { B}( 0))] {{\rm sgn}}\left[{\rm {Pf}}( { B}( {\text{π}}/q) ) \right] $ , 其中$ {\rm {Pf}}\left( { A}\right) = \dfrac{1}{2^{N}N!}\displaystyle\sum\nolimits_{P}{\rm {sgn}}( P) { A}_{P_{1}P_{2}}\cdots{ A}_{P_{2 N-1}P_{2 N}} $ 是反对称矩阵A的Pfaffian, P代表矩阵A中$ 2 N $ 个元素的置换,$ {\rm {sgn}}\left( P\right) $ 表示置换的符号.$ M = 1 $ 对应拓扑平庸相,$ M = -1 $ 对应拓扑非平庸相, 而拓扑相边界可以由$ M = 0 $ 来标记. 当$ \alpha = 1/2 $ 时,$\begin{aligned}& {\rm {Pf}}\left[ { B}\left( 0\right) \right] = \frac{-V^{2}\cos ^{2}\delta }{1-b^{2}\cos^{2}\delta }-4 t^{2},\\ & {\rm {Pf}}\left[ { B}\left( \frac{{\text{π}} }{2}\right) \right] = \frac{-V^{2}\cos ^{2}\delta}{1-b^{2}\cos ^{2}\delta }+4\varDelta ^{2} .\end{aligned}$ 显然,
$ {\rm{Pf}}\left[ { B}\left( 0\right) \right] < 0 $ , 系统的拓扑相边界由$ {\rm {Pf}}\left[ { B}\left( \dfrac{{\text{π}}}{2}\right) \right] = 0 $ 得出, 即$ \varDelta = \frac{V|\cos{\delta}|}{2\sqrt{1-b^{2}\cos^{2}\delta}}. $ 图2(a)展示了
$ b = 0.5 $ ,$ \alpha = 1/2 $ ,$ \delta = 0 $ 时, 系统拓扑相图. 图中的黑色实线对应方程(7)所示的解析结果, 红色三角表示的是通过数值求解BdG方程(4)得到的相变点. 可以看到数值结果与解析解得到的结果一致. 在区域Ⅰ, 系统处于拓扑非平庸相, 区域Ⅱ对应于系统处于拓扑平庸相. 我们可以看到, 当$ \delta = 0 $ 时, 系统会经历拓扑非平庸相到拓扑平庸相的转变. 由方程(7)可知, 当$ \delta $ 取值为$ {\text{π}}/2 $ 的奇数倍时, 任意小的$ \varDelta $ 将导致系统处于拓扑非平庸相, 而不依赖于周期调制的强度. 图3(a)—图3(c)展示了$ b = 0.5 $ ,$ \alpha = 1/2 $ ,$ \varDelta = 0.2 $ , 不同V时, 能谱随着相移$ \delta $ 变化的情况. 在V比较小的时候, 如图3(a)所示,$ V = 0.2 $ , 在整个相移参数空间中, Majorana零模一直存在. 随着V的增大, 能隙逐渐减小, 当它超过某个临界值时, 能隙将在某些$ \delta $ 的位置关闭, 随后再次打开, 而此时零模消失[图3(b),$ V = 0.5 $ ], 对应于系统从拓扑非平庸相到拓扑平庸相的转变. 然而当V足够大, 如图3(c),$ V = 3 $ , 除了在$ \delta = {\text{π}}/2 $ 和$ 3{\text{π}}/2 $ 处Majorana零模存在外, 几乎所有的$ \delta $ 区域都处于平庸相, 并且无论V值取多大, 这两点的零模始终存在, 这与我们的解析结果相一致.图 2 在
$b=0.5$ 时, 参数$\varDelta-V$ 平面的拓扑相图 (a)$\alpha= $ $1/2, L=120$ ; (b)$\alpha=1/3, L=120$ ; (c)$\alpha= ( \sqrt{5}-1) /2$ , L = 2584Figure 2. Topological phase diagram in
$\varDelta-V$ plane with$b=0.5$ . (a)$\alpha=1/2, L=120$ ; (b)$\alpha=1/3, L=120$ ; (c)$\alpha= ( \sqrt{5}-1) /2$ , L = 2584图 3 在开边界条件下, 本征能量随相移
$\delta$ 的变化.$b=0.5$ ,$\varDelta=0.2$ ,$L=2584$ Figure 3. Energy varying with phase shift
$\delta$ with$b=0.5$ ,$\varDelta=0.2$ and$L=2584$ under open boundary condition.当
$ \alpha\!=\! \dfrac13$ 时, 可以计算$ {\rm {Pf}}\left[ { B}\left( 0\right) \right] $ 和$ {\rm{Pf}}\left[ { B}\left( \dfrac{{\text{π}} }{3}\right) \right] $ , 由$ M = 0 $ , 可以解析地得到系统的拓扑相边界:$ \cos \left( 3\delta \right) = \frac{-2t^{2}\left[ \left( -4+3b^{2}\right) t+3b\left\vert V\right\vert \right] +6\varDelta ^{2}\left[ \left( 4-3b^{2}\right) t+b\left\vert V\right\vert \right] }{\left( 2bt-\left\vert V\right\vert \right) \big( b^{2}\left( 3\varDelta ^{2}+t^{2}\right) +2b\left\vert V\right\vert t+\left\vert V\right\vert ^{2}\big) }, $ 特别是, 当
$ b = 0 $ 时, 相边界可以写为一个简单的表达式[37]:$ V^3 |\cos{3\delta}| = 8 t(t^2+3\varDelta^2) $ . 在$ \cos{3\delta} = 0 $ 时, 系统始终处于拓扑非平庸相, 并且不依赖于V的取值. 图2(b)展示了$ b = 0.5 $ ,$ \alpha = 1/3 $ 和$ \delta = 0 $ 时的拓扑相图. 黑色实线为解析结果, 而红色三角为数值结果. 由图可知,$ \delta = 0 $ 时, 在某一特定的$ \varDelta $ 下, 随着周期调制强度V增强, 系统将出现一个拓扑相变. 图3(d)—图3(f)分别展示了$ b = 0.5 $ ,$ \alpha = 1/3 $ ,$ \varDelta = 0.2 $ ,$ V = 0.2 $ ,$ 2 $ 和$ 6 $ 时, 能量以相移$ \delta $ 为函数变化的情况. 在小V情况, 系统在不同的$ \delta $ 参数下, 始终出现Majorana零模[图3(d)], 而随着V的增加, 拓扑非平庸的区域逐渐减小[图3(e)], 当调制强度足够大时, 拓扑非平庸区域完全消失, 此时系统中并不存在某个特殊的$ \delta $ 使得Majorana零模一直存在[图3(f)], 这与$ b = 0 $ 的情况不符. 我们可以看到, 图3(f)中虽然某些$ \delta $ 下最低能量接近于零, 但它并不是Majorana零模, 其准粒子的最低能量不低于0.07. 由此可见, 对于$ b\ne 0 $ ,$ \alpha $ 为有理数的情况, 在某个固定的超导配对强度$ \varDelta $ 和调制强度V时, 系统的拓扑相变强烈地依赖于相移$ \delta $ . 然而在某些$ \alpha $ 值下, 并不存在与$ b = 0 $ 情况类似的特殊$ \delta $ 值, 使得拓扑相变不依赖于调制强度V. -
上一节中的结果表明
$ b\ne 0 $ 时, 周期调制的p波超导的拓扑性质依赖于$ \alpha $ 的取值, 一个有趣的问题是当系统的调制周期趋近于无穷时, 系统的相图将如何变化? 不失一般性, 在这一节中我们固定$ b = 0.5 $ ,$ \alpha = \left( \sqrt{5}-1\right) /2 $ 和$ \delta = 0 $ 作为参数来进行讨论. 图2(c)展示了$ b = 0.5 $ ,$ \delta = 0 $ 时准周期调制的p波超导线系统的量子相图, 其中红色三角为数值计算所得的相变界. 区域I为拓扑非平庸的超导相, 区域II为有能隙的平庸超导相, 而区域III为无能隙的平庸局域相. 当$ 0 < \varDelta < 0.33 $ 时, 系统经历I相$ \to $ II相$ \to $ III相的转变. 而当$ \varDelta > 0.33 $ 时, 系统经历I相$ \to $ III相的转变.为了得到图2(c)中所示的相图, 我们首先分别计算在开边界和周期边界条件下, 系统的准粒子最低激发能量, 如图4(a)所示. 以
$ \varDelta = 0.2 $ 为例, 图4(a)展示了最低激发能量$ E_{1} $ 随准周期调制强度V的变化. 图中黑色实线表示周期性边界的情况, 黑色方块表示开边界的情况. 当$ V < 1.5 $ 时, 开边界条件下展示了零能, 而周期边界条件下存在有限的能隙, 这表明在开边界条件下系统中存在零模. 在图4(b)和图4(c)中分别展示了在开边界条件下$ V = 1 $ 时, 最低激发态的空间分布$ \phi_1 $ 和$ \psi_1 $ , 这里$ \phi_{1, j} = \left( u_{1, j}+\nu_{1, j}\right) $ ,$ \psi_{1, j} = \left( u_{1, j}-\nu_{1, j}\right) $ [38]. 此时最低激发态$ \phi_1 $ 和$ \psi_1 $ 分别位于边界的左右两端,$ \phi_1 $ 和$ \psi_1 $ 的振幅不会重叠在一起, 而是分裂为两个在空间上独立的Majorana边缘态, 此时系统属于有Majorana零模的超导相. 当$ V\in(1.5, 2.5) $ 时, 开边界条件和周期边界条件下最低激发能量大于零, 展示了相同的能隙, 并没有展示边缘态, 并且在开边界条件下最低激发态$ \phi_1 $ 和$ \psi_1 $ 的振幅会重叠在一起, 且分布在整个空间, 此时系统属于超导相[如图4(b), 图4(c),$ V = 2 $ ]. 当$ V > 2.5 $ 时, 开边界和周期边界条件下, 能隙均消失, 其最低激发能量为零. 以$ V = 3 $ 为例, 其最低激发态$ \phi_1 $ 和$ \psi_1 $ 局域在空间某一点上, 并不局域在边界位置, 表明此区域的零能态不是Majorana零模[如图4(b),图4(c),$ V = 3 $ ]. 从准粒子的最低激发能量及其本征态的空间分布可以看出, 对于$ \alpha $ 为无理数的情况, 系统存在三种不同的相.图 4 (a)在开边界和周期性边界条件下最低激发态能量
$E_{1}$ 随准周期调制强度V的变化及其空间分布$\phi_1$ (b)和$\psi_1$ (c),$\alpha= $ $(\sqrt{5}-1)/2$ ,$b=0.5$ ,$\varDelta=0.2$ ,$L=2584$ Figure 4. (a) The lowest excitation energies,
$E_{1}$ , varying with the quasi-periodic modulation amplitude, V, under OBC and PBC, respectively. The spatial distribution of the lowest excited state$\phi_1$ (b),$\psi_1$ (c).$\alpha=(\sqrt{5}-1) /2$ ,$b=0.5$ ,$\varDelta=0.2$ ,$L=2584$ .为了进一步确定系统中三种不同相的拓扑特性, 我们用
$ Z_{2} $ 拓扑不变量来表征其拓扑性质. 在非公度势的情况下, 我们用散射矩阵$ {{S}} $ 来计算$ Z_{2} $ 拓扑不变量[52-54]. 散射矩阵$ {{S}} $ 与在费米能级$ E_{\rm F} = 0 $ 处的入射波和出射波的振幅有关,$ {{S}} = \left( \begin{array}{cc} {R} & {T}^{\prime } \\ {T} & {R}^{\prime } \end{array} \right), $ 这里, 2 × 2的子块
$ {R}, {R}^{\prime} $ 和$ {T}, {T}^{\prime} $ 分别为在超导线两端的反射和透射矩阵.$ Z_{2} $ 拓扑不变量定义为:$ M = {\rm {sgn}} [{\rm{Det}}({ R})] $ . 只有当$ M = -1 $ 时, 在超导量子线两端才会出现非平庸的Majorana费米子. 散射矩阵可以通过转移矩阵方法得到. 基于哈密顿量(4), 零能的薛定谔方程给出:$ \left( \begin{array}{c} \hat{t}_{i}{ {\varPhi}} _{i} \\ { {\varPhi}} _{i+1} \end{array} \right) = \tilde{ W}_{i}\left( \begin{array}{c} \hat{t}_{i-1}^{\dagger }{ {\varPhi}} _{i-1} \\ { {\varPhi}} _{i} \end{array} \right), $ 这里
$ \tilde{{ {W}}}_{i} = \left( \begin{array}{cc} 0 & \hat{t}_{i}^{\dagger } \\ -\hat{t}_{i}^{-1} & -\hat{t}_{i}^{-1}\hat{\lambda}_{i} \end{array} \right), $ $ \hat{t}_{i} = -t{ {\sigma}}_z+{\rm i}\Delta{ {\sigma}}_y $ ,$ \hat{\lambda}_{i} = V_i{ {\sigma}}_z $ ,$ { {\varPhi}} _{i} = \left( u_{i}, \nu _{i}\right) ^{\rm T} $ 是第i个格点上准粒子的波函数,$ { {\sigma}}_{y} $ 和$ { {\sigma}}_{z} $ 分别为y和z组分的泡利矩阵. 在超导量子线两端($ i = 1 $ 和L)总的转移矩阵为$ \tilde{ W} = \tilde{ W}_{L}\tilde{ W}_{L-1}\cdots \tilde{ W}_{2}\tilde{ W}_{1} $ . 通过相似变换, 把转移矩阵写在新的基矢下[55,56]:$ { W} = { U}^{\dagger }\tilde{ W}{ U},\quad { U} = \sqrt{\frac{1}{2}}\left( \begin{array}{cc} { I} & { I} \\ {\rm i}{ I} & -{\rm i}{ I} \end{array} \right), $ 这里I为2 × 2的单位阵. 在这个基矢下, 透射和反射矩阵的关系为
$ \left( \begin{aligned} { T} \\ 0 \end{aligned} \right) = { W}\left( \begin{aligned} { I} \\ { R} \end{aligned} \right),\left( \begin{aligned} { R}^{\prime } \\ { I}\; \end{aligned} \right) = { W}\left( \begin{aligned} 0\;\\ { T}^{\prime } \end{aligned} \right). $ 拓扑不变量M就可以通过计算转移矩阵W得到. 如图5(a)所示, 我们计算了
$ b = 0.5 $ ,$ \varDelta = 0.2 $ 时, 系统的拓扑不变量M随着调制强度V变化的情况. 从图中可以看出, 当$ V < 1.5 $ 时,$ M = -1 $ 对应于由Majorana零模所标记的拓扑非平庸的超导相. 而当$ V > 1.5 $ 时,$ M = 1 $ 对应为拓扑平庸相. 由此可以确定区域I为拓扑非平庸的超导相, 而在区域II和III中, 系统展现了拓扑平庸的特性.图 5 (a)
$Z_{2}$ 拓扑不变量随非公度势强度的变化; (b) V = 2时$\mathrm{IPR}_{1}$ 的标度分析; (c)$V=3$ 时$\mathrm{IPR}_{1}$ 的标度分析$b=0.5$ ,$\alpha= (\sqrt{5}-1) /2$ ,$\varDelta=0.2$ ,$L=2584$ Figure 5. (a)
$Z_{2}$ topological invariant varying with the strength of the potential V; (b) the scaling of$\mathrm{IPR}_{1}$ $V=2$ ; (c) the scaling of$\mathrm{IPR}_{1}$ $V=3$ . Here,$\alpha= (\sqrt{5}-1)/2$ ,$b=0.5$ ,$\varDelta=0.2$ ,$L=2584$ .区域II和III的最低激发态展现了不同的局域化特性, 通过计算逆参与率(inverse participation ratio, IPR)[57-61],
$ {\rm{IPR}}_{n} = \displaystyle\sum\nolimits_{j = 1}^{L}\left( u_{n, j}^{4}+\nu _{n, j}^{4}\right) $ [38,52], 区分系统最低激发态的局域和扩展性质. 这里n是能级指标,$ u_{n, j} $ 和$ \nu _{n, j} $ 是BdG方程(4)的本征态, 满足归一化条件,$ \displaystyle\sum\nolimits_{j = 1}^{L}\left( u_{n, j}^{2}+\nu _{n, j}^{2}\right) = 1 $ . 对于扩展态, IPR的值以$ 1/L $ 趋近零; 而对于局域态, 其IPR$ \propto (1/L)^{0} $ 趋于一个有限值. 图5(b)和图5(c)分别展示了$ \varDelta = 0.2 $ ,$ V = 2 $ 和$ 3 $ 时, 最低激发态$ {\rm{IPR}}_{1} $ 随着系统尺寸的标度行为.$ V = 3 $ 时, 最低激发态$ {\rm{IPR}}_{1} $ 不随尺寸L的变化而变化, 在L趋近于无穷时,$ {\rm{IPR}}_{1} $ 的值趋近于$ 0.45 $ , 表明此时其最低激发态为局域态. 而对于$ V = 2 $ 的情况, 最低激发态$ {\rm{IPR}}_{1} $ 随着$ 1/L $ 趋近于$ 0 $ , 展现扩展的特性. 由此可知, 区域II为拓扑平庸的超导相, 而区域III对应为拓扑平庸的局域相.当
$ \varDelta = 0 $ 时, 系统中存在迁移率边[45], 其解析表达式为$ E_{\rm c} = (2 t-V)/b $ . 图6(a)展示了$ \alpha = \left( \sqrt{5}-1\right) /2 $ ,$ b = 0.5, L = 144, \delta = 0 $ 和$ \varDelta = 0 $ 时不同能量$ E_{n} $ 的逆参与率随着调制强度V变化的情况, 其中蓝色实线表示迁移率边的解析解. 随着p波超导配对势的引入, 即$ \varDelta\neq0 $ , 系统中的迁移率边将如何改变? 首先考虑$ \varDelta\to 0 $ 的情况. 以$ \varDelta = 0.01 $ 为例[图6(b)], 可以看到原来能谱中间区域展现局域态特性的能态随着微小的超导配对项的引入开始变成扩展态, 而高能和低能部分并没有发生显著变化. 当$ \varDelta $ 为有限大时, 如图6(c)$ \varDelta = 0.5 $ 时, 可以看到高能部分的局域化特性并没有发生显著的变化, 中能部分局域化区域扩大, 而低能部分扩展区向局域化区域扩张. 随着$ \varDelta $ 值的进一步增加, 高能和中能部分的局域化区域进一步扩大, 而低能部分的局域化区不断缩小[如图6(d),$ \varDelta = 0.8 $ ]. 由此可见, 由于超导配对项的引入, 迁移率边将无法用一个解析的形式表示.图 6
${\rm {IPR}}$ 随准周期调制强度V和本征能量$E_{n}$ 的变化$\alpha=(\sqrt{5}-1) /2$ ,$b=0.5, L=144, \delta=0$ (a)$\varDelta=0$ ; (b)$\varDelta=0.01$ ; (c)$\varDelta=0.5$ ; (d)$\varDelta=0.8$ Figure 6.
${\rm {IPR}}$ varying with the amplitude of quasi-periodic modulation V and energy$E_{n}$ .$\alpha=(\sqrt{5}-1) /2$ ,$b=0.5, L=144$ ,$\delta=0$ : (a)$\varDelta=0$ ; (b)$\varDelta=0.01$ ; (c)$\varDelta=0.5$ ; (d)$\varDelta=0.8$ -
本文研究了一维调制的p波超导体的拓扑量子相变. 在公度势调制下, p波超导的拓扑性质强烈地依赖于
$ \alpha $ 和$ \delta $ 的取值. 当$ b = 0 $ 时, 系统中存在特殊的相移$ \delta $ 使得Majorana零模的存在不依赖于公度势调制强度V. 通过计算发现当$ b\neq 0 $ 时, 在公度势调制系统中, 存在特殊相移使得Majorana零模不受调制强度影响的结果并不是普适的. 在非公度势($ \alpha = \left( \sqrt{5}-1\right) /2 $ )调制下, 计算了相移$ \delta = 0 $ 时系统的低能激发谱、$ Z_{2} $ 拓扑不变量以及逆参与率(IPR)等, 发现当p波配对强度$ 0 < \varDelta < 0.33 $ 时, 随着非公度势强度V的增加, 系统将经历从拓扑非平庸超导相到拓扑平庸超导相到拓扑平庸局域相的转变. 而当$ \varDelta > 0.33 $ 时, 随着非公度势强度V的增加, 系统经历拓扑非平庸相到拓扑平庸局域相的转变, 这与$ b = 0 $ 的结果一致.
-
本文研究了一维公度势和非公度势调制下的p波超导量子线系统的拓扑相变. 在公度势调制下, 通过计算
$Z_2$ 拓扑不变量确定系统的相图, 指出系统的拓扑相变强烈地依赖于调制参数$\alpha$ 和相移$\delta$ . 在非公度势调制下, 以$\alpha=(\sqrt{5}-1)/2$ ,$\delta=0$ 为例, 计算系统的低能激发谱、$Z_2$ 拓扑不变量以及逆参与率等, 发现p波配对强度$\varDelta\in(0,0.33)$ 时, 系统存在拓扑非平庸超导相, 拓扑平庸超导相和拓扑平庸局域相的转变. 而当p波配对强度$\varDelta > 0.33$ 时, 系统存在拓扑非平庸超导相和拓扑平庸局域相的转变.-
关键词:
- p波超导体 /
- Majorana费米子 /
- Z2拓扑不变量
We consider a one-dimensional p-wave superconducting quantum wire with the modulated chemical potential, which is described by$\hat{H}= \displaystyle\sum\nolimits_{i}\left[ \left( -t\hat{c}_{i}^{\dagger }\hat{c}_{i+1}+\Delta \hat{c}_{i}\hat{c}_{i+1}+ h.c.\right) +V_{i}\hat{n}_{i}\right]$ ,$V_{i}=V\dfrac{\cos \left( 2{\text{π}} i\alpha + \delta \right) }{1-b\cos \left( 2{\text{π}} i\alpha+\delta \right) }$ and can be solved by the Bogoliubov-de Gennes method. When$b=0$ ,$\alpha$ is a rational number, the system undergoes a transition from topologically nontrivial phase to topologically trivial phase which is accompanied by the disappearance of the Majorana fermions and the changing of the$Z_2$ topological invariant of the bulk system. We find the phase transition strongly depends on the strength of potential V and the phase shift$\delta$ . For some certain special parameters$\alpha$ and$\delta$ , the critical strength of the phase transition is infinity. For the incommensurate case, i.e.$\alpha=(\sqrt{5}-1)/2$ , the phase diagram is identified by analyzing the low-energy spectrum, the amplitudes of the lowest excitation states, the$Z_2$ topological invariant and the inverse participation ratio (IPR) which characterizes the localization of the wave functions. Three phases emerge in such case for$\delta=0$ , topologically nontrivial superconductor, topologically trivial superconductor and topologically trivial Anderson insulator. For a topologically nontrivial superconductor, it displays zero-energy Majorana fermions with a$Z_2$ topological invariant. By calculating the IPR, we find the lowest excitation states of the topologically trivial superconductor and topologically trivial Anderson insulator show different scaling features. For a topologically trivial superconductor, the IPR of the lowest excitation state tends to zero with the increase of the size, while it keeps a finite value for different sizes in the trivial Anderson localization phase.-
Keywords:
- p-wave superconductor /
- Majorana fermions /
- Z2 topological invariant
[1] Alicea J 2012 Rep. Prog. Phys. 75 076501
Google Scholar
[2] Beenakker C W J 2013 Ann. Rev. Con. Mat. Phys. 4 113
Google Scholar
[3] Wilczek F 2009 Nat. Phys. 5 614
Google Scholar
[4] Elliott S R, Franz M 2015 Rev. Mod. Phys. 87 137
Google Scholar
[5] Kitaev A Y 2001 Phys. Usp. 44 131
Google Scholar
[6] Mourik V, Zuo K, Frolov S M, Plissard S R, Bakkers E P A M, Kouwenhoven L P 2012 Science 336 1003
Google Scholar
[7] Chen J, Yu P, Stenger J, Hocevar M, Car D, Plissard S R, Bakkers E P A M, Stanescu T D, Frolov S M 2017 Sci. Adv. 3 e1701476
Google Scholar
[8] Albrecht S M, Higginbotham A P, Madsen M, Kuemmeth F, Jespersen T S, Nygard J, Krogstrup P, Marcus C M 2016 Nature 531 206
Google Scholar
[9] Deng M T, Vaitiekenas S, Hansen E B, Danon J, Leijnse M, Flensberg K, Nygard J, Krogstrup P, Marcus C M 2016 Science 354 1557
Google Scholar
[10] Zhang H, Liu C X, Gazibegovic S, Xu D, Logan J A, Wang G Z, N Loo van, Bommer J D S, Moor M W A d, Car D, Veld R L M O H, Veldhoven P J, Koelling S, Verheijen M A, Pendharkar M, Pennachio D J, Shojaei B, Lee J S, Palmstrøm C J, Bakkers E P A M, Sarma S D, Kouwenhoven L P 2018 Nature 556 74
Google Scholar
[11] Nadj-Perge S, Drozdov I K, Bernevig B A, Yazdani A 2013 Phys. Rev. B 88 020407
Google Scholar
[12] Nadj-Perge S, Drozdov I K, Li J, Chen H, Jeon S, Seo J, MacDonald A H, Bernevig B A, Yazdani A 2014 Science 346 602
Google Scholar
[13] Jeon S, Xie Y L, Li Jian, Wang Z J, Bernevig B A, Yazdani A 2017 Science 358 772
[14] Hell M, Leijnse M, Flensberg K 2017 Phys. Rev. Lett. 118 107701
Google Scholar
[15] Pientka F, Keselman A, Berg E, Yacoby A, Stern A, Halperin B I 2017 Phys. Rev. X 7 021032
[16] Fornieri A, Whiticar A M, Setiawan F, Marín E P, Asbjórn C C D, Keselman A, Gronin S, Thomas C, Wang T, Kallaher R, Gardner G C, Berg E, Manfra M J, Stern A, Marcus C M, Nichele F 2019 Nature 569 89
Google Scholar
[17] Cook A, Franz M 2011 Phys. Rev. B 84 201105
Google Scholar
[18] Sun H H, Zhang K W, Hu L H, Li C, Wang G Y, Ma H Y, Xu Z A, Gao C L, Guan D D, Li Y Y, Liu CH, Qian D, Zhou Y, Fu L, Li S C, Zhang F C, Jia J F 2016 Phys. Rev. Lett. 116 257003
Google Scholar
[19] Fu L, Kane C L 2008 Phys. Rev. Lett. 100 096407
Google Scholar
[20] Ivanov D A 2001 Phys. Rev. Lett. 86 268
Google Scholar
[21] Zhu S L, Shao L B, Wang Z D, Duan L M 2011 Phys. Rev. Lett. 106 100404
Google Scholar
[22] Lindner N H, Berg E, Refael G, Stern A 2012 Phys. Rev. X 2 041002
[23] Nayak C, Simon S H, Stern A, Freedman M, Sarma S D 2008 Rev. Mod. Phys. 80 1083
Google Scholar
[24] Jiang L, Kitagawa T, Alicea J, Akhmerov A R, Pekker D, Refael G, Cirac J I, Demler E, Lukin M D, Zoller P 2011 Phys. Rev. Lett. 106 220402
Google Scholar
[25] Hubener H, Sentef M A, Giovannini U D, Kemper A F, Rubio A 2017 Nat. Commun. 8 13940
Google Scholar
[26] Cheng Q, Pan Y, Wang H, Zhang C, Yu D, Gover A, Zhang H, Li T, Zhou L, Zhu S 2019 Phys. Rev. Lett. 122 173901
Google Scholar
[27] Cadez T, Mondaini R, Sacramento P D 2019 Phys. Rev. B 99 014301
Google Scholar
[28] Wang H Y, Zhuang L, Liu W M 2019 arXiv: 1910.10911 [cond-mat.mes-hall]
[29] Takata K, Notomi M 2018 Phys. Rev. Lett. 121 213902
Google Scholar
[30] Zhou L 2019 arXiv: 1911.11978 [cond-mat.mes-hall]
[31] Zeng Q B, Yang Y B, Xu Y 2019 arXiv: 1901.08060 [cond-mat.mes-hall]
[32] Okuma N, Sato M 2019 Phys. Rev. Lett. 123 097701
Google Scholar
[33] Ezawa M 2019 Phys. Rev. B 100 045407
Google Scholar
[34] Wu Y J, Liu H W, Liu J, Jiang H, Xie X C https://doi.org/10.1093/nsr/nwz189 [2020-1-8]
[35] Amorim C S, Ebihara K, Yamakage A, Tanaka Y, Sato M 2015 Phys. Rev. B 91 174305
Google Scholar
[36] Chen C Z, Xie Y M, Liu J, Lee P A, Law K T 2018 Phys. Rev. B 97 104504
Google Scholar
[37] Lang L J, Chen S 2012 Phys. Rev. B. 86 205135
Google Scholar
[38] Cai X M, Lang L J, Chen S, Wang Y P 2013 Phys. Rev. Lett. 110 176403
Google Scholar
[39] Hegde S S, Vishveshwara S 2016 Phys. Rev. Lett. 94 115166
[40] DeGottardi W, Thakurathi M, Vishveshwara S, Sen D 2013 Phys. Rev. B 88 165111
Google Scholar
[41] Wakatsuki R, Ezawa M, Tanaka Y, Nagaosa N 2014 Phys. Rev. B 90 014505
Google Scholar
[42] Lang L J, Cai X M, Chen S 2012 Phys. Rev. Lett. 108 220401
Google Scholar
[43] Aubry S, André G 1980 Ann. Isr. Phys. Soc. 3 133
[44] Ganeshan S, Pixley J H, Sarma S D 2015 Phys. Rev. Lett 114 146601
Google Scholar
[45] Sen A, Damle K, Moessner R 2012 Phys. Rev. B 86 205134
Google Scholar
[46] Zhu J X 2016 Bogoliubov-de Gennes Method and Its Applications (New Mexico: Springer) p3
[47] Gennes P G d (translated by Pincus P A) 1999 Superconductivity of Metals and Alloys (Boulder: Westview Press) pp137–160
[48] Lieb E, Schultz T, Mattis D 1961 Ann. Phys. 16 407
Google Scholar
[49] Kraus Y E, Lahini Y, Ringel Z, Verbin M, Zliberberg O 2012 Phys. Rev. Lett. 109 106402
Google Scholar
[50] Hofstadter D R 1976 Phys. Rev. B 14 2239
Google Scholar
[51] Zhou B, Shen S Q 2011 Phys. Rev. B 84 054532
Google Scholar
[52] Liu T, Yan H Y, Guo H 2017 Phys. Rev. B 96 174207
Google Scholar
[53] Akhmerov A R, Dahlhaus J P, Hassler F, Wimmer M, Beenakker C W J 2011 Phys. Rev. Lett. 106 057001
Google Scholar
[54] Fulga I C, Hassler F, Akhmerov A R, Beenakker C W J 2011 Phys. Rev. B 83 155429
Google Scholar
[55] Snyman I, Tworzydlo J, Beenakker C W J 2008 Phys. Rev. B 78 045118
Google Scholar
[56] Choy T P, Edge J M, Akhmerov A R, Beenakker C W J 2011 Phys. Rev. B 84 195442
Google Scholar
[57] Thouless D J 1974 Phys. Rep. 13 93
Google Scholar
[58] Kohmoto M 1983 Phys. Rev. Lett. 51 1198
Google Scholar
[59] Schreiber M 1985 J. Phys. C 18 2493
Google Scholar
[60] Hashimoto Y, Niizeki K, Okabe Y 1992 J. Phys. A 25 5211
Google Scholar
[61] Ingolda G L, Wobst A, Aulbach Ch, Hanggi P 2002 Eur. Phys. J. B 30 175
Google Scholar
-
图 4 (a)在开边界和周期性边界条件下最低激发态能量
$E_{1}$ 随准周期调制强度V的变化及其空间分布$\phi_1$ (b)和$\psi_1$ (c),$\alpha= $ $(\sqrt{5}-1)/2$ ,$b=0.5$ ,$\varDelta=0.2$ ,$L=2584$ Fig. 4. (a) The lowest excitation energies,
$E_{1}$ , varying with the quasi-periodic modulation amplitude, V, under OBC and PBC, respectively. The spatial distribution of the lowest excited state$\phi_1$ (b),$\psi_1$ (c).$\alpha=(\sqrt{5}-1) /2$ ,$b=0.5$ ,$\varDelta=0.2$ ,$L=2584$ .图 5 (a)
$Z_{2}$ 拓扑不变量随非公度势强度的变化; (b) V = 2时$\mathrm{IPR}_{1}$ 的标度分析; (c)$V=3$ 时$\mathrm{IPR}_{1}$ 的标度分析$b=0.5$ ,$\alpha= (\sqrt{5}-1) /2$ ,$\varDelta=0.2$ ,$L=2584$ Fig. 5. (a)
$Z_{2}$ topological invariant varying with the strength of the potential V; (b) the scaling of$\mathrm{IPR}_{1}$ $V=2$ ; (c) the scaling of$\mathrm{IPR}_{1}$ $V=3$ . Here,$\alpha= (\sqrt{5}-1)/2$ ,$b=0.5$ ,$\varDelta=0.2$ ,$L=2584$ .图 6
${\rm {IPR}}$ 随准周期调制强度V和本征能量$E_{n}$ 的变化$\alpha=(\sqrt{5}-1) /2$ ,$b=0.5, L=144, \delta=0$ (a)$\varDelta=0$ ; (b)$\varDelta=0.01$ ; (c)$\varDelta=0.5$ ; (d)$\varDelta=0.8$ Fig. 6.
${\rm {IPR}}$ varying with the amplitude of quasi-periodic modulation V and energy$E_{n}$ .$\alpha=(\sqrt{5}-1) /2$ ,$b=0.5, L=144$ ,$\delta=0$ : (a)$\varDelta=0$ ; (b)$\varDelta=0.01$ ; (c)$\varDelta=0.5$ ; (d)$\varDelta=0.8$ -
[1] Alicea J 2012 Rep. Prog. Phys. 75 076501
Google Scholar
[2] Beenakker C W J 2013 Ann. Rev. Con. Mat. Phys. 4 113
Google Scholar
[3] Wilczek F 2009 Nat. Phys. 5 614
Google Scholar
[4] Elliott S R, Franz M 2015 Rev. Mod. Phys. 87 137
Google Scholar
[5] Kitaev A Y 2001 Phys. Usp. 44 131
Google Scholar
[6] Mourik V, Zuo K, Frolov S M, Plissard S R, Bakkers E P A M, Kouwenhoven L P 2012 Science 336 1003
Google Scholar
[7] Chen J, Yu P, Stenger J, Hocevar M, Car D, Plissard S R, Bakkers E P A M, Stanescu T D, Frolov S M 2017 Sci. Adv. 3 e1701476
Google Scholar
[8] Albrecht S M, Higginbotham A P, Madsen M, Kuemmeth F, Jespersen T S, Nygard J, Krogstrup P, Marcus C M 2016 Nature 531 206
Google Scholar
[9] Deng M T, Vaitiekenas S, Hansen E B, Danon J, Leijnse M, Flensberg K, Nygard J, Krogstrup P, Marcus C M 2016 Science 354 1557
Google Scholar
[10] Zhang H, Liu C X, Gazibegovic S, Xu D, Logan J A, Wang G Z, N Loo van, Bommer J D S, Moor M W A d, Car D, Veld R L M O H, Veldhoven P J, Koelling S, Verheijen M A, Pendharkar M, Pennachio D J, Shojaei B, Lee J S, Palmstrøm C J, Bakkers E P A M, Sarma S D, Kouwenhoven L P 2018 Nature 556 74
Google Scholar
[11] Nadj-Perge S, Drozdov I K, Bernevig B A, Yazdani A 2013 Phys. Rev. B 88 020407
Google Scholar
[12] Nadj-Perge S, Drozdov I K, Li J, Chen H, Jeon S, Seo J, MacDonald A H, Bernevig B A, Yazdani A 2014 Science 346 602
Google Scholar
[13] Jeon S, Xie Y L, Li Jian, Wang Z J, Bernevig B A, Yazdani A 2017 Science 358 772
[14] Hell M, Leijnse M, Flensberg K 2017 Phys. Rev. Lett. 118 107701
Google Scholar
[15] Pientka F, Keselman A, Berg E, Yacoby A, Stern A, Halperin B I 2017 Phys. Rev. X 7 021032
[16] Fornieri A, Whiticar A M, Setiawan F, Marín E P, Asbjórn C C D, Keselman A, Gronin S, Thomas C, Wang T, Kallaher R, Gardner G C, Berg E, Manfra M J, Stern A, Marcus C M, Nichele F 2019 Nature 569 89
Google Scholar
[17] Cook A, Franz M 2011 Phys. Rev. B 84 201105
Google Scholar
[18] Sun H H, Zhang K W, Hu L H, Li C, Wang G Y, Ma H Y, Xu Z A, Gao C L, Guan D D, Li Y Y, Liu CH, Qian D, Zhou Y, Fu L, Li S C, Zhang F C, Jia J F 2016 Phys. Rev. Lett. 116 257003
Google Scholar
[19] Fu L, Kane C L 2008 Phys. Rev. Lett. 100 096407
Google Scholar
[20] Ivanov D A 2001 Phys. Rev. Lett. 86 268
Google Scholar
[21] Zhu S L, Shao L B, Wang Z D, Duan L M 2011 Phys. Rev. Lett. 106 100404
Google Scholar
[22] Lindner N H, Berg E, Refael G, Stern A 2012 Phys. Rev. X 2 041002
[23] Nayak C, Simon S H, Stern A, Freedman M, Sarma S D 2008 Rev. Mod. Phys. 80 1083
Google Scholar
[24] Jiang L, Kitagawa T, Alicea J, Akhmerov A R, Pekker D, Refael G, Cirac J I, Demler E, Lukin M D, Zoller P 2011 Phys. Rev. Lett. 106 220402
Google Scholar
[25] Hubener H, Sentef M A, Giovannini U D, Kemper A F, Rubio A 2017 Nat. Commun. 8 13940
Google Scholar
[26] Cheng Q, Pan Y, Wang H, Zhang C, Yu D, Gover A, Zhang H, Li T, Zhou L, Zhu S 2019 Phys. Rev. Lett. 122 173901
Google Scholar
[27] Cadez T, Mondaini R, Sacramento P D 2019 Phys. Rev. B 99 014301
Google Scholar
[28] Wang H Y, Zhuang L, Liu W M 2019 arXiv: 1910.10911 [cond-mat.mes-hall]
[29] Takata K, Notomi M 2018 Phys. Rev. Lett. 121 213902
Google Scholar
[30] Zhou L 2019 arXiv: 1911.11978 [cond-mat.mes-hall]
[31] Zeng Q B, Yang Y B, Xu Y 2019 arXiv: 1901.08060 [cond-mat.mes-hall]
[32] Okuma N, Sato M 2019 Phys. Rev. Lett. 123 097701
Google Scholar
[33] Ezawa M 2019 Phys. Rev. B 100 045407
Google Scholar
[34] Wu Y J, Liu H W, Liu J, Jiang H, Xie X C https://doi.org/10.1093/nsr/nwz189 [2020-1-8]
[35] Amorim C S, Ebihara K, Yamakage A, Tanaka Y, Sato M 2015 Phys. Rev. B 91 174305
Google Scholar
[36] Chen C Z, Xie Y M, Liu J, Lee P A, Law K T 2018 Phys. Rev. B 97 104504
Google Scholar
[37] Lang L J, Chen S 2012 Phys. Rev. B. 86 205135
Google Scholar
[38] Cai X M, Lang L J, Chen S, Wang Y P 2013 Phys. Rev. Lett. 110 176403
Google Scholar
[39] Hegde S S, Vishveshwara S 2016 Phys. Rev. Lett. 94 115166
[40] DeGottardi W, Thakurathi M, Vishveshwara S, Sen D 2013 Phys. Rev. B 88 165111
Google Scholar
[41] Wakatsuki R, Ezawa M, Tanaka Y, Nagaosa N 2014 Phys. Rev. B 90 014505
Google Scholar
[42] Lang L J, Cai X M, Chen S 2012 Phys. Rev. Lett. 108 220401
Google Scholar
[43] Aubry S, André G 1980 Ann. Isr. Phys. Soc. 3 133
[44] Ganeshan S, Pixley J H, Sarma S D 2015 Phys. Rev. Lett 114 146601
Google Scholar
[45] Sen A, Damle K, Moessner R 2012 Phys. Rev. B 86 205134
Google Scholar
[46] Zhu J X 2016 Bogoliubov-de Gennes Method and Its Applications (New Mexico: Springer) p3
[47] Gennes P G d (translated by Pincus P A) 1999 Superconductivity of Metals and Alloys (Boulder: Westview Press) pp137–160
[48] Lieb E, Schultz T, Mattis D 1961 Ann. Phys. 16 407
Google Scholar
[49] Kraus Y E, Lahini Y, Ringel Z, Verbin M, Zliberberg O 2012 Phys. Rev. Lett. 109 106402
Google Scholar
[50] Hofstadter D R 1976 Phys. Rev. B 14 2239
Google Scholar
[51] Zhou B, Shen S Q 2011 Phys. Rev. B 84 054532
Google Scholar
[52] Liu T, Yan H Y, Guo H 2017 Phys. Rev. B 96 174207
Google Scholar
[53] Akhmerov A R, Dahlhaus J P, Hassler F, Wimmer M, Beenakker C W J 2011 Phys. Rev. Lett. 106 057001
Google Scholar
[54] Fulga I C, Hassler F, Akhmerov A R, Beenakker C W J 2011 Phys. Rev. B 83 155429
Google Scholar
[55] Snyman I, Tworzydlo J, Beenakker C W J 2008 Phys. Rev. B 78 045118
Google Scholar
[56] Choy T P, Edge J M, Akhmerov A R, Beenakker C W J 2011 Phys. Rev. B 84 195442
Google Scholar
[57] Thouless D J 1974 Phys. Rep. 13 93
Google Scholar
[58] Kohmoto M 1983 Phys. Rev. Lett. 51 1198
Google Scholar
[59] Schreiber M 1985 J. Phys. C 18 2493
Google Scholar
[60] Hashimoto Y, Niizeki K, Okabe Y 1992 J. Phys. A 25 5211
Google Scholar
[61] Ingolda G L, Wobst A, Aulbach Ch, Hanggi P 2002 Eur. Phys. J. B 30 175
Google Scholar
计量
- 文章访问数: 4233
- PDF下载量: 182
- 被引次数: 0