JP7420040B2 - 界面構造探索方法 - Google Patents
界面構造探索方法 Download PDFInfo
- Publication number
- JP7420040B2 JP7420040B2 JP2020166170A JP2020166170A JP7420040B2 JP 7420040 B2 JP7420040 B2 JP 7420040B2 JP 2020166170 A JP2020166170 A JP 2020166170A JP 2020166170 A JP2020166170 A JP 2020166170A JP 7420040 B2 JP7420040 B2 JP 7420040B2
- Authority
- JP
- Japan
- Prior art keywords
- stable
- liquid phase
- model
- descriptor
- initial
- 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
- 238000000034 method Methods 0.000 title claims description 50
- 239000007791 liquid phase Substances 0.000 claims description 43
- 238000004364 calculation method Methods 0.000 claims description 30
- 239000007790 solid phase Substances 0.000 claims description 23
- 238000012417 linear regression Methods 0.000 claims description 6
- 125000004429 atom Chemical group 0.000 description 13
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 12
- 229910052760 oxygen Inorganic materials 0.000 description 12
- 239000001301 oxygen Substances 0.000 description 12
- 239000007774 positive electrode material Substances 0.000 description 12
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 11
- 239000003792 electrolyte Substances 0.000 description 11
- 229910001416 lithium ion Inorganic materials 0.000 description 11
- 230000000052 comparative effect Effects 0.000 description 8
- 229910052744 lithium Inorganic materials 0.000 description 7
- 238000003775 Density Functional Theory Methods 0.000 description 6
- 229910013870 LiPF 6 Inorganic materials 0.000 description 5
- 239000007788 liquid Substances 0.000 description 5
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000006866 deterioration Effects 0.000 description 4
- 150000002500 ions Chemical class 0.000 description 4
- 150000003839 salts Chemical class 0.000 description 4
- KMTRUDSVKNLOMY-UHFFFAOYSA-N Ethylene carbonate Chemical compound O=C1OCCO1 KMTRUDSVKNLOMY-UHFFFAOYSA-N 0.000 description 3
- 229910013290 LiNiO 2 Inorganic materials 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000008151 electrolyte solution Substances 0.000 description 3
- JBTWLSYIZRCDFO-UHFFFAOYSA-N ethyl methyl carbonate Chemical compound CCOC(=O)OC JBTWLSYIZRCDFO-UHFFFAOYSA-N 0.000 description 3
- 239000002905 metal composite material Substances 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 2
- WYURNTSHIVDZCO-UHFFFAOYSA-N Tetrahydrofuran Chemical compound C1CCOC1 WYURNTSHIVDZCO-UHFFFAOYSA-N 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 239000002608 ionic liquid Substances 0.000 description 2
- 229910003002 lithium salt Inorganic materials 0.000 description 2
- 159000000002 lithium salts Chemical class 0.000 description 2
- 239000003960 organic solvent Substances 0.000 description 2
- 230000005501 phase interface Effects 0.000 description 2
- 238000007637 random forest analysis Methods 0.000 description 2
- ZZXUZKXVROWEIF-UHFFFAOYSA-N 1,2-butylene carbonate Chemical compound CCC1COC(=O)O1 ZZXUZKXVROWEIF-UHFFFAOYSA-N 0.000 description 1
- YBJCDTIWNDBNTM-UHFFFAOYSA-N 1-methylsulfonylethane Chemical compound CCS(C)(=O)=O YBJCDTIWNDBNTM-UHFFFAOYSA-N 0.000 description 1
- RNFJDJUURJAICM-UHFFFAOYSA-N 2,2,4,4,6,6-hexaphenoxy-1,3,5-triaza-2$l^{5},4$l^{5},6$l^{5}-triphosphacyclohexa-1,3,5-triene Chemical compound N=1P(OC=2C=CC=CC=2)(OC=2C=CC=CC=2)=NP(OC=2C=CC=CC=2)(OC=2C=CC=CC=2)=NP=1(OC=1C=CC=CC=1)OC1=CC=CC=C1 RNFJDJUURJAICM-UHFFFAOYSA-N 0.000 description 1
- YEVQZPWSVWZAOB-UHFFFAOYSA-N 2-(bromomethyl)-1-iodo-4-(trifluoromethyl)benzene Chemical compound FC(F)(F)C1=CC=C(I)C(CBr)=C1 YEVQZPWSVWZAOB-UHFFFAOYSA-N 0.000 description 1
- JWUJQDFVADABEY-UHFFFAOYSA-N 2-methyltetrahydrofuran Chemical compound CC1CCCO1 JWUJQDFVADABEY-UHFFFAOYSA-N 0.000 description 1
- GKZFQPGIDVGTLZ-UHFFFAOYSA-N 4-(trifluoromethyl)-1,3-dioxolan-2-one Chemical compound FC(F)(F)C1COC(=O)O1 GKZFQPGIDVGTLZ-UHFFFAOYSA-N 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- OIFBSDVPJOWBCH-UHFFFAOYSA-N Diethyl carbonate Chemical compound CCOC(=O)OCC OIFBSDVPJOWBCH-UHFFFAOYSA-N 0.000 description 1
- XTHFKEDIFFGKHM-UHFFFAOYSA-N Dimethoxyethane Chemical compound COCCOC XTHFKEDIFFGKHM-UHFFFAOYSA-N 0.000 description 1
- YCKRFDGAMUMZLT-UHFFFAOYSA-N Fluorine atom Chemical compound [F] YCKRFDGAMUMZLT-UHFFFAOYSA-N 0.000 description 1
- 229910015015 LiAsF 6 Inorganic materials 0.000 description 1
- 229910013063 LiBF 4 Inorganic materials 0.000 description 1
- 229910013684 LiClO 4 Inorganic materials 0.000 description 1
- 229910013131 LiN Inorganic materials 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 150000001450 anions Chemical class 0.000 description 1
- 229910052792 caesium Inorganic materials 0.000 description 1
- TVFDJXOCXUVLDH-UHFFFAOYSA-N caesium atom Chemical compound [Cs] TVFDJXOCXUVLDH-UHFFFAOYSA-N 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 1
- 150000001768 cations Chemical class 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 150000005676 cyclic carbonates Chemical class 0.000 description 1
- 238000003066 decision tree Methods 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- IEJIGPNLZYLLBP-UHFFFAOYSA-N dimethyl carbonate Chemical compound COC(=O)OC IEJIGPNLZYLLBP-UHFFFAOYSA-N 0.000 description 1
- VUPKGFBOKBGHFZ-UHFFFAOYSA-N dipropyl carbonate Chemical compound CCCOC(=O)OCCC VUPKGFBOKBGHFZ-UHFFFAOYSA-N 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 150000002170 ethers Chemical class 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000003063 flame retardant Substances 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 239000011737 fluorine Substances 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 229910052736 halogen Inorganic materials 0.000 description 1
- 150000002367 halogens Chemical class 0.000 description 1
- 230000012447 hatching Effects 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000329 molecular dynamics simulation Methods 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- 239000010450 olivine Substances 0.000 description 1
- 229910052609 olivine Inorganic materials 0.000 description 1
- MHYFEEDKONKGEB-UHFFFAOYSA-N oxathiane 2,2-dioxide Chemical compound O=S1(=O)CCCCO1 MHYFEEDKONKGEB-UHFFFAOYSA-N 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 150000003018 phosphorus compounds Chemical class 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- RUOJZAUFBMNUDX-UHFFFAOYSA-N propylene carbonate Chemical compound CC1COC(=O)O1 RUOJZAUFBMNUDX-UHFFFAOYSA-N 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000002516 radical scavenger Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 229910052596 spinel Inorganic materials 0.000 description 1
- 239000011029 spinel Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 150000003464 sulfur compounds Chemical class 0.000 description 1
- 239000004094 surface-active agent Substances 0.000 description 1
- YLQBMQCUIZJEEH-UHFFFAOYSA-N tetrahydrofuran Natural products C=1C=COC=1 YLQBMQCUIZJEEH-UHFFFAOYSA-N 0.000 description 1
- -1 that is Chemical compound 0.000 description 1
- DQWPFSLDHJDLRL-UHFFFAOYSA-N triethyl phosphate Chemical compound CCOP(=O)(OCC)OCC DQWPFSLDHJDLRL-UHFFFAOYSA-N 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy storage using batteries
Landscapes
- Secondary Cells (AREA)
- Battery Electrode And Active Subsutance (AREA)
Description
本発明は、界面構造探索方法に関する。
固相－液相間の界面挙動は、腐食、劣化などの現象を解明、制御する上でに非常に重要である。
固相－液相間の界面挙動の１つとして、例えばリチウムイオン二次電池においてリチウムイオンが正極材から電解液へと移動する反応も挙げられ、リチウムイオン二次電池の性能向上等の観点から、係る界面挙動を明らかにすることも求められている。このように、固相－液相間の界面挙動を明らかにすることは、腐食、劣化等の反応を抑制する上で、あるいはリチウムイオン二次電池のような各種デバイス等の性能を向上させる上で極めて重要である。
しかしながら、固相－液相の界面の状態は極めて微小なスケールの問題であり、界面の状態を分析的手法により明らかにすることは難易度が高く、分析的手法では断片的な情報しか得ることができない。
そのため、実験や分析によらず第一原理計算等の理論計算手法を用いて、固相－液相界面における劣化反応などを明らかにする手法が提案されている。例えば、日射遮蔽材料として用いられるＣｓ０．３３ＷＯ３の水分子との界面における劣化反応に着目して解析した具体的な例が特許文献１、非特許文献１などに報告されている。
ところで、例えばリチウムイオン二次電池の電解液は複数種の分子を含有している。このように液相が複数種の分子を含有する場合、安定な固相－液相界面を決定することが難しいという問題があった。このため、安定な固相－液相界面が推定しにくい構造においても、計算コストを抑えた上で安定な固相－液相界面を推定する手法が求められている。
S. Yoshio, K. Adachi, and M. Kubo, "Cesium desorption mechanism in Cs 0.33 WO 3 by first-principles molecular dynamics calculations," J. Appl. Phys., 2019, vol. 126, no. 7, p. 073101.
そこで上記従来技術が有する問題に鑑み、本発明の一側面では、固相－液相の安定な界面構造を容易に推定できる界面構造探索方法を提供することを目的とする。
上記課題を解決するため本発明の一態様によれば、
固相と、液相との界面構造を探索する界面構造探索方法であって、
前記固相の表面に、前記液相が含有する分子が配置され、前記液相が含有する前記分子の配置が異なる複数の初期構造を決定する初期構造決定工程と、
複数の前記初期構造について構造緩和計算を行い、それぞれ安定構造を探索する安定構造探索工程と、
前記安定構造探索工程で求めた安定構造を表す記述子を作成する記述子作成工程と、
前記記述子を用いて、前記界面構造のエネルギーを予測するモデルを作成するモデル作成工程と、
前記モデルを用いて、前記界面構造を決定する界面構造決定工程と、を有する界面構造探索方法を提供する。
固相と、液相との界面構造を探索する界面構造探索方法であって、
前記固相の表面に、前記液相が含有する分子が配置され、前記液相が含有する前記分子の配置が異なる複数の初期構造を決定する初期構造決定工程と、
複数の前記初期構造について構造緩和計算を行い、それぞれ安定構造を探索する安定構造探索工程と、
前記安定構造探索工程で求めた安定構造を表す記述子を作成する記述子作成工程と、
前記記述子を用いて、前記界面構造のエネルギーを予測するモデルを作成するモデル作成工程と、
前記モデルを用いて、前記界面構造を決定する界面構造決定工程と、を有する界面構造探索方法を提供する。
本発明の一態様によれば、固相-液相の安定な界面構造を容易に推定できる界面構造探索方法を提供することができる。
以下、本発明を実施するための形態について説明するが、本発明は、下記の実施形態に制限されることはなく、本発明の範囲を逸脱することなく、下記の実施形態に種々の変形および置換を加えることができる。
［界面構造探索方法］
本実施形態の界面構造探索方法について説明する。
［界面構造探索方法］
本実施形態の界面構造探索方法について説明する。
本実施形態の界面構造探索方法は、固相と、液相との界面構造を探索する界面構造探索方法であって、以下の工程を有することができる。
固相の表面に、液相が含有する分子が配置され、液相が含有する分子の配置が異なる複数の初期構造を決定する初期構造決定工程。
複数の初期構造について構造緩和計算を行い、それぞれ安定構造を探索する安定構造探索工程。
安定構造探索工程で求めた安定構造を表す記述子を作成する記述子作成工程。
記述子を用いて、界面構造のエネルギーを予測するモデルを作成するモデル作成工程。
モデル作成工程で用いたモデルを用いて、界面構造を決定する界面構造決定工程。
固相－液相の界面に限らず、構造候補が多数ある場合においては各構造のエネルギーを評価し、それらが分配関数に従って実現しているとする方法が一般的に用いられる。
低温を仮定することが妥当な状況下であれば、低エネルギー構造が高いウェイトを占めることになる。そのため、計算資源が限られている状況下においては、高いウェイトを占める低エネルギーの構造を多く計算対象とすることが望ましい。そのため、単純なランダム配置を用いた計算を行うより、ランダム配置した構造を初期構造として、構造緩和計算により安定構造を求めて、安定構造を採用して拡散などの計算を行い、物性値を評価することが好ましい。
しかしながら、液相側の分子が大きい場合や、分子種が複数ある場合、準安定構造が多数あるため、ある初期構造をもとにして得られた安定構造は最安定構造とは限らず、最安定構造から大きくかけ離れた構造である可能性も高い。そのため、多数の構造を初期構造として作成し、それら全てについて構造緩和計算により安定構造を探索する必要がある。そのため、安定構造を探索、選定する段階で非常に多くの計算コストがかかる。
そこで、本発明の発明者は、固相-液相の安定な界面構造を、計算量を抑制し、容易に推定できる界面構造探索方法について検討し、本発明を完成させた。本実施形態の界面構造探索方法について、工程毎に説明する。
（初期構造決定工程）
初期構造決定工程では、固相の表面に、液相が含有する分子が配置され、液相が含有する分子の配置が異なる複数の初期構造を決定できる。
（初期構造決定工程）
初期構造決定工程では、固相の表面に、液相が含有する分子が配置され、液相が含有する分子の配置が異なる複数の初期構造を決定できる。
本実施形態の界面構造探索方法では、後述するように記述子作成工程、モデル作成工程により、界面構造のエネルギーを予測するモデルである予測モデルを作成できる。係る予測モデルを用いることで、計算量を抑制して、固相－液相の安定な界面構造、すなわち低エネルギーとなる界面構造を容易に探索、推定できる。ただし、上述のように予測モデルを作成するためには、初期構造が１点では不足なため、複数の初期構造を作成できる。
初期構造決定工程において作成する初期構造の数は特に限定されず、予測モデルを作成するために必要となる数作成できる。初期構造を作成する際、固相の構造は、例えば固相の所定の結晶構造とすることができる。液相側については、例えば複数の分子種を含むように、液相に含まれる分子種を選択できる。また、液相に含まれる各分子種の物質量比に応じて、液相に含まれる分子数を決定できる。そして、液相側の分子の配置方法は特に限定されないが、例えばランダムに配置できる。
なお、初期構造決定工程を開始する前に、上述の固相の構造、具体的には液相と接する側の表面構造を決定できる。また、初期構造決定工程を開始する前に、液相が含有する分子種の種類を決定し、該液相が含有する分子数、すなわち初期構造内に配置する液相が含有する分子の個数を決定できる。上記固相の表面構造の決定や、液相が含有する分子種の種類や、液相が含有する分子数等は初期構造決定工程内で決定することもできる。
なお、例えばリチウムイオン二次電池の正極材料と電解液との界面構造を探索する場合、正極材料の具体的な組成は特に限定されず、リチウムイオンをドープおよび脱ドープすることができる、すなわちリチウムイオンを吸蔵・脱離できる、各種正極材料を用いることができる。正極材料としては、例えばスピネル型構造を有するリチウム金属複合酸化物や、層状構造を有するリチウム金属複合酸化物、オリビン型構造を有するリチウム金属複合酸化物等から選択された１種類以上を用いることができる。
また、液相の電解液としては、例えば支持塩としてのリチウム塩を有機溶媒に溶解したものを用いることができる。また、電解液として、イオン液体にリチウム塩が溶解したものを用いてもよい。なお、イオン液体とは、リチウムイオン以外のカチオンおよびアニオンから構成され、常温でも液体状の塩をいう。
有機溶媒としては、エチレンカーボネート、プロピレンカーボネート、ブチレンカーボネートおよびトリフルオロプロピレンカーボネートなどの環状カーボネートや、ジエチルカーボネート、ジメチルカーボネート、エチルメチルカーボネートおよびジプロピルカーボネートなどの鎖状カーボネート、さらにテトラヒドロフラン、２－メチルテトラヒドロフランおよびジメトキシエタンなどのエーテル化合物、エチルメチルスルホン、ブタンスルトンなどの硫黄化合物、リン酸トリエチル、リン酸トリオクチルなどのリン化合物等から選ばれる１種類を単独で用いてもよく、２種類以上を混合して用いることもできる。
支持塩としては、ＬｉＰＦ６、ＬｉＢＦ４、ＬｉＣｌＯ４、ＬｉＡｓＦ６、ＬｉＮ（ＣＦ３ＳＯ２）２、およびそれらの複合塩などを用いることができる。さらに、電解液は、ラジカル捕捉剤、界面活性剤および難燃剤などを含んでいてもよい。
ここでは、界面構造として、リチウムイオン二次電池の正極材料と電解液との界面構造の場合を例に説明したが、係る形態に限定されず、本実施形態の界面構造探索方法は各種界面構造に対して適用できる。
（安定構造探索工程）
安定構造探索工程では、複数の初期構造について構造緩和計算を行い、それぞれ安定構造を探索できる。
（安定構造探索工程）
安定構造探索工程では、複数の初期構造について構造緩和計算を行い、それぞれ安定構造を探索できる。
構造緩和計算の方法は特に限定されないが、密度汎関数理論（ＤＦＴ：Ｄｅｎｓｉｔｙ Ｆｕｎｃｔｉｏｎａｌ Ｔｈｅｏｒｙ）に基づく平面波基底第一原理計算を用いることが好ましい。
密度汎関数理論における汎関数としては、ＧＧＡ－ＰＢＥ（Ｇｅｎｅｒａｌｉｚｅｄ Ｇｒａｄｉｅｎｔ Ａｐｐｒｏｘｉｍａｔｉｏｎ－Ｐｅｒｄｅｗ，Ｂｕｒｋｅ，Ｅｒｎｚｅｒｈｏｆ）、またはより高精度な汎関数を用いることが好ましい。また、平面波のカットオフやｋ点のサンプリングはエネルギーが十分に収束するように、具体的には収束残差が１．０×１０－５ｅＶ以下になるように選択し、構造緩和は少なくとも原子に加わる力が０．０２ｅＶ／Å以下になるまで行うことが好ましい。そして、安定構造探索工程で求めた安定構造のエネルギーを求められる。
（記述子作成工程）
記述子作成工程では、安定構造探索工程で求めた安定構造を表す記述子を作成できる。具体的には、安定構造探索工程で得られた安定構造の位置情報を数値化する際に必要となる記述子を作成できる。記述子作成工程において作成する記述子としては、例えば、予め選択した原子と原子の結合距離や、予め選択した原子間の結合数、動径分布関数などを用いることができる。特に、本発明の発明者の検討によれば、安定構造には特定の原子間の結合が大きな影響を及ぼしているため、記述子作成工程において作成する記述子は、予め選択した原子間の結合数を含むことが好ましい。なお、予め選択した原子間の結合数とは、例えば固相、および液相に含まれる原子から選択した、原子間の結合数を意味する。本発明の発明者の検討によれば、例えばリチウムイオン二次電池の正極材料と電解液との界面構造においては、安定構造には、リチウムと、マイナスイオンとなる原子との間の結合数が大きな影響を有している。このため、リチウムイオン二次電池の正極材料と電解液との界面構造においては、記述子が、リチウムとマイナスイオンとなる原子との結合数を含むことが好ましい。ここでいうマイナスイオンとなる原子とは、酸素やハロゲンが挙げられ、特に酸素であることが好ましい。
（記述子作成工程）
記述子作成工程では、安定構造探索工程で求めた安定構造を表す記述子を作成できる。具体的には、安定構造探索工程で得られた安定構造の位置情報を数値化する際に必要となる記述子を作成できる。記述子作成工程において作成する記述子としては、例えば、予め選択した原子と原子の結合距離や、予め選択した原子間の結合数、動径分布関数などを用いることができる。特に、本発明の発明者の検討によれば、安定構造には特定の原子間の結合が大きな影響を及ぼしているため、記述子作成工程において作成する記述子は、予め選択した原子間の結合数を含むことが好ましい。なお、予め選択した原子間の結合数とは、例えば固相、および液相に含まれる原子から選択した、原子間の結合数を意味する。本発明の発明者の検討によれば、例えばリチウムイオン二次電池の正極材料と電解液との界面構造においては、安定構造には、リチウムと、マイナスイオンとなる原子との間の結合数が大きな影響を有している。このため、リチウムイオン二次電池の正極材料と電解液との界面構造においては、記述子が、リチウムとマイナスイオンとなる原子との結合数を含むことが好ましい。ここでいうマイナスイオンとなる原子とは、酸素やハロゲンが挙げられ、特に酸素であることが好ましい。
なお、予め選択した原子間の結合数を評価する場合、その結合形態は特に限定されず、例えば原子間距離が予め定めた所定の範囲内にある場合には結合しているとすることができる。係る距離は特に限定されないが、例えば２．６Å以下であることが好ましく、２．０Å以下であることがより好ましい。
（モデル作成工程）
モデル作成工程では、記述子作成工程で作成した記述子を用いて、界面構造のエネルギーを予測するモデルを作成できる。予測モデルには線形回帰や、決定木回帰、ランダムフォレスト、ガウス過程回帰等から選択された手法を用いることができる。
（モデル作成工程）
モデル作成工程では、記述子作成工程で作成した記述子を用いて、界面構造のエネルギーを予測するモデルを作成できる。予測モデルには線形回帰や、決定木回帰、ランダムフォレスト、ガウス過程回帰等から選択された手法を用いることができる。
ランダムフォレスト等の複雑な回帰式の場合、直感的に低エネルギーとなる界面構造を決めることが難しいため、例えば数値探索法により低エネルギーとなる界面構造を探索できる。ただし、モデル作成工程において作成するモデルとして、直感的に理解しやすい予測モデルを用いることが好ましく、例えばＬａｓｓｏ回帰や、Ｒｉｄｇｅ回帰などの線形回帰モデル等を用いることがより好ましい。
（界面構造決定工程）
そして、界面構造決定工程では、モデル作成工程で得られた予測モデルを用いて、エネルギーが低いことが予測される界面構造を作成し、構造緩和計算により界面構造を決定できる。
（界面構造決定工程）
そして、界面構造決定工程では、モデル作成工程で得られた予測モデルを用いて、エネルギーが低いことが予測される界面構造を作成し、構造緩和計算により界面構造を決定できる。
エネルギーが低いことが予測される界面構造は、記述子や予測モデルに線形回帰等の分かりやすいモデルを用いている場合には、予測モデルから手動で求めることもできる、また、複雑な記述子や予測モデルを用いている場合には、ベイズ最適化や粒子群最適化などの各種数値的な最適化手法により、低エネルギーが予測される界面構造を作成できる。
構造緩和計算については既に説明したため、ここでは説明を省略する。
なお、既述のモデル作成工程や、界面構造決定工程は、予測モデルを更新しながら複数回実施してもよい。
以下に具体的な実施例を挙げて説明するが、本発明はこれらの実施例に限定されるものではない。
［実施例１］
以下の手順により界面構造の探索を行った。
（初期構造決定工程）
初期構造決定工程では、固相である正極材料の表面に、液相である電解液の分子が配置され、分子の配置が異なる複数の初期構造を決定した。このため、初期構造は、固相として正極材料の表面構造と、液相として電解液の構造とを含む。
［実施例１］
以下の手順により界面構造の探索を行った。
（初期構造決定工程）
初期構造決定工程では、固相である正極材料の表面に、液相である電解液の分子が配置され、分子の配置が異なる複数の初期構造を決定した。このため、初期構造は、固相として正極材料の表面構造と、液相として電解液の構造とを含む。
なお、初期構造を作成する際、ユニットセル内に、正極材料としてはＬｉＮｉＯ２を、電解液としては、溶媒であるエチルメチルカーボネートを４分子と、エチレンカーボネート４分子と、電解質であるＬｉＰＦ６を２分子含む構造を作成した。なお、以下エチルメチルカーボネートはＥＭＣ、エチレンカーボネートはＥＣと記載する。
初期構造を作成する際、例えば図１に示すような構造を作成した。図１中正極材料１０と、電解液２０とが配置されている。図１中、同じハッチングは同じ種類の元素を意味しており、リチウム１１、酸素１２、ニッケル１３、リン１４、フッ素１５、炭素１６、水素１７の原子が配置されている。
初期構造は液相が含有する分子をランダムに配置して、全部で２４個作成した。
（安定構造探索工程）
安定構造探索工程では、初期構造作成工程で作成した２４個の初期構造について構造緩和計算を行い、それぞれに対応する２４個の安定構造を探索し、エネルギーを求めた。
（安定構造探索工程）
安定構造探索工程では、初期構造作成工程で作成した２４個の初期構造について構造緩和計算を行い、それぞれに対応する２４個の安定構造を探索し、エネルギーを求めた。
構造緩和計算、およびエネルギーの算出は、密度汎関数理論に基づく平面波基底第一原理計算を利用して行った。
第一原理計算は平面波基底第一原理計算ソフトであるＶＡＳＰ（Ｖｉｅｎｎａ Ａｂ ｉｎｉｔｉｏ Ｓｉｍｕｌａｔｉｏｎ Ｐａｃｋａｇｅ）を用いて行った。
また、第一原理計算は、密度汎関数理論（ＤＦＴ：Ｄｅｎｓｉｔｙ Ｆｕｎｃｔｉｏｎａｌ Ｔｈｅｏｒｙ）の範疇で行った。平面波基底のカットオフエネルギーは５００ｅＶとし、ｋ点を１×１×１とした。なお、計算は基底状態で行っている。
（記述子作成工程）
記述子作成工程では、記述子として、図２の横軸に示した様に、ＬｉＮｉＯ２のＬｉと以下の各酸素との結合数を用いた。具体的には、上記Ｌｉと、ＥＣの二重結合酸素（ｓｕｒｆＬｉ／ＥＣ＝Ｏ）、ＥＣの単結合酸素（ｓｕｒｆＬｉ／ＥＣ－Ｏ）、ＥＭＣの二重結合酸素（ｓｕｒｆＬｉ／ＥＭＣ＝Ｏ）、ＥＭＣの単結合酸素（ｓｕｒｆＬｉ／ＥＭＣ－Ｏ）の結合数を用いた。
（記述子作成工程）
記述子作成工程では、記述子として、図２の横軸に示した様に、ＬｉＮｉＯ２のＬｉと以下の各酸素との結合数を用いた。具体的には、上記Ｌｉと、ＥＣの二重結合酸素（ｓｕｒｆＬｉ／ＥＣ＝Ｏ）、ＥＣの単結合酸素（ｓｕｒｆＬｉ／ＥＣ－Ｏ）、ＥＭＣの二重結合酸素（ｓｕｒｆＬｉ／ＥＭＣ＝Ｏ）、ＥＭＣの単結合酸素（ｓｕｒｆＬｉ／ＥＭＣ－Ｏ）の結合数を用いた。
また、記述子として、ＬｉＰＦ６のＬｉと、以下の各酸素との結合数を用いた。具体的には、上記Ｌｉと、ＥＣの二重結合酸素（Ｌｉ／ＥＣ＝Ｏ）、ＥＣの単結合酸素（Ｌｉ／ＥＣ－Ｏ）、ＥＭＣの二重結合酸素（Ｌｉ／ＥＭＣ＝Ｏ）、ＥＭＣの単結合酸素（Ｌｉ／ＥＭＣ－Ｏ）)との結合数を用いた。
また、記述子としてさらにＬｉＰＦ６のＬｉと、ＬｉＮｉＯ２のＯの結合数（ｓｕｒｆＯ／Ｌｉ）、およびここまで説明した上記結合数の合計値（ｓｕｍ）を採用した。
図２は、各記述子による結合数当たりのエネルギーを示しており、縦軸はエネルギー／結合数であり、単位はｅＶになる。
（モデル作成工程）
モデル作成工程では、各安定構造が記述子作成工程で作成した１０個の記述子（パラメータ）で特徴づけられるとして、エネルギーを予測するモデルを作成した。予測モデルとしては線形回帰を用いた。
（構造決定工程）
モデル作成工程で作成した予測モデルの線形回帰の結果、図２からも明らかなようにＬｉＰＦ６のＬｉとＥＭＣの二重結合酸素の結合数が多い方が低エネルギーになることが明らかになった。このため、構造として可能な範囲で、ＬｉとＥＭＣの二重結合が多い構造を新たに初期構造として採用し、構造緩和計算により安定構造を求めたところ、ランダム探索を行った後述する比較例１の場合よりも大きく安定となる構造を得ることができた。すなわち、本実施例で行った界面構造探索方法によれば、固相－液相の安定な界面構造を容易に推定できることを確認できた。
［比較例１］
実施例１と同じ液相の構造のみが異なる初期構造をランダムに２４個作成した後、それぞれを構造緩和計算により２４個の安定構造を得た。比較例１ではさらに液相の構造のみが異なる２４個の構造をランダムに作成し、同様に構造緩和計算を行い２４個の安定構造を得た。なお、上記４８種類の初期構造は液相が含有する分子の配置のみが異なっており、含有する分子種や、その数は同じとなっている。
（モデル作成工程）
モデル作成工程では、各安定構造が記述子作成工程で作成した１０個の記述子（パラメータ）で特徴づけられるとして、エネルギーを予測するモデルを作成した。予測モデルとしては線形回帰を用いた。
（構造決定工程）
モデル作成工程で作成した予測モデルの線形回帰の結果、図２からも明らかなようにＬｉＰＦ６のＬｉとＥＭＣの二重結合酸素の結合数が多い方が低エネルギーになることが明らかになった。このため、構造として可能な範囲で、ＬｉとＥＭＣの二重結合が多い構造を新たに初期構造として採用し、構造緩和計算により安定構造を求めたところ、ランダム探索を行った後述する比較例１の場合よりも大きく安定となる構造を得ることができた。すなわち、本実施例で行った界面構造探索方法によれば、固相－液相の安定な界面構造を容易に推定できることを確認できた。
［比較例１］
実施例１と同じ液相の構造のみが異なる初期構造をランダムに２４個作成した後、それぞれを構造緩和計算により２４個の安定構造を得た。比較例１ではさらに液相の構造のみが異なる２４個の構造をランダムに作成し、同様に構造緩和計算を行い２４個の安定構造を得た。なお、上記４８種類の初期構造は液相が含有する分子の配置のみが異なっており、含有する分子種や、その数は同じとなっている。
以上により、比較例１では合計で４８個の安定構造を得た。結果を図３に示す。
図３は横軸が得られた安定構造のエネルギーを、縦軸がサンプル数を示している。図３中実施例１ランダムとして示しているのが、実施例１、および比較例１においてランダムに作成した２４個の初期構造から求めた安定構造のサンプルを示している。
図３中、比較例１ランダムとして示しているのが、比較例１においてさらにランダムに作成した２４個の初期構造から求めた安定構造のサンプルを示している。
しかしながら、比較例１で作成した合計４８個の安定構造のサンプルでは、図３中実施例１として示した、実施例１において作成したモデルを用いて予測、作成した安定構造よりも低いエネルギーの構造は得られなかった。
Claims (3)
- 固相と、液相との界面構造を探索する界面構造探索方法であって、
前記固相の表面に、前記液相が含有する分子が配置され、前記液相が含有する前記分子の配置が異なる複数の初期構造を決定する初期構造決定工程と、
複数の前記初期構造について構造緩和計算を行い、それぞれ安定構造を探索する安定構造探索工程と、
前記安定構造探索工程で求めた安定構造を表す記述子を作成する記述子作成工程と、
前記記述子を用いて、前記界面構造のエネルギーを予測するモデルを作成するモデル作成工程と、を有し、
前記モデルを用いて、前記界面構造を決定する界面構造決定工程と、を有する界面構造探索方法。 - 前記記述子作成工程において作成する前記記述子が、前記固相、および前記液相に含まれる原子から選択した、原子間の結合数を含む請求項１に記載の界面構造探索方法。
- 前記モデル作成工程において作成する前記モデルとして、線形回帰モデルを用いる請求項１または請求項２に記載の界面構造探索方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020166170A JP7420040B2 (ja) | 2020-09-30 | 2020-09-30 | 界面構造探索方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020166170A JP7420040B2 (ja) | 2020-09-30 | 2020-09-30 | 界面構造探索方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2022057759A JP2022057759A (ja) | 2022-04-11 |
| JP7420040B2 true JP7420040B2 (ja) | 2024-01-23 |
Family
ID=81111301
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020166170A Active JP7420040B2 (ja) | 2020-09-30 | 2020-09-30 | 界面構造探索方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7420040B2 (ja) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007072950A1 (ja) | 2005-12-22 | 2007-06-28 | Mitsui Mining & Smelting Co., Ltd. | 酸化亜鉛系透明導電膜のパターニング方法 |
| JP2009301437A (ja) | 2008-06-16 | 2009-12-24 | Toyota Central R&D Labs Inc | シミュレーションシステムおよびそのプログラム |
| JP2019210168A (ja) | 2018-05-31 | 2019-12-12 | 住友金属鉱山株式会社 | 置換元素の選択方法、電池正極材料の製造方法 |
-
2020
- 2020-09-30 JP JP2020166170A patent/JP7420040B2/ja active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007072950A1 (ja) | 2005-12-22 | 2007-06-28 | Mitsui Mining & Smelting Co., Ltd. | 酸化亜鉛系透明導電膜のパターニング方法 |
| JP2009301437A (ja) | 2008-06-16 | 2009-12-24 | Toyota Central R&D Labs Inc | シミュレーションシステムおよびそのプログラム |
| JP2019210168A (ja) | 2018-05-31 | 2019-12-12 | 住友金属鉱山株式会社 | 置換元素の選択方法、電池正極材料の製造方法 |
Non-Patent Citations (1)
| Title |
|---|
| 胡 春平 Chunping Hu，第一原理計算ソフトウェア Advance/PHASEの適用事例，アドバンスシミュレーション Ｖｏｌ．２３ Advance Simulation，日本，アドバンスソフト株式会社，2016年12月，第23巻，p.4-11 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2022057759A (ja) | 2022-04-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Review on modeling of the anode solid electrolyte interphase (SEI) for lithium-ion batteries | |
| Li et al. | Effects of imide–orthoborate dual-salt mixtures in organic carbonate electrolytes on the stability of lithium metal batteries | |
| Wu et al. | Recent progress in understanding solid electrolyte interphase on lithium metal anodes | |
| Li et al. | Low‐temperature electrolyte design for lithium‐ion batteries: prospect and challenges | |
| Freiberg et al. | Singlet oxygen reactivity with carbonate solvents used for Li-ion battery electrolytes | |
| Ma et al. | Critical factors dictating reversibility of the zinc metal anode | |
| Li et al. | Unraveling the new role of an ethylene carbonate solvation shell in rechargeable metal ion batteries | |
| Wu et al. | Development of cathode-electrolyte-interphase for safer lithium batteries | |
| Kim et al. | Stable electrode–electrolyte interfaces constructed by fluorine-and nitrogen-donating ionic additives for high-performance lithium metal batteries | |
| Guo et al. | Tailoring low-temperature performance of a lithium-ion battery via rational designing interphase on an anode | |
| Borodin et al. | Insights into the structure and transport of the lithium, sodium, magnesium, and zinc bis (trifluoromethansulfonyl) imide salts in ionic liquids | |
| Liu et al. | LiFSI and LiDFBOP dual-salt electrolyte reinforces the solid electrolyte interphase on a lithium metal anode | |
| Qian et al. | Mechanism study of unsaturated tripropargyl phosphate as an efficient electrolyte additive forming multifunctional interphases in lithium ion and lithium metal batteries | |
| Niehoff et al. | Interface investigations of a commercial lithium ion battery graphite anode material by sputter depth profile X-ray photoelectron spectroscopy | |
| Wang et al. | Reductive decomposition of solvents and additives toward solid-electrolyte interphase formation in lithium-ion battery | |
| Wu et al. | Designing anion-derived solid electrolyte interphase in a siloxane-based electrolyte for lithium-metal batteries | |
| Jorn et al. | Investigating the mechanism of lithium transport at solid electrolyte interphases | |
| Yang et al. | Dual-salt electrolytes to effectively reduce impedance rise of high-nickel lithium-ion batteries | |
| Han et al. | Theoretical study on fluoroethylene carbonate as an additive for the electrolyte of lithium ion batteries | |
| Hasanpoor et al. | Morphological evolution and solid–electrolyte interphase formation on LiNi0. 6Mn0. 2Co0. 2O2 cathodes using highly concentrated ionic liquid electrolytes | |
| Dhattarwal et al. | Mechanistic insight on the stability of ether and fluorinated ether solvent-based lithium bis (fluoromethanesulfonyl) electrolytes near Li metal surface | |
| Liu et al. | Distinct Differences in Li-Deposition/Dissolution Reversibility in Sulfolane-Based Electrolytes Depending on Li-Salt Species and Their Solvation Structures | |
| Li et al. | Superwettable Electrolyte Engineering for Fast Charging Li-Ion Batteries | |
| Lee et al. | Vulnerable solid electrolyte interphase deposition in sodium-ion batteries from insufficient overpotential development during formation | |
| JP7420040B2 (ja) | 界面構造探索方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230523 |
|
| 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: 20231212 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20231213 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20231225 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7420040 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |