(57)【要約】【課題】不均一系の輸送係数を解析することが可能な輸送係数の解析プログラムおよび輸送係数の解析方法を提供する。【解決手段】輸送係数の解析プログラムは、解析対象の系を複数の領域に分割する分割部Ｓ２と、系に含まれる原子の運動を計算する運動計算部Ｓ３と、運動計算部Ｓ３の計算結果に基づいて、複数の領域それぞれについて輸送係数に対応する物理量を算出する物理量算出部Ｓ４と、物理量の時間相関関数を算出する時間相関関数算出部Ｓ５と、時間相関関数に基づいて、複数の領域それぞれについて輸送係数を算出する輸送係数算出部Ｓ６と、を備える。【選択図】図３
2019148446.tif 000002
【特許請求の範囲】
【請求項１】
輸送係数の解析プログラムであって、
解析対象の系を複数の領域に分割する分割部（Ｓ２）と、
前記系に含まれる原子の運動を計算する運動計算部（Ｓ３）と、
前記運動計算部の計算結果に基づいて、前記複数の領域それぞれについて前記輸送係数に対応する物理量を算出する物理量算出部（Ｓ４）と、
前記物理量の時間相関関数を算出する時間相関関数算出部（Ｓ５）と、
前記時間相関関数に基づいて、前記複数の領域それぞれについて前記輸送係数を算出する輸送係数算出部（Ｓ６）と、を備える輸送係数の解析プログラム。
【請求項２】
前記輸送係数は、緩和弾性率である請求項１に記載の輸送係数の解析プログラム。
【請求項３】
前記時間相関関数算出部は、並列処理を用いて前記時間相関関数を算出する請求項１または２に記載の輸送係数の解析プログラム。
【請求項４】
輸送係数の解析方法であって、
解析対象の系を複数の領域に分割すること（Ｓ２）と、
前記系に含まれる原子の運動を計算すること（Ｓ３）と、
前記原子の運動を計算することの結果に基づいて、前記複数の領域それぞれについて前記輸送係数に対応する物理量を算出すること（Ｓ４）と、
前記物理量の時間相関関数を算出すること（Ｓ５）と、
前記時間相関関数に基づいて、前記複数の領域それぞれについて前記輸送係数を算出すること（Ｓ６）と、を備える輸送係数の解析方法。
【請求項５】
前記輸送係数は、緩和弾性率である請求項４に記載の輸送係数の解析方法。
【請求項６】
前記時間相関関数を算出することでは、並列処理を用いて前記時間相関関数を算出する請求項４または５に記載の輸送係数の解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、輸送係数の解析プログラムおよび輸送係数の解析方法に関するものである。
【背景技術】
【０００２】
従来、緩和弾性率、粘度、電気伝導度、熱伝導度、拡散係数等の輸送係数の解析については、測定による方法の他、シミュレーションによる方法が提案されている。例えば非特許文献１では、イオン間の距離に基づいて電気伝導度を解析する方法が提案されている。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】K.-M. Tu, R. Ishizuka, N. Matubayasi, "Spatial-decomposition analysis of electrical conductivity in concentrated electrolyte solution", The Journal of Chemical Physics 141, 044126 (2014)
【発明の概要】
【発明が解決しようとする課題】
【０００４】
接着面等の界面、特にナノオーダーの凹凸が形成された界面の特性が強調されるモデルでは、輸送係数の空間不均一性が系全体の特性に大きく影響すると予想され、このようなモデルの接着、電気伝導、熱伝導等の現象の解析は、産業上重要である。また、接着面、フィラー周りの電気特性、熱特性などの解析は、不具合の解明や新たな接着処理、フィラー処理等の実現、新規材料の創出等に役立つと考えられる。
【０００５】
しかしながら、非特許文献１に記載の方法では、このような不均一系の輸送係数を解析することができない。
【０００６】
本発明は上記点に鑑みて、不均一系の輸送係数を解析することが可能な輸送係数の解析プログラムおよび輸送係数の解析方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記目的を達成するため、請求項１に記載の発明では、輸送係数の解析プログラムであって、解析対象の系を複数の領域に分割する分割部（Ｓ２）と、系に含まれる原子の運動を計算する運動計算部（Ｓ３）と、運動計算部の計算結果に基づいて、複数の領域それぞれについて輸送係数に対応する物理量を算出する物理量算出部（Ｓ４）と、物理量の時間相関関数を算出する時間相関関数算出部（Ｓ５）と、時間相関関数に基づいて、複数の領域それぞれについて輸送係数を算出する輸送係数算出部（Ｓ６）と、を備える。
【０００８】
このように、系を複数の領域に分割し、各領域について輸送係数に対応する物理量を算出することにより、不均一系の輸送係数を解析することができる。
【０００９】
また、請求項４に記載の発明では、輸送係数の解析方法であって、解析対象の系を複数の領域に分割すること（Ｓ２）と、系に含まれる原子の運動を計算すること（Ｓ３）と、原子の運動を計算することの結果に基づいて、複数の領域それぞれについて輸送係数に対応する物理量を算出すること（Ｓ４）と、物理量の時間相関関数を算出すること（Ｓ５）と、時間相関関数に基づいて、複数の領域それぞれについて輸送係数を算出すること（Ｓ６）と、を備える。
【００１０】
このように、系を複数の領域に分割し、各領域について輸送係数に対応する物理量を算出することにより、不均一系の輸送係数を解析することができる。
【００１１】
なお、各構成要素等に付された括弧付きの参照符号は、その構成要素等と後述する実施形態に記載の具体的な構成要素等との対応関係の一例を示すものである。
【図面の簡単な説明】
【００１２】
【図１】第１実施形態が適用される不均一系の概略構成を示す図である。
【図２】図１のII部分の拡大図である。
【図３】輸送係数の解析処理のフローチャートである。
【図４】輸送係数の解析結果を示す図である。
【発明を実施するための形態】
【００１３】
以下、本発明の実施形態について図に基づいて説明する。なお、以下の各実施形態相互において、互いに同一もしくは均等である部分には、同一符号を付して説明を行う。
【００１４】
（第１実施形態）
第１実施形態について説明する。ここでは、図１に示す系について、輸送係数の１つである緩和弾性率を算出する場合について説明する。
【００１５】
図１に示す系は、金属層１の表面に接着剤２が塗布され、接着剤２に対して金属層１とは反対側の領域が真空とされたものである。金属層１は、ここではアルミニウムで構成されており、接着剤２は、ここでは炭素、酸素、水素で構成された樹脂とされている。
【００１６】
図２に示すように、金属層１と接着剤２との界面付近は空間的に不均一な構成とされており、このような構成の系については、従来の方法では輸送係数の解析が困難である。なお、接着剤２のうち金属層１との界面から離れた部分は均一な構成とされているが、この部分の輸送係数にも、金属層１との界面付近の特性が影響すると考えられる。
【００１７】
本実施形態では、図３に示すステップＳ１〜Ｓ６を順に実行することで、図１に示す系の緩和弾性率を算出し、界面付近の特性、および、界面から離れた部分の特性の解析を可能とする。なお、図３に示す処理は、例えばＬＡＭＭＰＳ、Ｐｙｔｈｏｎ等を用いた解析プログラムによって実行することができるが、他のソフトウェアおよびプログラミング言語を用いて解析を行ってもよい。
【００１８】
ステップＳ１では、解析対象の系のデータを取得する。このデータには、図２に示すような金属層１、接着剤２を構成する複数の原子の種類および初期位置の他、原子の初速度等が含まれている。なお、ステップＳ１では、このようなデータを解析プログラムの外部から取得してもよいし、解析プログラム内で設定してもよい。
【００１９】
ステップＳ２では、系を複数の領域に分割する。例えば、図１の紙面左右方向、奥行き方向、上下方向をｘ方向、ｙ方向、ｚ方向として、系をｘ方向およびｙ方向に２０分割し、ｚ方向に１００分割する。
【００２０】
ステップＳ３では、ＭＤ（分子動力学）法等に基づいて系に含まれる複数の原子の運動を計算し、計算の開始から終了までの複数の時刻における各原子の位置および速度を取得する。
【００２１】
ステップＳ４では、各時刻における各原子の位置および速度に基づいて、各時刻における各領域の輸送係数に対応する物理量を算出する。本実施形態では、各領域の応力テンソルを算出し、緩和弾性率に対応する物理量としてせん断応力を用いる。
【００２２】
ステップＳ５では、せん断応力のＴＣＦ（時間相関関数）を算出する。そして、ステップＳ６では、ステップＳ５で求めたＴＣＦに基づいて各領域の緩和弾性率を算出する。ステップＳ５およびステップＳ６では、具体的には、次のようにして緩和弾性率を求める。
【００２３】
すなわち、複数の領域の数をＮｓｕｂとし、時間をｔとし、複数の領域のうちｊ番目の領域の時間ｔにおける緩和弾性率をＧｊ（ｔ）とする。また、ｊ番目の領域の時間ｔにおける応力テンソルのａｂ成分をσｊａｂ（ｔ）とし、ｋ番目の領域の時間ｔにおける応力テンソルのａｂ成分をσｋａｂ（ｔ）とし、数式１を用いてＧｊ（ｔ）を求める。応力テンソルのａｂ成分は、例えば応力テンソルのｘｙ成分、ｙｚ成分、ｚｘ成分等である。なお、ＴＣＦの算出にはmulti−tau相関法などのスムージング手法を用いてもよい。
【００２４】
【数１】
2019148446.tif 000003
【００２５】
ステップＳ５では、ステップＳ４で算出された各時刻の応力テンソルについて順次計算を行ってもよいが、処理効率の向上のために、各時刻の応力テンソルについて並列処理を行うことが望ましい。
【００２６】
例えば、ステップＳ４において、応力テンソルを１ｆｓ毎に算出し、算出結果のデータを１００ｐｓ毎に出力する。そして、ステップＳ５にて、１００ｐｓ毎に分割された複数のデータを取得し、複数のデータそれぞれについて並列にＴＣＦを算出する。なお、ＴＣＦを算出する際には、遅れ時間に応じて、計算対象としているデータに、このデータに続く次の１００ｐｓのデータを追加して用いる。そして、１００ｐｓ毎のＴＣＦの計算結果を平均化する。
【００２７】
なお、Ｇｊ（ｔ）を用いて、数式２により、時間ｔにおける系全体の緩和弾性率Ｇ（ｔ）を算出することができる。ここで、Ｖは系全体の体積、ｋはボルツマン定数である。数式１、数式２の導出方法については後述する。
【００２８】
【数２】
2019148446.tif 000004
【００２９】
このようにして、系の輸送係数を求めることができる。なお、ステップＳ２は分割部に対応し、ステップＳ３は運動計算部に対応し、ステップＳ４は物理量算出部に対応し、ステップＳ５は時間相関関数算出部に対応し、ステップＳ６は輸送係数算出部に対応する。
【００３０】
図３に示す処理によって、例えば図４に示すような結果が得られる。なお、図４において、黒丸は、金属層１と接着剤２との界面付近の輸送係数の算出結果であり、実線は、この算出結果の近似曲線である。また、黒四角は、接着剤２のうち金属層１との界面から離れた部分の輸送係数の算出結果であり、破線は、この算出結果の近似曲線である。
【００３１】
図４に示す結果では、金属層１と接着剤２との界面付近が、接着剤２のうち金属層１との界面から離れた部分に比べて、Ｇｊ（ｔ）の減衰が遅く、弾性的になっている。金属層１の表面処理の方法や、接着剤２を構成する樹脂の種類等を様々に変化させて図３に示す処理を実行し、得られた結果を比較することで、用途に適した緩和弾性率が得られる材料等の条件を導き出すことができる。
【００３２】
数式１、数式２の導出方法について説明する。均一な系の緩和弾性率は、数式３で表され、せん断応力のＴＣＦに対応することが知られている。
【００３３】
【数３】
2019148446.tif 000005
【００３４】
また、原子にポテンシャルＵが働く場合、マクロな応力は、各原子からの寄与の和として、数式４のように表されることが知られている。ｖｉ１は、系に含まれる複数の原子のうちｉ１番目の原子の速度であり、ｒｉ１ｉ２はｉ１番目の原子とｉ２番目の原子との間の距離である。
【００３５】
【数４】
2019148446.tif 000006
【００３６】
ここで、系全体の応力のうちｉ１番目の原子からの寄与は、local atomic stress elementsと呼ばれる以下の量ｓｉ１ａｂで表される。
【００３７】
【数５】
2019148446.tif 000007
【００３８】
ｓｉ１ａｂを用いると、系全体の応力は数式６で表される。なお、ｓｉ１ａｂとしてはＬＡＭＭＰＳ等において具体的に定義されたものを用いればよい。
【００３９】
【数６】
2019148446.tif 000008
【００４０】
系を複数の領域に分割し、各領域についてｓｉ１ａｂの和を求めるようにすると、数式６は数式７のように変形される。Ｎａはｊ番目の領域に含まれる原子の個数、ｓｉ１，ｊａｂはｊ番目の領域に含まれる複数の原子のうちｉ１番目の原子のlocal atomic stress elementsである。
【００４１】
【数７】
2019148446.tif 000009
【００４２】
ここで、ｊ番目の領域の応力をσｊａｂとすると、σｊａｂは、数式８のようにｓｉ１、ｊａｂの和で表されるので、系全体の応力σａｂは、σｊａｂを用いて数式９のように表される。
【００４３】
【数８】
2019148446.tif 000010
【００４４】
【数９】
2019148446.tif 000011
【００４５】
数式９を数式３に代入して変形すると、緩和弾性率Ｇ（ｔ）は、数式１０のように表される。
【００４６】
【数１０】
2019148446.tif 000012
【００４７】
このように、系全体の緩和弾性率Ｇ（ｔ）は、領域ごとの相関関数の和で表される。そこで、数式１０のうち相関関数で表された部分を各領域の緩和弾性率Ｇｊ（ｔ）とし、数式１のように定義すれば、系全体の緩和弾性率Ｇ（ｔ）は、数式２に示すようになる。
【００４８】
以上説明したように、本実施形態では、系を複数の領域に分割し、各領域について輸送係数に対応する物理量を算出することで、各領域の輸送係数を算出することができる。これにより、不均一系について、界面付近の特性、および、界面から離れた部分の特性を解析することが可能となる。
【００４９】
（他の実施形態）
なお、本発明は上記した実施形態に限定されるものではなく、特許請求の範囲に記載した範囲内において適宜変更が可能である。
【００５０】
例えば、上記第１実施形態では緩和弾性率の解析方法について説明したが、粘度、電気伝導度、熱伝導度、拡散係数等の他の輸送係数についても、第１実施形態と同様に、ステップＳ１〜Ｓ６の処理により解析することができる。
【００５１】
すなわち、系全体の物理量を数式６のように各原子からの寄与の和で表すことができれば、寄与の和の部分を数式７のように領域ごとの和に分割して、系全体の物理量を数式９のように各領域の物理量の和を用いて表すことができる。また、輸送係数は、数式３と同様に、対応する物理量のＴＣＦを用いて表されるので、数式９に相当する式を、輸送係数を表す式のＴＣＦの部分に代入して変形すれば、各領域の物理量のＴＣＦの和を用いて輸送係数を表すことができる。そして、各領域についてのＴＣＦの部分を各領域の輸送係数とすれば、第１実施形態と同様に不均一系の輸送係数を解析することができる。
【符号の説明】
【００５２】
Ｓ２  分割部
Ｓ３  運動計算部
Ｓ４  物理量算出部
Ｓ５  時間相関関数算出部
Ｓ６  輸送係数算出部
【図１】
2019148446.tif 000013
【図２】
2019148446.tif 000014
【図３】
2019148446.tif 000015
【図４】
2019148446.tif 000016