(57)【要約】本開示による熱電変換材料は、化学式Ｌｉ3-aＢｉ1-bＧｅbにより表される組成を有する。ここで、ａおよびｂの値の範囲は：０≦ａ≦０．０００３、かつ−ａ＋０．０００３≦ｂ≦０．１０８；０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８；または０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［−０．１５７×（ｌｎ（ａ））2−２．２２×ｌｎ（ａ）−９．８１］であり、前記熱電変換材料はＢｉＦ3型の結晶構造を有し、かつ前記熱電変換材料はｐ型の極性を有する。
2020021783.tif 000002
【特許請求の範囲】
【請求項１】
化学式Ｌｉ3-aＢｉ1-bＧｅbにより表される組成を有する熱電変換材料。
ここで、
前記熱電変換材料はＢｉＦ3型の結晶構造を有し、
前記熱電変換材料はｐ型の極性を有し、かつ
以下の数式（Ｉ）〜（ＩＩＩ）のいずれか１つが充足される。
０≦ａ≦０．０００３、かつ−ａ＋０．０００３≦ｂ≦０．１０８  （Ｉ）；
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８  （ＩＩ）；または
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［−０．１５７×（ｌｎ（ａ））2−２．２２×ｌｎ（ａ）−９．８１］  （ＩＩＩ）
【請求項２】
請求項１に記載の熱電変換材料であって、
０．０００３≦ａ≦０．０８５、かつｂ＝０である。
【請求項３】
請求項１に記載の熱電変換材料であって、
ａ＝０、かつ０．０００３≦ｂ≦０．１０８である。
【請求項４】
熱電変換材料を用いて電力を得る方法であって、
前記方法は、前記熱電変換材料に温度差を印加して電力を生じる工程を具備し、
前記熱電変換材料は、    化学式Ｌｉ3-aＢｉ1-bＧｅbにより表される組成を有し、
ＢｉＦ3型の結晶構造を有し、かつ
ｐ型の極性を有し、かつ
以下の数式（Ｉ）〜（ＩＩＩ）のいずれか１つが充足される。    ０≦ａ≦０．０００３、かつ−ａ＋０．０００３≦ｂ≦０．１０８  （Ｉ）；
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８  （ＩＩ）；または
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［−０．１５７×（ｌｎ（ａ））2−２．２２×ｌｎ（ａ）−９．８１］  （ＩＩＩ）
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、熱電変換材料および熱電変換材料を用いて電力を得る方法に関する。
【背景技術】
【０００２】
熱電変換材料の両端に温度差が与えられると、温度差に比例した起電力が発生する。熱エネルギーが電気エネルギーに変換されるこの現象は、ゼーベック効果として知られている。熱電発電技術では、ゼーベック効果を利用し、熱エネルギーを直接電気エネルギーに変換する。【０００３】
熱電変換材料の技術分野においてよく知られているように、熱電変換デバイスに用いられる熱電変換材料の性能は、性能指数ＺＴにより評価される。ＺＴは、以下の数式（Ｉ）により表される。
ＺＴ＝Ｓ2σＴ／κ  （Ｉ）
ここで、
Ｓは、物質のゼーベック係数を表し、
σは、物質の電気伝導率、および
κは、熱伝導率κを表す。
ＺＴの値が高いほど、熱電変換効率は高い。【０００４】
非特許文献１は、Ｌｉ3Ｂｉ結晶物質の製造方法を開示している。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Gen-Tao Zhou et al., “Microwave-assisted solid-state synthesis and characterization of intermetallic compounds of Li3Bi and Li3Sb”, Journal of Materials Chemistry, 13, p.2607-2611, (2003)
【非特許文献２】C. G. Van de Walle et al., “First-principles calculations for defects and impurities: Applications to III-nitrides”, Journal of Applied Physics, 95, p.3851-3879, (2004).
【非特許文献３】Y. Koyama et al., “First principles study of dopant solubility and defect chemistry in LiCoO2”, Journal of Materials Chemistry A, 2, p.11235-11245, (2014).
【非特許文献４】G. K. H. Madsen et al., “BoltzTraP. A code for calculating band-structure dependent quantities”, Computer Physics Communications, Volume 175, p.67-71, (2006).
【非特許文献５】H. Wang et al., in Thermoelectric Nanomaterials, ed. K. Koumoto and T. Mori, Springer, Berlin Heidelberg, vol.182, ch. 1, p.3-32, (2013).
【非特許文献６】J. Chen et al., “First-Principles Predictions of Thermoelectric Figure of Merit for Organic Materials: Deformation Potential Approximation”, Journal of Chemical Theory and Computation, 8, p.3338-3347, (2012).
【非特許文献７】J. Yang et al., “Material descriptors for predicting thermoelectric performance”, Energy & Environmental Science, 8, p.983-994, (2015).
【非特許文献８】Sean M. McDeavitt, “Synthesis and casting of a lithium-bismuth compound for an ion-replacement electrorefiner”, Light Metals (Warrendale, Pennsylvania), pp. 1139-1142, (1999)
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本開示の目的は、新規な熱電変換材料の提供である。
【課題を解決するための手段】
【０００７】
本開示は、化学式Ｌｉ3-aＢｉ1-bＧｅbにより表される組成を有する熱電変換材料を提供する。
ここで、
前記熱電変換材料はＢｉＦ3型の結晶構造を有し、かつ
前記熱電変換材料はｐ型の極性を有し、かつ
以下の数式（Ｉ）〜（ＩＩＩ）のいずれか１つが充足される。
０≦ａ≦０．０００３、かつ−ａ＋０．０００３≦ｂ≦０．１０８  （Ｉ）；
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８  （ＩＩ）；または
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［−０．１５７×（ｌｎ（ａ））2−２．２２×ｌｎ（ａ）−９．８１］  （ＩＩＩ）
【発明の効果】
【０００８】
本開示は、新規な熱電変換材料を提供する。
【図面の簡単な説明】
【０００９】
【図１】図１は、Ｌｉ3Ｂｉの結晶構造の模式図を示す。
【図２】図２は、Ｌｉ3Ｂｉ結晶構造の回折Ｘ線強度分布を示すグラフである。
【図３】図３は、ａ−ｂ平面上に実施例１〜１６および比較例１〜４の点をプロットしたグラフである。
【発明を実施するための形態】
【００１０】
以下、本開示の実施形態について、図面を参照しながら説明する。
【００１１】
本開示の熱電変換材料は、以下の化学式（Ｉ）により表される組成を有する。
Ｌｉ3-aＢｉ1-bＧｅb    （Ｉ）
【００１２】
図１は、Ｌｉ3Ｂｉ結晶構造の模式図を示す。非特許文献１にも開示されているように、Ｌｉ3Ｂｉ結晶構造は、に、ＢｉＦ3型結晶構造またはＡｌＣｕ2Ｍｎ型結晶構造を有する。ＢｉＦ3型結晶構造およびＡｌＣｕ2Ｍｎ型結晶構造の両者は、空間群Ｆｍ−３ｍに属する。非特許文献１においては、Ｌｉ3Ｂｉ結晶物質が熱電変換材料として扱われていない。そのため、非特許文献１には、性能指数ＺＴは開示されていない。
【００１３】
本発明者らは、材料インフォマティクスと呼ばれる、データ科学に基づいた材料探索手法を利用することで、無機結晶構造データベース中の数万件の化合物について性能指数ＺＴの予測値を算出した。算出には、本発明者らが独自に確立した性能指数ＺＴの予測モデルが用いられた。この予測モデルは、従来の手法よりも精度が高い。このため、この予測モデルを用いることで、従来よりも信頼性の高い予測結果を得ることができる。そこで、本発明者らは、これまで熱電変換材料として扱われてこなかったＬｉ3Ｂｉ結晶物質が熱電変換材料として有望な材料であるかどうかを検討した。
【００１４】
Ｌｉ3Ｂｉ結晶物質が欠陥を含まない状態において、Ｌｉ3Ｂｉ結晶物質は、キャリアに乏しい。このため、欠陥のないＬｉ3Ｂｉ結晶物質が高い性能指数ＺＴを有することは望めない。そこで、本発明者らは、性能指数ＺＴを向上させるために、Ｌｉ3Ｂｉ結晶物質に欠陥を導入し、ｐ型のキャリアを生じさせることを検討した。その結果、本発明者らは、Ｌｉサイトを空孔で置換することによって生じる欠陥を有する物質、およびＢｉサイトを元素Ｇｅで置換することによって生じる欠陥を有する物質の２種類の物質を想到するに至った。
【００１５】
本発明者らは、計算により、Ｌｉ3-aＢｉ1-bＧｅb結晶物質のみが安定的に得られるａおよびｂの値の範囲を導き出した。さらに、その範囲において性能指数ＺＴを算出することで、本発明者らは、０．４以上の高い性能指数ＺＴが得られるａおよびｂの値の範囲を見出した。具体的には、後述の実施例１〜１６および比較例１〜４において実証されているように、以下の数式（Ｉ）、（ＩＩ）、または（ＩＩＩ）のうちいずれか１つを充足する場合に、熱電変換材料は０．４以上の高い性能指数ＺＴを有する。
０≦ａ≦０．０００３、かつ−ａ＋０．０００３≦ｂ≦０．１０８  （Ｉ）
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８  （ＩＩ）、または
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［−０．１５７×（ｌｎ（ａ））2−２．２２×ｌｎ（ａ）−９．８１］  （ＩＩＩ）
具体的には、ａおよびｂの値の上記の範囲において、性能指数ＺＴの値は０．４２以上である。表１を参照せよ。
【００１６】
（製造方法）
本開示の熱電変換材料の製造方法の一例が、非特許文献１および非特許文献８の開示内容に基づいて以下に説明される。
【００１７】
非特許文献１は、Ｌｉ３Ｂｉ結晶物質の製造方法を開示している。。非特許文献１に開示された製造方法においては、まず、アルゴンガスで満たされたグローブボックスの中で、Ｌｉ箔および粒状のＢｉがモル比４３：１２でアルミナ製の坩堝に入れられる。次に、坩堝の内部で、これらのＬｉ箔および粒状のＢｉが共に圧縮される。カーボンでコーティングされた石英管の中に、坩堝は入れられる。石英管の内部は、その真空度が６．６７×１０-3Ｐａに達するまで真空引きされる。石英管は、酸素ガスの炎を用いて密封される。密封された石英管はマイクロ波キャビティの中に入れられ、マイクロ波が石英管に対して照射される。このようにして、Ｌｉ3Ｂｉ結晶物質が得られる。グローブボックスの中で石英管を開封することにより、石英管から得られたＬｉ3Ｂｉ結晶物質が取り出される。
【００１８】
非特許文献８もまた、、Ｌｉ３Ｂｉ結晶物質の製造方法を開示している。非特許文献８に開示された製造方法においては、まず、不活性ガスで満たされたグローブボックスに接続された誘導鋳造炉に、タンタル製の坩堝が収容される。ＬｉおよびＢｉが、Ｌｉ：Ｂｉのモル比３：１で坩堝の中に導入される。炉の内部がゆっくりと加熱される。炉の内部の加熱にしたがってまず、Ｌｉが融解する。Ｌｉの融解の後にＢｉが融解する。このようにして、ＬｉのＢｉとの反応が進行し、反応物を得る。得られた反応物が、１２００℃まで加熱されて均質化される。次いで反応物は冷却されて、Ｌｉ3Ｂｉ結晶物質が得られる。
【００１９】
非特許文献１および非特許文献８に開示された両方の製法において、Ｌｉの欠損量は、出発物質のＬｉの量を、Ｂｉの量に対して変化させることで制御できる。
【００２０】
また、元素Ｂｉの原子半径（０．１５６ｎｍ）は、元素Ｇｅの原子半径（０．１２２ｎｍ）と大差がない。このため、出発物質としてＧｅを加えることで、一部の元素Ｂｉは元素Ｇｅによって置換されうると考えられる。Ｇｅによる置換量は、出発物質のＧｅの量を、Ｂｉの量に対して変化させることで制御できる。
【００２１】
以上から、本開示によるＬｉ3-aＢｉ1-bＧｅb結晶物質は、非特許文献１および非特許文献８に開示されている製造方法を参照することで製造できると考えられる。
【００２２】
（熱電変換材料を用いて電力を得る方法）
本実施形態において、本開示のＬｉ3-aＢｉ1-bＧｅb結晶物質に温度差を印加することで電力が得られる。
【００２３】
（実施例）
以下の実施例を参照しながら、本開示の熱電変換材料がさらに詳細に説明される。
【００２４】
（結晶構造の解析）
非特許文献１は、Ｌｉ3Ｂｉ結晶物質の結晶構造がＢｉＦ3型構造に属することを開示している。Ｘ線結晶回折法による測定に基づき、ＢｉＦ3型構造のＸ線回折ピークを確認する事ができる。図２は、ソフトウェア（ＲＩＥＴＡＮ、入手先ＵＲＬ：http://fujioizumi.verse.jp/download/download.html）を用いてＬｉ3Ｂｉの結晶構造因子Ｆおよび積分回折強度Ｉを計算することにより得られたＬｉ3Ｂｉ結晶構造の回折Ｘ線強度分布を示すグラフである。
【００２５】
結晶構造因子Ｆは、以下の関係式（１）によって求められた。
Ｆ＝Σｆiｅｘｐ（２πｉｒiΔｋ）・・・（１）
ここで、
ｒiは結晶中の原子の位置ベクトルを表し
ｆiはｒiの位置にある原子の原子散乱因子を表し、かつ
Δｋは散乱前後におけるＸ線の波数ベクトルの差を表す。
【００２６】
積分回折強度Ｉは、以下の関係式（２）によって求められた。
Ｉ＝ＩeＬ｜Ｆ｜2Ｎ2・・・（２）
ここで、
Ｉeは１個の電子の散乱強度を表し、
Ｎは結晶中の単位格子の数を表し、
Ｌは吸引因子を含み、かつ実験条件に依存する係数を表す。
【００２７】
本開示による熱電変換材料は、Ｌｉ3Ｂｉ結晶において、元素Ｌｉサイトの一部が空孔により置換されることにより生じた欠陥、または元素Ｂｉの一部が元素Ｇｅによって置換された欠陥を有している。これらの欠陥により、結晶構造に格子変形が生じうる。このため、本開示による熱電変換材料の回折Ｘ線強度分布のピークは、図２に示され、かつ非特許文献１に開示されたＬｉ3Ｂｉ結晶物質のピークに対して若干シフトしていると予測される。
【００２８】
（安定組成範囲の計算方法）
本発明者らは、密度汎関数法（以下、「ＤＦＴ」という）に基づいた計算により、Ｌｉ3Ｂｉ結晶物質のバンド構造を評価した。その結果、本発明者らは、Ｌｉ3Ｂｉ結晶物質は半導体であるという結果を得た。ＤＦＴによる計算には、本発明者らは、第一原理電子状態計算プログラム（Vienna Ab initio Simulation Package、入手先:https://www.vasp.at/）を使用した。このプログラムは、以下において、「ＶＡＳＰ」と称される。
【００２９】
性能指数ＺＴの最大化を図るためには、半導体に欠陥を導入し、ｐ型またはｎ型のキャリアを生じさせる必要がある。キャリアを生じさせることにより、熱電変換材料の両端に温度差を与えたときに大きな起電力が得られることが期待される。なお、どのような種類の欠陥を導入すれば高いキャリア濃度が得られるかは自明ではない。性能指数ＺＴを高めるためには、高いキャリア濃度が得られる欠陥の種類を見つけ出す必要がある。
【００３０】
非特許文献１に開示されているＬｉ3Ｂｉ結晶物質は欠陥を含まないので、Ｌｉ3Ｂｉ結晶物質のキャリアは乏しい。そのため、Ｌｉ3Ｂｉ結晶物質は高い性能指数ＺＴを有することは予想されない。性能指数ＺＴを高めるために、本発明者らは、例えば、Ｌｉ3Ｂｉの結晶構造において空孔を生じさせる欠陥、または、他の元素による置換を検討した。
【００３１】
本発明者らは、ｐ型のキャリアを生じさせるための欠陥として、（Ｉ）Ｌｉ3Ｂｉの結晶構造において空孔を生じさせる欠陥、および（ＩＩ）元素Ｂｉを他の元素で置換することによって生じる欠陥、の２種類の欠陥を検討した。この検討に当たって、本発明者らは、Ｌｉ3Ｂｉの結晶構造にこれらの欠陥を導入しても、ＢｉＦ3型構造が維持されることを前提とした。言い換えると、以下の計算では各結晶構造はＢｉＦ3型構造を有すると本発明者らは仮定した。
【００３２】
Ｌｉ3Ｂｉ結晶構造において空孔を生じさせる欠陥に関しては、Ｌｉサイトに空孔を生じさせる場合の欠陥形成エネルギーＥform、およびＢｉサイトに空孔を生じさせる場合の欠陥形成エネルギーＥformが算出された。
【００３３】
その結果、Ｌｉサイトに空孔を生じさせる場合はｐ型のキャリア濃度が高い一方、Ｂｉサイトに空孔を生じさせる場合はｐ型のキャリア濃度が低いことを本発明者らは見出した。欠陥形成エネルギーＥformの詳細は後述される。
【００３４】
元素Ｂｉを他の元素で置換することによって生じる欠陥に関しては、他の元素の候補として元素Ｇｅ、Ｇａ、およびＰｂを採用することを本発明者らは試した。本発明者らは、元素Ｇｅ、ＧａおよびＰｂを、元素Ｂｉと比較して価数が少なく、かつ元素Ｂｉと比較してイオン半径が大差ないという点から選択した。Ｂｉの一部をＧｅ、Ｇａ、またはＰｂのいずれか１つにより置換したＬｉ3Ｂｉ結晶構造の欠陥形成エネルギーＥformを本発明者らは算出した。
【００３５】
その結果、元素Ｇｅで置換する場合はｐ型のキャリア濃度が高い一方、元素Ｇａまたは元素Ｐｂで置換する場合はｐ型のキャリア濃度が低いことを本発明者らは見出した。
【００３６】
以上の考察により、本発明者らは、ｐ型キャリアを生じさせる欠陥として、Ｌｉサイトに空孔を生じさせる欠陥、および元素Ｂｉを元素Ｇｅで置換することで生じる欠陥の２種類を想到するに至った。
【００３７】
本発明者らは、本開示にかかるＬｉ3-aＢｉ1-bＧｅbの安定組成範囲を計算した。安定組成範囲とは、複数の結晶相が形成されていない単一の結晶相が得られる組成範囲を意味する。以下、「単一の結晶相」を「単相」という。、Ｌｉ3-aＢｉ1-bＧｅbの安定組成範囲とは、Ｌｉ3-aＢｉ1-bＧｅb結晶物質のみが得られる組成範囲を意味する。この範囲外では、Ｌｉ3-aＢｉ1-bＧｅb結晶物質以外にも、他の組成の結晶相が析出し、単相のＬｉ3-aＢｉ1-bＧｅb結晶物質を得ることができない。
【００３８】
安定組成範囲の計算は、非特許文献２および３に開示された内容に基づいて行われた。
【００３９】
具体的には、まず、Ｌｉ3-aＢｉ1-bＧｅbにおいて、ＢｉＦ3型の結晶構造が安定化するａおよびｂの値の範囲が、非特許文献２に開示されている半導体の欠陥形成理論に基づいた方法により計算された。
【００４０】
Ｌｉ3Ｂｉにおける欠陥形成エネルギーＥformが、以下の関係式（３）に基づいて評価された。
Ｅform（μi，ｑ，ＥF）＝Ｅdefect−Ｅpure−Σｎiμi＋ｑ（ＥVBM＋ＥF）・・・（３）
ここで、
Ｅdefectは欠陥が存在する場合の全エネルギーを表し、
Ｅpureは欠陥が存在しない場合の完全結晶の全エネルギーを表し、
ｎiは欠陥によるｉ番目の構成元素の増減量を表し、
μiはｉ番目の元素の化学ポテンシャルを表し、
ｑは欠陥が有する電荷量を表し、
ＥVBMは半導体であるＬｉ3Ｂｉの価電子帯上端の１電子エネルギーを表し、かつ
ＥFは電子のフェルミエネルギーを表す。
【００４１】
これらのエネルギー値は、一般化勾配近似の範囲内でＤＴＦ計算を適用することにより評価された。ただし、Ｌｉサイトに空孔が生じるパターン、およびＢｉがＧｅによって置換されるパターンを考慮して、Ｅformは計算された。
【００４２】
関係式（３）によって得られた欠陥形成エネルギーＥformを用いて、各欠陥の体積密度ＮDがボルツマン分布による以下の関係式（４）に基づいて評価された。
ＮD（μi，ｑ，ＥF）＝Ｎsite×ｅｘｐ［−Ｅform（μi，ｑ，ＥF）／ｋbＴ］・・・（４）
ここで、
Ｎsiteは考慮している欠陥が生じうるサイトの体積密度を表し、
ｋbはボルツマン定数を表し、かつ
Ｔは絶対温度を表す。
【００４３】
各欠陥の電荷ｑ×ＮDと、半導体にドープされたキャリアの電荷Ｑeとの総量が０になるという電荷中性条件は常に満たされている。この電荷中性条件から、フェルミエネルギーおよびキャリア濃度が各組成において決定された。このときの価電子帯のキャリア濃度ｐと、伝導体のキャリア濃度ｎとは、以下の関係式（５）〜（７）により求められた。
ｐ＝１−∫ＤVB（Ｅ）［１−ｆ（Ｅ；ＥF）］ｄＥ・・・（５）
ｎ＝∫ＤCB（Ｅ）ｆ（Ｅ；ＥF）ｄＥ・・・（６）
ｆ（Ｅ；ＥF）＝１／［ｅｘｐ（（Ｅ−ＥF）／ｋBＴ）＋１］・・・（７）
ここで、
ＤVB（Ｅ）およびＤCB（Ｅ）は、それぞれ、ＤＦＴ計算によって得られる価電子帯および伝導帯の電子状態密度を表し、
ｆ（Ｅ；ＥF）はフェルミ分布関数を表し、かつ
Ｔは熱電変換特性の評価を実施する絶対温度を表す。
全キャリアの電荷密度Ｑeは、数式Ｑe＝ｅ×（ｎ−ｐ）に従って算出された。ｅは、１電子の持つ電荷である。
【００４４】
非特許文献３と同等の手法を用いて、本発明者らは、各元素の化学ポテンシャルμiのとり得る許容範囲を計算した。さらに関係式（４）により、欠陥密度の範囲が計算された。金属Ｌｉ、金属Ｂｉ、またはＬｉＢｉが生じる範囲は結晶構造が安定化しない範囲であるとの仮定の下、単相のＬｉ3-aＢｉ1-bＧｅb結晶物質が得られる化学ポテンシャルの範囲は、上述の結晶構造が安定化しない範囲を除いた範囲として決定される。この取りうる化学ポテンシャルの範囲内で、本発明者らは、欠陥密度およびそれに関わる化学式組成を評価した。組成は、各元素の化学ポテンシャルの関数として決定された。例えば、ある化学ポテンシャルの条件において、Ｌｉ3Ｂｉ結晶内の全Ｌｉサイトのうち５％が欠損して空孔化した場合、組成は、Ｌｉ2.85Ｂｉである。
【００４５】
前記の手法に基づいて、Ｌｉ3-aＢｉ1-bＧｅbのＬｉ原子、Ｂｉ原子、およびＧｅ原子に対する化学ポテンシャルの許容範囲を本発明者らは評価した。その化学ポテンシャルの範囲でＬｉサイトに空孔を生じさせる欠陥の密度およびＢｉサイトを元素Ｇｅで置換することによって生じる欠陥の密度を本発明者ら評価し。
【００４６】
その結果、Ｌｉ3-aＢｉ結晶物質に関しては、０≦ａ≦０．６５４の範囲において、単相のＢｉＦ3型結晶が得られるという計算結果が得られた。
【００４７】
Ｌｉ3Ｂｉ1-bＧｅb結晶物質に関しては、０≦ｂ≦０．１４６の範囲において、単相のＢｉＦ3型結晶が得られるという計算結果が得られた。
【００４８】
Ｌｉ3-aＢｉ1-bＧｅb結晶物質に関しては、０≦ａ≦０．０８５かつ０≦ｂ≦０．１０８の範囲において単相のＢｉＦ3型結晶が得られる範囲が算出された。この結果、（ｉ）０≦ａ≦０．００３、かつ０≦ｂ≦０．１０８、および（ｉｉ）０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［−０．１５７×（ｌｎ（ａ））2−２．２２×ｌｎ（ａ）−９．８１］の範囲において、単相のＢｉＦ3型結晶が得られるという計算結果が得られた。
【００４９】
数式ｂ＝ｅｘｐ［−０．１５７×（ｌｎ（ａ））2−２．２２×ｌｎ（ａ）−９．８１］は、（ａ、ｂ）＝（０．００３、０．１０８）、（０．００９、０．０５８）、（０．０３０、０．０１８）、および（０．０８５、０．００５）の４点をつなぐなめらかなフィッティング曲線を算出することにより得られた。化学平衡状態下では、２種類の欠陥濃度ａとｂはａAｂB＝平衡定数で表される質量作用の法則で結ばれる。ｌｎ（ｂ）をｌｎ（ａ）に関する二次多項式で表した前記フィッティング曲線は、狭い組成区間内においては質量作用の法則、すなわちｌｎ（ａ）およびｌｎ（ｂ）の間に線形な関係が存在することを根拠としている。なお、これらの４点は、単相のＬｉ3-aＢｉ1-bＧｅb結晶物質が得られるａおよびｂの範囲において、固溶限界の組成に対応する。
【００５０】
後述されるように、上述のａおよびｂの範囲において単相のＬｉ3-aＢｉ1-bＧｅb結晶物質が得られるものの、その範囲に属するすべての物質が高い性能指数ＺＴを有するわけではない。
【００５１】
（熱電変換性能指数の計算方法）
熱電変換効率は、材料の性能指数ＺＴによって決定される。ＺＴは以下の関係式（８）により定義される。
ＺＴ＝Ｓ2σＴ／（κe＋κlat）・・・（８）
ここで、
Ｓはゼーベック係数を表し、
σは電気伝導率を表し、
Ｔは評価環境の絶対温度を表し、
κeは電子の熱伝導率を表し、かつ
κlatは格子熱伝導率を表す。
Ｓ、σ、およびκeについてはＶＡＳＰコードとＢｏｌｔｚＴｒａＰコード（非特許文献４を参照せよ）とを使用し、ボルツマン輸送理論に基づく評価が行われた。σを決定する際のパラメータである電子の緩和時間τは、以下の移動度μに関する関係式（９）と非特許文献５に開示されている理論式（１０）とを連立して解く事により計算された。
μ＝ｅτ／ｍ*・・・（９）
μ＝（８π）1/2（ｈ／２π）4ｅＢ／３ｍ*5/2（ｋbＴ）3/2ｇ2・・・（１０）
ここで、
ｅは電荷素量を表し、
ｍ*はキャリアの有効質量を表し、
Ｂは弾性定数を表し、かつ
ｇは変形ポテンシャルを表す。
ｍ*、Ｂ、およびｇの値は、ＶＡＳＰコードを用いた密度汎関数法により計算された。ｇの値は、非特許文献６に開示されている関係式ｇ＝−Δε／（Δｌ／ｌ）により計算された。Δεは、格子定数ｌをΔｌ変化させた時のバンド端エネルギー準位の変化量である。
【００５２】
格子熱伝導率は、非特許文献７に開示されているＤｅｂｙｅ−Ｃａｌｌａｗａｙモデルに基づいた以下の経験式（１１）を用いて計算された。
κL＝Ａ1Ｍｖ3／Ｖ2/3ｎ1/3＋Ａ2ｖ／Ｖ2/3（１−１／ｎ2/3）・・・（１１）
ここで、
Ｍは平均原子質量を表し、
ｖは縦波音波速度を表し、
Ｖは１原子辺りの体積を表し、かつ
ｎは単位胞内に含まれる原子の数を表す。
計算において、非特許文献５に開示されているＡ1およびＡ2の値が使用された。
【００５３】
（性能指数の計算結果）
上記手法により得たＬｉ3-aＢｉ1-bＧｅbのＢｉＦ3型結晶構造が最も安定である組成範囲において、熱電変換特性を本発明者らは評価した。表１は、実施例１〜１６および比較例１〜４の３００Ｋにおける熱電変換特性の予測結果を示している。
【００５４】
【表１】
2020021783.tif 000003
【００５５】
比較例１の材料は、欠陥を有しないＬｉ3Ｂｉ結晶物質である。
【００５６】
実施例１〜４および比較例２の材料は、Ｌｉを欠損させたＬｉ3-aＢｉ結晶物質である。
【００５７】
実施例５〜８および比較例３の材料は、ＢｉをＧｅによって置換させたＬｉ3Ｂｉ1-bＧｅb結晶物質である。
【００５８】
実施例９〜１６および比較例４の材料は、Ｌｉを欠損させ、かつＢｉをＧｅによって置換したＬｉ3-aＢｉ1-bＧｅb結晶物質である。
【００５９】
表１に示されるように、実施例１〜４に関しては、０．０００３≦ａ≦０．０８５の組成範囲において、Ｓが正値となるｐ型の特性が得られ、かつ比較例１〜４を大幅に上回る、０．４４以上の高いＺＴが得られた。
【００６０】
実施例５〜８に関しては、０．０００３≦ｂ≦１．０８の組成範囲においてＳが正値となるｐ型の特性が得られ、かつ比較例１〜４を大幅に上回る、０．４２以上の高いＺＴが得られた。
【００６１】
実施例９〜１６に関しては、０．０００１≦ａ≦０．０８５、及び０．０００２≦ｂ≦０．１０８の複数の組成条件においてＳが正値となるｐ型の特性が得られ、かつ比較例１〜４を大幅に上回る、０．４２以上の高いＺＴが得られた。
【００６２】
実施例１〜１６の材料のＺＴの値は、３００Ｋにおいていずれも０．４より大きかった。このため、本開示による熱電変換材料は、２００℃以下の低い温度域における熱発電において有用であると本発明者は考える。
【００６３】
図３は、ａ―ｂ平面上に実施例１〜１６および比較例１〜４の点をプロットしたグラフである。
【００６４】
図３に示される実施例１３、実施例１４、実施例１５、および実施例１６の４つの点（すなわち、（ａ、ｂ）＝（０．００３、０．１０８）、（０．００９、０．０５８）、（０．０３０、０．０１８）、および（０．０８５、０．００５））を通過する曲線は、数式：ｂ＝ｅｘｐ［−０．１５７×（ｌｎ（ａ））2−２．２２×ｌｎ（ａ）−９．８１］により表される。
【００６５】
実施例１、実施例５、および実施例９の３つの点（すなわち、（ａ、ｂ）＝（０．０００３、０）、（０、０．０００３）および（０．０００１、０．０００２））を通過する直線は、数式：ｂ＝−ａ＋０．０００３により表される。
【００６６】
図３に含まれる斜線部は、以下の３つの数式（Ｉ）、（ＩＩ）、および（ＩＩＩ）
０≦ａ≦０．０００３、かつ−ａ＋０．０００３≦ｂ≦０．１０８  （Ｉ）、
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８  （ＩＩ）、および
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［−０．１５７×（ｌｎ（ａ））2−２．２２×ｌｎ（ａ）−９．８１］  （ＩＩＩ）
により囲まれた領域である。
【産業上の利用可能性】
【００６７】
本開示の熱電変換材料は、熱エネルギーを電気エネルギーに変換する熱電変換装置に使用可能である。
【図１】
2020021783.tif 000004
【図２】
2020021783.tif 000005
【図３】
2020021783.tif 000006
【国際調査報告】
2020021783.tif 000007
2020021783.tif 000008
2020021783.tif 000009