尹傳磊， 朱學(xué)峰， 楊非凡， 李 蓉， 周毅堅(jiān)， 彭文博， 趙 洋

(河南科技大學(xué) 物理工程學(xué)院 河南省光電儲(chǔ)能材料與應(yīng)用重點(diǎn)實(shí)驗(yàn)室，洛陽 471023)

1 引 言

Su-Schrieffer-Heeger (SSH)模型是一種用來描述一維聚乙炔的緊束縛模型，具有十分優(yōu)美的數(shù)學(xué)結(jié)構(gòu)[1]. 近年來關(guān)于拓?fù)浣^緣體[2,3]的研究，引起了人們對(duì) SSH模型的廣泛關(guān)注. 更值得關(guān)注的是，由于 SSH模型具有多種物理現(xiàn)象[4-8]，特別是其不平庸的拓?fù)湎辔籟9,10]，人們對(duì) SSH模型的研究日益深入. SSH模型可以作為一維聚乙炔[1]、石墨烯絲帶[11]、非對(duì)角雙色光晶格[12]和p軌道光學(xué)階梯系統(tǒng)[13]等多種系統(tǒng)的數(shù)學(xué)投影. 然而上述文章大多是關(guān)于 SSH模型所表現(xiàn)出的量子狀態(tài)，很少有關(guān)于 SSH模型熱力學(xué)性質(zhì)的文章[14,15]. 此外，文獻(xiàn)[14]中熱容對(duì)溫度的依賴性與文獻(xiàn)[15]中的不同. 更值得注意的是， Zoli的計(jì)算沒有介紹二聚化[14]. 在計(jì)算擴(kuò)展后的 SSH模型的熱力學(xué)特性時(shí)，由于不含二聚化的 SSH模型不能用來描述聚乙炔，因此在計(jì)算擴(kuò)展SSH模型的熱力學(xué)特性時(shí)，二聚化不應(yīng)該被忽略.

本文利用統(tǒng)計(jì)力學(xué)方法和費(fèi)曼路徑積分方法，研究了含三個(gè)過渡項(xiàng)的擴(kuò)展 SSH模型在周期性邊界條件下的熱力學(xué)性質(zhì). 在這項(xiàng)工作中不涉及定壓熱容，所以文中所用熱容均為指定熱容. 根據(jù)統(tǒng)計(jì)力學(xué)的一般表述，一個(gè)系統(tǒng)的熱力學(xué)性質(zhì)可以由它的配分函數(shù)得出[16]. 而該配分函數(shù)又可以通過費(fèi)曼路徑積分表述[17]. 此外，在zoli的研究中[14]，配分函數(shù)由費(fèi)曼路徑積分得到. 又因?yàn)樵撃Ｐ偷呐浞趾瘮?shù)可以用一個(gè)連乘來表示，所以其亥姆霍茲自由能、熵、內(nèi)能、熱容均可用求和的方式進(jìn)行表述. 文章用圖解法，解釋了模型ⅰ的熱容與溫度的比率(HCOTR)對(duì)溫度的依賴性，并解釋了擴(kuò)展的 SSH模型(模型ⅱ)只包含兩種躍遷項(xiàng)(NN, NNN躍遷項(xiàng))和只包含一種躍遷項(xiàng)(NN躍遷項(xiàng))的 SSH模型(模式識(shí)別)[1]對(duì)溫度的依賴性.

2 理論模型

一個(gè)原胞包含a和b兩個(gè)格點(diǎn)的一維晶格的哈密頓量是[1, 9, 10]

H=HSSH+HNNN+HNNNN

(1)

其中

(2)

(3)

(4)

這里t1=t+ε和t2=t-ε均是原胞中的躍遷幅度，ta是沿著a子格的躍遷幅度，tb是沿著b子格的躍遷幅度，ca,n(或者cb,n)是a(或者b)子格中的湮滅算符，t3和t4是NNNN躍遷幅度. 為了簡化模型，將t=1作為能量的單元[9, 10].

(5)

(6)

根據(jù)矩陣?yán)碚摰姆椒ǎ覀兊玫搅司仃噃(k)的特征值

(7)

其中

(8)

3 結(jié)果與討論

模型ⅰ的配分函數(shù)可以表示為費(fèi)曼路徑積分的形式

(9)

其中歐幾里德作用量

(10)

泛函積分的測度

(11)

且

(12)

(13)

且該場可以用傅里葉級(jí)數(shù)表示

其中松原頻率為

(14)

l屬于所有整數(shù)的集合. (9)式表示的配分函數(shù)可以寫為

(15)

其中測度

(16)

且作用量

(17)

通過計(jì)算路徑積分(15)，得到了模型ⅰ的配分函數(shù)

(18)

這里Z±(k)=1+e-βE±(k). 模型ⅰ配分函數(shù)的自然對(duì)數(shù)乘-T，得到它的亥姆霍茲自由能

(19)

利用熱力學(xué)關(guān)系，得到它的熵

(20)

和內(nèi)能

(21)

通過算其內(nèi)能對(duì)溫度的偏導(dǎo)數(shù)，我們得到它的熱容

(22)

圖1(a)和(b)顯示了模型ⅰ(由實(shí)線表示)、模型ⅱ(由點(diǎn)線表示)、模型ⅲ(由短劃線表示)的HCOTR的溫度依賴性.

為了獲得圖1,我們設(shè)置ε=0.065、ta=0.5、tb=0.6、t3=0.08、t4=0.08. 當(dāng)溫度由絕對(duì)零度逐漸升高時(shí)：模型ⅱ和模型ⅲ的HCOTR均呈現(xiàn)出先增大后減小再增大的趨勢，而模型ⅰ的HCOTR先減小，在達(dá)到最小值后增大. 在低溫(0-0.3)下每個(gè)模型都有一個(gè)HCOTR的極值，如圖1 (a)所示. 模型ⅱ的極值低于模型ⅲ的極值，但高于模型ⅰ的極值.

由于HCOTR在低溫下隨溫度增加而上升是玻璃態(tài)的特殊性質(zhì)，所以這三個(gè)極值點(diǎn)是三個(gè)相變點(diǎn)[14,18]. 如圖1(a)所示，模型ⅰ的相變溫度低于模型ⅱ的相變溫度，且模型ⅱ的相變溫度低于模型ⅲ的相變溫度. 這表明躍遷項(xiàng)可以降低相變溫度.

在極低溫度下，模型ⅱ和模型ⅲ的HCOTR都隨著溫度升高而增大，而模型ⅰ的HCOTR隨溫度升高而降低. 模型ⅰ和其他兩種模型HCOTR的單調(diào)性之間的這種差異，根本原因是NNNN躍遷項(xiàng).

圖1 模型ⅰ(由實(shí)線表示)，模型ⅱ(由點(diǎn)線表示)與模型ⅲ(由短劃線表示)的HCOTR隨溫度增加而變化的曲線. 每條曲線的ε值都設(shè)定為0.065，ta值都設(shè)定為0.5，tb值都設(shè)定為0.6，t3值都設(shè)定為0.08，t4值都設(shè)定為0.08.(a)和(b)除溫度域不同外，其它參數(shù)均相同.Fig. 1 HCOTR curves of model ⅰ (represented by the solid line), model ⅱ (represented by the dotted line) and model iii (represented by the dashed line) with increasing temperature. The value of each curve is set to 0.065, the value of ta is set to 0.5, the value of tb is set to 0.6, the value of t3 is set to 0.08, and the value of t4 is set to 0.08. (a) and (b) except for the difference in temperature range, other parameters are the same ones.

如圖1(b)所示，在高溫(0.3～1)下這三種模型的HCOTR的溫度依賴性相似，當(dāng)溫度由0.3逐漸升高時(shí)：這三種模型的HCOTR均先增大，在達(dá)到極大值后減小. 在高溫(0.3～1)下模型ⅱ的最大HCOTR高于模型ⅲ的HCOTR，而低于模型ⅰ的HCOTR.

在圖2(a)和(b)描繪的是在ε取四個(gè)不同值時(shí)模型ⅰ的HCOTR的溫度依賴性，其中實(shí)線對(duì)應(yīng)ε=0、點(diǎn)線對(duì)應(yīng)ε=0.15、點(diǎn)劃線對(duì)應(yīng)ε=0.065、短劃線對(duì)應(yīng)ε=0.3 . 為了獲得圖2，我們設(shè)置ta=0.5、tb=0.6、t3=0.08、t4=0.08 . 如圖2(b)所示，在高溫(0.3～1)下模型ⅰ中的ε取四個(gè)不同值時(shí)HCOTR有類似的溫度依賴性且HCOTR的最大值隨ε減少而增加；當(dāng)溫度由0.3逐漸升高到1時(shí)，這四個(gè)ε值對(duì)應(yīng)的HCOTR都先增大,到最大值后減小.

然而如圖2(a)所示，在極低溫(0-0.05)下不同ε值對(duì)應(yīng)的HCOTR的單調(diào)性不同，這表明NN躍遷項(xiàng)會(huì)影響模型ⅰ的HCOTR的單調(diào)性.ε=0.3時(shí)模型ⅰ有兩個(gè)相變點(diǎn)，這表明NN躍遷項(xiàng)會(huì)影響模型ⅰ的相變.

圖2 在ε取四個(gè)不同值時(shí)模型ⅰ的HCOTR的溫度依賴性. 實(shí)線對(duì)應(yīng)ε=0；點(diǎn)線對(duì)應(yīng)ε=0.15；點(diǎn)劃線對(duì)應(yīng)ε=0.065；短劃線對(duì)應(yīng)ε=0.3. 對(duì)于每條曲線，ta值都設(shè)定為0.5，tb值都設(shè)定為0.6，t3值都設(shè)定為0.08，t4值都設(shè)定為0.08. (a)和(b)除溫度域不同外，其它參數(shù)均相同.Fig. 2 The temperature dependence of the HCOTR of model i when four different values ε are used. The solid line corresponds to ε=0; the dotted line corresponds to ε=0.15; the dash dot line corresponds to ε=0.065; the dashed line corresponds to ε=0.3. For each curve, the ta value is set to 0.5, the tb value is set to 0.6, the t3 value is set to 0.08, and the t4 value is set to 0.08. (a) and (b) except for the difference in temperature range, other parameters are the same ones.

