用具有負介電常數的模擬光子晶體同軸系統獲得負群速

姜榮，黃志洵

(中國傳媒大學信息工程學院，北京 100024)

1 引言

Sommerfeld和Brillouin[1]在其經典的波速理論中指出在反常色散條件下，群速可能大于真空中光速，甚至為負。這引起了科學家們對反常色散媒質中的群速研究產生了極大的興趣。1970年Garrett[2]等人的研究表明：當頻率為f的電磁波通過非吸收但色散的媒質傳播時，當媒質具有接近頻率f的吸收線時就會出現反常色散的情況，此時群速可以大于光速，甚至為負，并且通過計算還表明對于吸收媒質來講群速的概念同樣具有意義。在1982年S.Chu and S.Wong[3]利用Garrett[2]提出的近吸收線的反常色散媒質首次測量到負群速。除此以外，科學家還利用增益線的方法實現負群速傳播，其中尤其是2000年王力軍[4]應用輔助增益線的方法實現反常色散首次在波形不失真的條件下觀測到負群速，這使得觀測到的群速更加具有意義。除了在這些有源的媒質實現電磁波的反常色散外，科學家們同樣對無源媒質中電磁波的反常色散傳播進行研究[5，6]，例如利用同軸晶體結構實現電磁波的反常色散傳播[7-10]，其中文獻[7，8]觀測到超光速的群速，[9，10]觀測到負的群速；利用波導系統實現類Fabry-Perot干涉儀達到電磁波的反常色散傳播并觀測負群速等[11]。

本文通過使用不同阻抗的同軸線周期相連形成模擬光子同軸晶體的同軸系統，由于阻抗不匹配，在某些頻率形成阻帶，在這些區域內出現反常色散的現象，通過計算得出反常色散區域處的超光速群速。然后假設其具有負折射率，計算電磁波的反常色散傳播，獲得了電磁波的負群速傳播。

2 理論分析

使用圖1所示不同阻抗的同軸電纜周期相連，得到模擬CPC的結構，來實現電磁波的反常色散性傳播。為了計算色散需要確定折射率實部nr與相移φ的關系。由于在CPC中，不同阻抗的同軸電纜周期相連，在交界處不連續，存在入射波和反射波，因此電磁波的相位和振幅都發生變化。當電磁波通過總長度為L的CPC傳播時相移可表示為：

(1)

其中k媒質中的波數，k0真空中的波數，c真空中的光速，因此CPC的折射率為：

(2)

為了計算相移需要計算CPC的反射系數R和傳輸系數T。

圖1 同軸晶體結構。有兩種不同阻抗的同軸電纜交替周期相連組成，共2m+1節。其中一種同軸電纜的阻抗為z1，長度為l1，波數為k1；另一種同軸電纜的阻抗為z2，長度為l2，波數為k2。

首先計算第2n節處的反射和透射系數。根據同軸線的電壓電流公式以及連續的邊界條件，可得：

a2nexp(-jφ2n)+Rznaznexp(jφ2n)

=T2na2nexp(-jφ2n)+R2n+1T2na2nexp(jφ2n)

(3)

(4)

整理后得第2n節的反射系數R2n和透射系數T2n分別為：

(5)

(6)

同理，可得第2n-1節的反射系數和透射系數

(7)

(8)

其中φ2n-1=n(l1k1+l2k2)。

令初始相移為零，得第2n-1節相移φ2n-1和第2n節的相移φ2n分別為：

(9)

(10)

通過迭代法可得總相移，最后算出有效折射率nr，再利用色散媒質中的群速公式進行計算

(11)

從而根據公式可得負群速傳播的條件為：

(12)

3 CPC色散曲線和群速計算

首先，令圖1所示的同軸晶體系統中z1=50，z2=75，l1=5m，lz=5m，m=7，電磁波在兩種同軸線中的傳播速度為0.66c，利用上述方法進行電腦計算，得到該系統的色散曲線和群速，見圖2、圖3，

圖2 CPC的色散曲線

圖3 CPC中電磁波的群速曲線

可見該CPC在8～12MHz出現反常色散，而在9～11MHz群速大于光速，最大群速可達 1.66c。

圖4 不同參數對色散的影響

圖5 不同參數對群速的影響

此后，我們對不同參數的CPC進行電腦計算，得到圖4、圖5。可以看出，改變系統參數會對色散曲線和群速都產生影響，當阻抗z2增大時反常色散區變寬，程度加深，群速變大(vg=2.66c)。當節數m增加時反常色散區基本不變，程度加深，群速變大，此時獲得最大群速vg=3.25c，并且群速超光速的頻帶變寬。當增大每節長度l1或l2時，反常色散區改變，在兩個區發生反常色散，但程度基本不變，群速大小的與圖3相同(vg=1.66c)，只是出現兩個超光速頻帶。通過公式(11)我們可以看出減小有效折射率nr并且大大的加深反常色散的程度可能會使群速變為負。根據上述研究，我們增大阻抗z2增加節數m，來加深反常色散，并使電磁波在同軸線中傳播的速度達到c，即使如此群速也難變為負值。

4 獲得負群速的一種方法

根據公式(11)若有效折射率nr為負，那么在該系統中群速就更容易變為負。假設電磁波在兩種同軸線中的相速變為-0.66c，即折射率為nr=-1.5，其他參數不變，進行電腦計算結果見圖6、圖7，

圖6 折射率為-1.5的CPC的色散曲線

圖7 折射率為-1.5的CPC中電磁波的群速曲線

[1]Sommerfeld A，Brillouin L.Wave Propagation and Group Velocity[M].New York：Academic Press Inc，1960.

[2]Garrett C Q B，McCumber D E.Propagation of a Gaussian Light Pulse through an Anomalous Dispersion Medium[J].Phys Rev A，1970，1：305-313.

[3]Chu S，Wong S.Linear Pulse Propagation in an Absorbing Medium [J].Phy Rev Lett，1982，48：738-741.

[4]Wang L J，Kuzmich A，Dogarlu A.Gain-assisted Superluminal Light Propagation [J].Nature，2000，406：277-279.

[5]黃志洵.無源媒質中電磁波的異常傳播[J].中國傳媒大學學報(自然科學版)，2013，20：4-20。

[6]Huang Z X，Lu G Z，Guan J.Superluminal and Negative Group Velocity in the Electromagnetic Wave Propagation [J].Eng Sci，2003，1：35-39.

[7] Hache A，Poirier L.Long-range superluminal pulse propagation in a coaxial photonic crystal [J].Appl Phys Lett，2002，80：518-520.

[8]Lu G Z，Huang Z X，Guan J.Study on the Superluminal Group Velocity in a Coaxial Photonic Crystal [J].Eng Sci，2004，2：67-69.

[9]Munday J N，Robertson W M.Negative group velocity pulse tunneling through a coaxial photonic crystal [J].Appl Phys Lett，2002，81：2127-2129.

[10]周渭，李智奇.電領域群速超光速的特性實驗[J].北京石油化工學院學報，2009，3：48-53.

[11]Carot A，Nimtz G.Giant Negative Group Time Delay by Microwave Adaptors [J].Euro Phys Lett，2012，98：64002 1-4.