~マルチスケール模型での計算 現場数学(その6)

更新日

投稿日

 

1.  計算の限界

 如何(いか)にスーパーコンピューターを活用しても、計算出来る範囲はやはり限られています。例えば、タンパク質の3次元構造を電子状態から計算しようとすれば、百万個どころではない電子数を扱わなければならなくなります。現在の標準第一原理計算である密度汎関数法のレベルでは1万電子程度が限度であり、計算量は電子数の3乗に比例するので、現状より百万倍速い計算機が必要ということになります。

 機械工学の基礎であるクラックの問題でも同様の規模の計算量が必要とされています。一方、スーパーコンピューターは、アムダールの法則が継続すると仮定して5年で1桁(けた)速くなるので、これらの問題を電子レベルで計算することが30年後には可能になると予想されます。そうなれば、実験事実を使わずに、物理の第一原理だけでタンパク質や機械用材料の原子構造、物性値を予言することができるため、新物質・新材料・新生体材料の精密で迅速な理論設計が実現するはず…なのでしょうか?

2.  マルチスケール模型の登場

 30年待っている?まさか、そうはいかないので、現在使える計算機を使って何とか必要な解を求めなければなりません。そのためには、精度を落とした計算方法を採用することが考えられます。例えば、量子式学の方程式を解くのではありますが、相互作用をパラメーター化したタイトバインディング模型を使う、古典分子動力学を使う、有限要素法を適用した連続体模型を使う…などの方法があります。しかし、電子状態は原子位置が変われば大幅に変わるという具合に極めて微妙です。それに実際、反応などに重要な寄与をする原子の数は限られているのです。そこで、肝心な部分は精密に計算し、その周りは少し雑に計算、さらにもっと遠方はもっと雑に扱う、という方策が考案されました。

 下の図のように、メッシュサイズを大幅に変え、中心近くのメッシュ内では第一原理計算、その次のメッシュ内では古典分子動力学、さらに外側では連続体模型を適用するのです。もちろん、その境界部分での接続を如何に上手に扱うかが一番の問題になります。もちろん、図の様な異なるサイズのメッシュを使う方法は昔から用いられ、古典分子動力学法では標準的な手法です。しかし、そこでは、メッシュサイズごとに全く異なる計算模型を適用していた訳ではありません。

図. メッシュを使った古典分子動力学の標準的手法

 鉄道模型の世界でも、昔からマルチスケールという言葉が使われてきました。まず、単位であるゲージに9mm、16.5mm、32mm、45mmとあります(それぞれ、Nゲージ、16番、Oゲージ、Nゲージと呼ばれます)。さらに、16.5mmの線路幅には、HO、16番、O O、Sn3-1/2、On30とあり、それぞれ、1/87、1/80、1/76、1/64、1/48模型になります。それは本物の鉄道の方に広軌と狭軌などと異なる規格があるからなのです。さらに、遊園地で見るような大きい模型もあり、それらは、16.5mmの軌道の何倍かになっていることが多いのです。本当は単一ゲージでその何倍かだけがあるのが良かったのでしょうけども、そうなっていないのは、国の事情、人間の歴史の複雑さの表れです。

3.  今後の発展 ~ マルチスケール模型は現場数学に不可欠

 最初に述べたように、計算機は5年で一桁速くなっていますので、そのうち、現実に存在する物質をそのまま計算機の中に設定して、それを電子状態から全て計算できるようになると考える人が多いと思います。しかし、それはほぼ不可能なのです。なぜならば、上述の話は、我々が対象とする物質は一様であると想定していました。しかし、本当の物質はそんなに単純ではありません。結晶も実は多結晶と呼ばれていて、単位となる結晶粒から構成されているのです。その各粒は1020くらいの原子から構成されているのが普通なのです。つまり、単位として解かなければならない結晶粒の精密計算が可能となるだけでも今後、80年もかかることになるのです。さらに、引っ張り強度などは、主として粒界(結晶粒間)の相互作用によるので…と…どんどん、気が遠くなるように必要な計算量は増えます。タンパク質の計算でも、現在は単一のタンパク質の3次元構造を問題...

 

1.  計算の限界

 如何(いか)にスーパーコンピューターを活用しても、計算出来る範囲はやはり限られています。例えば、タンパク質の3次元構造を電子状態から計算しようとすれば、百万個どころではない電子数を扱わなければならなくなります。現在の標準第一原理計算である密度汎関数法のレベルでは1万電子程度が限度であり、計算量は電子数の3乗に比例するので、現状より百万倍速い計算機が必要ということになります。

 機械工学の基礎であるクラックの問題でも同様の規模の計算量が必要とされています。一方、スーパーコンピューターは、アムダールの法則が継続すると仮定して5年で1桁(けた)速くなるので、これらの問題を電子レベルで計算することが30年後には可能になると予想されます。そうなれば、実験事実を使わずに、物理の第一原理だけでタンパク質や機械用材料の原子構造、物性値を予言することができるため、新物質・新材料・新生体材料の精密で迅速な理論設計が実現するはず…なのでしょうか?

2.  マルチスケール模型の登場

 30年待っている?まさか、そうはいかないので、現在使える計算機を使って何とか必要な解を求めなければなりません。そのためには、精度を落とした計算方法を採用することが考えられます。例えば、量子式学の方程式を解くのではありますが、相互作用をパラメーター化したタイトバインディング模型を使う、古典分子動力学を使う、有限要素法を適用した連続体模型を使う…などの方法があります。しかし、電子状態は原子位置が変われば大幅に変わるという具合に極めて微妙です。それに実際、反応などに重要な寄与をする原子の数は限られているのです。そこで、肝心な部分は精密に計算し、その周りは少し雑に計算、さらにもっと遠方はもっと雑に扱う、という方策が考案されました。

 下の図のように、メッシュサイズを大幅に変え、中心近くのメッシュ内では第一原理計算、その次のメッシュ内では古典分子動力学、さらに外側では連続体模型を適用するのです。もちろん、その境界部分での接続を如何に上手に扱うかが一番の問題になります。もちろん、図の様な異なるサイズのメッシュを使う方法は昔から用いられ、古典分子動力学法では標準的な手法です。しかし、そこでは、メッシュサイズごとに全く異なる計算模型を適用していた訳ではありません。

図. メッシュを使った古典分子動力学の標準的手法

 鉄道模型の世界でも、昔からマルチスケールという言葉が使われてきました。まず、単位であるゲージに9mm、16.5mm、32mm、45mmとあります(それぞれ、Nゲージ、16番、Oゲージ、Nゲージと呼ばれます)。さらに、16.5mmの線路幅には、HO、16番、O O、Sn3-1/2、On30とあり、それぞれ、1/87、1/80、1/76、1/64、1/48模型になります。それは本物の鉄道の方に広軌と狭軌などと異なる規格があるからなのです。さらに、遊園地で見るような大きい模型もあり、それらは、16.5mmの軌道の何倍かになっていることが多いのです。本当は単一ゲージでその何倍かだけがあるのが良かったのでしょうけども、そうなっていないのは、国の事情、人間の歴史の複雑さの表れです。

3.  今後の発展 ~ マルチスケール模型は現場数学に不可欠

 最初に述べたように、計算機は5年で一桁速くなっていますので、そのうち、現実に存在する物質をそのまま計算機の中に設定して、それを電子状態から全て計算できるようになると考える人が多いと思います。しかし、それはほぼ不可能なのです。なぜならば、上述の話は、我々が対象とする物質は一様であると想定していました。しかし、本当の物質はそんなに単純ではありません。結晶も実は多結晶と呼ばれていて、単位となる結晶粒から構成されているのです。その各粒は1020くらいの原子から構成されているのが普通なのです。つまり、単位として解かなければならない結晶粒の精密計算が可能となるだけでも今後、80年もかかることになるのです。さらに、引っ張り強度などは、主として粒界(結晶粒間)の相互作用によるので…と…どんどん、気が遠くなるように必要な計算量は増えます。タンパク質の計算でも、現在は単一のタンパク質の3次元構造を問題としていますが、生体内ではタンパク質が単独で存在していることはなく、さらに水の中にあるのです。

 生体全体のシミュレーションという課題が最近流行っていますが、とても全てを第一原理計算で扱うような対象物ではありません。もし百年後に、全て第一原理で数値計算できたとしても、今度は単に膨大な数値が計算機から吐き出されるだけで、何を計算しているのかを人間が理解することが不可能になります。グラフィック処理して眺めれば良い、というサイズの問題さえも遙(はる)かに超えてしまうのです。つまり、対象物(シミュレーション計算結果)を理解するためにも、マルチスケール模型が必須なのです。これは、日常、人間を扱うときに、頭、手、足、胴体と分けて扱い、さらに、頭には目、口、鼻、耳と分離していることと同等なのです。そうでないと、各人の識別方法を伝えることさえも不可能になります。こうして、マルチスケール模型は将来とも必須な現場数学のシミュレーション技法なのです。

   続きを読むには・・・


この記事の著者

川添 良幸

市販の材料設計シミュレーションプログラムでは満足できない御社技術者に、本当に意味のある物理の基本に基づいた設計法を伝授します。我々の計算と実験結果が合わない場合は実験の方に問題があると言えるレベルを達成しています。

市販の材料設計シミュレーションプログラムでは満足できない御社技術者に、本当に意味のある物理の基本に基づいた設計法を伝授します。我々の計算と実験結果が合わな...


「SQC一般」の他のキーワード解説記事

もっと見る
~長方形の面積S=a×b?-誤差の概念、誤差の伝播― 現場数学(その3)

  1. 誤差の概念について  子どもの頃、面積を求める問題は分かりやすく、誰でも算数が出来た気になりました。その時、悩みもせず、例えば...

  1. 誤差の概念について  子どもの頃、面積を求める問題は分かりやすく、誰でも算数が出来た気になりました。その時、悩みもせず、例えば...


~空間と時間の刻みは独立ではない ~ 流体計算 現場数学(その8)

  ♦ 解の不安定性には常に注意! 1. 我々の住んでいる空間  次元という言葉は魅力的です。「ステレオでは味わえない5....

  ♦ 解の不安定性には常に注意! 1. 我々の住んでいる空間  次元という言葉は魅力的です。「ステレオでは味わえない5....


層別比較に便利なグラフ

 層別比較に便利なツールとしては、『箱ひげ図』があります。箱ひげ図は少数データの扱いも容易で作成も簡便なので便利ですが、分布の中心が中央値で、分布の幅がパ...

 層別比較に便利なツールとしては、『箱ひげ図』があります。箱ひげ図は少数データの扱いも容易で作成も簡便なので便利ですが、分布の中心が中央値で、分布の幅がパ...


「SQC一般」の活用事例

もっと見る
ビックデータ時代と米国大統領戦

 ビックデータ時代を考える事例として、今回の第45代米国大統領選挙を見てみます。第45代の米国大統領に、ドナルド・トランプ氏が決まりましたが、トランプ氏当...

 ビックデータ時代を考える事例として、今回の第45代米国大統領選挙を見てみます。第45代の米国大統領に、ドナルド・トランプ氏が決まりましたが、トランプ氏当...


統計教育が必須な検査担当者

 測定や検査結果の値で合格判定を行い、製品性能の品質保証をする事は製造業に取っては当たり前と言える日常的作業です。一方で検査や測定の正確さや信頼性を担保す...

 測定や検査結果の値で合格判定を行い、製品性能の品質保証をする事は製造業に取っては当たり前と言える日常的作業です。一方で検査や測定の正確さや信頼性を担保す...


統計手法による 2018年LPGA公式記録データの解析事例(その1)

【統計手法による 2018年LPGA公式記録データの解析 連載目次】 統計手法による 2018年LPGA公式記録データの解析事例(その1) 統...

【統計手法による 2018年LPGA公式記録データの解析 連載目次】 統計手法による 2018年LPGA公式記録データの解析事例(その1) 統...