(57)【要約】【課題】粗視化粒子の粒子間の相互作用を高精度で求めることができる高分子材料の解析方法を提供する。【解決手段】粗視化粒子の粒子間の相互作用を求める高分子材料の解析方法では、粒子間の相互作用を、χｉｊ＝０．２３８２×（εｉｊ−εｋｋ）＋１.５８８２、又は、（εｉｊ−εｋｋ）＝２．６８５８χｉｊを用いて求める。【選択図】図１
2013185822.tif 000002
【特許請求の範囲】
【請求項１】
粗視化粒子の粒子間の相互作用を求める高分子材料の解析方法であって、
前記粒子間の相互作用を、
χｉｊ＝０．２３８２×（εｉｊ−εｋｋ）＋１.５８８２
を用いて求めることを特徴とする高分子材料の解析方法。
【請求項２】
粗視化粒子の粒子間の相互作用を求める高分子材料の解析方法であって、
前記粒子間の相互作用を、
（εｉｊ−εｋｋ）＝２．６８５８χｉｊ
を用いて求めることを特徴とする高分子材料の解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、粗視化粒子の粒子間の相互作用を求める高分子材料の解析方法に関する。
【背景技術】
【０００２】
電池では、高分子中のイオンの輸送現象が電池性能を決定する要因となる。例えば、燃料電池における高分子膜中の水素イオンの輸送や、触媒層高分子電解質中の水素イオンの輸送、リチウムポリマー電池における高分子膜中のリチウムイオンの輸送が上記の現象に該当し、いずれも数ｎｍ〜数μｍのスケールで生じる。このよう数ｎｍ〜数μｍの現象について、それを直接計測して解析することは困難である。
【０００３】
そこで、スケールの非常に小さい現象の解析には、複数の原子を一つの粒子として扱う粗視化という方法が用いられている。この粗視化による解析方法では、粗視化粒子の粒子間の相互作用を求めることが必要とされる。この相互作用を求める方法としては、例えば２成分からなる混合溶液を生成し、その生成時の溶解性から相互作用を求める方法が知られている（例えば、下記非特許文献１参照）。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Gary S.Grest at al.、「Efficient continuum model for simulating polymer blends andcopolymers」、J. Chem. Phys., Vol. 105, No.23、15 December 1996、pp10583-10594
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上記非特許文献１に記載の方法では、２成分からなる混合溶液を生成し、そのときの溶解性に基づいて粒子間の相互作用を求めているため、混ざり易い成分同士であれば短時間で相互作用を求めることができるが、混ざり難い成分同士についてはどうしても誤差が大きくなるといった問題がある。
【０００６】
本発明は、上記課題を解決するためになされたものであり、粗視化粒子の粒子間の相互作用を高精度で求めることができる高分子材料の解析方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、本発明に係る高分子材料の解析方法は、粗視化粒子の粒子間の相互作用を求める高分子材料の解析方法であって、粒子間の相互作用を、
χｉｊ＝０．２３８２×（εｉｊ−εｋｋ）＋１.５８８２
を用いて求めることを特徴とする。
【０００８】
この高分子材料の解析方法では、
χｉｊ＝０．２３８２×（εｉｊ−εｋｋ）＋１.５８８２
を用いて、粗視化粒子の粒子間の相互作用を求めている。これにより、高分子などの混ざり難いχパラメータを有するものであっても、求められる相互作用の誤差を小さくすることができる。したがって、粗視化粒子の粒子間の相互作用を高精度で求めることができる。
【０００９】
本発明に係る高分子材料の解析方法は、粗視化粒子の粒子間の相互作用を求める高分子材料の解析方法であって、粒子間の相互作用を、
（εｉｊ−εｋｋ）＝２．６８５８χｉｊ
を用いて求めることを特徴とする。
【００１０】
この高分子材料の解析方法では、
（εｉｊ−εｋｋ）＝２．６８５８χｉｊ
を用いて、粗視化粒子の粒子間の相互作用を求めている。これにより、高分子などの混ざり難いχパラメータを有するものであっても、求められる相互作用の誤差を小さくすることができる。したがって、粗視化粒子の粒子間の相互作用を高精度で求めることができる。
【発明の効果】
【００１１】
本発明によれば、粗視化粒子の粒子間の相互作用を高精度で求めることができる。
【図面の簡単な説明】
【００１２】
【図１】高分子材料の解析方法を示すフローチャートである。
【図２】ポリクロロトリフルオロエチレンを粗視化した粒子の一例を示す図である。
【図３】χパラメータと相互作用パラメータεとの関係を示すグラフである。
【発明を実施するための形態】
【００１３】
以下、添付図面を参照して、本発明の好適な実施形態について詳細に説明する。なお、図面の説明において同一又は相当要素には同一符号を付し、重複する説明は省略する。
【００１４】
［第１実施形態］
図１は、第１実施形態に係る高分子材料の解析方法を示すフローチャートである。以下、図１を参照しながら、高分子材料の解析方法について説明する。本実施形態の高分子材料の解析方法では、高分子の原子を複数含んだ粗視化粒子の相互作用を、フローリー・ハギンズ（Flory-Huggins）理論においてパラメータが２以上である領域で定めたχパラメータと相互作用パラメータεとの間の単調関数（１価関数）により求める。
【００１５】
図１に示すように、まず高分子を複数の原子を含むユニットに分け、ユニット毎に粗視化粒子を作成する（ステップＳ０１）。図２は、ポリクロロトリフルオロエチレンを粗視化した粒子の一例を示す図である。図２に示す例では、破線で囲んだ部分を１つのユニットＵとしている（図２では、３つのユニットＵを示している）。
【００１６】
図１に戻って、次に、χパラメータを算出する（ステップＳ０２）。χパラメータは、各ユニットＵの溶解度パラメータδから算出する。溶解度パラメータδには、実験値、推算値及び計算値などを用いることができる。なお、実験値は高分子ハンドブック、推算値は原子団寄与法、計算値は分子動力学法などにより設定できる。χパラメータは、溶解度パラメータδの差の二乗に比例するように設定する。具体的には、χパラメータは、以下の式（１）から求める。
【数１】
2013185822.tif 000003
上記式（１）において、δ：溶解度パラメータ、Ｖ：モル体積、ｋＢ：ボルツマン係数、Ｔ：絶対温度である。χパラメータは、式（１）に上記の溶解度パラメータδの値を代入して求める。
【００１７】
続いて、粗視化粒子の粒子間の相互作用を導出する（ステップＳ０３）。粗視化粒子の粒子間の相互作用は、χパラメータと相互作用パラメータεとの間の単調関数に基づいて求める。図３は、χパラメータと相互作用パラメータεとの関係を示すグラフである。図３では、実線はχパラメータが２以上の領域で定めた場合を示し、破線はχパラメータが２以下の領域で定めた場合を示している。本実施形態では、フローリー・ハギンズ理論においてパラメータが２以上の領域で定めたχパラメータを用いる。
【００１８】
レナード-ジョーンズ・ポテンシャル（Lennard-Jonespotential）の相互作用パラメータεには、下記式（２）を用いる。
χｉｊ＝０．２３８２×（εｉｊ−εｋｋ）＋１.５８８２  …（２）
上記式（２）において、「ｉ」及び「ｊ」は、それぞれ原子の種類を示している。このようにして、粗視化粒子の粒子間の相互作用を導出する
【００１９】
以上説明したように、本実施形態では、
χｉｊ＝０．２３８２×（εｉｊ−εｋｋ）＋１.５８８２
を用いて、粗視化粒子の粒子間の相互作用を求めている。これにより、高分子などの混ざり難いχパラメータを有するものであっても、求められる相互作用の誤差を小さくすることができる。すなわち、図３に示すχパラメータが２以下の領域で定めた関数を用いることにより、混合溶液が混ざり易いχパラメータが２以下で定められる関数により相互作用パラメータεを設定した場合に比べて、電池などで多く使用される混合溶液が混ざり難い場合において誤差を小さくできる。
【００２０】
詳細には、例えば、燃料電池の高分子膜における水素イオンの輸送経路である水形状について、従来の解析方法（非特許文献１に記載の方法）と本実施形態の解析方法とで比較した結果、従来の解析方法では実験誤差が７１％であるのに対し、本実施形態では実験誤差が１１％と大幅に誤差が低減した。
【００２１】
［第２実施形態］
続いて、第２実施形態について説明する。第２実施形態に係る高分子材料の解析方法は、相互作用の導出方法が第１実施形態と異なっている。本実施形態の高分子材料の解析方法では、２成分系の自由エネルギー密度の式と相互作用の式とから設定したχパラメータと相互作用パラメータεの間の単調関数により求める。
【００２２】
本実施形態の高分子材料の解析方法では、第１実施形態と同様に、図１に示すように、まず高分子を複数の原子を含むユニットＵに分け、ユニットＵ毎に粗視化粒子を作成する（ステップＳ０１）。
【００２３】
続いて、χパラメータを算出し（ステップＳ０２）、その後に、粗視化粒子の粒子間の相互作用を導出する（ステップＳ０３）。粗視化粒子の粒子間の相互作用は、自由エネルギーの式と相互作用の式とから求める。
【００２４】
自由エネルギー密度の式は、下記式（３）により表される。
【数３】
2013185822.tif 000004
上記式（３）において、ρ：密度である。
【００２５】
また、２成分系であることから、ａＡＡ＝ａＢＢ，ρＡ＋ρＢ≒ｃｏｎｓｔ.と条件を設定すると、χパラメータは、以下の式（４）で表される。
【数４】
2013185822.tif 000005
上記式（４）において、σ＝ｃｏｎｓｔ．の条件とし、式（４）を以下の式（５）に変換する。
【数５】
2013185822.tif 000006
上記式（５）において、σ：排除体積である。そして、「Ａ」は以下の式（６）で表される。
【数６】
2013185822.tif 000007
上記式（６）において、ｒｃｕｔ：カットオフ（原子間）距離である。また、上記式（５）において、Δε＝εｉｊ−εｋｋである。
【００２６】
上記のように、χパラメータと相互作用パラメータεとの間の単調関数（１価関数）から相互作用パラメータεを設定するときの単調関数を、２成分系の自由エネルギー密度の式（３）と相互作用の式（５）とから求める。レナード-ジョーンズ・ポテンシャル（Lennard-Jones potential）の相互作用パラメータεには、下記式（７）を用いる。
（εｉｊ−εｋｋ）＝２．６８５８χｉｊ  …（７）
以上により、粗視化粒子の粒子間の相互作用を導出する。
【００２７】
以上説明したように、本実施形態では、
（εｉｊ−εｋｋ）＝２．６８５８χｉｊ
により、粒子間の相互作用を求めている。これにより、高分子などの混ざり難いχパラメータをもつものであっても、相互作用の誤差を小さくすることができる。したがって、粗視化粒子の粒子間の相互作用を高精度で求めることができる。すなわち、混合溶液においてχパラメータが混ざり難い場合に、混ざり難いといった条件下で得られる理論式から設定した関数を用いることにより、高分子材料を粗視化した粒子間の相互作用を精度よく求めることができる。
【００２８】
詳細には、例えば、燃料電池の高分子膜における水素イオンの輸送経路である水形状について、従来の解析方法（非特許文献１に記載の方法）と本実施形態の解析方法とで比較した結果、従来の解析方法では実験誤差が７４％であるのに対し、本実施形態では実験誤差が２９％と大幅に誤差が低減した。
【００２９】
本発明は、上記実施形態に限定されるものではない。例えば、上記実施形態に加えて、粗視化した粒子間の相互作用を求めるときに、上記の式（２）及び式（７）の両方を用いてもよい。
【図１】
2013185822.tif 000008
【図２】
2013185822.tif 000009
【図３】
2013185822.tif 000010