姜雅金，姜潤澤，張禧征

（天津師范大學(xué)物理與材料科學(xué)學(xué)院，天津 300387）

近年來，隨著量子力學(xué)應(yīng)用研究的迅速發(fā)展，人們已經(jīng)能夠通過量子模擬手段，即利用人工可調(diào)控的量子系統(tǒng)，去模擬待研究的量子系統(tǒng)，探究當(dāng)前實(shí)驗(yàn)條件下難以觀測和解釋的物理現(xiàn)象.例如用光晶格與耦合腔陣列[1]等系統(tǒng)來模擬傳統(tǒng)凝聚態(tài)物理中的強(qiáng)關(guān)聯(lián)系統(tǒng)，最終在這些可控系統(tǒng)中達(dá)到制備、傳輸和控制量子態(tài)的目的.在制備、傳輸和控制量子態(tài)的過程中，人們已經(jīng)完成了對厄米體系中粒子在化學(xué)勢作用下動力學(xué)行為的基本運(yùn)算和分析.

當(dāng)Bender[2]提出狄拉克厄米性并非實(shí)數(shù)能譜和幺正的時間演化的必要條件時，研究者們發(fā)現(xiàn)系統(tǒng)哈密頓量中的相互作用勢還可以由實(shí)數(shù)域延擴(kuò)到復(fù)數(shù)域，且不改變其能譜的實(shí)數(shù)性.因此，非厄米系統(tǒng)引起了人們的廣泛關(guān)注.非厄米系統(tǒng)是用來描述萬物間相互聯(lián)系、相互作用的真實(shí)物理系統(tǒng)，將等效的非厄米哈密頓量引入后，能夠更加簡單直接地描述出外部環(huán)境對研究者所關(guān)心系統(tǒng)的影響，但非厄米系統(tǒng)中的單粒子波包在緊束縛模型上發(fā)生散射時存在共振穿透現(xiàn)象.所以，研究人員對“是否可以將這種共振穿透現(xiàn)象作為理論基礎(chǔ)，對非厄米系統(tǒng)中的緊束縛模型進(jìn)行特殊的結(jié)構(gòu)設(shè)計(jì)來實(shí)現(xiàn)單通性系統(tǒng)的制備”這一問題進(jìn)行了深入思考.

本研究主要以無限長的單鏈為基礎(chǔ)，通過調(diào)控非對稱耦合腔內(nèi)虛數(shù)勢[3]及腔內(nèi)格點(diǎn)間不等幅跳躍，觀察并分析單粒子波包入射至非厄米散射中心后在特定參數(shù)下實(shí)現(xiàn)的信息單向輸運(yùn)[4]，旨在為量子信息的輸運(yùn)提供理論基礎(chǔ).

1 方法

1.1 計(jì)算方法

量子力學(xué)以波函數(shù)公設(shè)、微觀粒子動力學(xué)公設(shè)和算符公設(shè)等幾個基本假設(shè)為起點(diǎn).由于厄米算符的本征值為實(shí)數(shù)時即為可觀測量，因此，算符公設(shè)必須用厄米算符來表示量子力學(xué)中的力學(xué)量.描述真實(shí)物理系統(tǒng)時，由于厄米哈密頓量用于表示系統(tǒng)的能量，其本征值為實(shí)數(shù)，因此必須滿足厄米性.從物理角度分析，滿足厄米性的哈密頓量描述的是封閉的物理系統(tǒng)，但無論生活中還是研究中，完全封閉的物理系統(tǒng)趨于理想化.因此，本研究僅對非厄米系統(tǒng)中的一維緊束縛模型展開分析.此外，經(jīng)電子源產(chǎn)生單電子脈沖，帶電粒子入射到一維緊束縛鏈上，將在每個格點(diǎn)能量最低的能帶上發(fā)生躍遷.根據(jù)緊束縛近似原理，總有單電子在近似情況下可通過能帶間的帶隙發(fā)生隧穿效應(yīng).同時，由于帶隙足夠大，處于低能級的帶電粒子很難躍遷至上一能級，因此可近似認(rèn)為此過程中沒有能量損耗.本研究將在該理論基礎(chǔ)上對單粒子波包進(jìn)行探究.

首先，假設(shè)一條長度為無窮大的一維緊束縛鏈，晶格均勻排列分布于該鏈上，在N=0的位置上存在虛數(shù)勢

對于非厄米系統(tǒng)，虛數(shù)勢是一種由于系統(tǒng)與外界環(huán)境耦合而產(chǎn)生的在位勢，一般用iγ來表示，同時格點(diǎn)間跳躍常數(shù)為-J，如圖1所示.

圖1 存在虛數(shù)勢的一維緊束縛鏈Fig.1 One-dimensional tight binding chain with imaginary potential

圖1為具有缺陷的一維緊束縛鏈，可采用Bethe Ansatz方法求解，得到系統(tǒng)的哈密頓量

設(shè)公式（2）中哈密頓量滿足本征值方程

公式（3）中：Ψ為哈密頓量H的本征態(tài)，滿足

由于勢的影響，位于N=0處的嘗試波函數(shù)可由相鄰格點(diǎn)間的公式給出.

將公式（4）代入公式（3）可得

令嘗試波函數(shù)

滿足公式（2）哈密頓量的非厄米系統(tǒng)中，相鄰格點(diǎn)間存在相互作用，且在虛數(shù)勢的影響下，得到波包由左向右入射后（規(guī)定向右為正方向），行進(jìn)過程中遇到N=0處的非厄米散射中心時，會發(fā)生奇特的散射現(xiàn)象，因此將嘗試波函數(shù)fj代入公式（9）—（11），得到左入射時單粒子波包在虛數(shù)勢影響下的透射系數(shù)tL和反射系數(shù)rL分別為

則透射幾率T和反射幾率R可以表示為

最后，在公式（2）哈密頓量的作用下，單粒子波包的能級光譜可以表示為

1.2 圖像分析

由公式（12）和公式（13）可知，反射系數(shù)與透射系數(shù)只依賴于波矢k、跳躍常數(shù)J與虛數(shù)勢實(shí)部γ，與系統(tǒng)的長度N以及勢所處的位置Nc無關(guān).為了能夠觀察單粒子波包的反射幾率和透射幾率受γ影響的變化規(guī)律，作數(shù)值解，如圖2所示，其中波矢k=π/2，跳躍常數(shù)J=1.由圖2可知，對于滿足公式（2）中哈密頓量的非厄米系統(tǒng)而言，系統(tǒng)中的散射中心帶給波包的影響非常顯著.傳播過程中，單粒子波包遇到虛數(shù)勢時會發(fā)生一定幾率的透射和反射，并且隨著γ的增大，反射幾率R與透射幾率T也增大，遞增程度如圖2所示.這說明滿足公式（1）的非厄米散射中心帶給系統(tǒng)的影響可以類比于從外界接入一個通道并注入能量，使散射中心處粒子的透射行為和反射行為增強(qiáng)（γ＞0），本研究將這樣的系統(tǒng)認(rèn)為是增益的.

圖2 γ＞0時反射幾率和透射幾率受γ影響的變化規(guī)律Fig.2 Changing law of reflection probability and transmission probability affected by imaginary potential γ when γ＞0

取系統(tǒng)長度N=100，虛數(shù)勢位于系統(tǒng)中央，即Nc=50處，利用數(shù)值法模擬寬度ω=1/3的單粒子波包在空間內(nèi)的演化規(guī)律，結(jié)果如圖3所示.從圖3可以看出，初始位置N0=10處的單粒子波包在Nc處受虛數(shù)勢的影響分裂為反射波與透射波，其中，波矢k=π/2，跳躍常數(shù)J=1，γ=0.8.

圖3 單粒子波包的演化示意圖（Nc=50，N0=10）Fig.3 Evolutionary schematic of the single-particle wave packet（Nc=50，N0=10）

此外，單粒子波包之所以能夠不散傳播，主要是因?yàn)檫\(yùn)動的高斯波包[5]在任意大小的φ值條件下，都能夠保形并且勻速地進(jìn)行傳輸，其中φ決定了群速的大小即波包演化速率正比于sin φ.然而φ值大小不同的高斯波包對波包的半寬也有不同的要求.從高效、局域傳輸?shù)慕嵌瓤紤]，取φ=π/2最為理想[6].

同理，利用Bethe Ansatz方法精確求解哈密頓量本征值及粒子的散射，可以得到單粒子波包在γ＜0時的透射系數(shù)tL和反射系數(shù)rL

由此可得反射幾率R和透射幾率T受γ影響的變化規(guī)律，結(jié)果如圖4所示，其中，波矢k=π/2，跳躍常數(shù)J=1.

圖4 γ＜0時反射幾率和透射幾率受γ影響的變化規(guī)律Fig.4 Changing law of reflection probability and transmission probability affected by imaginary potential γ when γ＜0

滿足公式（1）的非厄米散射中心帶給系統(tǒng)的影響可以類比于系統(tǒng)向外界搭建一個通道輸出能量，即向外界耗散掉一部分系統(tǒng)能量[7]，使得散射中心處粒子的透射行為和反射行為被削弱（γ＜0）.然而由圖4可以看出，隨著γ的遞增，雖然透射幾率呈衰減狀態(tài)，但反射幾率卻不減反增，這是因?yàn)椴ò趥鬏斶^程中產(chǎn)生了干涉效應(yīng)[8]，導(dǎo)致系統(tǒng)的反射幾率由于干涉相長而呈微弱的遞增現(xiàn)象，但整個系統(tǒng)還是可以看成是一個耗散系統(tǒng).

綜上所述，量子散射過程是測量儀器與被測系統(tǒng)形成理想糾纏的過程[9]，被測系統(tǒng)所處的狀態(tài)直接影響緊束縛模型中電子空間波函數(shù)的反射和透射.當(dāng)系統(tǒng)中接入一個γ＞0的通道使其轉(zhuǎn)化為增益系統(tǒng)時[10]，其反射幾率R和透射幾率T均隨γ的增大而增大，散射行為增強(qiáng)，可將其視為一個“源”[11]；同理，當(dāng)接入一個γ＜0的通道使其轉(zhuǎn)化為耗散系統(tǒng)時，散射行為被削弱，可將其視為一個“漏”[12].然而，當(dāng)離散量子系統(tǒng)中同時引入兩個正負(fù)相反的虛數(shù)勢，且格點(diǎn)之間被調(diào)制為不等幅跳躍結(jié)構(gòu)時，單粒子波包就產(chǎn)生了奇特的散射現(xiàn)象.

2 單通性系統(tǒng)

2.1 非對稱耦合散射中心

在研究單通性系統(tǒng)的構(gòu)建前，首先要對散射中心內(nèi)引入虛數(shù)勢和不等幅跳躍的一維緊束縛鏈進(jìn)行討論.假設(shè)一條長度為無窮大的單鏈，晶格均勻排列分布于鏈上，N=0和N=1位置上的虛數(shù)勢分別為iγ和-iγ（γ＞0），且雙虛數(shù)勢所在的格點(diǎn)間存在不等幅跳躍-α和-β，由此可視為將孤立的緊束縛系統(tǒng)引入了含有虛數(shù)勢的非對稱耦合散射中心，如圖5所示.

圖5 存在虛數(shù)勢和不等幅跳躍的一維緊束縛鏈Fig.5 One-dimensional tight binding chain with imaginary potential and unequal amplitude jumping

圖5中，系統(tǒng)的哈密頓量為

由于雙虛數(shù)勢的引入以及非對稱耦合結(jié)構(gòu)對系統(tǒng)的影響，同樣利用Bethe Ansatz方法可以在非對稱耦合散射中心處得到與公式（9）—（11）略有差別的下列公式

在滿足公式（19）哈密頓量的非厄米系統(tǒng)中，單粒子波包由左向右入射后（規(guī)定向右為正方向），遇到含有雙虛數(shù)勢的非對稱耦合散射中心時，波包會發(fā)生散射，將嘗試波函數(shù)fj代入公式（20）—（23），最終可以得到透射系數(shù)tL和反射系數(shù)rL

經(jīng)計(jì)算可知，反射系數(shù)與透射系數(shù)不僅依賴于波矢k、跳躍常數(shù)J與虛數(shù)勢實(shí)部γ三者間的比值，還與系統(tǒng)中非對稱耦合散射中心處不等幅跳躍α和β的大小密切相關(guān).圖6為反射幾率R與透射幾率T受參數(shù)α、β和γ影響后的變化規(guī)律.

由圖6（a）可知，隨著虛數(shù)勢中γ的增大，系統(tǒng)的透射幾率T逐漸增大，反射幾率R逐漸減小，透射行為更加明顯.說明單粒子波包從左向右入射至緊束縛鏈上，由于不等幅跳躍α=2和β=3，即α＜β，粒子透射的能力強(qiáng)于反射能力，所以表現(xiàn)為透射幾率隨著勢的增大而增大，反射幾率隨之減小.因此，在雙虛數(shù)勢±iγ作用下，系統(tǒng)中粒子的透射現(xiàn)象更明顯.

由圖6（b）可知，隨著不等幅跳躍α的增大，系統(tǒng)的透射幾率T逐漸減小，反射幾率R逐漸增大.對于從左入射的單粒子波包而言，系統(tǒng)中不等幅跳躍α表示著波包經(jīng)過散射中心后反射回到原路徑的困難程度，α的數(shù)值越大，單粒子波包反射的能力越強(qiáng)，粒子的反射現(xiàn)象更為明顯.

圖6 單粒子波包由左入射時，反射幾率和透射幾率在不同參數(shù)作用下的變化規(guī)律Fig.6 Changing law of reflection probability and transmission probability under the action of different parameters when the single-partical wave packet incident from left

對于從左入射的單粒子波包而言，系統(tǒng)中不等幅跳躍β表示波包經(jīng)過散射中心后按原路徑透射的困難程度，β的數(shù)值越大，單粒子波包透射的能力越強(qiáng)，透射幾率越大，反射幾率越小.然而由圖6（c）可以看出，隨著不等幅跳躍β的增大，透射幾率逐漸增大，反射幾率非但沒有隨之減小，反呈遞增的狀態(tài)，這是由于干涉效應(yīng)，即相位的干涉相長導(dǎo)致了反射振幅的異常變化，表現(xiàn)為反射行為和透射行為均增強(qiáng).

取系統(tǒng)長度N=100，且假設(shè)虛數(shù)勢位于系統(tǒng)中Nc=50和Nc=51處，利用數(shù)值方法模擬寬度為ω=1/2的單粒子高斯波包在空間內(nèi)的演化規(guī)律，其中，波矢k=π/2，γ=0.2，跳躍常數(shù)J=1，不等幅跳躍分別為α=2和β=3，結(jié)果如圖7所示.

圖7 單粒子波包的演化示意圖（Nc=50、51，N0=10）Fig.7 Evolutionary schematic of the single-particle wave packet（Nc=50、51，N0=10）

由圖7可知，隨著時間推移，波包從格點(diǎn)為N0=10的初始位置開始行進(jìn)，行進(jìn)過程中在Nc=50和Nc=51處遇到散射中心±iγ.此后，波包分裂為反射波和透射波，一部分按原定方向繼續(xù)前進(jìn)，另一部分按與入射方向相反的路徑返回.

當(dāng)單粒子波包從鏈端的反方向入射，即由右向左入射時，利用Bethe Ansatz方法繼續(xù)精確求解公式（19）哈密頓量的本征值及粒子散射.

首先，設(shè)嘗試波函數(shù)為

將嘗試波函數(shù)fj分別代入公式（20）—（23），最終得到單粒子波包右入射時在該非對稱耦合散射中心影響下的透射系數(shù)tR和反射系數(shù)rR

其反射幾率R和透射幾率T受不同參數(shù)α、β和γ影響的變化規(guī)律如圖8所示.

由圖8（a）可知，隨著虛數(shù)勢中γ的增大，系統(tǒng)的透射幾率T和反射幾率R均逐漸增大，且反射行為更明顯.說明單粒子波包從右向左入射至緊束縛鏈上，由于不等幅跳躍α=2和β=3，即α＜β，粒子的反射能力強(qiáng)于透射能力，表現(xiàn)為反射幾率處處大于透射幾率.因此，在雙虛數(shù)勢±iγ作用下的系統(tǒng)中，粒子反射現(xiàn)象更為明顯.

由圖8（b）可知，隨著不等幅跳躍α的增大，系統(tǒng)的透射幾率T和反射幾率R均逐漸增大，且透射行為更明顯.對于右入射的單粒子波包，系統(tǒng)中不等幅跳躍α表示波包經(jīng)過散射中心后按原路徑透射的困難程度，α的數(shù)值越大，單粒子波包透射的能力越強(qiáng)，透射幾率越大.然而圖8（b）中粒子的反射幾率非但沒有減小，反而越來越大，這是由于干涉效應(yīng)，即相位的干涉相長導(dǎo)致了反射振幅的異常變化，表現(xiàn)為反射和透射行為均增強(qiáng).

圖8 單粒子波包由右入射時，反射幾率和透射幾率在不同參數(shù)作用下的變化規(guī)律Fig.8 Changing law of reflection probability and transmission probability under the action of different parameters when the single-partical wave packet incident from right

由圖8（c）可知，隨著不等幅跳躍β的增大，透射幾率T減小，反射幾率R增大，反射行為更明顯.說明對于從右向左入射的單粒子波包，系統(tǒng)中不等幅跳躍β表示波包經(jīng)過散射中心后反射回到原路徑的困難程度，β的數(shù)值越大，單粒子波包反射的能力越強(qiáng)，反射幾率越大.

取系統(tǒng)長度N=100，假設(shè)虛數(shù)勢位于系統(tǒng)中Nc=50和Nc=51處，利用數(shù)值方法模擬寬度為ω=1/2的單粒子波包在空間內(nèi)的演化規(guī)律，其中波矢k=π/2一定，γ=0.2，跳躍常數(shù)J=1，不等幅跳躍分別為α=2和β=3，如圖9所示.

圖9 單粒子波包的演化示意圖（Nc=50、51，N0=90）Fig.9 Evolutionary schematic of the single-particle wave packet（Nc=50、51，N0=90）

圖9展示了波包隨著時間的演化過程，從格點(diǎn)為N0=90的初始位置開始行進(jìn)，行進(jìn)過程中在Nc=50和Nc=51處遇到了散射中心±iγ.此后，波包分裂為反射波和透射波，一部分按原定方向繼續(xù)前進(jìn)，另一部分按與入射方向相反的路徑返回.值得注意的是，對于一個相反方向入射的單粒子波包，其隨時間演化的行進(jìn)軌跡應(yīng)該與原定入射方向的波包行進(jìn)軌跡基本完全相反.

結(jié)合上述兩個相反方向入射的粒子波包在該系統(tǒng)下的散射行為，提出假設(shè)：對于一個反向入射粒子波包散射行為幾乎與原入射方向完全相反的系統(tǒng)，存在一種情況使其某一方向的反射幾率與相反方向的透射幾率均等于0，即rL=tR=0，rR=tL=0.兩組關(guān)系式中任意一種能夠得以驗(yàn)證，均可為實(shí)驗(yàn)中制備具有單向輸運(yùn)性即單通性的量子開關(guān)元件[13]提供理論基礎(chǔ).為驗(yàn)證這種類似二極管性質(zhì)的系統(tǒng)是否成立，理論上可以通過調(diào)制系統(tǒng)參數(shù)加以計(jì)算.

2.2 單通性系統(tǒng)的構(gòu)建

為了找到某一臨界點(diǎn)使系統(tǒng)具備單向輸運(yùn)性，將分兩種情況進(jìn)行討論：

（1）考慮滿足公式（19）哈密頓量的非厄米系統(tǒng).單粒子波包由任意方向入射至一維無限長的緊束縛鏈上，令系統(tǒng)內(nèi)α=0，β≠0，rL=tR=0，可得虛數(shù)勢γ與波矢k的關(guān)系為γ=ieik.接下來令波矢k=π/6，根據(jù)關(guān)系式γ=ieik得到利用數(shù)值模擬求得單粒子波包在特定參數(shù)下的演化規(guī)律，規(guī)定系統(tǒng)長度為N=200，跳躍常數(shù)J=1，波包寬度ω=1/20，單粒子波包雙向入射時的初始位置分別為N0=10和N0=190，演化結(jié)果如圖10所示.

