删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

一维化学势调制的p波超导体中的拓扑量子相变

本站小编 Free考研考试/2021-12-29
摘要:本文研究了一维公度势和非公度势调制下的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拓扑不变量

English Abstract


全文HTML

--> --> -->
1.引　言
早在上世纪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波超导量子线系统中的拓扑量子相变.
2.理论模型与方法
考虑一维具有(准)周期调制的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 $.
3.结果分析与讨论
2
3.1.周期调制的p波超导线
-->

3.1.周期调制的p波超导线

这一小节讨论$ \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, δ = 0
Figure1. 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 = 2584
Figure2. 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$
Figure3. 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.
2
3.2.准周期调制的p波超导线
-->

3.2.准周期调制的p波超导线

上一节中的结果表明$ 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$
Figure4. (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$
Figure5. (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$
Figure6. ${\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$

4.结　论
本文研究了一维调制的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 $的结果一致.
相关话题/系统 计算 空间 拓扑 化学
  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效，10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材，产品类型包含电子书、题库、全套资料以及视频，无论您是考研复习、考证刷题，还是考前冲刺等，不同类型的产品可满足您学习上的不同需求。 ...
    考试优惠券 本站小编 Free壹佰分学习网 2022-09-19
  • 弱Soret效应混合流体对流系统的分岔与非线性演化
    摘要:混合流体Rayleigh-Bénard(RB)对流是研究非平衡耗散系统的自组织斑图及非线性动力学特性的典型模型.本文利用高精度数值方法模拟了底部均匀加热的矩形腔体中混合流体RB对流,研究了具有极微弱Soret效应(分离比$\psi=-0.02$)的混合流体对流的分岔特性及斑图的形成和演化,给出 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 量子点-Su-Schrieffer-Heeger原子链系统的电子输运特性
    摘要:Su-Schrieffer-Heeger(SSH)原子链是典型的具有拓扑边缘态的一维系统,并且已在光子和冷原子系统中实验实现.本文在紧束缚近似下,利用传输矩阵方法研究了量子点-SSH原子链系统的电子输运特性,这里,量子点的作用是调节SSH原子链与电极的隧穿耦合强度.当量子点与SSH原子链弱耦合 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 一种新型的液闪阵列成像屏空间分辨特性
    摘要:针对低强度射线成像,自主研制了一种像元为0.1mm高探测效率的液闪阵列屏.为此,基于倾斜刀口边缘响应的测量原理,建立了理论模拟方法和实验研究方法,对该液闪阵列屏开展了空间分辨性能研究.通过理论模拟,给出了液闪阵列屏在14MeV中子和1.25MeV伽马射线激发下的调制传递函数,并与像元为0.1, ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 气力提升系统气液两相流数值模拟分析
    摘要:污水处理、油田采油、液态金属冷却反应堆和磁流体动力转换器等领域采用气力提升系统有其显著优势.由于不同液体介质与气体介质密度对气力提升系统性能影响较大,因此本文基于Fluent仿真软件,采用欧拉模型、k-ω剪切应力输运湍流模型数值模拟了氮气-水、氮气-煤油、氮气-水银及空气-水、氩气-水、氮气- ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 一种结合图像复原技术的自适应光学系统控制方法
    摘要:在天文高分辨成像领域,自适应光学校正和事后图像复原都必不可少,但传统的自适应光学系统控制方法以提升光学成像质量为目的,并未考虑图像复原环节,因此,研究一种结合两者以获得高质量复原图像为目标的控制方法具有重要意义.本文对传统自适应光学技术结合事后图像解卷积的方法进行了分析,阐述了其存在的缺陷.首 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 腔光力系统制备微波非经典态研究进展
    摘要:腔光力系统作为一种新型的混合量子系统,因其超强耦合度、低温超导条件下极低的噪声、较长的相干时间等优势而成为被广受关注的量子实验平台.本文简要介绍腔光力学及腔光力系统基本原理,对常见腔光力系统进行分类,详细介绍利用广义腔光力系统进行微波非经典量子态制备的相关进展,对其性能优势和待解决问题进行分析 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • Ar原子序列双光双电离产生光电子角分布的理论计算
    摘要:基于多组态Dirc-Fock方法和密度矩阵理论,给出了原子序列双光双电离光电子角分布的计算表达式,发展了相应的计算程序.利用该程序对Ar原子3p壳层序列双光双电离过程进行了理论研究,给出了光电离的总截面、磁截面、剩余离子取向以及光电子角分布的各向异性参数与入射光子能量的函数关系.结果显示在光电 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 脉冲强磁场下的电极化测量系统
    摘要:多铁性材料是当前物质科学研究的热点,具有重要的科学研究意义和应用前景.低温和强磁场实验环境为研究多铁性材料提供了一种有效途径.脉冲强磁场下的电极化测量系统能实现最高磁场强度60T、最低温度0.5K的铁电特性测量.该系统采用热释电方法,具有磁场强度高、控温范围广、转角测量等特点,可用于强磁场下的 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于新的五维多环多翼超混沌系统的图像加密算法
    摘要:本文提出了一种基于新的五维多环多翼超混沌系统的数字图像加密方法.首先,将明文图像矩阵和五条混沌序列分别通过QR分解法分解成一个正交矩阵和一个上三角矩阵,将混沌系统产生的五条混沌序列分别通过LU分解法分解成一个上三角矩阵和一个下三角矩阵,分别将两个上三角矩阵和一个下三角矩阵相加,得到五个离散后的 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于量子点接触的开放双量子点系统电子转移特性
    摘要:基于量子点接触探测器(QPC)理论上研究了双量子点(DQD)系统在耗散环境和纯退相环境影响下的电子转移特性.结果表明,耗散环境中探测器导致的退相干会增大平均电流和Fanofactor随时间演化的值,并观察到量子芝诺效应的存在.在对称的DQD情况下,弛豫减小了平均电流随时间演化的震荡振幅.在非对 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29

Free考研考试FreeKaoYan.Com
欢迎来到Free考研考试，"为实现人生的Free而奋斗"

© 2020 FreeKaoYan! . All rights reserved.