JP7523027B2 - 熱電変換材料および熱電変換材料を用いて電力を得る方法 - Google Patents
熱電変換材料および熱電変換材料を用いて電力を得る方法 Download PDFInfo
- Publication number
- JP7523027B2 JP7523027B2 JP2020532157A JP2020532157A JP7523027B2 JP 7523027 B2 JP7523027 B2 JP 7523027B2 JP 2020532157 A JP2020532157 A JP 2020532157A JP 2020532157 A JP2020532157 A JP 2020532157A JP 7523027 B2 JP7523027 B2 JP 7523027B2
- Authority
- JP
- Japan
- Prior art keywords
- thermoelectric conversion
- conversion material
- crystal structure
- type
- patent document
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H10—SEMICONDUCTOR DEVICES; ELECTRIC SOLID-STATE DEVICES NOT OTHERWISE PROVIDED FOR
- H10N—ELECTRIC SOLID-STATE DEVICES NOT OTHERWISE PROVIDED FOR
- H10N10/00—Thermoelectric devices comprising a junction of dissimilar materials, i.e. devices exhibiting Seebeck or Peltier effects
- H10N10/80—Constructional details
- H10N10/85—Thermoelectric active materials
- H10N10/851—Thermoelectric active materials comprising inorganic compositions
- H10N10/853—Thermoelectric active materials comprising inorganic compositions comprising arsenic, antimony or bismuth
-
- C—CHEMISTRY; METALLURGY
- C01—INORGANIC CHEMISTRY
- C01G—COMPOUNDS CONTAINING METALS NOT COVERED BY SUBCLASSES C01D OR C01F
- C01G29/00—Compounds of bismuth
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02N—ELECTRIC MACHINES NOT OTHERWISE PROVIDED FOR
- H02N11/00—Generators or motors not provided for elsewhere; Alleged perpetua mobilia obtained by electric or magnetic means
Landscapes
- Chemical & Material Sciences (AREA)
- Inorganic Chemistry (AREA)
- Organic Chemistry (AREA)
- Inorganic Compounds Of Heavy Metals (AREA)
- Crystals, And After-Treatments Of Crystals (AREA)
Description
本開示は、熱電変換材料および熱電変換材料を用いて電力を得る方法に関する。
熱電変換材料の両端に温度差が与えられると、温度差に比例した起電力が発生する。熱エネルギーが電気エネルギーに変換されるこの現象は、ゼーベック効果として知られている。熱電発電技術では、ゼーベック効果を利用し、熱エネルギーを直接電気エネルギーに変換する。
熱電変換材料の技術分野においてよく知られているように、熱電変換デバイスに用いられる熱電変換材料の性能は、性能指数ＺＴにより評価される。ＺＴは、以下の数式（Ｉ）により表される。
ＺＴ＝Ｓ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型の結晶構造を有し、かつ
前記熱電変換材料はｐ型の極性を有し、かつ
以下の数式（Ｉ）～（ＩＩＩ）のいずれか１つが充足される。
０≦ａ≦０．０００３、かつ－ａ＋０．０００３≦ｂ≦０．１０８ （Ｉ）；
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８ （ＩＩ）；または
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［－０．１５７×（ｌｎ（ａ））2－２．２２×ｌｎ（ａ）－９．８１］ （ＩＩＩ）
本開示は、新規な熱電変換材料を提供する。
以下、本開示の実施形態について、図面を参照しながら説明する。
本開示の熱電変換材料は、以下の化学式（Ｉ）により表される組成を有する。
Ｌｉ3-aＢｉ1-bＧｅb （Ｉ）
Ｌｉ3-aＢｉ1-bＧｅb （Ｉ）
図１は、Ｌｉ3Ｂｉ結晶構造の模式図を示す。非特許文献１にも開示されているように、Ｌｉ3Ｂｉ結晶構造は、に、ＢｉＦ3型結晶構造またはＡｌＣｕ2Ｍｎ型結晶構造を有する。ＢｉＦ3型結晶構造およびＡｌＣｕ2Ｍｎ型結晶構造の両者は、空間群Ｆｍ－３ｍに属する。非特許文献１においては、Ｌｉ3Ｂｉ結晶物質が熱電変換材料として扱われていない。そのため、非特許文献１には、性能指数ＺＴは開示されていない。
本発明者らは、材料インフォマティクスと呼ばれる、データ科学に基づいた材料探索手法を利用することで、無機結晶構造データベース中の数万件の化合物について性能指数ＺＴの予測値を算出した。算出には、本発明者らが独自に確立した性能指数ＺＴの予測モデルが用いられた。この予測モデルは、従来の手法よりも精度が高い。このため、この予測モデルを用いることで、従来よりも信頼性の高い予測結果を得ることができる。そこで、本発明者らは、これまで熱電変換材料として扱われてこなかったＬｉ3Ｂｉ結晶物質が熱電変換材料として有望な材料であるかどうかを検討した。
Ｌｉ3Ｂｉ結晶物質が欠陥を含まない状態において、Ｌｉ3Ｂｉ結晶物質は、キャリアに乏しい。このため、欠陥のないＬｉ3Ｂｉ結晶物質が高い性能指数ＺＴを有することは望めない。そこで、本発明者らは、性能指数ＺＴを向上させるために、Ｌｉ3Ｂｉ結晶物質に欠陥を導入し、ｐ型のキャリアを生じさせることを検討した。その結果、本発明者らは、Ｌｉサイトを空孔で置換することによって生じる欠陥を有する物質、およびＢｉサイトを元素Ｇｅで置換することによって生じる欠陥を有する物質の２種類の物質を想到するに至った。
本発明者らは、計算により、Ｌｉ3-aＢｉ1-bＧｅb結晶物質のみが安定的に得られるａおよびｂの値の範囲を導き出した。さらに、その範囲において性能指数ＺＴを算出することで、本発明者らは、０．４以上の高い性能指数ＺＴが得られるａおよびｂの値の範囲を見出した。具体的には、後述の実施例１～１６および比較例１～４において実証されているように、以下の数式（Ｉ）、（ＩＩ）、または（ＩＩＩ）のうちいずれか１つを充足する場合に、熱電変換材料は０．４以上の高い性能指数ＺＴを有する。
０≦ａ≦０．０００３、かつ－ａ＋０．０００３≦ｂ≦０．１０８ （Ｉ）
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８ （ＩＩ）、または
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［－０．１５７×（ｌｎ（ａ））2－２．２２×ｌｎ（ａ）－９．８１］ （ＩＩＩ）
具体的には、ａおよびｂの値の上記の範囲において、性能指数ＺＴの値は０．４２以上である。表１を参照せよ。
０≦ａ≦０．０００３、かつ－ａ＋０．０００３≦ｂ≦０．１０８ （Ｉ）
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８ （ＩＩ）、または
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［－０．１５７×（ｌｎ（ａ））2－２．２２×ｌｎ（ａ）－９．８１］ （ＩＩＩ）
具体的には、ａおよびｂの値の上記の範囲において、性能指数ＺＴの値は０．４２以上である。表１を参照せよ。
（製造方法）
本開示の熱電変換材料の製造方法の一例が、非特許文献１および非特許文献８の開示内容に基づいて以下に説明される。
本開示の熱電変換材料の製造方法の一例が、非特許文献１および非特許文献８の開示内容に基づいて以下に説明される。
非特許文献１は、Ｌｉ３Ｂｉ結晶物質の製造方法を開示している。。非特許文献１に開示された製造方法においては、まず、アルゴンガスで満たされたグローブボックスの中で、Ｌｉ箔および粒状のＢｉがモル比４３：１２でアルミナ製の坩堝に入れられる。次に、坩堝の内部で、これらのＬｉ箔および粒状のＢｉが共に圧縮される。カーボンでコーティングされた石英管の中に、坩堝は入れられる。石英管の内部は、その真空度が６．６７×１０-3Ｐａに達するまで真空引きされる。石英管は、酸素ガスの炎を用いて密封される。密封された石英管はマイクロ波キャビティの中に入れられ、マイクロ波が石英管に対して照射される。このようにして、Ｌｉ3Ｂｉ結晶物質が得られる。グローブボックスの中で石英管を開封することにより、石英管から得られたＬｉ3Ｂｉ結晶物質が取り出される。
非特許文献８もまた、、Ｌｉ３Ｂｉ結晶物質の製造方法を開示している。非特許文献８に開示された製造方法においては、まず、不活性ガスで満たされたグローブボックスに接続された誘導鋳造炉に、タンタル製の坩堝が収容される。ＬｉおよびＢｉが、Ｌｉ：Ｂｉのモル比３：１で坩堝の中に導入される。炉の内部がゆっくりと加熱される。炉の内部の加熱にしたがってまず、Ｌｉが融解する。Ｌｉの融解の後にＢｉが融解する。このようにして、ＬｉのＢｉとの反応が進行し、反応物を得る。得られた反応物が、１２００℃まで加熱されて均質化される。次いで反応物は冷却されて、Ｌｉ3Ｂｉ結晶物質が得られる。
非特許文献１および非特許文献８に開示された両方の製法において、Ｌｉの欠損量は、出発物質のＬｉの量を、Ｂｉの量に対して変化させることで制御できる。
また、元素Ｂｉの原子半径（０．１５６ｎｍ）は、元素Ｇｅの原子半径（０．１２２ｎｍ）と大差がない。このため、出発物質としてＧｅを加えることで、一部の元素Ｂｉは元素Ｇｅによって置換されうると考えられる。Ｇｅによる置換量は、出発物質のＧｅの量を、Ｂｉの量に対して変化させることで制御できる。
以上から、本開示によるＬｉ3-aＢｉ1-bＧｅb結晶物質は、非特許文献１および非特許文献８に開示されている製造方法を参照することで製造できると考えられる。
（熱電変換材料を用いて電力を得る方法）
本実施形態において、本開示のＬｉ3-aＢｉ1-bＧｅb結晶物質に温度差を印加することで電力が得られる。
本実施形態において、本開示のＬｉ3-aＢｉ1-bＧｅb結晶物質に温度差を印加することで電力が得られる。
（実施例）
以下の実施例を参照しながら、本開示の熱電変換材料がさらに詳細に説明される。
以下の実施例を参照しながら、本開示の熱電変換材料がさらに詳細に説明される。
（結晶構造の解析）
非特許文献１は、Ｌｉ3Ｂｉ結晶物質の結晶構造がＢｉＦ3型構造に属することを開示している。Ｘ線結晶回折法による測定に基づき、ＢｉＦ3型構造のＸ線回折ピークを確認する事ができる。図２は、ソフトウェア（ＲＩＥＴＡＮ、入手先ＵＲＬ：http://fujioizumi.verse.jp/download/download.html）を用いてＬｉ3Ｂｉの結晶構造因子Ｆおよび積分回折強度Ｉを計算することにより得られたＬｉ3Ｂｉ結晶構造の回折Ｘ線強度分布を示すグラフである。
非特許文献１は、Ｌｉ3Ｂｉ結晶物質の結晶構造がＢｉＦ3型構造に属することを開示している。Ｘ線結晶回折法による測定に基づき、ＢｉＦ3型構造のＸ線回折ピークを確認する事ができる。図２は、ソフトウェア（ＲＩＥＴＡＮ、入手先ＵＲＬ：http://fujioizumi.verse.jp/download/download.html）を用いてＬｉ3Ｂｉの結晶構造因子Ｆおよび積分回折強度Ｉを計算することにより得られたＬｉ3Ｂｉ結晶構造の回折Ｘ線強度分布を示すグラフである。
結晶構造因子Ｆは、以下の関係式（１）によって求められた。
Ｆ＝Σｆiｅｘｐ（２πｉｒiΔｋ）・・・（１）
ここで、
ｒiは結晶中の原子の位置ベクトルを表し
ｆiはｒiの位置にある原子の原子散乱因子を表し、かつ
Δｋは散乱前後におけるＸ線の波数ベクトルの差を表す。
Ｆ＝Σｆiｅｘｐ（２πｉｒiΔｋ）・・・（１）
ここで、
ｒiは結晶中の原子の位置ベクトルを表し
ｆiはｒiの位置にある原子の原子散乱因子を表し、かつ
Δｋは散乱前後におけるＸ線の波数ベクトルの差を表す。
積分回折強度Ｉは、以下の関係式（２）によって求められた。
Ｉ＝ＩeＬ｜Ｆ｜2Ｎ2・・・（２）
ここで、
Ｉeは１個の電子の散乱強度を表し、
Ｎは結晶中の単位格子の数を表し、
Ｌは吸引因子を含み、かつ実験条件に依存する係数を表す。
Ｉ＝ＩeＬ｜Ｆ｜2Ｎ2・・・（２）
ここで、
Ｉeは１個の電子の散乱強度を表し、
Ｎは結晶中の単位格子の数を表し、
Ｌは吸引因子を含み、かつ実験条件に依存する係数を表す。
本開示による熱電変換材料は、Ｌｉ3Ｂｉ結晶において、元素Ｌｉサイトの一部が空孔により置換されることにより生じた欠陥、または元素Ｂｉの一部が元素Ｇｅによって置換された欠陥を有している。これらの欠陥により、結晶構造に格子変形が生じうる。このため、本開示による熱電変換材料の回折Ｘ線強度分布のピークは、図２に示され、かつ非特許文献１に開示されたＬｉ3Ｂｉ結晶物質のピークに対して若干シフトしていると予測される。
（安定組成範囲の計算方法）
本発明者らは、密度汎関数法（以下、「ＤＦＴ」という）に基づいた計算により、Ｌｉ3Ｂｉ結晶物質のバンド構造を評価した。その結果、本発明者らは、Ｌｉ3Ｂｉ結晶物質は半導体であるという結果を得た。ＤＦＴによる計算には、本発明者らは、第一原理電子状態計算プログラム（Vienna Ab initio Simulation Package、入手先:https://www.vasp.at/）を使用した。このプログラムは、以下において、「ＶＡＳＰ」と称される。
本発明者らは、密度汎関数法（以下、「ＤＦＴ」という）に基づいた計算により、Ｌｉ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（μ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（μi，ｑ，ＥF）＝Ｎsite×ｅｘｐ［－Ｅform（μi，ｑ，ＥF）／ｋbＴ］・・・（４）
ここで、
Ｎsiteは考慮している欠陥が生じうるサイトの体積密度を表し、
ｋbはボルツマン定数を表し、かつ
Ｔは絶対温度を表す。
各欠陥の電荷ｑ×ＮDと、半導体にドープされたキャリアの電荷Ｑeとの総量が０になるという電荷中性条件は常に満たされている。この電荷中性条件から、フェルミエネルギーおよびキャリア濃度が各組成において決定された。このときの価電子帯のキャリア濃度ｐと、伝導体のキャリア濃度ｎとは、以下の関係式（５）～（７）により求められた。
ｐ＝１－∫ＤVB（Ｅ）［１－ｆ（Ｅ；ＥF）］ｄＥ・・・（５）
ｎ＝∫ＤCB（Ｅ）ｆ（Ｅ；ＥF）ｄＥ・・・（６）
ｆ（Ｅ；ＥF）＝１／［ｅｘｐ（（Ｅ－ＥF）／ｋBＴ）＋１］・・・（７）
ここで、
ＤVB（Ｅ）およびＤCB（Ｅ）は、それぞれ、ＤＦＴ計算によって得られる価電子帯および伝導帯の電子状態密度を表し、
ｆ（Ｅ；ＥF）はフェルミ分布関数を表し、かつ
Ｔは熱電変換特性の評価を実施する絶対温度を表す。
全キャリアの電荷密度Ｑeは、数式Ｑ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・・・（１０）
ここで、
ｅは電荷素量を表し、
ｍ*はキャリアの有効質量を表し、
Ｂは弾性定数を表し、かつ
ｇは変形ポテンシャルを表す。
ｍ*、Ｂ、およびｇの値は、ＶＡＳＰコードを用いた密度汎関数法により計算された。ｇの値は、非特許文献６に開示されている関係式ｇ＝－Δε／（Δｌ／ｌ）により計算された。Δεは、格子定数ｌをΔｌ変化させた時のバンド端エネルギー準位の変化量である。
熱電変換効率は、材料の性能指数ＺＴによって決定される。ＺＴは以下の関係式（８）により定義される。
ＺＴ＝Ｓ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の値が使用された。
κL＝Ａ1Ｍｖ3／Ｖ2/3ｎ1/3＋Ａ2ｖ／Ｖ2/3（１－１／ｎ2/3）・・・（１１）
ここで、
Ｍは平均原子質量を表し、
ｖは縦波音波速度を表し、
Ｖは１原子辺りの体積を表し、かつ
ｎは単位胞内に含まれる原子の数を表す。
計算において、非特許文献５に開示されているＡ1およびＡ2の値が使用された。
（性能指数の計算結果）
上記手法により得たＬｉ3-aＢｉ1-bＧｅbのＢｉＦ3型結晶構造が最も安定である組成範囲において、熱電変換特性を本発明者らは評価した。表１は、実施例１～１６および比較例１～４の３００Ｋにおける熱電変換特性の予測結果を示している。
上記手法により得たＬｉ3-aＢｉ1-bＧｅbのＢｉＦ3型結晶構造が最も安定である組成範囲において、熱電変換特性を本発明者らは評価した。表１は、実施例１～１６および比較例１～４の３００Ｋにおける熱電変換特性の予測結果を示している。
比較例１の材料は、欠陥を有しないＬｉ3Ｂｉ結晶物質である。
実施例１～４および比較例２の材料は、Ｌｉを欠損させたＬｉ3-aＢｉ結晶物質である。
実施例５～８および比較例３の材料は、ＢｉをＧｅによって置換させたＬｉ3Ｂｉ1-bＧｅb結晶物質である。
実施例９～１６および比較例４の材料は、Ｌｉを欠損させ、かつＢｉをＧｅによって置換したＬｉ3-aＢｉ1-bＧｅb結晶物質である。
表１に示されるように、実施例１～４に関しては、０．０００３≦ａ≦０．０８５の組成範囲において、Ｓが正値となるｐ型の特性が得られ、かつ比較例１～４を大幅に上回る、０．４４以上の高いＺＴが得られた。
実施例５～８に関しては、０．０００３≦ｂ≦１．０８の組成範囲においてＳが正値となるｐ型の特性が得られ、かつ比較例１～４を大幅に上回る、０．４２以上の高いＺＴが得られた。
実施例９～１６に関しては、０．０００１≦ａ≦０．０８５、及び０．０００２≦ｂ≦０．１０８の複数の組成条件においてＳが正値となるｐ型の特性が得られ、かつ比較例１～４を大幅に上回る、０．４２以上の高いＺＴが得られた。
実施例１～１６の材料のＺＴの値は、３００Ｋにおいていずれも０．４より大きかった。このため、本開示による熱電変換材料は、２００℃以下の低い温度域における熱発電において有用であると本発明者は考える。
図３は、ａ―ｂ平面上に実施例１～１６および比較例１～４の点をプロットしたグラフである。
図３に示される実施例１３、実施例１４、実施例１５、および実施例１６の４つの点（すなわち、（ａ、ｂ）＝（０．００３、０．１０８）、（０．００９、０．０５８）、（０．０３０、０．０１８）、および（０．０８５、０．００５））を通過する曲線は、数式：ｂ＝ｅｘｐ［－０．１５７×（ｌｎ（ａ））2－２．２２×ｌｎ（ａ）－９．８１］により表される。
実施例１、実施例５、および実施例９の３つの点（すなわち、（ａ、ｂ）＝（０．０００３、０）、（０、０．０００３）および（０．０００１、０．０００２））を通過する直線は、数式：ｂ＝－ａ＋０．０００３により表される。
図３に含まれる斜線部は、以下の３つの数式（Ｉ）、（ＩＩ）、および（ＩＩＩ）
０≦ａ≦０．０００３、かつ－ａ＋０．０００３≦ｂ≦０．１０８ （Ｉ）、
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８ （ＩＩ）、および
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［－０．１５７×（ｌｎ（ａ））2－２．２２×ｌｎ（ａ）－９．８１］ （ＩＩＩ）
により囲まれた領域である。
０≦ａ≦０．０００３、かつ－ａ＋０．０００３≦ｂ≦０．１０８ （Ｉ）、
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８ （ＩＩ）、および
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［－０．１５７×（ｌｎ（ａ））2－２．２２×ｌｎ（ａ）－９．８１］ （ＩＩＩ）
により囲まれた領域である。
本開示の熱電変換材料は、熱エネルギーを電気エネルギーに変換する熱電変換装置に使用可能である。
Claims (4)
- 化学式Ｌｉ3-aＢｉ1-bＧｅbにより表される組成を有する熱電変換材料。
ここで、
前記熱電変換材料はＢｉＦ3型の結晶構造を有し、
前記熱電変換材料はｐ型の極性を有し、かつ
以下の数式（Ｉ）～（ＩＩＩ）のいずれか１つが充足される。
０≦ａ≦０．０００３、かつ－ａ＋０．０００３≦ｂ≦０．１０８ （Ｉ）；
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８ （ＩＩ）；または
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［－０．１５７×（ｌｎ（ａ））2－２．２２×ｌｎ（ａ）－９．８１］ （ＩＩＩ） - 請求項１に記載の熱電変換材料であって、
０．０００３≦ａ≦０．０８５、かつｂ＝０である。 - 請求項１に記載の熱電変換材料であって、
ａ＝０、かつ０．０００３≦ｂ≦０．１０８である。 - 熱電変換材料を用いて電力を得る方法であって、
前記方法は、前記熱電変換材料に温度差を印加して電力を生じる工程を具備し、
前記熱電変換材料は、 化学式Ｌｉ3-aＢｉ1-bＧｅbにより表される組成を有し、
ＢｉＦ3型の結晶構造を有し、かつ
ｐ型の極性を有し、かつ
以下の数式（Ｉ）～（ＩＩＩ）のいずれか１つが充足される。 ０≦ａ≦０．０００３、かつ－ａ＋０．０００３≦ｂ≦０．１０８ （Ｉ）；
０．０００３≦ａ≦０．００３、かつ０≦ｂ≦０．１０８ （ＩＩ）；または
０．００３≦ａ≦０．０８５、かつ０≦ｂ≦ｅｘｐ［－０．１５７×（ｌｎ（ａ））2－２．２２×ｌｎ（ａ）－９．８１］ （ＩＩＩ）
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018137695 | 2018-07-23 | ||
| JP2018137695 | 2018-07-23 | ||
| PCT/JP2019/015881 WO2020021783A1 (ja) | 2018-07-23 | 2019-04-12 | 熱電変換材料および熱電変換材料を用いて電力を得る方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2020021783A1 JPWO2020021783A1 (ja) | 2021-08-02 |
| JP7523027B2 true JP7523027B2 (ja) | 2024-07-26 |
Family
ID=69180397
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020532157A Active JP7523027B2 (ja) | 2018-07-23 | 2019-04-12 | 熱電変換材料および熱電変換材料を用いて電力を得る方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US11437558B2 (ja) |
| JP (1) | JP7523027B2 (ja) |
| CN (1) | CN111819704A (ja) |
| WO (1) | WO2020021783A1 (ja) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016535392A (ja) | 2013-10-29 | 2016-11-10 | マサチューセッツ インスティテュート オブ テクノロジー | 自己回復液体／固体状態バッテリ |
| US20170162864A1 (en) | 2014-02-03 | 2017-06-08 | Ramot At Tel-Aviv University Ltd. | Electrode compositions and alkali metal batteries comprising same |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4171795B2 (ja) * | 1998-05-22 | 2008-10-29 | 独立行政法人産業技術総合研究所 | 熱電変換素子用酸化物部材 |
| US6312617B1 (en) * | 1998-10-13 | 2001-11-06 | Board Of Trustees Operating Michigan State University | Conductive isostructural compounds |
| CN101523627A (zh) * | 2006-08-24 | 2009-09-02 | 住友化学株式会社 | 热电转换材料、其制造方法以及热电转换元件 |
| DE102007039060B4 (de) * | 2007-08-17 | 2019-04-25 | Evonik Degussa Gmbh | Thermokraftelement oder Peltier-Elemente aus gesinterten Nanokristallen aus Silicium, Germanium oder Silicium-Germanium-Legierungen |
| JP6127281B1 (ja) * | 2015-10-27 | 2017-05-17 | パナソニックＩｐマネジメント株式会社 | 熱電変換材料 |
-
2019
- 2019-04-12 WO PCT/JP2019/015881 patent/WO2020021783A1/ja not_active Ceased
- 2019-04-12 JP JP2020532157A patent/JP7523027B2/ja active Active
- 2019-04-12 CN CN201980017532.7A patent/CN111819704A/zh active Pending
-
2020
- 2020-10-23 US US17/079,403 patent/US11437558B2/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016535392A (ja) | 2013-10-29 | 2016-11-10 | マサチューセッツ インスティテュート オブ テクノロジー | 自己回復液体／固体状態バッテリ |
| US20170162864A1 (en) | 2014-02-03 | 2017-06-08 | Ramot At Tel-Aviv University Ltd. | Electrode compositions and alkali metal batteries comprising same |
Non-Patent Citations (1)
| Title |
|---|
| Xiuxian Yang，外4名，Low lattice thermal conductivity and excellent thermoelectric behavior in Li3Sb and Li3Bi，Journal of Physics: Condensed Matter，英国，IOP Publishing Ltd，2018年10月03日，Vol. 30，425401-1 - 425401-8 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPWO2020021783A1 (ja) | 2021-08-02 |
| US20210043822A1 (en) | 2021-02-11 |
| CN111819704A (zh) | 2020-10-23 |
| WO2020021783A1 (ja) | 2020-01-30 |
| US11437558B2 (en) | 2022-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sajjad et al. | Ultralow lattice thermal conductivity in double perovskite Cs2PtI6: a promising thermoelectric material | |
| Bafekry et al. | Semiconducting chalcogenide alloys based on the (Ge, Sn, Pb)(S, Se, Te) formula with outstanding properties: a first-principles calculation study | |
| Gharsallah et al. | Giant seebeck effect in Ge-doped SnSe | |
| Wei et al. | Thermoelectric figure-of-merit of fully dense single-crystalline SnSe | |
| Pei et al. | Convergence of electronic bands for high performance bulk thermoelectrics | |
| Hua et al. | Phase Stability and Ordering in Rock Salt-Based Thermoelectrics: NaSbX2, AgSbX2, and Their Alloys with PbX and SnX (X= S, Se, Te) | |
| Even et al. | Density Functional Theory Simulations of Semiconductors for Photovoltaic Applications: Hybrid Organic‐Inorganic Perovskites and III/V Heterostructures | |
| Yan et al. | Achieving high average power factor in tetrahedrite Cu12Sb4S13 via regulating electron-phonon coupling strength | |
| Wen et al. | Thermoelectric properties of p-Type Cu3VSe4 with high seebeck coefficients | |
| Pan et al. | A2Cu3In3Te8 (A= Cd, Zn, Mn, Mg): a type of thermoelectric material with complex diamond-like structure and low lattice thermal conductivities | |
| Hariharan et al. | Comprehensive study on structural, electronic, optical, elastic, and transport properties of natural mercury sulphohalides via DFT computation | |
| Wafula et al. | DFT+ U investigation of electronic, optical, and thermoelectric properties of YAuX (X= Si or Ge or Sn) half-Heusler alloys | |
| Goumri-Said | A comprehensive study of electronic, optical, and thermoelectric characteristics of Cs2PbI2Br2 inorganic layered ruddlesden–popper mixed halide perovskite through systematic first-principles analysis | |
| Iyer et al. | Li3Sb7Se12: A Pavonite Homologue with Ultralow Thermal Conductivity | |
| Bouchrit et al. | First principle investigation on electronic and thermoelectric properties of CoTiP half-Heusler with Co, Ti, and P vacancies | |
| Tassanov et al. | Layered NaBa2M3Q3 (Q2)(M= Cu or Ag; Q= S or Se) chalcogenides and local ordering in their mixed-anion compositions | |
| JP7523027B2 (ja) | 熱電変換材料および熱電変換材料を用いて電力を得る方法 | |
| JP7523028B2 (ja) | 熱電変換材料および熱電変換材料を用いて電力を得る方法 | |
| JP7523029B2 (ja) | 熱電変換材料および熱電変換材料を用いて電力を得る方法 | |
| Wang et al. | Ternary compounds Cu3 R Te3 (R= Y, Sm, and Dy): A family of new thermoelectric materials with trigonal structures | |
| Lakshan et al. | Disorder-Mediated Structural Transformation in the Cu4TiSe4–x S x (0≤ x≤ 4) System and Its Effects on the Thermal Transport Property | |
| Ali et al. | Exploration of physical aspects of Li2AgAsZ6 (Z= F, Cl, Br, I) double perovskites for energy harvesting perspectives | |
| Velu et al. | Effects of Te-doping on the thermoelectric properties of InGaSb crystals | |
| Tengå et al. | Sphalerite− Chalcopyrite Polymorphism in Semimetallic ZnSnSb2 | |
| Chan et al. | First-principles computational design and synthesis of hybrid carbon–silicon clathrates |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220208 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230502 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20231010 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20231205 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20240118 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20240604 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240701 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7523027 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |