Anharmonic Phonon Renormalization and Defect Tolerance of the Thermoelectric Power Factor in Monolayer SnSe

技術要約：単層SnSeにおける非調和フォノン繰り込みと欠陥耐性に関する熱電出力因子

問題提起
単層のスズセレン化物（SnSe）は有望な熱電材料であるが、その性能は、相に依存する非調和格子力学および固有の点欠陥の影響に関する未解決の課題によって制約されている。バルクのSnSeは、強いフォノン非調和性に起因する高い熱電効率（$ZT$）を示すが、単層形態における熱電出力因子（PF）を制御する具体的なメカニズムは依然として不明である。主な課題は以下の通りである：

1. 非調和効果： 有限温度におけるフォノン繰り込みが、低温度の極性$\alpha$-SnSe（$Pnma$）および高温の非極性$\beta$-SnSe（$Cmcm$）の両相において、電子－フォノン（e-ph）散乱率および輸送係数をどの程度変化させるか。
2. 欠陥耐性： 固有の空孔（SnおよびSe）がどのように弾性散乱チャネルを導入してキャリア移動度とPFを低下させるのか、また、性能を維持するための欠陥濃度の定量的な限界が存在するかどうか。
3. 安定性： 調和近似は高温の$\beta$-SnSe相において動的不安定性（虚数周波数）を予測するため、動作温度（800–1000 K）での輸送計算のために格子を安定化させる手法が必要となる。

手法
著者らは、密度汎関数理論（DFT）、確率的自己一貫調和近似（SSCHA）、およびボルツマン輸送方程式（BTE）計算を組み合わせた包括的な第一原理フレームワークを用いている。

• 構造および電子構造モデリング： Quantum ESPRESSO（PBE汎関数）を用いたDFT計算を、単層の$\alpha$-および$\beta$-SnSeに対して行った。輸送解析のために、緻密なバンド構造を生成するためにWannier補間が使用された。
• フォノン繰り込み： SSCHAを用いて、有限温度における繰り込まれたフォノン周波数を計算した。$\alpha$-SnSeに対しては、これは動的に安定な調和スペクトルに対する補正として機能した。$\beta$-Seに対しては、虚数のソフトモードを除去し、800、900、1000 Kにおける$Cmcm$構造を安定化させるためにSSCHAが不可欠であった。
• 散乱計算：
  • 電子－フォノン（e-ph）： 散乱率は、2D極性補正を含むPerturboコード内のフェルミの黄金律を用いて計算された。散乱率はフォノンモード（LO, TO, ZO）ごとに分解された。
  • 電子－欠陥（e-d）： Sn空孔（$V_{Sn}$）およびSe空孔（$V_{Se}$）による散乱は、欠陥を希薄限界における弾性散乱体として扱うLuらによる定式化を用いてモデル化した。
• 輸送特性： 電気伝導度（$\sigma$）およびゼーベック係数（$S$）は、輸送分布関数から導出された。全伝導度は、部分抵抗率を加算するマティッセンの則によって得られ、ゼーベック係数はノルトハイム・ゴーターの則によって結合された。
• 欠陥耐性指標： 欠陥のないフォノン制限値に対してPFが15%低下する濃度として、運用上の臨界欠陥濃度（$C_{crit}$）を定義した。

主要な貢献および結果

1. $\alpha$-SnSeにおける非調和繰り込み：

   • $\alpha$-SnSeの調和スペクトルは動的に安定である。SSCHAによる繰り込みは、音響枝を低周波数へ、光学枝を高周波数へとシフトさせるが、n型輸送の全体像を質的に変えるものではない。
   • 散乱チャネル： 電子散乱は光学フォノンによって支配される。100 Kでは、横波光学モード（TO-2）が最大の寄与を与える。300 Kでは、縦波光学モード（LO-1）が支配的になる。
   • 出力因子： フォノン繰り込みは数値的なPF値を変化させるが（例：300 Kにおけるn型PFは8.96から10.01 $\mu$W/(K$^2$·cm)へ増加）、全体的な大きさや温度傾向は維持される。p型PFはこれらの繰り込みに対してより高い感受性を示す。

2. $\beta$-SnSeにおける安定化と輸送：

   • SSCHAは、調和ソフトモードの不安定性を除去することで高温の$Cmcm$相を正常に安定化させ、800–1000 Kでの輸送計算を可能にした。
   • 支配的な散乱： 調査された温度範囲全体を通じて、LO/TO-2モードが主要な電子散乱チャネルであり、その散乱率は1000 Kで著しく増加する。
   • 出力因子： 低密度窓（$\sim 10^{12}$ cm$^{-2}$）において、n型PFは800–900 Kで**15–19 $\mu$W/(K$^2$·cm)**に達し、p型値を大幅に上回る。この優位性は、ゼーベック係数が同等である一方で、n型ドーピングにおける電気伝導度が高いことに主に起因する。
   • 高密度限界： 全体的なPF最大値（~43–50 $\mu$W/(K$^2$·cm)）は、キャリア密度が$10^{14}$ cm$^{-2}$付近で発生する。著者らは、これらはフォノン制限による上限である可能性が高いと指摘している。なぜなら、このような密度を高温で維持することは実験的に困難であり、リジッドバンド近似に抵触する可能性があるためである。

3. 欠陥耐性と臨界濃度：

   • 空孔の影響： Se空孔（$V_{Se}$）は、Sn空孔（$V_{Sn}$）よりも弱い電子－欠陥散乱（より長い緩和時間）を誘起する。したがって、両方の相において、p型輸送はn型輸送よりも欠陥耐性が低い。
   • $C_{crit}$の値：
     • 最も低い耐性限界は、100 Kにおけるp型$\alpha$-SnSeと$V_{Sn}$で見られ、その値は$C_{crit} = 8.841 \times 10^{-5}$（約88 ppm）である。
     • 最も高い耐性は、n型$\beta$-SnSeと$V_{Se}$において観察され、15%の劣化閾値は調査された範囲の上限（$5 \times 10^{-3}$ または 5000 ppm）においても到達しない。
   • 温度依存性： 一般に、e-ph抵抗率の割合的な寄与が増大するため、温度とともに$C_{crit}$は増加する。これは、e-d散乱の相対的な影響を希釈するためである。

意義および主張
本論文は、2D SnSeにおける非調和格子力学と欠陥散乱の相互作用を理解するための定量的な枠組みを提供すると主張している。その主な貢献は以下の通りである：

• メカニズムの区別： 安定な$\alpha$-SnSeにおける有限温度のフォノン繰り込みと、$\beta$-SnSeにおける非調和安定化を区別し、後者が高温輸送を計算するための前提条件であることを示した。
• 欠陥限界： 出力を維持するための具体的な運用上の欠陥濃度限界（$C_{crit}$）を確立し、Sn空孔がSe空孔よりも厳格な制約を課すこと、特にp型輸送において顕著であることを特定した。
• 材料最適化： Se空孔は、電子散乱を弱める一方で格子熱伝導率を抑制する可能性がある（将来の検証のための仮説として言及）ため、熱電応用において有利なトレードオフを提供することを示唆している。
• 性能ベンチマーク： 単層SnSeの具体的なPFベンチマークを提供し、キャリア密度を制御できれば、高温における$\beta$相のn型ドーピングが高性能のポテンシャルを持つことを強調している。

著者らは、自身の計算がフォノン制限による上限と欠陥耐性の限界を提供するものである一方で、これらの特定の単層特性および空孔が単層形態の格子熱伝導率（$\kappa_L$）に与える影響についての直接的な実験的検証が依然として必要であると述べ、控えめな姿勢を維持している。