圖10 單粒子波包在特定參數(shù)下演化的平面示意圖（ω=1/20）Fig.10 Plane schematic diagram of the evolution of a singleparticle wave packet under specific parameters（ω=1/20）

（2）繼續(xù)觀察滿足公式（19）哈密頓量的非厄米系統(tǒng).與第一種情況不同的是，單粒子波包由任意方向入射至一維無限長的緊束縛鏈上，此時令α≠0，β=0，rR=tL=0，可得γ與波矢k的關(guān)系為γ=ieik.接下來令波矢k=π/3，根據(jù)關(guān)系式γ=-ieik得到γ=-i/2+利用數(shù)值模擬求得高斯波包在特定參數(shù)下的演化規(guī)律，規(guī)定調(diào)整波包寬度為ω=1/5，其余條件不變，演化結(jié)果如圖11所示.

對比圖10和圖11可知，兩種不同情況下分別設(shè)置不等幅跳躍中的某一參數(shù)為0，均可以很好地找到滿足單向輸運(yùn)性的參數(shù)值及關(guān)系式.此外，由于非對稱耦合散射中心的誘導(dǎo)，單粒子波包發(fā)生共振穿透等行為說明了其單向輸運(yùn)性特點(diǎn)的存在.解析與數(shù)值模擬的吻合驗(yàn)證了前文提出關(guān)于單通性系統(tǒng)存在的假設(shè)，說明理論上確實(shí)存在這樣兩個臨界點(diǎn)，它們可以有效保證單粒子波包的單向輸運(yùn)性，實(shí)驗(yàn)中可以通過調(diào)節(jié)參數(shù)的大小來實(shí)現(xiàn)這種裝置的制備.

圖11 單粒子波包在特定參數(shù)下演化的平面示意圖（ω=1/5）Fig.11 Plane schematic diagram of the evolution of a singleparticle wave packet under specific parameters（ω=1/5）

對比現(xiàn)有成熟的理論基礎(chǔ)，研究人員已經(jīng)完成了對厄米體系中粒子在化學(xué)勢作用下動力學(xué)行為的基本運(yùn)算和分析，且可以通過操縱光晶格等系統(tǒng)來完成對量子態(tài)的調(diào)制.但在非厄米體系中，僅僅研究在位勢影響下非厄米哈密頓量作用在一維緊束縛鏈上所得波包的動力學(xué)行為還遠(yuǎn)遠(yuǎn)不夠.

3 結(jié)論

本文主要研究了在非對稱耦合散射中心作用下，單粒子波包在一維緊束縛鏈上的動力學(xué)行為，得到以下結(jié)論：

（1）對于給定的左向散射增益系統(tǒng)（代表外界一直在輸入，如+iγ視為一個增益，-iγ視為一個耗散），可以通過滿足公式（2）哈密頓量來計(jì)算模擬并觀察該系統(tǒng)內(nèi)單粒子波包的動力學(xué)行為.

（2）利用解析及數(shù)值模擬的方法，得到一維緊束縛鏈模型上單粒子波包在非對稱耦合散射中心的作用下隨時間演化的規(guī)律.由模擬結(jié)果可知，在虛數(shù)勢和不等幅跳躍的雙作用下，可以找到符合條件的臨界參數(shù)來制備具有單向輸運(yùn)性即單通性的半導(dǎo)體.

（3）非厄米系統(tǒng)中非對稱耦合結(jié)構(gòu)的散射中心可以誘導(dǎo)波包發(fā)生共振穿透等行為，使得信息的單向性輸運(yùn)在系統(tǒng)中成為可能，這在量子通信[14]領(lǐng)域有很好的應(yīng)用價值，如制備具有單通性的半導(dǎo)體裝置.

單粒子波包作為量子通信中的信號單元，在傳播過程中受到非對稱耦合散射中心的影響會發(fā)生不同的動力學(xué)行為變化，利用人工調(diào)制系統(tǒng)參數(shù)操縱粒子傳播的不同效果可以為人們研究量子固體器件提供理論基礎(chǔ)，并為今后信號傳輸?shù)亩鄻有匝芯刻峁┍憷?