楊蕭玥

（西南大學(xué) 物理科學(xué)與技術(shù)學(xué)院，重慶400715）

1 概述

近年來，由于能源的短缺及環(huán)境問題的加重，一種可以將熱能和電能相互轉(zhuǎn)換的功能材料——熱電材料引起了科學(xué)家們的廣泛關(guān)注[1-3]。熱電材料也叫溫差電材料，主要是利用固體內(nèi)部載流子和聲子的輸運(yùn)及其相互作用實(shí)現(xiàn)熱電轉(zhuǎn)換，具有眾多優(yōu)點(diǎn)，如：體積小、重量輕且運(yùn)行安靜，不會(huì)造成任何環(huán)境污染，節(jié)約能源，使用壽命長，易于控制等。由此，在溫差發(fā)電和制冷領(lǐng)域其都具有重要的應(yīng)用價(jià)值和廣泛的應(yīng)用前景。同時(shí)，利用熱電材料制備的微型元件如微型電源、微區(qū)冷卻、光通信激光二極管和紅外線傳感器的調(diào)溫系統(tǒng)也大大拓展了熱電材料的應(yīng)用領(lǐng)域[4]。熱電器件的工作效率取決于熱電材料的性能，材料的熱電性能可通過熱電優(yōu)值ZT 來評(píng)估，ZT=S2σT/（κe+κl），其中S 代表塞貝克系數(shù)，σ 表示電導(dǎo)率，T 表示絕對(duì)溫度，κe和κl分別是電子熱導(dǎo)率和晶格熱導(dǎo)率。好的熱電材料需要較高的ZT 值，因此需要較大的功率因數(shù)（S2σ） 和較低的熱導(dǎo)率（κe+κl）[5]。但是，這些輸運(yùn)參數(shù)（S，σ和κe）是相互依賴的[6]，因此我們需要尋找低的晶格熱導(dǎo)率κl。

由于量子限制效應(yīng)，二維半導(dǎo)體可以減少S 和σ 之間的耦合，從而提高熱電性能。隨著剝離和合成方法的發(fā)展，許多二維范德華材料已經(jīng)成功合成，包括石墨烯[7]，磷烯[8]和過渡金屬二鹵化物[9-11]。理論上已經(jīng)研究了許多具有高ZT 值的IV-VI 族化合物，例如二維SnSe 和GeSe，它們?cè)跍囟葹?00K 時(shí)的最大ZT 值分別為3.27[12]和1.99[13]。因此，二維鉛硫?qū)倩衔锟赡芫哂泻玫臒犭娦阅堋５侥壳盀橹?，?duì)一些已成功合成或預(yù)測(cè)的二維熱電材料的理論計(jì)算都是有限的，二維材料可能具有許多不同的結(jié)構(gòu)，因此，理論上可以研究所有可能存在的化合物形式，并探索其潛在的熱電性能，為實(shí)驗(yàn)上成功合成該材料奠定基礎(chǔ)。C. Kamal 等人系統(tǒng)研究了IV-VI 族單層（XY）材料（X=C，Si，Ge，Sn 和Y=O，S，Se，Te）的可能結(jié)構(gòu)，計(jì)算了這些材料的能帶并證明了其穩(wěn)定性[9]，我們選擇二維CSe，采用第一性原理計(jì)算和玻耳茲曼輸運(yùn)理論研究了其晶格熱導(dǎo)率。

2 計(jì)算方法

本文基于密度泛函理論（DFT），使用投影綴加平面波方法（PAW）與VASP 軟件包，來優(yōu)化CSe 的幾何結(jié)構(gòu)，其中，平面波截?cái)嗄茉O(shè)置為400 eV。我們使用20×20×1 的k 點(diǎn)網(wǎng)格，能量和力的收斂標(biāo)準(zhǔn)分別為10-6eV 和0.01 eV/魡，應(yīng)用了15魡的真空區(qū)域，以消除相鄰CSe 單層之間的相互作用。二階簡諧原子間力常數(shù)是由PHONOPY 軟件包使用3×3×1 的超胞和5×5×1 的k 點(diǎn)采樣獲得的。CSe 的三階非諧力常數(shù)通過使用3×3×1 超胞的有限差分方法獲得。以二階簡諧力常數(shù)和三階非諧力常數(shù)為輸入文件，可以通過求解聲子玻耳茲曼變換來獲得晶格熱導(dǎo)率κl。

3 結(jié)果和討論

優(yōu)化后的二維CSe 結(jié)構(gòu)的俯視圖和側(cè)視圖如圖1 所示，該結(jié)構(gòu)和藍(lán)磷的結(jié)構(gòu)類似。CSe 在平面上是各向同性的，屬于P3m1 的空間群，每個(gè)原胞中有兩個(gè)原子。優(yōu)化后CSe 單層的晶格常數(shù)為3.07魡，與別人的計(jì)算結(jié)果基本一致[14]。這表明我們?cè)O(shè)置的參數(shù)是合理的，得到的結(jié)果比較可靠，可進(jìn)行下一步的計(jì)算。

圖1 單層CSe 晶體結(jié)構(gòu)示意圖

接著，我們探索CSe 的聲子輸運(yùn)性質(zhì)，由此來獲得晶格熱導(dǎo)率。計(jì)算的二維CSe 聲子色散關(guān)系如圖2 所示。CSe 的原胞中有兩個(gè)原子，因此聲子譜具有三個(gè)聲學(xué)支（ZA，TA 和LA）和三個(gè)光學(xué)支，最下面的三條曲線是聲學(xué)支，包括兩個(gè)橫向振動(dòng)的分支（ZA，TA）和一個(gè)縱向振動(dòng)的分支（LA），其他的曲線是光學(xué)支，晶格的熱導(dǎo)率主要由聲學(xué)支貢獻(xiàn)，光學(xué)支對(duì)晶格熱傳導(dǎo)的影響較小。其中，ZA 分支在Г 點(diǎn)處有一個(gè)約為0.07 THz 的虛頻，這是因?yàn)橛?jì)算所使用的超胞大小有限，是可以忽略的系統(tǒng)誤差。從聲子譜的圖像可以看出，CSe 的結(jié)構(gòu)是穩(wěn)定的，可以用于進(jìn)行其他性質(zhì)的計(jì)算。

圖2 單層CSe 的聲子譜

根據(jù)聲子譜的數(shù)據(jù)，我們可以得到單層CSe 的晶格熱導(dǎo)率。通過分子動(dòng)力學(xué)的計(jì)算，CSe 單層的相變溫度約為750 K，因此我們選擇計(jì)算的溫度區(qū)間為300 K 到700 K，晶格熱導(dǎo)率κl與溫度的關(guān)系于圖3 所示。室溫下晶格熱導(dǎo)率κl為9.92 Wm-1K-1，它隨著溫度的升高而降低，這是因?yàn)楫?dāng)溫度升高時(shí)，晶格的振動(dòng)會(huì)更加強(qiáng)烈，增強(qiáng)了聲子的散射從而使晶格熱傳導(dǎo)降低，當(dāng)溫度為700K 時(shí)，κl的大小為4.21Wm-1K-1，低于部分磷烯和砷烯[15，16]。這表明CSe 單層可能具有良好的熱電性能。

圖3 晶格熱導(dǎo)率隨溫度的變化關(guān)系

4 結(jié)論

綜上所述，我們利用第一性原理計(jì)算和半經(jīng)典玻爾茲曼輸運(yùn)理論研究了二維CSe 單層的晶格結(jié)構(gòu)和晶格熱導(dǎo)率。首先，我們通過二維CSe 的幾何優(yōu)化確定了最穩(wěn)定的晶體結(jié)構(gòu)，其晶格常數(shù)為3.07魡。聲子譜證明了該結(jié)構(gòu)的動(dòng)態(tài)穩(wěn)定性。此外，通過求解聲子玻爾茲曼輸運(yùn)理論，得到二階簡諧和三階非諧原子間力常數(shù)，從而計(jì)算出了晶格熱導(dǎo)率，當(dāng)溫度范圍在300K 到700K 之間時(shí)，晶格熱導(dǎo)率κl的取值范圍為9.92Wm-1K-1到4.21Wm-1K-1，是一個(gè)較低的值，這表明，二維CSe 單層可能是潛在的熱電材料。