全原子古典分子動力学法による高分子物性計算を全自動化するソフトウェアRadonPyをリリース 高分子材料物性大地図の作成に向けた第一歩

材料データとデータ科学・計算科学を融合した材料研究の新しい形態“マテリアルズインフォマティクス”（Materials Informatics: 以下、MIという。）がいま大きな注目を集めています。MIでは、機械学習を適用して、データのパターンを読み解き、広大なデザインスペースから革新的な特性を有する新材料とその作製方法を予測します。ここ数年間でMIの技術は材料研究の様々な領域に広まり、実際に多くの新材料が発見されてきました。しかしながら、高分子材料^[用語1]のMIの普及は、他の材料系に比べると大幅に遅れています。MIにおける最も重要な学術資源はいうまでもなくデータです。しかしながら、高分子物性のデータベースに関しては、国立研究開発法人物質・材料研究機構（以下、NIMSという。）が開発しているPoLyInfo^[用語2]（論文等の実験値を網羅的に収集したデータベース）を除き、データ駆動型研究に資する体系的なオープンデータを創出しようという動きは極めて低調です。このままでは高分子材料のMIに携わる多くの研究者にとって、大学の研究室や一企業で生産可能なデータが標準になってしまう可能性があります。

このような現状を打破するために、東京工業大学 物質理工学院 材料系の森川淳子教授、統計数理研究所（以下、統数研という。）の林慶浩助教と吉田亮教授、東京大学の塩見淳一郎教授、NIMSのPoLyInfo開発チームは、全原子古典分子動力学法^[用語3]（以下、MD計算あるいはMDシミュレーションという。）による高分子物性計算を全自動化するソフトウェアRadonPyを開発しました。高分子物性のMD計算は計算条件の設定により大きくゆらぎ、計算量も膨大であるため、自動計算による大規模データベースの構築は技術的に困難であると考えられてきました。RadonPyは、高分子材料の繰り返し単位（モノマー）の化学構造と重合度^[用語4]、温度等の計算条件を入力とし、アモルファスポリマー^[用語5a]や高分子溶液等の系に対し、熱物性、機械特性、光学特性を含む様々な物性を自動計算します。今回リリースしたバージョンは、15種類の物性の自動計算機能を実装しています。RadonPyには、多数の煩雑な手続きからなるMD計算の全工程を自動化し、パイプラインを簡便に構築するための諸機能が実装されています。RadonPyの特色の一つは、様々な骨格の高分子材料に適用可能な検証済みパラメータセットと計算条件（プリセット）を標準搭載していることです。プリセットの設定では、NIMSの協力を得てPoLyInfoの実験データをベンチマークとして使用しました。また、MD計算の物性値と複雑な現実系（実験値）の間のバイアスやばらつきを補正するために、転移学習^[用語6]という機械学習の解析技術を用いました。

同グループは、RadonPyを用いて10万種類以上の分子骨格を包含する高分子物性オープンデータベースを創出するプロジェクトを推進しています。現在、統数研に加えて3大学と19企業がプロジェクトに参画し、多くの研究者が産学の垣根を越えてRadonPyとデータベースの共同開発に取り組んでいます。この取り組みは文部科学省「富岳」成果創出加速プログラムの支援を受け、参画者は、スーパーコンピュータ「富岳」^[用語7]の計算資源を最大限に活用して日々膨大なデータを生産・蓄積しています。プロジェクトの目標は、高分子材料物性の大地図を作成することです。データを生産していく過程で、複数物性の同時分布や物性間のトレードオフが生み出すパレートフロンティア^[用語8]、さらにフロンティアを形成する新材料が明らかになっていきます。材料合成等の膨大なコストを伴う実験的アプローチのみでは、このような網羅的な観測は実現不可能です。世界最大級の二つの高分子物性データベース、PoLyInfoとRadonPyデータベースを統合的に活用することで、不確かで複雑な現実系と不完全な計算モデルの壁を乗り越えられるかもしれません。

本研究成果は 2022年11月8日10時（日本時間8日19時）に npj Computational Materials 誌（Nature Publishing Group）にて発表されました。

データ駆動型研究における最も重要な学術資源はいうまでもなくデータです。しかしながら、材料研究において利用できるデータ資源は、量的にも質的にも非常に乏しいという難しさがあります。特に、高分子材料のデータ資源の乏しさは際立っています。PoLyInfoを除く既存の高分子物性データベースは、いずれもデータ量が非常に少なく、データを自動抽出・成形・加工するためのツールが整備されていません（表1）。また、試料の作製条件や高次構造等の情報も十分には記録されていません。データ駆動型高分子材料研究に資する体系的かつ包括的なオープンデータの創出が喫緊の課題となっています。

表1.

代表的な高分子物性オープンデータの一覧。データ科学の他の応用分野に比べるとデータ量は非常に少なく、データを自動抽出するためのツールも整備されていないケースがほとんどである。

オープンデータの整備に遅れが生じている原因として、次の3点が考えられます。

このような背景からコミュニティが協調してオープンデータを創出していこうという動きが低調になっています。これらの要因の多くは文化的な問題であるため、短中期的には、大学の研究室や一企業で生産可能なレベルのデータが、MIの標準的な解析対象になっていくことが予想されます。

そこでMIでは、データの不足を補うためにシミュレーションから得られるデータを統合的に解析することが重要になってきます。現在、世界各国にて様々な材料を対象に大規模な計算物性データベースが開発されています。特に無機化合物や低分子化合物については、数万から数百万種類以上の物質を包含する第一原理計算データベース（Materials Project、AFLOW、OQMD、QM9等）が開発されたことで、MIの技術的進歩と実践展開が飛躍的に進みました。一方、高分子材料については、物性計算の自動化の技術的な難しさと分子シミュレーションの膨大な計算コストが壁となり、計算物性データベースの開発はほとんど進んでいません。

RadonPyは、全原子MDシミュレーションによる高分子物性計算を全自動化する世界初のオープンソースソフトウェアです。RadonPyは、高分子の繰り返し単位とその化学構造、重合度、温度等の計算条件を入力とし、分子モデリング、電荷計算、平衡・非平衡MD計算、平衡化完了の自動判定、収束に失敗した場合の再スタートのスケジューリング、ポストプロセス段階での物性計算等、MDシミュレーションの全工程を完全に自動実行します（図1）。ファーストリリースで公開されたバージョンは、熱物性、機械特性、光学特性を含む15種類の物性を自動計算するアルゴリズムを実装しています。計算可能な系は、線状高分子^[用語9]（ホモポリマー^[用語10]、コポリマー^[用語11]）のアモルファス状態^[用語5b]や溶液中の高分子等を含みます。

図1.

高分子物性自動計算ソフトウェアRadonPyの概要。高分子の繰り返し単位とその化学構造、重合度、温度等の計算条件を入力とし、分子モデリング、電荷計算、平衡・非平衡MD計算、平衡化完了の自動判定、収束に失敗した場合の再スタートのスケジューリング、ポストプロセス段階での物性計算等、MDシミュレーションの全工程を完全に自動実行する。ファーストリリースで公開されたバージョンは、熱物性や光学特性を含む15種類の物性を自動計算できる。

今回の研究は、大規模データベース構築に向けたパイロットスタディとして実施されました。本研究では、1,000種類以上のアモルファスポリマーを対象に15種類の物性を計算し、PoLyInfoの実験値を用いた検証実験が実施されました。重合度等の計算条件や分子骨格の種類が計算値に与える影響や、物性毎の予測精度や性能限界が徹底的に調べられました。このような規模で高分子物性MD計算の性能が系統的に検証された事例はありませんでした。

さらに、転移学習という機械学習の解析技術が、計算物性の系統バイアスやばらつきを補正する上で強力な性能を発揮することが明らかになりました（図2）。様々な高分子材料に対して普遍的に適用できる計算条件は存在しません。したがって、全自動計算で大量生産されるデータには、必ずバイアスやばらつきが生じます。また、実験条件や試料に関する非観測因子や計測系の特性により、実験値にもバイアスとばらつきが生じます。このような複雑な現実系と不完全な計算モデルの間のギャップを機械学習で埋めることができます。

図2.

転移学習によるMD計算物性（定圧比熱容量（C_p）、線膨張係数、体積膨張率）のキャリブレーション。MD計算値と実験値の間に存在するバイアスとばらつき（上）が転移学習を用いたキャリブレーションによって大幅に改善された（下）。これは実験と計算両方のバイアスとばらつきを転移学習は予測可能であることを意味する。

さらに今回の研究により、高分子材料の複数物性の同時分布やパレートフロンティアの位置、そこに存在する高分子材料の構造的特徴が明らかになってきました（図3）。特に本研究は高分子材料の熱伝導率に着目しました。近年のモバイル機器の小型化・高性能化に伴う発熱量の増加により、モバイル機器の絶縁樹脂、成形樹脂、接着剤、コーティング剤等に応用できる材料へのニーズが高まり、成形性に優れた高分子材料の高熱伝導化の研究に注目が集まっています。今回の計算では、熱伝導率が0.4 W/(m・K)を超える8種類のアモルファスポリマーが同定されました。通常のアモルファスポリマーで到達できる熱伝導率は高々0.2～0.3 W/(m・K) 程度といわれていますが、同定された材料はこの水準を大きく上回ります。また、RadonPyに実装されている熱伝導成分分割の解析手法を適用し、アモルファスポリマーの高熱伝導化のメカニズムと設計指針を明らかにしました。その結果、水素結合ユニットを高密度に含む高分子では、分子間の水素結合や双極子相互作用を介して熱伝導性が向上すること、あるいは高い硬直性と線形性を持つ高分子では共有結合を介して熱伝導率が向上することが明らかになりました。

図3.

RadonPyによるハイスループット自動計算で明らかとなった高分子材料の複数物性（熱伝導率、密度、定圧比熱容量 (C_p)、体積膨張率、線膨張係数、屈折率）の同時分布とパレートフロンティア

現在、同グループはRadonPyを用いて10万種類以上の分子骨格を包含する世界最大級の高分子物性データベースを創出するプロジェクトを推進しています。2022年10月には、RadonPyとデータベースの共同開発を目的とする産学協働コンソーシアムを正式に発足しました（非公式の活動は2021年4月に始動）。コンソーシアムには現在、統数研以外に3大学と企業19社が参画しており、約90名の研究者が産学の垣根を越えてRadonPyとデータベースの共同開発を推進しています。この取り組みは文部科学省「富岳」成果創出加速プログラムの支援を受け、産学の多くの研究者が「富岳」の計算資源を最大限に活用し、日々膨大なデータを生産・蓄積しています。プロジェクトが掲げるスローガンは、高分子材料物性大地図の作成です。データを生産・蓄積する過程で、複数物性の同時分布とパレートフロンティア、さらにフロンティアを形成する新材料を明らかにしていきます。特に、生分解性プラスチックや高熱伝導性高分子材料に関する体系的なデータを創出し、脱炭素・循環型社会やサーマルマネジメントに資する新材料の創製を目指します。また、産学の組織の垣根を超えたデータの共創というモデルケースを社会に発信していくこともコンソーシアムの活動目的の1つです。

本事業がマイルストーンに到達した暁には、PoLyInfoとRadonPyデータベースという実験・計算の世界最大級の高分子物性データベースが利用可能になります。全ての計算モデルは不完全です。計算資源の制約により、実材料開発における混錬や加熱、フィラー導入等のプロセス情報をMD計算に取り組むことは難しく、必ず計算結果には現実系との乖離が生じます。同様に、実験データも不完全であり不確かです。実験と計算という車の両輪は2つが揃ってこそ意味があり、データ科学が両輪を回すエンジンの役割を担います。

RadonPy 公開ウェブサイト
RadonPy project RadonPy｜GitHub