ホーム > 特許ランキング > 国立大学法人 東京大学
▶ 国立大学法人 東京大学の特許一覧 ▶ 学校法人北里研究所の特許一覧
特開2024-4296結晶構造予測方法、結晶構造予測装置、及び結晶構造予測プログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024004296
(43)【公開日】2024-01-16
(54)【発明の名称】結晶構造予測方法、結晶構造予測装置、及び結晶構造予測プログラム
(51)【国際特許分類】
G16C 20/20 20190101AFI20240109BHJP
【ＦＩ】
G16C20/20
【審査請求】未請求
【請求項の数】8
【出願形態】ＯＬ
(21)【出願番号】P 2022103891
(22)【出願日】2022-06-28
(71)【出願人】
【識別番号】504137912
【氏名又は名称】国立大学法人 東京大学
(71)【出願人】
【識別番号】598041566
【氏名又は名称】学校法人北里研究所
(74)【代理人】
【識別番号】100106909
【弁理士】
【氏名又は名称】棚井 澄雄
(74)【代理人】
【識別番号】100188558
【弁理士】
【氏名又は名称】飯田 雅人
(72)【発明者】
【氏名】渡辺 豪
(72)【発明者】
【氏名】岡本 敏宏
(72)【発明者】
【氏名】竹谷 純一
(72)【発明者】
【氏名】伊藤 良将
(72)【発明者】
【氏名】關 拓和
(72)【発明者】
【氏名】佐藤 俊輔
(57)【要約】
【課題】対象分子の分子構造から、当該対象分子により形成される結晶構造を、結晶構造における分子の熱的揺動の影響を考慮して、より高精度に予測できる。
【解決手段】結晶構造予測方法は、対象分子の分子構造に基づいて、前記対象分子の結晶構造の候補となる候補構造を特定する候補構造特定工程と、前記候補構造を用いた分子動力学計算を行うことにより前記候補構造の評価を行い、前記評価に基づいて前記対象分子の結晶構造を予測する結晶構造予測工程と、を含む。
【選択図】図１０
【解決手段】結晶構造予測方法は、対象分子の分子構造に基づいて、前記対象分子の結晶構造の候補となる候補構造を特定する候補構造特定工程と、前記候補構造を用いた分子動力学計算を行うことにより前記候補構造の評価を行い、前記評価に基づいて前記対象分子の結晶構造を予測する結晶構造予測工程と、を含む。
【選択図】図１０
【特許請求の範囲】
【請求項1】
対象分子の分子構造に基づいて、前記対象分子の結晶構造の候補となる候補構造を特定する候補構造特定工程と、
前記候補構造を用いた分子動力学計算を行うことにより前記候補構造の評価を行い、前記評価に基づいて前記対象分子の結晶構造を予測する結晶構造予測工程と、
を含むことを特徴とする結晶構造予測方法。
前記候補構造を用いた分子動力学計算を行うことにより前記候補構造の評価を行い、前記評価に基づいて前記対象分子の結晶構造を予測する結晶構造予測工程と、
を含むことを特徴とする結晶構造予測方法。
【請求項2】
前記結晶構造予測工程において、前記分子動力学計算から得た前記候補構造のエネルギーに基づいて前記評価を行う、
請求項１に記載の結晶構造予測方法。
請求項１に記載の結晶構造予測方法。
【請求項3】
前記候補構造を、そのエネルギーが低い順に順位付けし、前記候補構造の前記順位を、前記結晶構造の予測の優先順位とする、
請求項２に記載の結晶構造予測方法。
請求項２に記載の結晶構造予測方法。
【請求項4】
前記候補構造のエネルギーが、前記候補構造におけるポテンシャルエネルギーと運動エネルギーの和を含む、
請求項２から３のいずれかに記載の結晶構造予測方法。
請求項２から３のいずれかに記載の結晶構造予測方法。
【請求項5】
前記結晶構造予測工程において、所定の温度を設定して前記分子動力学計算を行う、
請求項１から３のいずれかに記載の結晶構造予測方法。
請求項１から３のいずれかに記載の結晶構造予測方法。
【請求項6】
前記対象分子が有機半導体である、
請求項１から３のいずれかに記載の結晶構造予測方法。
請求項１から３のいずれかに記載の結晶構造予測方法。
【請求項7】
対象分子の分子構造に基づいて、前記対象分子の結晶構造の候補となる候補構造を特定する候補構造特定手段と、
前記候補構造を用いた分子動力学計算を行うことにより前記候補構造の評価を行い、前記評価に基づいて前記対象分子の結晶構造を予測する結晶構造予測手段と、
を有することを特徴とする結晶構造予測装置。
前記候補構造を用いた分子動力学計算を行うことにより前記候補構造の評価を行い、前記評価に基づいて前記対象分子の結晶構造を予測する結晶構造予測手段と、
を有することを特徴とする結晶構造予測装置。
【請求項8】
対象分子の分子構造に基づいて、前記対象分子の結晶構造の候補となる候補構造を特定する候補構造特定処理と、
前記候補構造を用いた分子動力学計算を行うことにより前記候補構造の評価を行い、前記評価に基づいて前記対象分子の結晶構造を予測する結晶構造予測処理と、
をコンピュータに行わせることを特徴とする結晶構造予測プログラム。
前記候補構造を用いた分子動力学計算を行うことにより前記候補構造の評価を行い、前記評価に基づいて前記対象分子の結晶構造を予測する結晶構造予測処理と、
をコンピュータに行わせることを特徴とする結晶構造予測プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本件は、結晶構造予測方法、結晶構造予測装置、及び結晶構造予測プログラムに関する。
【背景技術】
【0002】
今日、有機半導体はフレキシブルな次世代デバイスの材料として注目されており、より高性能の革新的な有機半導体が必要とされている。有機半導体の性能には、その結晶構造が大きく影響するが、有機半導体の設計を実験的な手法に基づいて行う場合、有機半導体の結晶化等に膨大な時間的・金銭的コストが必要となり、新規な有機半導体を効率的に設計することができない。
【0003】
また、有機半導体などの分子の結晶構造を探索するための計算科学的な手法（プログラム）も提案されている。具体的には、例えば、分子（対象分子）の分子構造情報（１分子の配座の情報）を用いて結晶構造の候補を探索し、探索した結晶構造候補に対して理論的に得られた粉末Ｘ線回折パターンと、実際に計測で得られた粉末Ｘ線回折パターンとの類似性を評価することで、結晶構造の予測精度を高める技術が提案されている（例えば、非特許文献１参照）。
【0004】
しかしながら、上述した解析手法を用いることで多くの結晶構造候補が生成される。この解析手法によって探索された多くの結晶構造候補には、正しい構造である可能性が高い順（例えば、エネルギーが低い順）で順位付けがなされても、正しい結晶構造に対応する候補が下位に順位付けされてしまうなど、予測される結晶構造の信頼性が十分ではなかった。その解析結果から得られた多くの結晶構造候補のなかからより適した結晶構造を予測するためには、その結晶構造候補の分子の合成や粉末Ｘ線回折といった検証のための実験が夫々必要になり、これに係る時間的・金銭的コストが大きくなることがあった。
このように、対象分子の分子構造情報に基づいて当該分子の結晶構造を予測する解析手法を単独で利用しても、適した結晶構造を高精度に予測できないことがあった。
このように、対象分子の分子構造情報に基づいて当該分子の結晶構造を予測する解析手法を単独で利用しても、適した結晶構造を高精度に予測できないことがあった。
【先行技術文献】
【非特許文献】
【0005】
【非特許文献1】Ｉｓｈｉｉ， Ｈ．， Ｏｂａｔａ， Ｓ．， Ｎｉｉｔｓｕ， Ｎ． ｅｔ ａｌ．， Ｃｈａｒｇｅ ｍｏｂｉｌｉｔｙ ｃａｌｃｕｌａｔｉｏｎ ｏｆ ｏｒｇａｎｉｃ ｓｅｍｉｃｏｎｄｕｃｔｏｒｓ ｗｉｔｈｏｕｔ ｕｓｅ ｏｆ ｅｘｐｅｒｉｍｅｎｔａｌ ｓｉｎｇｌｅ－ｃｒｙｓｔａｌ ｄａｔａ， Ｓｃｉ． Ｒｅｐ． １０， ２５２４ （２０２０）．
【発明の概要】
【発明が解決しようとする課題】
【0006】
本発明は、一つの側面では、対象分子の分子構造から、当該対象分子により形成される結晶構造を、結晶構造における分子の熱的揺動の影響を考慮して、より高精度に予測できる結晶構造予測方法、結晶構造予測装置、及び結晶構造予測プログラムを提供することを目的とする。
【課題を解決するための手段】
【0007】
上記の課題を解決するための手段の一つの実施態様は、以下のとおりである。
すなわち、一つの実施態様では、結晶構造予測方法は、対象分子の分子構造に基づいて、対象分子の結晶構造の候補となる候補構造を特定する候補構造特定工程と、候補構造を用いた分子動力学計算を行うことにより候補構造の評価を行い、評価に基づいて対象分子の結晶構造を予測する結晶構造予測工程と、を含む。
すなわち、一つの実施態様では、結晶構造予測方法は、対象分子の分子構造に基づいて、対象分子の結晶構造の候補となる候補構造を特定する候補構造特定工程と、候補構造を用いた分子動力学計算を行うことにより候補構造の評価を行い、評価に基づいて対象分子の結晶構造を予測する結晶構造予測工程と、を含む。
【0008】
また、一つの実施態様では、結晶構造予測装置は、対象分子の分子構造に基づいて、対象分子の結晶構造の候補となる候補構造を特定する候補構造特定手段と、候補構造を用いた分子動力学計算を行うことにより候補構造の評価を行い、評価に基づいて対象分子の結晶構造を予測する結晶構造予測手段と、を有する。
【0009】
さらに、一つの実施態様では、結晶構造予測プログラムは、対象分子の分子構造に基づいて、対象分子の結晶構造の候補となる候補構造を特定する候補構造特定処理と、候補構造を用いた分子動力学計算を行うことにより候補構造の評価を行い、評価に基づいて対象分子の結晶構造を予測する結晶構造予測処理と、をコンピュータに行わせる。
【発明の効果】
【0010】
本件は、一つの側面では、分子構造（１分子の配座）の情報から、当該分子により形成される結晶構造を、結晶構造における分子の熱的揺動の影響を考慮して高精度に予測できる結晶構造予測方法、結晶構造予測装置、及び結晶構造予測プログラムを提供できる。
【図面の簡単な説明】
【0011】
【図1】第１の実施形態の結晶構造予測装置の概略構成図である。
【図2】実施形態の処理部の構成図である。
【図3】実施形態の対象分子とその結晶構造について説明するための図である。
【図4A】実施形態の分子構造が異なる対象分子の結晶構造の一例を示す図である。
【図4B】実施形態の分子構造が異なる対象分子の結晶構造の一例を示す図である。
【図5】実施形態の対象分子の構造式を示す図である。
【図6】実施形態のＣＯＮＦＬＥＸによる抽出結果と、分子動力学計算における全エネルギーの対応を説明するための図である。
【図7】実施形態のＣＯＮＦＬＥＸによる抽出結果の例に基づいた２つの分子を空間内に配置する例を説明するための図である。
【図8】実施形態のＣＯＮＦＬＥＸによる抽出結果の例に基づいた分子を結晶構造にして空間内に配置する例を説明するための図である。
【図9】実施形態のＣＯＮＦＬＥＸによる抽出結果と、分子動力学計算の解析結果を説明するための図である。
【図10】実施形態における、対象分子の結晶構造を予測する際の流れの一例を示すフローチャートである。
【図11】第２の実施形態の対象分子の構造式を示す図である。
【図12】実施形態のＣＯＮＦＬＥＸによる抽出結果と、分子動力学計算における全エネルギーの対応を説明するための図である。
【図13】実施形態のＣＯＮＦＬＥＸによる抽出結果と、分子動力学計算の解析結果を説明するための図である。
【発明を実施するための形態】
【0012】
＜第１の実施形態＞
以下、本発明の第１の実施形態に係る結晶構造予測装置について、図を参照しながら説明する。
以下、本発明の第１の実施形態に係る結晶構造予測装置について、図を参照しながら説明する。
【0013】
本件は、「分子構造（１分子の配座）の情報から、当該分子により形成される結晶構造を予測する際の予測精度が十分でなかった」という、本発明者らの知見に基づくものである。より具体的には、本件は、分子により形成される結晶構造の安定性には当該分子の熱的揺動が影響することから、結晶構造を高精度に予測するためには、結晶構造における分子の熱的揺動の影響を考慮した評価を行う必要があるという、本発明者らの知見に基づくものである。
【0014】
本件においては、上述した結晶構造における分子の熱的揺動の影響を考慮した評価を行うために、まず、対象分子の分子構造に基づいて、対象分子の結晶構造の候補となる候補構造を特定した後、候補構造を用いた分子動力学計算による候補構造の評価を解析的に行い、評価に基づいて対象分子の結晶構造を予測する。
つまり、本件では、特定した候補構造を用いて、熱的揺動を含めた分子のダイナミクスをシミュレート可能な分子動力学計算を行うことにより、結晶構造における分子の熱的揺動の影響を考慮した評価を行うことができるため、従来技術と比べて非常に高精度に結晶構造を予測することができる。
つまり、本件では、特定した候補構造を用いて、熱的揺動を含めた分子のダイナミクスをシミュレート可能な分子動力学計算を行うことにより、結晶構造における分子の熱的揺動の影響を考慮した評価を行うことができるため、従来技術と比べて非常に高精度に結晶構造を予測することができる。
【0015】
【0016】
結晶構造予測装置１０は、表示部１と、入力部２と、処理部３と、主記憶部４と、出力部５と、外部インターフェース（Ｉ／Ｏ）部６と、記憶部７とを備える。結晶構造予測装置１０の各部は例えばバス８に接続されている。このような結晶構造予測装置１０は、例えば、結晶構造予測のためのプログラムを実行するコンピュータ等により構築することができる。
【0017】
表示部１は、例えば液晶ディスプレイ等の表示デバイスを含み、結晶構造予測装置１０の操作、結晶構造予測装置１０の解析結果などを表示デバイスに表示させる。入力部２は、キーボード等の各種入力手段を含む。例えば、入力部２に対するユーザの操作によって、処理に必要なデータや情報が入力される。処理部３は、ＣＰＵなどのプロセッサを含み、結晶構造予測のための各種制御、演算等を行う。主記憶部４は結晶構造予測のためのプログラム等を格納している。出力部５は、プリンタ等の外部装置に処理結果を出力する。外部インターフェース（Ｉ／Ｏ）部６は、ＬＡＮやインターネット等のネットワークＮＷに接続され、ネットワークを介して他の端末装置２０等と通信する。記憶部７は、少なくとも、候補構造特定処理と、結晶構造予測処理を含む各種プログラム７１と、各種データ等を書き換え可能に記憶する。分子構造情報７２と、候補構造情報７３と、初期構造情報７４と、分子動力学解析結果７５は、上記の各種データ等の一例である。結晶構造予測装置１０は、１つの装置として構成することに制限はなく、適宜分割することができる。
例えば、結晶構造予測装置１０を計算サーバにして、作業用の端末装置２０を用いるサーバクライアント型に構成することができる。この場合、解析用の主たるプログラムを計算サーバの結晶構造予測装置１０に配置して、端末装置２０からの要求に応じて、計算サーバが要求に応じた処理を実行するとよい。
また、結晶構造予測装置１０を処理の種類に応じて、結晶構造予測装置１０Ａ、１０Ｂに分割してもよい。結晶構造予測装置１０Ａ、１０Ｂは、互いに異なる種類の処理を実施するように構成される。この場合、結晶構造予測装置１０Ａ、１０Ｂと端末装置２０が、互いに通信可能に構成されていて、各装置間のデータが通信を利用して、又は共有メモリを利用して相互にデータを送受することができる。
上記の結晶構造予測装置１０Ａ、１０Ｂと端末装置２０は、図１に示す結晶構造予測装置１０と同様な構成を備えるものであってよい。
なお、以下の説明では、結晶構造予測装置１０Ａ、１０Ｂと端末装置２０などを、１つの装置（結晶構造予測装置１０）に纏めて構成した場合を中心に説明する。
例えば、結晶構造予測装置１０を計算サーバにして、作業用の端末装置２０を用いるサーバクライアント型に構成することができる。この場合、解析用の主たるプログラムを計算サーバの結晶構造予測装置１０に配置して、端末装置２０からの要求に応じて、計算サーバが要求に応じた処理を実行するとよい。
また、結晶構造予測装置１０を処理の種類に応じて、結晶構造予測装置１０Ａ、１０Ｂに分割してもよい。結晶構造予測装置１０Ａ、１０Ｂは、互いに異なる種類の処理を実施するように構成される。この場合、結晶構造予測装置１０Ａ、１０Ｂと端末装置２０が、互いに通信可能に構成されていて、各装置間のデータが通信を利用して、又は共有メモリを利用して相互にデータを送受することができる。
上記の結晶構造予測装置１０Ａ、１０Ｂと端末装置２０は、図１に示す結晶構造予測装置１０と同様な構成を備えるものであってよい。
なお、以下の説明では、結晶構造予測装置１０Ａ、１０Ｂと端末装置２０などを、１つの装置（結晶構造予測装置１０）に纏めて構成した場合を中心に説明する。
【0018】
結晶構造予測装置１０の処理部３は、例えば、分子構造作成部３１と、結晶構造生成部３２と、初期結晶構造生成部３３と、分子動力学計算部３４と、エネルギー算出部３５と、結晶構造予測部３６と、を備える。処理部３内の各部の詳細については、後述する。
例えば、結晶構造生成部３２と、初期結晶構造生成部３３は、候補構造特定手段３１０の一例である。候補構造特定手段３１０は、分子構造作成部３１を含めて構成されてもよい。分子動力学計算部３４と、エネルギー算出部３５と、結晶構造予測部３６は、結晶構造予測手段３２０の一例である。上記の機能分割は一例であり、これに制限されない。例えば、初期結晶構造生成部３３は、候補構造特定手段３１０に含まれるものではなく結晶構造予測手段３２０に含めてもよい。
例えば、結晶構造生成部３２と、初期結晶構造生成部３３は、候補構造特定手段３１０の一例である。候補構造特定手段３１０は、分子構造作成部３１を含めて構成されてもよい。分子動力学計算部３４と、エネルギー算出部３５と、結晶構造予測部３６は、結晶構造予測手段３２０の一例である。上記の機能分割は一例であり、これに制限されない。例えば、初期結晶構造生成部３３は、候補構造特定手段３１０に含まれるものではなく結晶構造予測手段３２０に含めてもよい。
【0019】
上記の候補構造特定手段３１０は、対象分子の分子構造に基づいて、対象分子の結晶構造の候補となる候補構造を特定する。結晶構造予測手段３２０は、候補構造を用いた分子動力学計算を行うことにより候補構造の評価を行い、評価に基づいて対象分子の結晶構造を予測するとよい。
本件で開示する技術の一例としての結晶構造予測装置１０は、上述したように、候補構造特定手段３１０と、結晶構造予測手段３２０とを少なくとも有し、必要に応じてその他の手段を有してよい。
本件で開示する技術の一例としての結晶構造予測装置１０は、上述したように、候補構造特定手段３１０と、結晶構造予測手段３２０とを少なくとも有し、必要に応じてその他の手段を有してよい。
【0020】
結晶構造予測装置１０は、例えば、後述する結晶構造予測方法に基づく各種処理を実行する装置として形成されているとよい。また、本件の結晶構造予測装置における好適な態様は、例えば、本件の結晶構造予測方法における好適な態様と同様とすることができる。
【0021】
（結晶構造予測方法）
本件で開示する技術の一例としての結晶構造予測方法は、候補構造特定工程（例えば、後述の図１１のステップＳ１０２）と、結晶構造予測工程（例えば、後述の図１０のステップＳ１０３－Ｓ１０５）とを少なくとも含み、必要に応じてその他の工程を含む。また、本件の結晶構造予測方法は、例えば、これらの工程を、コンピュータが行うことにより実施することができる。さらに、本件の結晶構造予測方法は、例えば、本件の結晶構造予測装置１０、本件の結晶構造予測プログラムの実行などにより行うことができる。
本件で開示する技術の一例としての結晶構造予測方法は、候補構造特定工程（例えば、後述の図１１のステップＳ１０２）と、結晶構造予測工程（例えば、後述の図１０のステップＳ１０３－Ｓ１０５）とを少なくとも含み、必要に応じてその他の工程を含む。また、本件の結晶構造予測方法は、例えば、これらの工程を、コンピュータが行うことにより実施することができる。さらに、本件の結晶構造予測方法は、例えば、本件の結晶構造予測装置１０、本件の結晶構造予測プログラムの実行などにより行うことができる。
【0022】
＜候補構造特定工程＞
まず、候補構造特定工程について説明する。
結晶構造予測装置１０の処理部３は、候補構造特定工程のなかで、対象分子の分子構造に基づいて、対象分子の結晶構造の候補となる候補構造を特定する。
まず、候補構造特定工程について説明する。
結晶構造予測装置１０の処理部３は、候補構造特定工程のなかで、対象分子の分子構造に基づいて、対象分子の結晶構造の候補となる候補構造を特定する。
【0023】
図３は、実施形態の対象分子とその結晶構造について説明するための図である。結晶構造を形成する単位格子内に対象分子はエネルギー的に安定な状態になるような位置に置かれている。この結晶構造における対象分子の位置は、対象分子の分子配座（立体構造における幾何的な原子の位置関係）により異なる場合がある。
【0024】
図４Ａと図４Ｂは、実施形態の分子配座が異なる対象分子の結晶構造の一例を示す図である。図４Ａと図４Ｂに示す対象分子は、同一の構造式で表されるが、配座やトポロジカルな構造が異なる。このように、同一の対象分子であっても、分子配座が異なる場合などには、単位格子内における分子の位置が大きく異なる場合がある。
【0025】
例えば、結晶構造生成部３２と初期結晶構造生成部３３は、候補構造特定手段３１０として、候補構造特定工程のなかで、対象分子の分子構造（１分子の配座）の情報を用いて、当該対象分子の結晶構造の候補となる候補構造を特定する。対象分子の候補構造を特定することができれば、具体的な手法に特に制限はなく目的に応じて適宜選択することができる。なお、対象分子の分子構造の情報には、対象分子の配座の情報が含まれる。結晶構造生成部３２と初期結晶構造生成部３３は、候補構造特定工程において、少なくとも対象分子の１分子の配座の情報を用いて特定される対象分子に適した候補構造を識別する。対象分子の分子構造（１分子の配座）の情報は、記憶部７の分子構造情報７２に格納されているとよい。
【0026】
＜＜対象分子＞＞
図５に、実施形態の対象分子の一例を構造図に示す。図５は、実施形態の対象分子の構造式を示す図である。なお、以下の説明に示す具体例（実施例）の対象分子は、例えばＰｈＣ２－ＢＱＱＤＩである。対象分子は、これに制限されず、以下に示すように適宜変更できる。
図５に、実施形態の対象分子の一例を構造図に示す。図５は、実施形態の対象分子の構造式を示す図である。なお、以下の説明に示す具体例（実施例）の対象分子は、例えばＰｈＣ２－ＢＱＱＤＩである。対象分子は、これに制限されず、以下に示すように適宜変更できる。
【0027】
例えば、対象分子としては、集合して結晶を形成しうる分子であれば、特に制限はなく目的に応じて適宜選択することができる。例えば、分子結晶を形成しうる分子とすることができる。分子結晶を形成しうる分子としては、例えば、結晶性の有機化合物などが挙げられる。結晶性の有機化合物としては、例えば、π－π相互作用によりスタッキング可能な低分子化合物を好適に用いることができ、具体的には、π電子共役構造（環状構造など）を有する芳香族化合物などを用いることができる。
本件の技術では、π－π相互作用が見込める剛直な中心骨格（コア）を有する分子に対して、特に高精度で結晶構造を予測することができる。
例えば、中心骨格（コア）は、benzo[de]isoquinolino[1,8-gh]quinoline diimide （ＢＱＱＤＩ)骨格をなす。この具体例のＰｈＣ２－ＢＱＱＤＩの中心骨格（コア）において、上記の中心骨格の両端部に、分子が延伸する方向にイミド基が配置されている。各イミド基には、コアに対する側鎖として、芳香環を含む官能基（例えば、フェネチル基：－Ｃ２Ｈ４Ｐｈ）、アルキル鎖（アルキル基：－Ｒ）などがそれぞれ配置されている。このような構造を有する結晶性の有機化合物に対して、本件の技術を好適に適用可能である。別の具体例については後述する。
本件の技術では、π－π相互作用が見込める剛直な中心骨格（コア）を有する分子に対して、特に高精度で結晶構造を予測することができる。
例えば、中心骨格（コア）は、benzo[de]isoquinolino[1,8-gh]quinoline diimide （ＢＱＱＤＩ)骨格をなす。この具体例のＰｈＣ２－ＢＱＱＤＩの中心骨格（コア）において、上記の中心骨格の両端部に、分子が延伸する方向にイミド基が配置されている。各イミド基には、コアに対する側鎖として、芳香環を含む官能基（例えば、フェネチル基：－Ｃ２Ｈ４Ｐｈ）、アルキル鎖（アルキル基：－Ｒ）などがそれぞれ配置されている。このような構造を有する結晶性の有機化合物に対して、本件の技術を好適に適用可能である。別の具体例については後述する。
【0028】
また、本件において、結晶性の有機化合物としては、有機半導体として用いることができるものであることが好ましい。言い換えると、本件では、対象分子が有機半導体であることが好ましい。なお、「有機半導体」は、有機半導体の分子を示すことがある。
ところで、軽量かつフレキシブルな次世代電子デバイスの実現に向けて有機半導体の性能の向上が求められており、有機半導体の性能に直結する結晶構造を予測することの重要性は極めて高い。本件の技術では、上述したように、一般にπ電子共役構造（環状構造など）を有する有機半導体の結晶構造について、特に高精度に予測することができる。
ところで、軽量かつフレキシブルな次世代電子デバイスの実現に向けて有機半導体の性能の向上が求められており、有機半導体の性能に直結する結晶構造を予測することの重要性は極めて高い。本件の技術では、上述したように、一般にπ電子共役構造（環状構造など）を有する有機半導体の結晶構造について、特に高精度に予測することができる。
【0029】
また、結晶構造生成部３２による、対象分子の分子構造に基づいて候補構造を特定する手法としては、結晶構造の候補となる構造を分子構造の情報から生成できる手法であれば、特に制限はなく目的に応じて適宜選択することができる。候補構造の特定に、例えば、公知の結晶構造探索プログラム（ソフトウェア）を用いることもでき、具体的には、ＣＯＮＦＬＥＸ（登録商標）を用いることができる。
ＣＯＮＦＬＥＸを用いて候補構造を特定する場合、ＣＯＮＦＬＥＸにおける計算条件等は適宜選択することができる。ＣＯＮＦＬＥＸにおける計算条件としては、例えば、力場パラメータ、結晶構造を探索する空間群、構造最適化アルゴリズムなどが挙げられる。例えば、ＣＯＮＦＬＥＸによって、分子の配座空間を探索するなどの方法で、基準にする対象分子の配座異性体の最適化構造を検索して、その検索によって抽出したなかから対象分子（モデル分子）の候補構造を特定してもよい。
ＣＯＮＦＬＥＸを用いて候補構造を特定する場合、ＣＯＮＦＬＥＸにおける計算条件等は適宜選択することができる。ＣＯＮＦＬＥＸにおける計算条件としては、例えば、力場パラメータ、結晶構造を探索する空間群、構造最適化アルゴリズムなどが挙げられる。例えば、ＣＯＮＦＬＥＸによって、分子の配座空間を探索するなどの方法で、基準にする対象分子の配座異性体の最適化構造を検索して、その検索によって抽出したなかから対象分子（モデル分子）の候補構造を特定してもよい。
【0030】
なお、特定する候補構造の数（候補構造の生成数）としては、特に制限はなく目的に応じて適宜選択することができる。なお、特定する候補構造の数を、後述する結晶構造予測工程における分子動力学計算の計算コストや想定される候補構造の精度に合わせて選択することが好ましい。
結晶構造生成部３２による候補構造の特定には、結晶エネルギー（Ｅｃｒｙｓｔ）を利用するとよい。結晶エネルギー（Ｅｃｒｙｓｔ）は、次の式（１）によって規定される。
結晶構造生成部３２による候補構造の特定には、結晶エネルギー（Ｅｃｒｙｓｔ）を利用するとよい。結晶エネルギー（Ｅｃｒｙｓｔ）は、次の式（１）によって規定される。
【0031】
Ｅｃｒｙｓｔ＝Ｅｉｎｔｒａ＋Ｅｉｎｔｅｒ （１）
【0032】
上記の式（１）におけるＥｉｎｔｒａが分子内相互作用（非結合相互作用）に係るエネルギーであり、Ｅｉｎｔｅｒが分子間相互作用（結合相互作用）に係るエネルギーである。このように結晶エネルギー（Ｅｃｒｙｓｔ）が、分子内相互作用と分子間相互作用とに基づくものとして規定できる。
【0033】
結晶エネルギー（Ｅｃｒｙｓｔ）はＣＯＮＦＬＥＸを用いて算出できる。ＣＯＮＦＬＥＸによる抽出結果を用いて候補構造を特定する場合には、例えば、ＣＯＮＦＬＥＸにより出力される結晶エネルギー（Ｅｃｒｙｓｔ）が低い順に所定数の候補構造を特定することが好ましい。この場合、本件の技術においては、ＣＯＮＦＬＥＸにより出力された候補構造の全てについて、後述する分子動力学計算を実行して評価を行うことは必須ではなく、適宜選択した候補構造に対して、結晶構造予測工程の分子動力学計算による評価を行ってもよい。
【0034】
図６から図８を参照して、実施形態のＣＯＮＦＬＥＸによる抽出結果について説明する。
図６は、実施形態のＣＯＮＦＬＥＸによる抽出結果と、分子動力学計算における全エネルギーの対応を説明するための図である。図７は、実施形態のＣＯＮＦＬＥＸによる抽出結果の例に基づいた２つの分子を空間内に配置する例を説明するための図である。図８は、実施形態のＣＯＮＦＬＥＸによる抽出結果の例に基づいた分子を結晶構造にして空間内に配置する例を説明するための図である。このように分子配置は、結晶構造に関係する。
図６は、実施形態のＣＯＮＦＬＥＸによる抽出結果と、分子動力学計算における全エネルギーの対応を説明するための図である。図７は、実施形態のＣＯＮＦＬＥＸによる抽出結果の例に基づいた２つの分子を空間内に配置する例を説明するための図である。図８は、実施形態のＣＯＮＦＬＥＸによる抽出結果の例に基づいた分子を結晶構造にして空間内に配置する例を説明するための図である。このように分子配置は、結晶構造に関係する。
【0035】
例えば、ＣＯＮＦＬＥＸによって、対象分子の複数の候補構造が抽出され、抽出結果が候補構造情報７３に格納される。図６に示すグラフの横軸は、ＣＯＮＦＬＥＸの解析結果によって順位付けされた上位から３５個の候補構造に対応付けられている。この図の場合、横軸の左端の原点に近いほど結晶エネルギーが低く、実際の結晶構造に近い候補構造である可能性が高いものであると、ＣＯＮＦＬＥＸが判定したものになる。
なお、この順位は、候補構造を識別するため、その名称（ＩＤ）の最後にこの順位を付与し、「ＰｈＣ２－ＢＱＱＤＩ－１」や「候補１」のような名称で候補構造を呼称する。図６に示すグラフの縦軸は、各候補構造について検証を行った結果に基づく分子動力学計算における全エネルギーの大きさを示す。図６に示すようにＣＯＮＦＬＥＸの解析結果による順位に対して、全エネルギーの大きさとの間に相関性が低くなる場合があることが明らかになった。なお、分子動力学計算における全エネルギーの詳細については後述する。
なお、この順位は、候補構造を識別するため、その名称（ＩＤ）の最後にこの順位を付与し、「ＰｈＣ２－ＢＱＱＤＩ－１」や「候補１」のような名称で候補構造を呼称する。図６に示すグラフの縦軸は、各候補構造について検証を行った結果に基づく分子動力学計算における全エネルギーの大きさを示す。図６に示すようにＣＯＮＦＬＥＸの解析結果による順位に対して、全エネルギーの大きさとの間に相関性が低くなる場合があることが明らかになった。なお、分子動力学計算における全エネルギーの詳細については後述する。
【0036】
以下、ＣＯＮＦＬＥＸの結果よりも解析精度を高める分子動力学計算の手法の詳細について説明する。
なお、後述する分子動力学計算のために、上記の結晶構造生成部３２（ＣＯＮＦＬＥＸ）によって抽出する候補構造の個数は、予め設定されたパラメータによって規定するとよい。
初期結晶構造生成部３３は、それらの候補構造のなかで、結晶エネルギー（Ｅｃｒｙｓｔ）が比較的低い所定数の候補構造を選択する。以下の説明では、例えば７個を例示する。
なお、後述する分子動力学計算のために、上記の結晶構造生成部３２（ＣＯＮＦＬＥＸ）によって抽出する候補構造の個数は、予め設定されたパラメータによって規定するとよい。
初期結晶構造生成部３３は、それらの候補構造のなかで、結晶エネルギー（Ｅｃｒｙｓｔ）が比較的低い所定数の候補構造を選択する。以下の説明では、例えば７個を例示する。
【0037】
図７の（ａ）から（ｇ）に、前述の図６に示した候補構造について、その空間内の分子配置をモデル化して夫々示す。このように、同じ分子式で識別される分子であっても、分子内の配座の違いにより、個々の分子が単位格子内に配置される位置が異なるものになる。
【0038】
図８の（ａ）から（ｇ）に、図７と同様に、前述の図６で示したＣＯＮＦＬＥＸでの順位の番号を示した候補構造の結晶構造をモデル化して夫々示す。このように、個々の分子の構造によって、これらが結晶構造で空間内に配置されると、互いに異なる配置になり、幾何学的に異なっていることが読み取れる。このように、同じ分子式で識別される分子であっても、分子内の配座の違いにより、これが結晶構造で空間内に配置されると、互いに異なる配置で安定することがある。本件の技術は、上記の夫々の配置例の安定性を、以下に示す「結晶構造予測工程」の手法を用いて評価する。
【0039】
＜結晶構造予測工程＞
次に、結晶構造予測工程について説明する。
結晶構造予測装置１０の処理部３は、結晶構造予測工程のなかで、候補構造を用いた分子動力学計算を行うことにより候補構造の評価を行い、評価に基づいて対象分子の結晶構造を予測する。例えば、分子動力学計算部３４と、エネルギー算出部３５と、結晶構造予測部３６は、結晶構造予測工程のなかで、候補構造特定工程において特定された候補構造を利用した分子動力学計算に基づく評価を行い結晶構造を予測する。結晶構造を予測することができれば、具体的な手法に特に制限はなく目的に応じて適宜選択が可能である。
次に、結晶構造予測工程について説明する。
結晶構造予測装置１０の処理部３は、結晶構造予測工程のなかで、候補構造を用いた分子動力学計算を行うことにより候補構造の評価を行い、評価に基づいて対象分子の結晶構造を予測する。例えば、分子動力学計算部３４と、エネルギー算出部３５と、結晶構造予測部３６は、結晶構造予測工程のなかで、候補構造特定工程において特定された候補構造を利用した分子動力学計算に基づく評価を行い結晶構造を予測する。結晶構造を予測することができれば、具体的な手法に特に制限はなく目的に応じて適宜選択が可能である。
【0040】
ここで、候補構造を用いた分子動力学計算としては、例えば、候補構造に基づいた初期構造に対する分子動力学計算とすることができる。候補構造に基づいた初期構造としては、特に制限はなく目的に応じて適宜選択することができ、例えば、特定された候補構造自体であってもよく、候補構造をさらに積層した構造であってもよい。前述の初期結晶構造生成部３３により生成された候補構造に基づいた初期構造を適用できる。
なお、ＣＯＮＦＬＥＸにより候補構造を特定する場合、通常は、候補構造として結晶構造における単位構造（ユニットセル）が出力されるため、分子動力学計算を行う際には、出力された単位構造を３次元的に積層して数百分子を含む構造を生成して、生成した大きな構造に対して分子動力学計算を行うことにより、結晶構造の予測精度をより向上させることができる。
なお、ＣＯＮＦＬＥＸにより候補構造を特定する場合、通常は、候補構造として結晶構造における単位構造（ユニットセル）が出力されるため、分子動力学計算を行う際には、出力された単位構造を３次元的に積層して数百分子を含む構造を生成して、生成した大きな構造に対して分子動力学計算を行うことにより、結晶構造の予測精度をより向上させることができる。
【0041】
＜＜分子動力学計算＞＞
分子動力学計算部３４は、結晶構造予測工程において分子動力学計算を実施する。分子動力学計算部３４は、この分子動力学計算において、候補構造（候補構造を積層した構造を含む）における、熱的揺動（熱揺らぎ）を含めた原子・分子のダイナミクスをシミュレートする。この分子動力学計算として、上記のダイナミクスをシミュレート可能なものであれば、具体的な手法に特に制限はなく目的に応じて適宜選択することができる。
ここで、分子動力学（Ｍｏｌｅｃｕｌａｒ Ｄｙｎａｍｉｃｓ、ＭＤ）とは、ニュートンの運動方程式を数値的に解くことにより、原子などの粒子（質点）の運動をシミュレーションする方法を意味する。
分子動力学法による分子動力学計算（シミュレーション）は、例えば、分子動力学計算プログラムを用いて行うことができる。分子動力学計算プログラムとしては、例えば、ＧＲＯＭＡＣＳ（ｈｔｔｐ：／／ｗｗｗ．ｇｒｏｍａｃｓ．ｏｒｇ／）、ＡＭＢＥＲ（ｈｔｔｐ：／／ａｍｂｅｒｍｄ．ｏｒｇ／）、ＣＨＡＲＭＭ（ｈｔｔｐ：／／ｗｗｗ．ｃｈａｒｍｍ．ｏｒｇ／ｃｈａｒｍｍ／）、ＮＡＭＤ（ｈｔｔｐ：／／ｗｗｗ．ｋｓ．ｕｉｕｃ．ｅｄｕ／Ｒｅｓｅａｒｃｈ／ｎａｍｄ／）、ＬＡＭＭＰＳ（ｈｔｔｐｓ：／／ｗｗｗ．ｌａｍｍｐｓ．ｏｒｇ／）、ＧＥＮＥＳＩＳ（ｈｔｔｐｓ：／／ｗｗｗ．ｒ－ｃｃｓ．ｒｉｋｅｎ．ｊｐ／ｌａｂｓ／ｃｂｒｔ／）などが挙げられる。
分子動力学計算部３４は、結晶構造予測工程において分子動力学計算を実施する。分子動力学計算部３４は、この分子動力学計算において、候補構造（候補構造を積層した構造を含む）における、熱的揺動（熱揺らぎ）を含めた原子・分子のダイナミクスをシミュレートする。この分子動力学計算として、上記のダイナミクスをシミュレート可能なものであれば、具体的な手法に特に制限はなく目的に応じて適宜選択することができる。
ここで、分子動力学（Ｍｏｌｅｃｕｌａｒ Ｄｙｎａｍｉｃｓ、ＭＤ）とは、ニュートンの運動方程式を数値的に解くことにより、原子などの粒子（質点）の運動をシミュレーションする方法を意味する。
分子動力学法による分子動力学計算（シミュレーション）は、例えば、分子動力学計算プログラムを用いて行うことができる。分子動力学計算プログラムとしては、例えば、ＧＲＯＭＡＣＳ（ｈｔｔｐ：／／ｗｗｗ．ｇｒｏｍａｃｓ．ｏｒｇ／）、ＡＭＢＥＲ（ｈｔｔｐ：／／ａｍｂｅｒｍｄ．ｏｒｇ／）、ＣＨＡＲＭＭ（ｈｔｔｐ：／／ｗｗｗ．ｃｈａｒｍｍ．ｏｒｇ／ｃｈａｒｍｍ／）、ＮＡＭＤ（ｈｔｔｐ：／／ｗｗｗ．ｋｓ．ｕｉｕｃ．ｅｄｕ／Ｒｅｓｅａｒｃｈ／ｎａｍｄ／）、ＬＡＭＭＰＳ（ｈｔｔｐｓ：／／ｗｗｗ．ｌａｍｍｐｓ．ｏｒｇ／）、ＧＥＮＥＳＩＳ（ｈｔｔｐｓ：／／ｗｗｗ．ｒ－ｃｃｓ．ｒｉｋｅｎ．ｊｐ／ｌａｂｓ／ｃｂｒｔ／）などが挙げられる。
【0042】
ここで、分子動力学計算においては、分子中に存在する各原子が、分子動力学計算のなかでどのような力を受けているのかを関数として数式化した分子力場を用いることができる。分子力場に基づく分子動力学計算では、原子間にはたらく力を、原子間の相互作用を表すパラメータ（結合距離や結合角など）を変数とし、原子の種類や結合様式によって決まるポテンシャル関数で数値として表す。
本件の技術において、適用することができる分子力場としては、特に制限なく目的に応じて適宜選択することができ、例えば、有機化合物についての汎用分子力場であるＧＡＦＦ（ｇｅｎｅｒａｌ ＡＭＢＥＲ ｆｏｒｃｅ ｆｉｅｌｄ）などを用いることができる。
本件の技術において、適用することができる分子力場としては、特に制限なく目的に応じて適宜選択することができ、例えば、有機化合物についての汎用分子力場であるＧＡＦＦ（ｇｅｎｅｒａｌ ＡＭＢＥＲ ｆｏｒｃｅ ｆｉｅｌｄ）などを用いることができる。
【0043】
それぞれの原子にはたらく力によるエネルギーの種類としては、例えば、結合長のエネルギー、結合角のエネルギー、二面角のエネルギー、ファンデルワールス相互作用エネルギー、静電相互作用エネルギーなどが挙げられる。なお、分子動力学計算においては、通常、分子を構成する全ての原子にはたらく力の総和が「ポテンシャルエネルギー」となる。
また、分子動力学計算において静電相互作用を算出する際には、通常、分子を形成する各原子に点電荷を設定し、その原子における点電荷に基づいて静電相互作用を算出する。分子を形成する各原子に点電荷としては、例えば、ＲＥＳＰ（Ｒｅｓｔｒａｉｎｅｄ Ｅｌｅｃｔｒｏｓｔａｔｉｃ Ｐｏｔｅｎｔｉａｌ）電荷を用いることが好ましい。
上記のように、エネルギー算出部３５は、上記の各種エネルギーに基づく「ポテンシャルエネルギー」を算出する。
また、エネルギー算出部３５は、分子動力学計算における運動エネルギーを、各原子の速度とその質量とに基づいて算出する。
本件の技術においては、候補構造のエネルギーとして、ポテンシャルエネルギーと運動エネルギーの両方を考慮することが好ましい。言い換えると、本件の技術では、候補構造のエネルギーが、候補構造におけるポテンシャルエネルギーと運動エネルギーの和を含むことが好ましい。
こうすることにより、本件の技術では、分子にかかる力（ポテンシャル）と分子のダイナミクス（運動）をより適切に評価することができるため、さらに高精度な評価を行うことができる。例えば、エネルギー算出部３５は、算出された各種エネルギーに基づく「ポテンシャルエネルギー」と「運動エネルギー」を夫々算出し、さらにその和である「全エネルギー」を算出して、これを利用するとよい。
なお、以下では、ポテンシャルエネルギーと運動エネルギーの和を全エネルギーと称することがある。
また、分子動力学計算において静電相互作用を算出する際には、通常、分子を形成する各原子に点電荷を設定し、その原子における点電荷に基づいて静電相互作用を算出する。分子を形成する各原子に点電荷としては、例えば、ＲＥＳＰ（Ｒｅｓｔｒａｉｎｅｄ Ｅｌｅｃｔｒｏｓｔａｔｉｃ Ｐｏｔｅｎｔｉａｌ）電荷を用いることが好ましい。
上記のように、エネルギー算出部３５は、上記の各種エネルギーに基づく「ポテンシャルエネルギー」を算出する。
また、エネルギー算出部３５は、分子動力学計算における運動エネルギーを、各原子の速度とその質量とに基づいて算出する。
本件の技術においては、候補構造のエネルギーとして、ポテンシャルエネルギーと運動エネルギーの両方を考慮することが好ましい。言い換えると、本件の技術では、候補構造のエネルギーが、候補構造におけるポテンシャルエネルギーと運動エネルギーの和を含むことが好ましい。
こうすることにより、本件の技術では、分子にかかる力（ポテンシャル）と分子のダイナミクス（運動）をより適切に評価することができるため、さらに高精度な評価を行うことができる。例えば、エネルギー算出部３５は、算出された各種エネルギーに基づく「ポテンシャルエネルギー」と「運動エネルギー」を夫々算出し、さらにその和である「全エネルギー」を算出して、これを利用するとよい。
なお、以下では、ポテンシャルエネルギーと運動エネルギーの和を全エネルギーと称することがある。
【0044】
本分子動力学計算における「シミュレーション時間」とは、一般に、ニュートンの運動方程式に基づいて離散時間解析することによって、短い刻みの時間ごとに原子の位置、又はその位置の変化を繰り返し計算することにより、その原子を含む分子の動きをシミュレートする時間を意味する。
本件の技術においては、上記のシミュレーション時間における所定の時間範囲を解析の対象として選択し、エネルギー等を算出することが好ましい。つまり、解析の時間範囲として、分子動力学計算によるシミュレーション時間の一部を選択して、選択した時間範囲のデータを用いてエネルギー等を算出することが好ましい。解析の時間範囲としては、例えば、シミュレーション時間における後半（例えば、シミュレーション時間が２０ｎｓ（ナノ秒）であれば、１５ｎｓ～２０ｎｓや１９ｎｓ～２０ｎｓなど）とすることが好ましい。
また、上記の短い刻みの時間（時間刻み幅）は、例えば、０．１ｆｓ（フェムト秒）以上１０ｆｓ以下であることが好ましく、０．５ｆｓ以上２．０ｆｓ以下であることがより好ましい。ここで、時間刻み幅での原子の位置の変化を繰り返し計算する際における繰り返し回数を「ループ回数」とすると、シミュレーション時間は、時間刻み幅とループ回数の積で表される。
本件の技術においては、上記のシミュレーション時間における所定の時間範囲を解析の対象として選択し、エネルギー等を算出することが好ましい。つまり、解析の時間範囲として、分子動力学計算によるシミュレーション時間の一部を選択して、選択した時間範囲のデータを用いてエネルギー等を算出することが好ましい。解析の時間範囲としては、例えば、シミュレーション時間における後半（例えば、シミュレーション時間が２０ｎｓ（ナノ秒）であれば、１５ｎｓ～２０ｎｓや１９ｎｓ～２０ｎｓなど）とすることが好ましい。
また、上記の短い刻みの時間（時間刻み幅）は、例えば、０．１ｆｓ（フェムト秒）以上１０ｆｓ以下であることが好ましく、０．５ｆｓ以上２．０ｆｓ以下であることがより好ましい。ここで、時間刻み幅での原子の位置の変化を繰り返し計算する際における繰り返し回数を「ループ回数」とすると、シミュレーション時間は、時間刻み幅とループ回数の積で表される。
【0045】
【0046】
本件の技術において、シミュレーション時間としては、候補構造の評価を行うために必要となるデータを取得できれば、特に制限はなく目的に応じて適宜選択できる。例えば、実施例に示す結果は、分子動力学計算によるシミュレーション時間を２０ｎｓにして、解析の時間範囲を、１５ｎｓから２０ｎｓにした事例を示す。
【0047】
図９に、抽出された７個の候補構造についてＣＯＮＦＬＥＸの解析結果による順位と分子動力学計算の結果の全エネルギーによる順位を対比して示し、さらに実際に実験的に確認された結果と比較して、有効性が認められた候補構造の「評価」の欄に「〇」印を付与する。図９に示す順位は、図６に示した全エネルギーの大きさに基づくものである。この７個の候補構造を、全エネルギーが小さい順に並べると、候補２５、候補３５、候補９、候補５、候補３、候補７、候補１７の順になる。
なお、Ｘ線結晶構造解析により得られた実際の結晶構造は、ＣＯＮＦＬＥＸにより抽出された候補３と一致する構造であり、本実施形態では候補３が正しい結晶構造（正解の構造）である。実験的に確認された結果と、分子動力学計算の結果から有効性が認められた候補構造とが一致することが確認できた。なお、候補５は候補３とは一部の分子の配向が反転しただけの、互いに類似した結晶構造であり、実験ではこれまでに確認されていないが、実際の結晶構造としては存在しうる可能性がある。
上記のとおり、分子動力学計算の解析結果に基づいて抽出された候補構造には、ＣＯＮＦＬＥＸによって抽出された候補１、候補２が含まれていない。換言すれば、ＣＯＮＦＬＥＸによって抽出された候補１、候補２は、分子動力学計算の結果に基づいて抽出された７個の候補構造よりも注目度が低いものと判断された。ＣＯＮＦＬＥＸの結果によって抽出された結果をそのまま利用する場合、上記の候補１、候補２を、候補３の前に検証することが一般的である。本実施形態の分子動力学計算の解析結果には、上記の候補１、候補２などが含まれていない。ＣＯＮＦＬＥＸの結果の全３５個の候補を、順に検証する場合に比べて解析対象を絞ることができる。
なお、Ｘ線結晶構造解析により得られた実際の結晶構造は、ＣＯＮＦＬＥＸにより抽出された候補３と一致する構造であり、本実施形態では候補３が正しい結晶構造（正解の構造）である。実験的に確認された結果と、分子動力学計算の結果から有効性が認められた候補構造とが一致することが確認できた。なお、候補５は候補３とは一部の分子の配向が反転しただけの、互いに類似した結晶構造であり、実験ではこれまでに確認されていないが、実際の結晶構造としては存在しうる可能性がある。
上記のとおり、分子動力学計算の解析結果に基づいて抽出された候補構造には、ＣＯＮＦＬＥＸによって抽出された候補１、候補２が含まれていない。換言すれば、ＣＯＮＦＬＥＸによって抽出された候補１、候補２は、分子動力学計算の結果に基づいて抽出された７個の候補構造よりも注目度が低いものと判断された。ＣＯＮＦＬＥＸの結果によって抽出された結果をそのまま利用する場合、上記の候補１、候補２を、候補３の前に検証することが一般的である。本実施形態の分子動力学計算の解析結果には、上記の候補１、候補２などが含まれていない。ＣＯＮＦＬＥＸの結果の全３５個の候補を、順に検証する場合に比べて解析対象を絞ることができる。
【0048】
ここで、分子動力学計算で設定するアンサンブルとしては、特に制限はなく目的に応じて適宜選択できる。
本件の技術においては、より安定な結果が得られるように分子動力学計算のシミュレーションを実施するために、例えば、粒子数、体積、温度一定（ＮＶＴアンサンブル）の分子動力学計算による構造緩和を行った後に、粒子数、圧力、温度一定（ＮＰＴアンサンブル）の分子動力学計算による計算系の平衡化を行うことが好ましい。
本件の技術においては、より安定な結果が得られるように分子動力学計算のシミュレーションを実施するために、例えば、粒子数、体積、温度一定（ＮＶＴアンサンブル）の分子動力学計算による構造緩和を行った後に、粒子数、圧力、温度一定（ＮＰＴアンサンブル）の分子動力学計算による計算系の平衡化を行うことが好ましい。
【0049】
また、分子動力学計算における設定温度としては、特に制限はなく目的に応じて適宜選択でき、例えば、対象分子が結晶構造をとりうる温度とすることが好ましい。言い換えると、本件の技術では、結晶化可能な対象分子について解析する。その解析において、対象分子が結晶化しうる温度を設定して分子動力学計算を行うことが好ましい。
ここで、分子動力学計算で設定する所定の温度としては、上述した対象分子が結晶構造をとりうる温度のほかに、例えば、室温（３００Ｋなど）、対象分子と類似する分子が結晶化する温度、対象分子を材料として用いるときに想定される温度（デバイスを動作させるときの温度）などとすることができる。
こうすることにより、本件の技術では、対象分子による結晶構造におけるリアリティの高い熱的揺動をシミュレートすることができるため、より適切に熱的揺動の影響を考慮して、さらに高精度に結晶構造を予測することができる。
ここで、分子動力学計算で設定する所定の温度としては、上述した対象分子が結晶構造をとりうる温度のほかに、例えば、室温（３００Ｋなど）、対象分子と類似する分子が結晶化する温度、対象分子を材料として用いるときに想定される温度（デバイスを動作させるときの温度）などとすることができる。
こうすることにより、本件の技術では、対象分子による結晶構造におけるリアリティの高い熱的揺動をシミュレートすることができるため、より適切に熱的揺動の影響を考慮して、さらに高精度に結晶構造を予測することができる。
【0050】
＜＜候補構造の評価＞＞
結晶構造予測部３６は、結晶構造予測工程のなかで、分子動力学計算に基づく評価を行い結晶構造を予測する。
本件の技術において、分子動力学計算に基づいて候補構造の評価を行う手法としては、分子動力学計算によって得られたデータ（トラジェクトリ、エネルギーなど）に基づく評価であれば、特に制限はなく目的に応じて適宜選択できる。例えば、候補構造の評価は、候補構造（候補構造を積層した構造を含む）におけるエネルギーに基づく評価であることが好ましい。言い換えると、本件の技術の候補構造の評価の手法では、分子動力学計算から得た候補構造のエネルギーに基づいて評価を行うことが好ましい。
こうすることにより、本件の技術では、候補構造における分子の熱的揺動の影響を考慮して、候補構造のエネルギー的な安定性を評価できるため、より高精度な評価を行うことができる。候補構造のエネルギーを、単位量当たりの熱量（ｋＪ／ｍｏｌ）を単位にして示すことにより、エネルギー的な安定性を評価することが可能になる。
結晶構造予測部３６は、結晶構造予測工程のなかで、分子動力学計算に基づく評価を行い結晶構造を予測する。
本件の技術において、分子動力学計算に基づいて候補構造の評価を行う手法としては、分子動力学計算によって得られたデータ（トラジェクトリ、エネルギーなど）に基づく評価であれば、特に制限はなく目的に応じて適宜選択できる。例えば、候補構造の評価は、候補構造（候補構造を積層した構造を含む）におけるエネルギーに基づく評価であることが好ましい。言い換えると、本件の技術の候補構造の評価の手法では、分子動力学計算から得た候補構造のエネルギーに基づいて評価を行うことが好ましい。
こうすることにより、本件の技術では、候補構造における分子の熱的揺動の影響を考慮して、候補構造のエネルギー的な安定性を評価できるため、より高精度な評価を行うことができる。候補構造のエネルギーを、単位量当たりの熱量（ｋＪ／ｍｏｌ）を単位にして示すことにより、エネルギー的な安定性を評価することが可能になる。
【0051】
また、候補構造のエネルギーを算出する手法としては、特に制限はなく目的に応じて適宜選択でき、例えば、分子動力学計算プログラムに付随した解析ツールなどを用いて行うことができる。
【0052】
本件の技術において、対象分子の結晶構造を予測する手法としては、候補構造についての分子動力学計算による評価に基づいた手法であれば、特に制限はなく目的に応じて適宜選択できる。例えば、所定の評価基準に従い、候補構造のエネルギー（全エネルギー）の大きさ（評価値）の順に順位付け（ランキング）する手法などを用いることができる。この順位付けの手法によれば、正しい結晶構造である可能性が高い順に順位付ける（ランキング）ことができる。
評価値に応じて、正しい結晶構造である可能性が高い順に順位付けする手法としては、特に制限はなく目的に応じて適宜選択できる。例えば、結晶構造の評価値を低い順に順位付けして整理すれば、その評価値が低いほど正しい結晶構造が抽出されたとみなすことができる。この候補構造の順位を、結晶構造の予測の優先順位とすれば、より正しい結晶構造を抽出できる。
評価値に応じて、正しい結晶構造である可能性が高い順に順位付けする手法としては、特に制限はなく目的に応じて適宜選択できる。例えば、結晶構造の評価値を低い順に順位付けして整理すれば、その評価値が低いほど正しい結晶構造が抽出されたとみなすことができる。この候補構造の順位を、結晶構造の予測の優先順位とすれば、より正しい結晶構造を抽出できる。
【0053】
例えば、一般にエネルギーが低い程安定な構造であると考えられる。本件の技術において、候補構造のエネルギーに基づいて評価を行う場合に、より安定な構造を求めるならば、エネルギーが低い候補構造の評価が高くなることが好ましい。言い換えると、本件の技術においては、候補構造のエネルギーが低い程、評価の結果が良いものとして結晶構造を予測することが好ましい。
より具体的には、本件の技術では、候補構造についての分子動力学計算から算出したエネルギーが低い順に、正しい結晶構造である可能性が高いものとして順位付ける。
この順位付けに従えば、エネルギー的に安定である順に、上位に順位付けされることになるため、正しい結晶構造を高精度に予測することができる。これにより、この順位に従って選択された結晶構造が、より安定な結晶構造を有していると予測できる。
より具体的には、本件の技術では、候補構造についての分子動力学計算から算出したエネルギーが低い順に、正しい結晶構造である可能性が高いものとして順位付ける。
この順位付けに従えば、エネルギー的に安定である順に、上位に順位付けされることになるため、正しい結晶構造を高精度に予測することができる。これにより、この順位に従って選択された結晶構造が、より安定な結晶構造を有していると予測できる。
【0054】
また、本件の技術においては、エネルギーによる候補構造の評価だけではなく、他の指標を用いて候補構造を評価してよい。他の指標としては、候補構造の安定性を反映した指標であれば、特に制限はなく目的に応じて適宜選択することができ、例えば、この指標に温度因子（Ｂ－ｆａｃｔｏｒ）を用いることができる。
温度因子は、分子における熱的揺動の大きさを定量的に表した指標である。結晶構造においては、各分子の温度因子の値が小さい方が熱的安定性が高いと考えることができる。このため、本件の技術では、例えば、分子動力学計算の結果を解析して算出した温度因子の値が小さい候補構造を、正しい結晶構造である可能性が高い候補構造として予測してもよい。さらに、本件の技術では、エネルギーによる評価と温度因子による評価を両方とも行ってもよい。
なお、温度因子は、例えば、分子動力学計算を行うプログラムに付随した解析プログラムなどの公知のプログラムを用いて求めることができる。
温度因子は、分子における熱的揺動の大きさを定量的に表した指標である。結晶構造においては、各分子の温度因子の値が小さい方が熱的安定性が高いと考えることができる。このため、本件の技術では、例えば、分子動力学計算の結果を解析して算出した温度因子の値が小さい候補構造を、正しい結晶構造である可能性が高い候補構造として予測してもよい。さらに、本件の技術では、エネルギーによる評価と温度因子による評価を両方とも行ってもよい。
なお、温度因子は、例えば、分子動力学計算を行うプログラムに付随した解析プログラムなどの公知のプログラムを用いて求めることができる。
【0055】
＜その他の工程＞
その他の工程としては、特に制限はなく、目的に応じて適宜選択することができる。
その他の工程としては、特に制限はなく、目的に応じて適宜選択することができる。
【0056】
（結晶構造予測プログラム）
本件で開示する技術の一例としての結晶構造予測プログラムは、上述したように、候補構造特定処理と、結晶構造予測処理とをコンピュータに行わせ、必要に応じてその他の処理をコンピュータに行わせるためのプログラムである。
本件で開示する技術の一例としての結晶構造予測プログラムは、上述したように、候補構造特定処理と、結晶構造予測処理とをコンピュータに行わせ、必要に応じてその他の処理をコンピュータに行わせるためのプログラムである。
【0057】
本件の結晶構造予測プログラムは、例えば、本件の結晶構造予測方法をコンピュータに実行させるプログラムとすることができる。また、本件の結晶構造予測プログラムにおける好適な態様は、例えば、本件の結晶構造予測方法における好適な態様と同様にすることができる。
【0058】
本件の結晶構造予測プログラムは、使用するコンピュータシステムの構成及びオペレーティングシステムの種類・バージョンなどに応じて、公知の各種のプログラム言語を用いて作成することができる。
【0059】
本件の結晶構造予測プログラムは、内蔵ハードディスク、外付けハードディスクなどの記録媒体に記録しておいてもよいし、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭ、ＭＯディスク、ＵＳＢメモリなどの記録媒体に記録しておいてもよい。
さらに、本件の結晶構造予測プログラムを、上記の記録媒体に記録する場合には、必要に応じて、コンピュータシステムが有する記録媒体読取装置を通じて、これを直接又はハードディスクにインストールして使用することができる。また、コンピュータシステムから情報通信ネットワークを通じてアクセス可能な外部記憶領域（他のコンピュータ、クラウド・コンピューティング（クラウド）における記憶装置など）に本件の結晶構造予測プログラムを記録しておいてもよい。この場合、外部記憶領域に記録された本件の結晶構造予測プログラムは、必要に応じて、外部記憶領域から情報通信ネットワークを通じてこれを直接、又はハードディスクにインストールして使用することができる。
なお、本件の結晶構造予測プログラムは、複数の記録媒体に、任意の処理ごとに分割されて記録されていてもよい。
さらに、本件の結晶構造予測プログラムを、上記の記録媒体に記録する場合には、必要に応じて、コンピュータシステムが有する記録媒体読取装置を通じて、これを直接又はハードディスクにインストールして使用することができる。また、コンピュータシステムから情報通信ネットワークを通じてアクセス可能な外部記憶領域（他のコンピュータ、クラウド・コンピューティング（クラウド）における記憶装置など）に本件の結晶構造予測プログラムを記録しておいてもよい。この場合、外部記憶領域に記録された本件の結晶構造予測プログラムは、必要に応じて、外部記憶領域から情報通信ネットワークを通じてこれを直接、又はハードディスクにインストールして使用することができる。
なお、本件の結晶構造予測プログラムは、複数の記録媒体に、任意の処理ごとに分割されて記録されていてもよい。
【0060】
【0061】
まず、Ｓ１０１では、処理部３の分子構造作成部３１は、分子モデリングソフトにより、対象分子の分子構造を作成する。言い換えると、Ｓ１０１においては、例えば、処理部３の分子構造作成部３１は、分子モデリングソフトを用いて、結晶構造を予測する対象となる対象分子の分子構造の構造データを作成する。
ここで、分子構造の構造データ作成に用いる分子モデリングソフトとしては、特に制限はなく目的に応じて適宜選択することができ、例えば、Ｗｉｎｍｏｓｔａｒ（登録商標）などを用いることができる。また、作成する構造データのフォーマットとしては、特に制限はなく目的に応じて適宜選択することができ、例えば、ｍｏｌファイルなどを用いることができる。
ここで、分子構造の構造データ作成に用いる分子モデリングソフトとしては、特に制限はなく目的に応じて適宜選択することができ、例えば、Ｗｉｎｍｏｓｔａｒ（登録商標）などを用いることができる。また、作成する構造データのフォーマットとしては、特に制限はなく目的に応じて適宜選択することができ、例えば、ｍｏｌファイルなどを用いることができる。
【0062】
続いて、Ｓ１０２では、処理部３の候補構造特定手段３１０（分子構造作成部３１）は、対象分子の分子構造に基づいて結晶構造の候補を探索して候補構造を特定する。より具体的には、Ｓ１０２においては、例えば、Ｓ１０１で作成した対象分子の分子構造の構造データを用いて、結晶構造探索プログラムによる結晶構造探索を行うことにより、結晶構造の候補となる候補構造を複数生成する。ここで、結晶構造探索プログラムとしては、上述したように、例えば、ＣＯＮＦＬＥＸ（登録商標）を用いることができる。
また、本件の技術においては、Ｓ１０２での結晶構造探索プログラムによる結晶構造探索の前に、結晶構造探索プログラム等を用いて、Ｓ１０１で作成した分子構造に対する、１分子での構造最適化を行ってもよい。
また、本件の技術においては、Ｓ１０２での結晶構造探索プログラムによる結晶構造探索の前に、結晶構造探索プログラム等を用いて、Ｓ１０１で作成した分子構造に対する、１分子での構造最適化を行ってもよい。
【0063】
次に、Ｓ１０３では、処理部３の候補構造特定手段３１０（初期結晶構造生成部３３）は、特定した候補構造を積層して、分子動力学計算の初期構造を作成する。より具体的には、Ｓ１０３においては、例えば、候補構造特定手段３１０は、Ｓ１０２で特定した候補構造（結晶構造における単位構造）の構造データのそれぞれに対して、３次元的に積層する処理を行うことにより、数百分子を含む初期構造を作成して、初期構造情報７４に追加する。なお、Ｓ１０３で行う処理は、適宜公知のプログラムを用いて行うことができる。
【0064】
そして、Ｓ１０４では、処理部３の結晶構造予測手段３２０（分子動力学計算部３４）は、作成した初期構造を用いて、分子動力学計算を実行する。より具体的には、Ｓ１０４においては、例えば、結晶構造予測手段３２０は、Ｓ１０３で作成され初期構造情報７４に格納されている初期構造を用いて、対象分子が結晶化しうる温度（例えば、室温）を設定した条件で、比較的短時間の構造緩和計算（ＮＶＴアンサンブル）を行った後、比較的長時間の平衡化計算（ＮＰＴアンサンブル）を実行する。結晶構造予測手段３２０は、分子動力学計算の結果を分子動力学解析結果７５に追加する。
【0065】
なお、分子動力学計算を実行する際には、結晶構造予測手段３２０は、例えば、力場パラメータの設定、分子を形成する各原子における点電荷の設定、エネルギー極小化計算などの事前の処理を行うことが通常であるため、必要に応じてこれらの処理を行うことが好ましい。
【0066】
続いて、Ｓ１０５では、処理部３の結晶構造予測手段３２０（エネルギー算出部３５）は、分子動力学計算の結果を用いて、候補構造におけるエネルギーを算出する。より具体的には、Ｓ１０５においては、例えば、結晶構造予測手段３２０は、Ｓ１０４で実行した分子動力学計算の計算結果（分子動力学解析結果７５）のデータを解析することにより、候補構造から作成したそれぞれの計算系についての全エネルギー（ポテンシャルエネルギーと運動エネルギーの和）を算出する。
なお、全エネルギーを算出するために用いる計算結果のデータの範囲としては、例えば、上述した平衡化計算におけるシミュレーション時間の後半（平衡状態に達していると考えられる時間範囲）のデータを用いることが好ましい。より好ましくは、シミュレーション時間の終了に比較的近い所定の期間を含むように上記の範囲を予め決定するとよい。
なお、全エネルギーを算出するために用いる計算結果のデータの範囲としては、例えば、上述した平衡化計算におけるシミュレーション時間の後半（平衡状態に達していると考えられる時間範囲）のデータを用いることが好ましい。より好ましくは、シミュレーション時間の終了に比較的近い所定の期間を含むように上記の範囲を予め決定するとよい。
【0067】
次に、Ｓ１０６では、処理部３の結晶構造予測手段３２０（結晶構造予測部３６）は、算出したエネルギーの値に基づいて、候補構造を順位付けする。より具体的には、Ｓ１０６においては、例えば、結晶構造予測手段３２０は、Ｓ１０５で算出したエネルギーの値が低い順に、評価の結果が良いものとして候補構造を順位付けることにより、候補構造を評価するとよい。
【0068】
そして、Ｓ１０７では、結晶構造予測手段３２０（結晶構造予測部３６）は、エネルギーが低い候補構造を正しい結晶構造である可能性が高い構造として予測すると、処理を終了させる。より具体的には、Ｓ１０７においては、例えば、結晶構造予測手段３２０は、Ｓ１０６で行った順位付けの結果に基づいて、エネルギーが低い順に所定数の候補構造を選択し、選択された所定数の候補構造が正しい結晶構造である可能性が高い構造であるとして予測（特定）すると、処理を終了させる。
【0069】
また、図１０においては、本件の技術の一例における処理の流れについて、特定の順序に従って説明したが、本件の技術においては、技術的に可能な範囲で、適宜各ステップの順序を入れ替えることができる。また、本件の技術においては、技術的に可能な範囲で、複数のステップを一括して行ってもよい。
【0070】
以上、説明したように、本件の結晶構造予測方法は、対象分子の分子構造に基づいて、対象分子の結晶構造の候補となる候補構造を特定し、候補構造を用いた分子動力学計算を行うことにより候補構造の評価を行い、評価に基づいて対象分子の結晶構造を予測する。こうすることにより、本件の結晶構造予測方法は、対象分子の分子構造の情報から、当該対象分子により形成される結晶構造を、結晶構造における分子の熱的揺動の影響を考慮して、より高精度に予測できる。
【0071】
（第２の実施形態）
図１１から図１３を参照して、第２の実施形態について説明する。第１の実施形態では、具体的な対象分子としてＰｈＣ２－ＢＱＱＤＩを適用した事例について説明した。対象分子の種類はＰｈＣ２－ＢＱＱＤＩに制限されない。第２の実施形態では、これに代わる別の対象分子としてＰｈＣ２－ＰＤＩを適用する事例について説明する。対象分子の中心骨格（コア）がperylenetetracarboxylic diimide （ＰＤＩ）に代わる。
図１１から図１３を参照して、第２の実施形態について説明する。第１の実施形態では、具体的な対象分子としてＰｈＣ２－ＢＱＱＤＩを適用した事例について説明した。対象分子の種類はＰｈＣ２－ＢＱＱＤＩに制限されない。第２の実施形態では、これに代わる別の対象分子としてＰｈＣ２－ＰＤＩを適用する事例について説明する。対象分子の中心骨格（コア）がperylenetetracarboxylic diimide （ＰＤＩ）に代わる。
【0072】
図１１は、実施形態の対象分子の構造式を示す図である。図１２は、実施形態のＣＯＮＦＬＥＸによる抽出結果と、分子動力学計算における全エネルギーの対応を説明するための図である。図１３は、実施形態のＣＯＮＦＬＥＸによる抽出結果と、分子動力学計算の解析結果を説明するための図である。なお、図１１が前述の図５に対応し、図１２が前述の図６に対応し、図１３が前述の図９に対応する。
【0073】
第２の実施形態における対象分子（ＰｈＣ２－ＰＤＩ）の場合も第１の実施形態と同様に、候補構造特定工程において抽出する候補構造の個数を、例えば、３５個にする。
【0074】
図１３に、抽出された８個の候補構造についてＣＯＮＦＬＥＸの解析結果による順位と分子動力学計算の結果の全エネルギーによる順位を対比して示し、さらに実際に実験的に確認された結果と比較して、有効性が認められた候補構造の「評価」の欄に「〇」印を付与する。図１３に示す順位は、図１２に示した全エネルギーの大きさに基づくものである。この８個の候補構造を、全エネルギーが小さい順に並べると、候補１、候補３４、候補１０、候補２、候補９、候補１９、候補２５、候補５の順になる。
なお、実験による検証の結果から確認された結晶構造は、候補３４と一致する構造であり、本実施形態では候補３４が正しい結晶構造（正解の構造）である。このように、候補３４は、分子動力学計算の結果が示す候補構造の比較的上位にあり、実際に実験的に確認された結果と一致した。
比較例としてのＣＯＮＦＬＥＸによって抽出された候補の順に検証を行うと、３４番目にようやくその対象と判断される場合が想定される。
これに対して、分子動力学計算の結果を用いれば、候補３４は、候補１に次いで２番目になる。このため、分子動力学計算の結果の順に従えば、高精度に結晶構造を予測することができる。本実施形態の分子動力学計算の解析結果を用いることで、ＣＯＮＦＬＥＸの結果の全３５個の候補を、順に検証する場合に比べて解析対象を絞ることのほか、解析を実施する順序も最適化できる。
なお、実験による検証の結果から確認された結晶構造は、候補３４と一致する構造であり、本実施形態では候補３４が正しい結晶構造（正解の構造）である。このように、候補３４は、分子動力学計算の結果が示す候補構造の比較的上位にあり、実際に実験的に確認された結果と一致した。
比較例としてのＣＯＮＦＬＥＸによって抽出された候補の順に検証を行うと、３４番目にようやくその対象と判断される場合が想定される。
これに対して、分子動力学計算の結果を用いれば、候補３４は、候補１に次いで２番目になる。このため、分子動力学計算の結果の順に従えば、高精度に結晶構造を予測することができる。本実施形態の分子動力学計算の解析結果を用いることで、ＣＯＮＦＬＥＸの結果の全３５個の候補を、順に検証する場合に比べて解析対象を絞ることのほか、解析を実施する順序も最適化できる。
【0075】
第２の実施形態によれば、第１の実施形態と同様に、分子動力学計算の解析結果を用いることで、候補構造のなかからより適した対象分子を抽出できる。
【0076】
上記のように本実施形態によれば、対象分子の種類が代わっても、第１の実施形態と同様の効果を奏する。
【0077】
以上のとおり、本開示に係るいくつかの実施形態を説明したが、これら全ての実施形態は、例として提示したものであり、発明の範囲を限定することを意図していない。これらの実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で種々の省略、置き換え、変更を行うことができる。これらの実施形態及びその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【0078】
３ 処理部
７ 記憶部
１０ 結晶構造予測装置
３１ 分子構造作成部
３２ 結晶構造生成部
３３ 初期結晶構造生成部
３４ 分子動力学計算部
３５ エネルギー算出部
３６ 結晶構造予測部
７２ 分子構造情報
７３ 候補構造情報
７４ 初期構造情報
７５ 分子動力学解析結果
３１０ 候補構造特定手段
３２０ 結晶構造予測手段
７ 記憶部
１０ 結晶構造予測装置
３１ 分子構造作成部
３２ 結晶構造生成部
３３ 初期結晶構造生成部
３４ 分子動力学計算部
３５ エネルギー算出部
３６ 結晶構造予測部
７２ 分子構造情報
７３ 候補構造情報
７４ 初期構造情報
７５ 分子動力学解析結果
３１０ 候補構造特定手段
３２０ 結晶構造予測手段
【図1】
【図2】
【図3】
【図4A】
【図4B】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【手続補正書】
【提出日】2022-09-22
【手続補正1】
【補正対象書類名】図面
【補正対象項目名】図３
【補正方法】変更
【補正の内容】
【図3】
【手続補正2】
【補正対象書類名】図面
【補正対象項目名】図４Ａ
【補正方法】変更
【補正の内容】
【図4A】
【手続補正3】
【補正対象書類名】図面
【補正対象項目名】図４Ｂ
【補正方法】変更
【補正の内容】
【図4B】
【手続補正4】
【補正対象書類名】図面
【補正対象項目名】図７
【補正方法】変更
【補正の内容】
【図7】
【手続補正5】
【補正対象書類名】図面
【補正対象項目名】図８
【補正方法】変更
【補正の内容】
【図8】
