李蔭銘，孔 鵬，畢仁貴，何兆劍，鄧 科

(吉首大學(xué)物理與機(jī)電工程學(xué)院，湖南 吉首 416000)

近幾十年來(lái)，隨著拓?fù)浣^緣體[1-4]概念的提出，新奇的物理效應(yīng)使得其成為凝聚態(tài)物理備受關(guān)注的熱點(diǎn)課題.在電子系統(tǒng)中，人們已經(jīng)發(fā)現(xiàn)了描述波函數(shù)在相關(guān)色散帶上的量子化行為的拓?fù)涮匦訹5].受到電子系統(tǒng)中拓?fù)湎嚓P(guān)概念的啟發(fā)，研究人員通過(guò)類比經(jīng)典體系和量子體系，發(fā)現(xiàn)體態(tài)拓?fù)湎喾诸惡拖鄳?yīng)的拓?fù)浔砻鎽B(tài)(或邊界態(tài))也存在于光子晶體、聲子晶體等人工晶體中，其更加容易制備的優(yōu)勢(shì)使得在經(jīng)典波中實(shí)現(xiàn)拓?fù)湫?yīng)成為研究熱點(diǎn)[5-9].近年來(lái)，人們?cè)谔剿鞴庾泳w及聲子晶體的拓?fù)湫再|(zhì)中，出現(xiàn)了許多有趣的現(xiàn)象，如單向傳輸、定向無(wú)衍射傳輸?shù)萚10-13]，這為操縱聲波、電磁波尋找到更多新的自由度.

在聲學(xué)系統(tǒng)實(shí)現(xiàn)拓?fù)湫?yīng)的研究中，最早出現(xiàn)的是通過(guò)打破體系的時(shí)間反演對(duì)稱性，從而實(shí)現(xiàn)聲量子霍爾效應(yīng)，但此實(shí)現(xiàn)方式往往需要引入外部場(chǎng)，如引入循環(huán)流體等[2].隨后有學(xué)者提出了聲學(xué)量子自旋霍爾效應(yīng)，理論的提出已成功證明了在時(shí)間反演不變中系統(tǒng)也能實(shí)現(xiàn)聲學(xué)拓?fù)湫?yīng)[14]，相較于之前研究需要打破時(shí)間反演對(duì)稱性，這無(wú)疑是提供了一種更為簡(jiǎn)單的解決方案.2016年，Lu等[15]利用打破鏡面對(duì)稱性的方法實(shí)現(xiàn)了量子谷霍爾效應(yīng)，提出了聲學(xué)谷態(tài)的概念，這也是在聲學(xué)系統(tǒng)中實(shí)現(xiàn)拓?fù)鋺B(tài)最簡(jiǎn)單的一種辦法.目前，大部分研究較多關(guān)注在流體介質(zhì)中實(shí)現(xiàn)聲學(xué)拓?fù)鋺B(tài)，而固體材料中彈性波的拓?fù)湫再|(zhì)相關(guān)研究卻相對(duì)較少，這是由于聲波在固體材料中是以橫向極化和縱向極化兩種模式共同存在，在彈性波中尋找拓?fù)鋺B(tài)，用于實(shí)現(xiàn)類似于空氣聲中的邊界輸運(yùn)、選擇性激發(fā)等還存在一定困難，不過(guò)這也意味著此研究方向蘊(yùn)含更大的研究空間.

筆者基于彈性材料設(shè)計(jì)了二維蜂窩聲子晶體，在多頻段實(shí)現(xiàn)了SV模式波的多個(gè)狄拉克簡(jiǎn)并點(diǎn).當(dāng)相鄰兩個(gè)六邊形完全相同時(shí)，整個(gè)體系受時(shí)間反演和空間反演對(duì)稱性的保護(hù)，在布里淵區(qū)角點(diǎn)處出現(xiàn)3個(gè)線性簡(jiǎn)并的狄拉克點(diǎn)，分別對(duì)應(yīng)在不同的頻段中.通過(guò)調(diào)整蜂窩晶格的幾何結(jié)構(gòu)，破壞體系的空間對(duì)稱性，使其在K(K′)點(diǎn)處的多個(gè)狄拉克點(diǎn)同時(shí)打開(kāi)，在多個(gè)頻段形成谷態(tài).

1 二維彈性聲子晶體的設(shè)計(jì)

使用細(xì)長(zhǎng)管道將相鄰的六邊形A和B兩兩連接起來(lái)，組成規(guī)則的蜂窩六角晶格聲子晶體，見(jiàn)圖1.紅色實(shí)線包圍的是單胞結(jié)構(gòu)，該聲子晶體的晶格常數(shù)a=4 mm，其中六邊形A和B的晶格常數(shù)用a1、a2表示，均為0.55 mm，連接通道的寬度b=0.2 mm.使用不銹鋼作為整個(gè)模型的材料，參數(shù)如下：彈性模量E=210.6 GPa，剪切模量μ=81 GPa，質(zhì)量密度ρ=7 780 kg/m3.

圖1 聲子晶體的幾何結(jié)構(gòu)

圖2 聲子晶體的能帶結(jié)構(gòu)

圖3 不同的微擾參數(shù)Δr對(duì)應(yīng)的單胞結(jié)構(gòu)

2 體模態(tài)的渦旋場(chǎng)及拓?fù)湫再|(zhì)分析

圖4 體模態(tài)的渦旋場(chǎng)和拓?fù)湎嘧?/p>

圖5 能帶反轉(zhuǎn)時(shí)上、下能谷頻率隨Δr的變化

HK=vD(δkxσx+δkyσy)+Δgσz，

(1)

其中：vD是狄拉克點(diǎn)附近的有效群速度；δkx和δky則分別表示動(dòng)量空間中波矢在K點(diǎn)附近的偏移量；σi是泡利矩陣.當(dāng)Δr=0.06時(shí)，上、下能谷之間的帶隙頻率可以表示為Δg=f+-f-，對(duì)于較低一條帶，貝利曲率

(2)

3 總結(jié)

筆者基于彈性材料設(shè)計(jì)了二維蜂窩聲子晶體，在系統(tǒng)具有C3v對(duì)稱性時(shí)，在布里淵區(qū)角點(diǎn)處K(K′)點(diǎn)出現(xiàn)了3個(gè)線性簡(jiǎn)并的狄拉克點(diǎn).引入微擾參數(shù)打破體系的空間對(duì)稱性，使得原本出現(xiàn)的狄拉克點(diǎn)打開(kāi)，帶隙形成能谷.并通過(guò)對(duì)低頻段的上、下能谷的體模態(tài)渦旋場(chǎng)進(jìn)行數(shù)值模擬分析，發(fā)現(xiàn)該聲子晶體隨著結(jié)構(gòu)參數(shù)的改變發(fā)生了拓?fù)湎嘧儯瑫r(shí)揭示了谷態(tài)出現(xiàn)的原因.最后利用k·p微擾理論對(duì)其在兩種微擾參數(shù)下能谷的拓?fù)浜蛇M(jìn)行分析，發(fā)現(xiàn)其具有非零的谷陳數(shù).