張曉霞,李志堅(jiān)
(山西大學(xué) 理論物理研究所,山西 太原 030006)
一維分離時(shí)間量子行走的拓?fù)涮匦匝芯?/p>
張曉霞,李志堅(jiān)*
(山西大學(xué) 理論物理研究所,山西 太原 030006)
分離時(shí)間量子行走是由硬幣旋轉(zhuǎn)算符和依賴于硬幣態(tài)的條件平移算符相繼重復(fù)作用于系統(tǒng)初態(tài)的一種量子動(dòng)力學(xué)方案。文章在兩種作用算符中各引入一個(gè)可控參數(shù)研究系統(tǒng)的拓?fù)渥兓匦浴Mㄟ^(guò)計(jì)算系統(tǒng)的幾何相位,給出系統(tǒng)在可調(diào)參數(shù)空間的拓?fù)湎鄨D。進(jìn)而考慮空間不均勻的一維分離時(shí)間量子行走,分析其在邊界上的占有概率隨時(shí)間的變化,驗(yàn)證了拓?fù)浔Wo(hù)邊界態(tài)的存在。對(duì)比拓?fù)湎鄨D與擴(kuò)散系數(shù)在參數(shù)空間的變化特征,發(fā)現(xiàn)可通過(guò)擴(kuò)散系數(shù)揭露量子行走的拓?fù)涮匦浴?/p>
量子行走;拓?fù)湎嘧?幾何相位
量子行走是經(jīng)典隨機(jī)行走的量子對(duì)應(yīng)[1],描述了量子粒子在離散格點(diǎn)上的動(dòng)力學(xué)行為。在量子行走中由于量子態(tài)的相干疊加和干涉效應(yīng),使其與經(jīng)典隨機(jī)行走相比有顯著的不同,例如量子行走有比經(jīng)典隨機(jī)行走更快的傳播速度。量子行走包括連續(xù)時(shí)間量子行走和分離時(shí)間量子行走兩種形式。分離時(shí)間量子行走沒(méi)有勢(shì)能的概念,但可以通過(guò)硬幣算符的參數(shù)來(lái)等效量子勢(shì)[2-3]。因此,可以通過(guò)分離時(shí)間量子行走模擬各種各樣的物理過(guò)程,如光合作用[4],量子擴(kuò)散[5]和電子擊穿[6]等。實(shí)驗(yàn)上,人們已經(jīng)利用光格子中的超冷原子[6]、離子阱中的原子[7]、光學(xué)平臺(tái)中的光子[8]等一些物理系統(tǒng)實(shí)現(xiàn)了單個(gè)或多個(gè)粒子的量子行走。
另一方面,從整數(shù)量子霍爾效應(yīng)[9]的發(fā)現(xiàn)開(kāi)始,物理系統(tǒng)的拓?fù)浣Y(jié)構(gòu)對(duì)系統(tǒng)性質(zhì)的影響受到了人們的廣泛關(guān)注。近些年,理論和實(shí)驗(yàn)相繼發(fā)現(xiàn)了拓?fù)涑瑢?dǎo)體和拓?fù)浣^緣體等新型材料。然而具有這種拓?fù)湫再|(zhì)的天然材料很少,分離時(shí)間量子行走則為人工制造這種材料提供了可能[10]。通常定義的分離時(shí)間量子行走只具有一種非平庸的拓?fù)浣Y(jié)構(gòu)。后來(lái)人們提出了劈裂步分離時(shí)間量子行走的設(shè)計(jì)方案[11-12],這種方案引入了兩個(gè)硬幣算符,通過(guò)調(diào)節(jié)這兩個(gè)硬幣算符中的參數(shù)可以實(shí)現(xiàn)平庸拓?fù)湎嗪头瞧接雇負(fù)湎嗟霓D(zhuǎn)變。分離時(shí)間量子行走可以給出更為豐富的拓?fù)浣Y(jié)構(gòu),證實(shí)在具有不同拓?fù)湫再|(zhì)的兩種相的邊界處存在拓?fù)浔Wo(hù)的邊界態(tài)[13]。本文中,我們?cè)跅l件平移算符中引入一個(gè)可控參數(shù),和硬幣算符的一個(gè)參數(shù)形成二維的參數(shù)空間,在該空間中給出分離時(shí)間量子行走的拓?fù)湎鄨D,進(jìn)而討論拓?fù)湎噢D(zhuǎn)變和拓?fù)浔Wo(hù)的邊界態(tài)等物理特性。
一維分離時(shí)間量子行走包含位置和硬幣兩種自由度,其態(tài)矢量定義在位置空間Hp={|m〉,m∈Z}和硬幣空間Hc={|c〉,c=↑,↓}的直積希爾伯特空間中。分離時(shí)間量子行走的每一步時(shí)間演化U都由硬幣算符Rw和條件平移算符S相繼作用產(chǎn)生,即
(1)
其中Im為位置空間的單位算符。在標(biāo)準(zhǔn)的二維硬幣量子行走模型中,量子行走在一維位置空間中的兩個(gè)行走方向分別由兩個(gè)正交硬幣態(tài)決定。本文中,我們?cè)跅l件平移算符中引入一個(gè)參數(shù)φ,當(dāng)二維硬幣態(tài)發(fā)生翻轉(zhuǎn)時(shí),量子行走只以概率幅isinφ進(jìn)行部分行走,翻轉(zhuǎn)方向決定行走方向,另一部分則以cosφ為概率幅保持硬幣態(tài)和位置不變。相應(yīng)的條件平移算符可表示為
(2)
容易證明方程(2)滿足幺正性條件SS+=1。硬幣算符Rw也是一個(gè)作用于二維硬幣空間上的幺正算符,一般地,它可由繞任意軸w旋轉(zhuǎn)的旋轉(zhuǎn)算符Rw=e-iθ w·σ表示,其中σ為泡利矩陣,θ為旋轉(zhuǎn)角。本文中我們選定繞y軸旋轉(zhuǎn),則其可以表示為
(3)
該模型在實(shí)驗(yàn)上可以通過(guò)光學(xué)平臺(tái)加以實(shí)現(xiàn)[14-15],硬幣自由度和位置自由度分別編碼在光子的極化態(tài)和軌道角動(dòng)量態(tài)中。利用1/4波片調(diào)節(jié)光束的極化方向,實(shí)現(xiàn)量子硬幣的拋擲過(guò)程,此時(shí)的量子硬幣可以處于|↓〉態(tài)和|↑〉態(tài)的任意疊加態(tài)上。條件平移算符方程(2)可以通過(guò)一個(gè)稱為q-plate的光學(xué)元件來(lái)實(shí)現(xiàn)。它是一個(gè)結(jié)構(gòu)特殊的可以產(chǎn)生相位延遲的液晶雙折射光學(xué)鏡片,并具有拓?fù)潆姾蓂。在一個(gè)橫截面內(nèi),它的光軸不是平行指向。由于雙折射效應(yīng),它的功能可以分為兩部分。一部分是它以一定的概率使得入射光子的每個(gè)圓極化分量發(fā)生螺旋性翻轉(zhuǎn),使得自旋角動(dòng)量和軌道角動(dòng)量相互轉(zhuǎn)變,從而改變了入射光子的角動(dòng)量態(tài)。由此可以利用q-plate設(shè)計(jì)產(chǎn)生光子的自旋軌道耦合,實(shí)現(xiàn)在光子極化態(tài)控制下的軌道角動(dòng)量的改變,這在數(shù)學(xué)上的表示對(duì)應(yīng)方程(2)的前兩項(xiàng)。另一部分是它使剩下概率的光以尋常光直接通過(guò),既不改變?nèi)肷涔庾拥淖孕莿?dòng)量也不改變光子的軌道角動(dòng)量,這在數(shù)學(xué)上的表示對(duì)應(yīng)方程(2)的后兩項(xiàng)。入射光子的角動(dòng)量態(tài)改變的概率與q-plate產(chǎn)生的延遲相位有關(guān),這一相位由方程(2)中引入的參數(shù)φ描述,它可以通過(guò)電場(chǎng)、溫度等外加方式來(lái)進(jìn)行調(diào)控。
在量子行走中,雖然演化算符是由幺正算符S和Ry直接定義而來(lái),但其可以等效為有效哈密頓量Heff的形式,即一步演化算符可表示為U=e-iHeff。經(jīng)過(guò)t步,量子行走演化為
因?yàn)閠取分立的值,因此分離時(shí)間量子行走可以看作一個(gè)由有效哈密頓量描述的系統(tǒng)進(jìn)行時(shí)間演化時(shí)的閃頻模擬,而且有效哈密頓量Heff的本征值以2π為周期。將方程(2)作傅里葉變換,S在動(dòng)量空間可表示為
(4)
相應(yīng)地,動(dòng)量空間中的有效哈密頓量為
(5)
其中nφ,θ(k)為單位矢量,表示動(dòng)量為k的量子行走在硬幣空間的本征態(tài)的極化方向;Eφ,θ(k)是能量本征值,與動(dòng)量k滿足色散關(guān)系
cos[Eφ,θ(k)]=cosθcosφ-sinksinθsinφ.
(6)
至此,我們可以清楚地看到量子行走可以等效為具有自旋軌道耦合的有效哈密頓量形式,這為實(shí)現(xiàn)非平庸的拓?fù)浣Y(jié)構(gòu)提供了必要的條件。
由于量子行走的坐標(biāo)和時(shí)間都為離散變量,而且具有平移不變性,其動(dòng)量和能量都是具有周期性結(jié)構(gòu)的準(zhǔn)動(dòng)量和準(zhǔn)能量。不失一般性,我們將二者都定義在-π到π區(qū)間之內(nèi)。給定一個(gè)動(dòng)量k,兩個(gè)硬幣自由度使得方程(6)包含正負(fù)兩支能帶,調(diào)節(jié)參數(shù)θ和φ可以使兩支能帶之間的能隙打開(kāi)或閉合。圖1(a)給出當(dāng)θ和φ取不同值時(shí),量子行走的能量Eφ,θ(k)隨動(dòng)量k的變化圖。當(dāng)|sinφ|≠|(zhì)sinθ|時(shí),正負(fù)兩支能量曲線關(guān)于Eφ,θ(k)=0對(duì)稱,而且不會(huì)相交,即兩支能帶之間存在一定的能隙;當(dāng)|sinφ|=|sinθ|≠1時(shí),兩支能帶在k=-π/2或π/2處閉合,簡(jiǎn)并能量Eφ,θ(k)=0或π。特別地,如果|θ|=|φ|=π/2,量子能量曲線隨能量k成線性變化,而且同時(shí)在k=±π/2處閉合。對(duì)應(yīng)每支能帶,量子行走的群速分別為
(7)
從方程(6)和(7)中可以看出,當(dāng)θ,φ=0或π時(shí),Eφ,θ(k)=0或π,此時(shí)群速?zèng)]有定義。取不同的θ和φ值,對(duì)應(yīng)能量正支的群速隨動(dòng)量k的變化如圖1(b)所示。從圖可以看出,如果兩支能帶之間存在間隙,則群速會(huì)連續(xù)變化,形成一條平滑的曲線;但是如果能帶發(fā)生閉合,則群速會(huì)出現(xiàn)不連續(xù)的情況,特別是當(dāng)|θ|=|φ|=π/2時(shí),群速除了在k=±π/2處沒(méi)有定義外,其絕對(duì)值為常數(shù)1。
Fig.1 Variations of (a) energy and (b) group velocity with momentum are shown for different values of θ and φ in discrete time quantum walk圖1 參數(shù)θ和φ取不同的值時(shí),分離時(shí)間量子行走的(a)能量和(b)群速度隨動(dòng)量的變化情況
當(dāng)|sinφ|≠|(zhì)sinθ|時(shí),隨著k的變化,系統(tǒng)的兩支能帶之間總是存在能隙,在不同的取值區(qū)域,雖然系統(tǒng)的能帶結(jié)構(gòu)相似,但相應(yīng)態(tài)矢量的極化方向不同,會(huì)表現(xiàn)出不同拓?fù)涮匦?。下面我們用幾何相?/p>
(8)
作為拓?fù)洳蛔兞坎⒔o出系統(tǒng)的相圖,以表征系統(tǒng)的拓?fù)涮匦浴榱朔奖阌?jì)算,我們將演化算符U繞y軸旋轉(zhuǎn)角度π-θ,使得單位向量nφ,θ(k)旋轉(zhuǎn)至x-y平面。這時(shí)演化算符變?yōu)?/p>
U′=ei(π-θ)σy/2Ue-i(π-θ)σy/2=(cosθcosφ-sinθsinφsink)I-isinφcoskσx-i(sinθcosφ+cosθsinφsink)σy,
(9)
其本征矢可解析地表示為
(10)
其中δ(k)滿足tanδ(k)=(sinφcosk)/(sinθcosφ+cosθsinφsink)。將方程(10)代入(8)可得幾何相位
(11)
上式表明,在參數(shù)θ和φ的不同取值區(qū)域,該量子行走的幾何相位或者為±π,或者為0,它們分別對(duì)應(yīng)著系統(tǒng)具有非平庸和平庸的拓?fù)湎?。在圖2中,我們給出幾何相位,即拓?fù)洳蛔兞侩S著θ和φ變化的拓?fù)湎鄨D,實(shí)線和虛線分別表示兩支能帶在Eφ,θ(k)=0和Eφ,θ(k)=π處閉合。當(dāng)能帶間隙閉合時(shí),有效哈密頓量的本征態(tài)是極化方向相反的兩個(gè)硬幣態(tài)的任意疊加態(tài),此時(shí)拓?fù)湎鄾](méi)有定義。在兩支能帶閉合的直線兩側(cè),系統(tǒng)具有不同的拓?fù)湎唷?/p>
Fig.2 Diagram is the topological phase diagram of discrete-time quantum walk in parameter space (θ,φ),±π and 0 in figure are the values of Zak phase in corresponding parameter region圖2 分離時(shí)間量子行走在參數(shù)(θ,φ)空間中的拓?fù)湎鄨D,圖中標(biāo)記的±π和0是相應(yīng)區(qū)域的幾何相的大小
根據(jù)體-邊界對(duì)應(yīng)原理,在具有不同拓?fù)浣Y(jié)構(gòu)的兩種相的邊界處會(huì)出現(xiàn)拓?fù)浔Wo(hù)的邊界態(tài)。為了研究這種邊界態(tài),下面我們考慮一個(gè)空間不均勻的量子行走,也就是說(shuō),在方程(1)-(3)中,我們固定平移參數(shù)φ,當(dāng)位置坐標(biāo)m≤0時(shí),選取硬幣參數(shù)θ=θ-,而當(dāng)位置坐標(biāo)m>0時(shí),選取θ=θ+。當(dāng)(θ-,φ)和(θ+,φ)在圖2中所對(duì)應(yīng)的點(diǎn)具有不同的幾何相位時(shí),m=0兩側(cè)區(qū)域具有不同的拓?fù)湎?。容易證明這種空間不均勻的量子行走,其演化算符U仍然滿足幺正性。邊界態(tài)是否存在可以通過(guò)粒子在邊界處的占有概率來(lái)判斷,如果粒子在邊界處的占有概率隨著演化時(shí)間的增加不再發(fā)生變化,并保持恒定的占有率,則說(shuō)明邊界處存在邊界態(tài)。這是因?yàn)殡S著演化步數(shù)的增加,量子行走一般會(huì)彈道擴(kuò)散,但如果邊界處存在束縛態(tài),那么量子行走就會(huì)以一定的概率留在邊界處,而且這一概率不會(huì)隨時(shí)間改變,其大小由初態(tài)在束縛態(tài)上的投影決定。如果邊界處不存在束縛態(tài),則量子行走在邊界處的概率會(huì)隨著演化時(shí)間的增加趨于0.
(a)(b) θ-=-π/4,θ+=5π/16;(c)(d) θ-=-3π/5,θ+=5π/16;(e)(f) θ-=π/8,θ+=3π/4Fig.3 Subfigures (a), (c) and (e) show the probability distribution in position space when t=50 and the subfigures (b), (d)and (f) show the probability distribution in position space as a function of time t. The other parameters are chosen as φ=5π/12圖3 (a),(c)和(e)為當(dāng)t=50時(shí),空間不均勻的量子行走在位置空間中的概率分布圖;(b),(d)和(f)為空間不均勻的量子行走的概率分布隨時(shí)間和空間變化的二維投影圖,其中參數(shù)φ=5π/12
從上文已經(jīng)看到,通過(guò)調(diào)節(jié)方程(1)-(3)中的平移參數(shù)φ和硬幣參數(shù)θ,可以使得量子行走的拓?fù)涮匦园l(fā)生改變。研究這種拓?fù)滢D(zhuǎn)變與可觀察效應(yīng)的關(guān)系具有重要意義。下面我們研究擴(kuò)散系數(shù)對(duì)這些參數(shù)的依賴關(guān)系,從而來(lái)判斷系統(tǒng)所具有的拓?fù)涮匦?。擴(kuò)散系數(shù)是物質(zhì)的物理性質(zhì)之一,表示粒子的擴(kuò)散能力,某一時(shí)刻的擴(kuò)散系數(shù)D[16]定義為
(12)
Fig.4 Diffusion coefficient as a function of parameter angle φ are plotted for different values of θ圖4 選取不同的參數(shù)θ,擴(kuò)散系數(shù)隨參數(shù)φ的變化圖
分離時(shí)間量子行走是操控一個(gè)粒子的量子動(dòng)力學(xué)方案,我們通過(guò)在分離時(shí)間量子行走中的旋轉(zhuǎn)算符和條件平移算符中各引入一個(gè)可控參數(shù),研究其隨著這些參數(shù)變化的拓?fù)渥兓匦?。由于該分離時(shí)間量子行走的準(zhǔn)能量與準(zhǔn)動(dòng)量間具有獨(dú)特的色散關(guān)系,可以通過(guò)改變引進(jìn)的參數(shù)使得兩支能帶間的能隙打開(kāi)或閉合,相應(yīng)地,群速在能隙閉合時(shí)會(huì)發(fā)生突變。通過(guò)計(jì)算系統(tǒng)的幾何相位,發(fā)現(xiàn)該系統(tǒng)在不同參數(shù)區(qū)域有幾何相位為0的平庸拓?fù)湎嗪蛶缀蜗辔粸椤捆械姆瞧接雇負(fù)湎?由此給出系統(tǒng)在可調(diào)參數(shù)空間的拓?fù)湎鄨D。在此基礎(chǔ)上,我們考慮了一個(gè)空間不均勻的一維分離時(shí)間量子行走,也就是改變可調(diào)參數(shù)使得具有相同或不同拓?fù)湎嗟膬煞N結(jié)構(gòu)相連接,通過(guò)計(jì)算邊界上量子行走的占有概率隨著時(shí)間的變化情況,驗(yàn)證了拓?fù)浔Wo(hù)邊界態(tài)的存在。對(duì)比拓?fù)湎鄨D與擴(kuò)散系數(shù)在參數(shù)空間的變化特征,發(fā)現(xiàn)可觀察的擴(kuò)散系數(shù)在不同的拓?fù)湎鄥^(qū)域會(huì)表現(xiàn)出兩種完全不同的變化形式,一種是拋物型變化,一種是平直不變,二者的連接處是一個(gè)拐點(diǎn),剛好對(duì)應(yīng)兩支能帶的閉合點(diǎn),由此可以通過(guò)擴(kuò)散系數(shù)揭露量子行走的拓?fù)涮匦浴?/p>
[1]AharonvY,DavidovichL,ZaguryN.QuantumRandomWalks[J].Phys Rev A,1993,48:1687-1690.DOI:10.1103/PhysRevA. 48.1687.
[2]LiZJ,IzaacJA,WangJB.Position-Defect-InducedReflection,Trapping,Transmission,andResonanceinQuantumWalks[J].Phys Rev A,2013,87:012314.DOI:10.1103/PhysRevA.87.012314.
[3] 梁樹(shù)琴,李志堅(jiān).粒子在非均勻環(huán)上的量子行走[J].山西大學(xué)學(xué)報(bào):自然科學(xué)版,2014,37:222-226.DOI:10.1341/j.cnki.shanxi.univ(nat.sci.).2014.02.011.
[4]SensionRJ.QuantumPathtoPhotosynthesis[J].Nature,2007,446:740-741.DOI:10.1038/446740a.
[5]GodoyS,FujitaS.AQuantumRandom-WalkModelforTunnelingDiffusionina1DLattice.AQuantumCorrectiontoFick’sLaw[J].J Chem Phys,1992,97:5148-5154.DOI:10.1063/1.463812.
[6]KarskiM,F?rsterL,ChoiJM,et al.QuantumWalkinPositionSpacewithSingleOpticallyTrappedAtoms[J].Science,2009,325:174-177.DOI:10.1126/science.1174436.
[7]Z?hringerF,KirchmairG,GerritsmaR,et al.RealizationofaQuantumWalkwithOneandTwoTrappedIons[J].Phys Rev Lett,2010,104:100503.DOI:10.1103/PhysRevLett.104.100503.
[8]SchreiberA,CassemiroKN,PotocekV,et al.PhotonsWalkingtheLine:AQuantumWalkwithAdjustableCoinOperations[J].Phys Rev Lett,2010,104:050502.DOI:10.1103/PhysRevLett.104.050502.
[9]KlitzingKV,DordaG,PepperM.NewMethodforHigh-AccuracyDeterminationoftheFine-StructureConstantBasedonQuantizedHallResistance[J].Phys Rev Lett,1980,45:494-497.DOI:10.1103/PhysRevLett.45.494.
[10]KitagawaT.TopologicalPhenomenainQuantumWalks:ElementaryIntroductiontothePhysicsofTopologicalPhases[J].Quantum Inf Process,2012,11:1107-1148.DOI:10.1007/s11128-012-0425-4.
[11]ObuseH,AsbthJK,NishimuraY,et al.UnveilingHiddenTopologicalPhasesofaOne-DimensionalHadamardQuantumWalk[J].Phys Rev B,2015,92:045424.DOI:10.1103/PhysRevB.92.045424.
[12]RakovszkyT,AsbothJK.Localization,Delocalization,andTopologicalPhaseTransitionsintheOne-DimensionalSplit-StepQuantumWalk[J].Phys Rev A,2015,92:052311.DOI:10.1103/PhysRevA.92.052311.
[13]GrohT,BrakhaneS,AltW,et al.RobustnessofTopologicallyProtectedEdgeStatesinQuantumWalkExperimentswithNeutralAtoms[J].Phys Rev A,2016,94:013620.DOI:10.1103/PhysRevA.94.013620.
[14]CardanoF,MaffeiM,MassaF,et al.StatisticalMomentsofQuantum-WalkDynamicsRevealTopologicalQuantumTransitions[J].Nat Commun,2016,7:11439.DOI: 10.1038/ncomms11439.
[15]GoyalSK,RouxFS,ForbesA,et al.ImplementingQuantumWalksUsingOrbitalAngularMomentumofClassicalLight[J].Phys Rev Lett,2013,110:263602.DOI:10.1103/PhysRevLett.110.263602.
[16]RomanelliA,SiriR,AbalG,et al.DecoherenceintheQuantumWalkontheLine:I-Measurements[J].Phys A,2004,347:137-152.DOI:10.1016/j.physa.2004.08.070.
Topological Property of One-Dimensional Discrete-Time Quantum Walk
ZHANG Xiaoxia,LI Zhijian*
(InstituteofTheoreticalPhysics,ShanxiUniversity,Taiyuan030006,China)
Discrete-time quantum walk is a quantum dynamic scheme, which involves the coin flipping and conditional position shift based on the outcome of the coin state.We study the topological properties of one-dimensional discrete-time quantum walk by introducing two controllable parameters in coin operator and conditional shift operator respectively. The topological phase diagram in this parameter space is present by calculating the geometric phase. Furthermore, we consider an inhomogeneous quantum walk to demonstrate the existing of the topologically protected edge states. When two parts with different topological phase joint together, the occupancy probability at the boundary keeps a constant after lone time evolution. By comparing the topological phase diagram with the change of the diffusion coefficient in the parameter space, we find that the diffusion coefficient can reveal the topological properties of quantum walk.
quantum walk;topological phase transition;geometric phase
10.13451/j.cnki.shanxi.univ(nat.sci.).2017.01.014
2016-12-05;
2016-12-26
山西省回國(guó)留學(xué)人員科研資助項(xiàng)目(2015-012);山西省面上自然科學(xué)基金(201601D011009)
張曉霞(1990-),女,山西大同人,碩士研究生, E-mail :928970149@qq.com
*通信作者:李志堅(jiān)(LI Zhijian),E-mail:zjli@sxu.edu.cn
O413
A
0253-2395(2017)01-0100-06