圖3(a)和(b)顯示了在ta、tb取四組不同值時(shí)模型ⅰ的HCOTR的溫度依賴性. 實(shí)線對(duì)應(yīng)ta=0.3、tb=0.5，點(diǎn)線對(duì)應(yīng)ta=0.2、tb=0.7，點(diǎn)劃線對(duì)應(yīng)ta=0.6、tb=0.3，短劃線對(duì)應(yīng)ta=0.8、tb=0.3. 為了獲得圖3，我們設(shè)置ε=0.065、t3=0.08、t4=0.08. 如圖3(b)所示，在高溫(0.3～1)下模型ⅰ中的ta、tb取四組不同值時(shí)HCOTR有類似的溫度依賴性. 當(dāng)溫度由0.3逐漸升高到1時(shí)，這四組ta、tb值對(duì)應(yīng)的HCOTR都先增大，到最大值后減小.

然而如圖3(a)所示，在低溫(0.3-1)下不同ta、tb值對(duì)應(yīng)的HCOTR的單調(diào)性不同，這表明NNN躍遷項(xiàng)會(huì)影響模型ⅰ的HCOTR的單調(diào)性.ta=0.2、tb=0.7時(shí)模型ⅰ有兩個(gè)相變點(diǎn)，這表明NNN躍遷項(xiàng)可以影響模型ⅰ的相變.

圖4(a)和(b)描繪了在t3、t4取四組不同值時(shí)模型ⅰ的HCOTR的溫度依賴性. 實(shí)線對(duì)應(yīng)t3=0.08、t4=0.08，點(diǎn)線對(duì)應(yīng)t3=0.01、t4=0.07，點(diǎn)劃線對(duì)應(yīng)t3=0.06、t4=0.01，短劃線對(duì)應(yīng)t3=0.08、t4=0.03 . 為了獲得圖4，我們設(shè)置ε=0.065、ta=0.5、tb=0.6 . 在高溫(0.1～1)下模型ⅰ中的t3、t4取四組不同值時(shí)HCOTR有類似的溫度依賴性. 當(dāng)溫度由0.1逐漸升高到1時(shí)，這四組t3、t4值對(duì)應(yīng)的HCOTR均先增大到最大值后減小.

然而如圖4(a)所示，在低溫(0-0.1)下不同t3、t4值對(duì)應(yīng)的HCOTR的單調(diào)性不同，這表明NNNN躍遷項(xiàng)會(huì)影響模型ⅰ的HCOTR的單調(diào)性.

4 總 結(jié)

利用統(tǒng)計(jì)力學(xué)的一般表述和費(fèi)曼路徑積分，計(jì)算了模型ⅰ在周期性邊界條件下的熱力學(xué)性質(zhì). 對(duì)模型ⅱ和模型ⅲ也進(jìn)行了類似的計(jì)算并與模型ⅰ作對(duì)比. HCOTR對(duì)溫度的依賴曲線有最小值點(diǎn)表明這三種模型均存在相變.在低溫下，每個(gè)模型都有相變且躍遷項(xiàng)可以降低相變溫度. 在極低的溫度下，NNNN躍遷項(xiàng)可以引起模型ⅰ與其他兩種模型在HCOTR的單調(diào)性方面的差異. 在低溫下，NN、NNN、NNNN三種躍遷項(xiàng)都可能影響模型ⅰ的HCOTR的單調(diào)性. NN或NNN躍遷項(xiàng)可以影響模型ⅰ的相變.

圖3 在ta、tb取四組不同值時(shí)模型ⅰ中的HCOTR的溫度依賴性. 實(shí)線對(duì)應(yīng)ta=0.3、tb=0.5；點(diǎn)線對(duì)應(yīng)ta=0.2、tb=0.7；點(diǎn)劃線對(duì)應(yīng)ta=0.6、tb=0.3；短劃線對(duì)應(yīng)ta=0.8、tb=0.3. 對(duì)于每條曲線，ε=0.065、t3=0.08、t4=0.08.(a)和(b)除溫度域不同外，其它參數(shù)均相同.Fig. 3 The temperature dependence of HCOTR in model i when ta and tb take four different values. The solid line corresponds to ta=0.3, tb =0.5; the dotted line corresponds to ta=0.2, tb=0.7; dash dot line corresponds to ta=0.6, tb=0.3; the dashed line corresponds to ta=0.8, tb=0.3. For each curve, ε=0.065, t3=0.08, t4=0.08. (a) and (b) except for the difference in temperature range, other parameters are the same ones.