JP2006282929A - 分子構造予測方法 - Google Patents
分子構造予測方法 Download PDFInfo
- Publication number
- JP2006282929A JP2006282929A JP2005107361A JP2005107361A JP2006282929A JP 2006282929 A JP2006282929 A JP 2006282929A JP 2005107361 A JP2005107361 A JP 2005107361A JP 2005107361 A JP2005107361 A JP 2005107361A JP 2006282929 A JP2006282929 A JP 2006282929A
- Authority
- JP
- Japan
- Prior art keywords
- molecular
- molecular structure
- calculation
- trial
- fragments
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 40
- 238000000329 molecular dynamics simulation Methods 0.000 claims abstract description 30
- 239000000126 substance Substances 0.000 claims abstract description 5
- 229920001187 thermosetting polymer Polymers 0.000 claims description 33
- 230000000704 physical effect Effects 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 11
- 239000000470 constituent Substances 0.000 claims description 9
- 239000012634 fragment Substances 0.000 abstract description 39
- 229920000642 polymer Polymers 0.000 abstract description 17
- 230000027455 binding Effects 0.000 abstract 1
- 238000009739 binding Methods 0.000 abstract 1
- 238000004519 manufacturing process Methods 0.000 abstract 1
- 125000004429 atom Chemical group 0.000 description 14
- KPUWHANPEXNPJT-UHFFFAOYSA-N disiloxane Chemical class [SiH3]O[SiH3] KPUWHANPEXNPJT-UHFFFAOYSA-N 0.000 description 7
- 239000000203 mixture Substances 0.000 description 7
- 238000000862 absorption spectrum Methods 0.000 description 4
- 238000012900 molecular simulation Methods 0.000 description 4
- 238000000342 Monte Carlo simulation Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 125000005375 organosiloxane group Chemical group 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 3
- LIVNPJMFVYWSIS-UHFFFAOYSA-N silicon monoxide Inorganic materials [Si-]#[O+] LIVNPJMFVYWSIS-UHFFFAOYSA-N 0.000 description 3
- 238000002910 structure generation Methods 0.000 description 3
- 229910018540 Si C Inorganic materials 0.000 description 2
- 229910018557 Si O Inorganic materials 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 125000004430 oxygen atom Chemical group O* 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 229910010271 silicon carbide Inorganic materials 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- WZJUBBHODHNQPW-UHFFFAOYSA-N 2,4,6,8-tetramethyl-1,3,5,7,2$l^{3},4$l^{3},6$l^{3},8$l^{3}-tetraoxatetrasilocane Chemical compound C[Si]1O[Si](C)O[Si](C)O[Si](C)O1 WZJUBBHODHNQPW-UHFFFAOYSA-N 0.000 description 1
- 229910004283 SiO 4 Inorganic materials 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 238000004132 cross linking Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229920002521 macromolecule Polymers 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000005268 plasma chemical vapour deposition Methods 0.000 description 1
- 239000002861 polymer material Substances 0.000 description 1
- -1 polysiloxane Polymers 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 238000003077 quantum chemistry computational method Methods 0.000 description 1
- 230000005610 quantum mechanics Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Other Resins Obtained By Reactions Not Involving Carbon-To-Carbon Unsaturated Bonds (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
【課題】試行錯誤的探索によりアモルファス網状高分子モデリングを行う際、適切な分子構造モデルが迅速に得られるよう試行構造を生成する方法を提供する
【解決手段】高分子のフラグメントの個々の構造最適化を行うステップＳ２と分子動力学法によりフラグメントの乱雑な空間的配置を生成するステップＳ３とこの空間的配置を構成するフラグメントを暫定的に化学結合で連結するステップＳ４−１とフラグメントの連結結合のうちの２組の間で連結原子対を交換して絡み合いを持たないようにネットワークを変更するステップＳ４−２とネットワークに対して構造最適化と分子動力学計算を行うステップＳ５と上記ステップで得られる試行構造から、微視的な構造パラメータが適切な値を取っているものを分子構造モデルとして選択するステップとＳ６と実験的物性値を再現することを条件に適切な分子構造モデルを選択するステップを備える。
【選択図】図１
【解決手段】高分子のフラグメントの個々の構造最適化を行うステップＳ２と分子動力学法によりフラグメントの乱雑な空間的配置を生成するステップＳ３とこの空間的配置を構成するフラグメントを暫定的に化学結合で連結するステップＳ４−１とフラグメントの連結結合のうちの２組の間で連結原子対を交換して絡み合いを持たないようにネットワークを変更するステップＳ４−２とネットワークに対して構造最適化と分子動力学計算を行うステップＳ５と上記ステップで得られる試行構造から、微視的な構造パラメータが適切な値を取っているものを分子構造モデルとして選択するステップとＳ６と実験的物性値を再現することを条件に適切な分子構造モデルを選択するステップを備える。
【選択図】図１
Description
本発明は、アモルファス網状高分子のモデリングのための分子構造予測方法に関する。ここでモデリングとは、現実の網状高分子の特徴を有する分子構造モデルを予測することを意味する。
近年、高分子の物性研究の手段として、分子動力学法、モンテカルロ法、量子化学計算法等の分子シミュレーションの手法がしばしば用いられるようになっている。
これらの手法は原子、電子の運動を古典力学、量子力学の方程式に従って計算する方法で、物質の性質を分子レベルに遡って求めるために使用される手法である。
これらの手法は原子、電子の運動を古典力学、量子力学の方程式に従って計算する方法で、物質の性質を分子レベルに遡って求めるために使用される手法である。
シミュレーションによって、着目している物性値が正しく算出されるためには、適切な分子構造モデルに対して計算がなされることが必要である。
ここで、適切な分子構造モデルとは、現実の高分子の特徴を有する分子構造のことを意味し、具体的には原子や原子団の組成および密度等の巨視的な構造パラメータが実物とほぼ等しく、結合長、結合角、二面体角、非結合原子間距離等の微視的な構造パラメータが実物のものに近い値をとる構造のことを意味する。以下では、この意味で適切な分子構造モデルを予測することをモデリングと言う。
ここで、適切な分子構造モデルとは、現実の高分子の特徴を有する分子構造のことを意味し、具体的には原子や原子団の組成および密度等の巨視的な構造パラメータが実物とほぼ等しく、結合長、結合角、二面体角、非結合原子間距離等の微視的な構造パラメータが実物のものに近い値をとる構造のことを意味する。以下では、この意味で適切な分子構造モデルを予測することをモデリングと言う。
高分子の興味ある物性の多くはアモルファス相の物性であるが、線状高分子については、アモルファス構造のモデリングは比較的簡単で、これまで多くシミュレーションが行われてきた。
一方、網状高分子に関してはアモルファス構造の適当なモデリングの方法が無いために、これまで分子シミュレーションが実施されていない。
一方、網状高分子に関してはアモルファス構造の適当なモデリングの方法が無いために、これまで分子シミュレーションが実施されていない。
網状高分子のアモルファス構造のモデリングができないのは、網状高分子が剛直であるために、雑に作成した出発構造を分子動力学シミュレーション等を用いて適切な分子構造へ構造緩和させる方法（この方法は、線状高分子をはじめとする複雑分子や、分子集合体のモデリングに常套的に使用される方法である。）が効果的でないためである。
したがって、人為的にある程度適切なアモルファス網状高分子の試行構造を作成することが必要で、このための手法が今のところ開発されていない。
なお、高分子材料の形態予測に関しては、以下の文献がある。
特開２００３−１０５０９０号公報
なお、高分子材料の形態予測に関しては、以下の文献がある。
網状高分子においては、図９に示すように、網目を構成する直鎖部分が他の直鎖部分と多重的に交差して端点Ａを共有しているので、網状高分子の適当な分子構造モデルを作成する際には、端点共有の制約の下で結合や結合角が適切な値を取る分子構造を求めることが必要となる。従って、アモルファス網状高分子のモデリングは、試行構造の生成と取捨選択のステップを経るものとなる。
また、アモルファス網状高分子では、原子や原子団の組成だけが判明していて、それらの結合の仕方(トポロジー)が具体的に判明していないのが通常であるため、モデリングにおいては可能なトポロジーをも試行錯誤的に探索する必要がある。
このような試行探索において、モデリングの対象とする高分子の架橋の度合いが高い場合、失敗試行の割合が多くなるために相当数の試行構造を作成する必要がある。
このような試行探索において、モデリングの対象とする高分子の架橋の度合いが高い場合、失敗試行の割合が多くなるために相当数の試行構造を作成する必要がある。
本発明の課題は、上記の試行錯誤的探索によりアモルファス網状高分子モデリングを行う際、適切な分子構造モデルが迅速に得られるように試行構造を生成する方法を提供することである。また、この試行構造生成法を利用したアモルファス網状高分子のモデリング方法を提供することである。
前記課題を達成するために、
請求項１記載の発明は、次に挙げる７個のステップを含む方法である。
１)アモルファス網状高分子の構成原子団（以下、フラグメントと言うことがある。）の個々の構造最適化を行うステップと、
２)分子動力学法あるいはモンテカルロ法を使用してフラグメントの乱雑な空間的配置を生成するステップと、
３)与えられた空間的配置を構成するフラグメントを暫定的に化学結合により連結するステップと、
４)フラグメントの連結結合にうちの２組の間で連結原子対を交換することによって、絡み合いを持たないように連結結合のネットワークを変更するステップと、
５)出来上がったネットワークの分子構造に対して構造最適化計算と分子動力学計算およびモンテカルロ計算を行うステップと、
６）上記５つのステップによって得られるアモルファス網状高分子の試行構造から、結合長、結合角などの微視的な構造パラメータが適切な値を取っているものを分子構造モデルとして選択するステップと、
７）実験的に得られる物性値を再現することを条件に、アモルファス網状高分子の分子構造モデルとしてより適切な分子構造モデルを選択するステップを備えるものである。
請求項１記載の発明は、次に挙げる７個のステップを含む方法である。
１)アモルファス網状高分子の構成原子団（以下、フラグメントと言うことがある。）の個々の構造最適化を行うステップと、
２)分子動力学法あるいはモンテカルロ法を使用してフラグメントの乱雑な空間的配置を生成するステップと、
３)与えられた空間的配置を構成するフラグメントを暫定的に化学結合により連結するステップと、
４)フラグメントの連結結合にうちの２組の間で連結原子対を交換することによって、絡み合いを持たないように連結結合のネットワークを変更するステップと、
５)出来上がったネットワークの分子構造に対して構造最適化計算と分子動力学計算およびモンテカルロ計算を行うステップと、
６）上記５つのステップによって得られるアモルファス網状高分子の試行構造から、結合長、結合角などの微視的な構造パラメータが適切な値を取っているものを分子構造モデルとして選択するステップと、
７）実験的に得られる物性値を再現することを条件に、アモルファス網状高分子の分子構造モデルとしてより適切な分子構造モデルを選択するステップを備えるものである。
本発明の特徴は、生成する試行構造が高分子の分子シミュレーションに使用されることを想定して、周期境界条件を満足するような構造を生成できることである。これに関係して、ステップ４) の分子動力学法と構造最適化計算は、格子定数も変数とする分子動力学法と構造最適化計算である。
本発明では、考察対象のアモルファス網状高分子をフラグメントへ分割するステップを含んでおらず、適当に分割されたフラグメントが与えられていることを前提とする。
網状高分子をフラグメントへ分割する仕方は、モデリング対象の網状高分子の可能なネットワークトポロジーが網羅されるという条件で、これ以上分割する必要のない原子団のまとまりとすれば良い。
網状高分子をフラグメントへ分割する仕方は、モデリング対象の網状高分子の可能なネットワークトポロジーが網羅されるという条件で、これ以上分割する必要のない原子団のまとまりとすれば良い。
本発明では、与えられたネットワークからネットワークに適合する分子構造を生成する試行を回避するために、はじめにフラグメントの空間配置を適当に与え、次にこれに対して適当なネットワークを作成する手順を取っている。
以下に、ステップ１）ないしステップ７）の意味、目的を簡単に説明する。
ステップ１）では、アモルファス網状高分子のフラグメントの個々の構造最適化を行うが、これは高分子中における構成原子団の構造が高分子の種類によらず
ほぼ一定で、孤立の構成原子団の構造が大部分保持さているからである。
ステップ１）では、アモルファス網状高分子のフラグメントの個々の構造最適化を行うが、これは高分子中における構成原子団の構造が高分子の種類によらず
ほぼ一定で、孤立の構成原子団の構造が大部分保持さているからである。
ステップ２）では、分子動力学法またはモンテカルロ法を使用してフラグメントのユニットセルにおける乱雑な空間的配置を生成するが、ここで分子動力学法やモンテカルロ法を使用するのは、構成原子団が相互に大きな重なりを持たない空間的配置を生成するためで、現実の高分子においては確かに構成原子団が相互に大きく重なっていない。
本発明では、ステップ３）とステップ４）によって与えられたフラグメントの配置に対して絡み合を持たないネットワークを決定している。これは、ネットワークで連結した構造が、後の構造緩和のステップによって適切な分子構造へと緩和されやすくなるようにするためである。与えられたフラグメントの配置に対して絡み合を持たないネットワークを決める為に、フラグメントの結合価不飽和原子で互いに近いもの同士を連結するがその方法は、詳細に後述する。
ステップ５）では、出来上がったネットワークの分子構造に対して構造最適化計算と分子動力学計算およびモンテカルロ計算を行うが、このことは現実の高分子がエネルギー的に安定な構造を取っているという事実に基づいている。
このように、ステップ５）までの手法により、多数のアモルファス網状高分子の試行構造が得られる。
このように、ステップ５）までの手法により、多数のアモルファス網状高分子の試行構造が得られる。
ステップ６）では、ステップ５）までの手順によって得られた多数の試行構造から、微視的な構造パラメータが適切な値を取っているものを分子構造モデルの候補として選択するが、このことは、微視的な構造パラメータには結合の種類毎に標準的な値が存在するという事実に基づいており、高分子においても構造パラメータが標準値のごく近傍の値を取ることを前提にしている。
ステップ７）では、実験的に物性値が測定されている場合に、この測定値が計算によってよく再現されるような構造を選択するが、これにより、アモルファス網状高分子の特徴をよりよく取り込んだ分子構造モデルを得ることができる。
なお、ステップ６）、ステップ７）の手順において、将来、物性値を予測することができる理論的な方法や経験的な法則が開発された場合には、その物性値をこれらのステップ６）、７）の選択に利用することが可能である。
なお、ステップ６）、ステップ７）の手順において、将来、物性値を予測することができる理論的な方法や経験的な法則が開発された場合には、その物性値をこれらのステップ６）、７）の選択に利用することが可能である。
本発明の方法を使用することにより、指定された原子団組成を持つアモルファス網状高分子の試行構造を迅速に得ることができ、この試行構造から分子構造モデルを速やかに得ることができる。こうして生成される分子構造モデルを使用すれば、これまで出来なかった網状高分子の分子シミュレーションが可能になる。
以下、図面を参照して、本発明の実施形態について説明する。
図１は、与えられた網状高分子のフラグメントから分子構造モデルを得るまでの一連のデータ処理のフローチャートである。上述のステップ１）ないし７）はそれぞれフローチャートのＳ２，Ｓ３、Ｓ４−１，Ｓ４−２、Ｓ５、Ｓ６−１およびＳ６−２に該当する。
図１は、与えられた網状高分子のフラグメントから分子構造モデルを得るまでの一連のデータ処理のフローチャートである。上述のステップ１）ないし７）はそれぞれフローチャートのＳ２，Ｓ３、Ｓ４−１，Ｓ４−２、Ｓ５、Ｓ６−１およびＳ６−２に該当する。
図２は、このフローチャートによる一連の処理を視覚的に表現した図であり、この図では、図２−１）に示した５種のフラグメントをＳｉ…Ｏ間で結合してＳｉ５Ｏ８Ｃ２Ｈ８をユニットとするポリシロキサンを作る場合を例に取っている。
図２−１）ないし−５）はそれぞれ、図１のフローチャートの各ステップＳ１，Ｓ３，Ｓ４−１，Ｓ４−２，Ｓ５の処理を視覚的に表現するものである。
図２−１）ないし−５）はそれぞれ、図１のフローチャートの各ステップＳ１，Ｓ３，Ｓ４−１，Ｓ４−２，Ｓ５の処理を視覚的に表現するものである。
・フローチャートのステップＳ１では、網状高分子のフラグメントの種類と個数が与えられる（図２−１））。
・ステップＳ２は、フラグメントの種類ごとに単独フラグメントの構造最適化を行う部分である。
・ステップＳ２−１では、フラグメントの構造最適化のための初期構造（仮構造）を与える。構造最適化に使用するソフトウェアの入力書式に従って初期構造を指定する。
・ステップＳ２は、フラグメントの種類ごとに単独フラグメントの構造最適化を行う部分である。
・ステップＳ２−１では、フラグメントの構造最適化のための初期構造（仮構造）を与える。構造最適化に使用するソフトウェアの入力書式に従って初期構造を指定する。
・ステップＳ２−２では、ステップＳ２−１で与えた初期構造を出発にフラグメントの構造を最適化する。
・ステップＳ３は、フラグメントの乱雑な空間的配置を複数個発生させる部分である（図２−２））。
・ステップＳ３−１では、ステップＳ−２で決定したフラグメントの最適化構造を、必要な個数だけ複製し適当に空間的に配置する。
・ステップＳ３は、フラグメントの乱雑な空間的配置を複数個発生させる部分である（図２−２））。
・ステップＳ３−１では、ステップＳ−２で決定したフラグメントの最適化構造を、必要な個数だけ複製し適当に空間的に配置する。
・ステップＳ３−２では、ステップＳ３−１で決定したフラグメントの空間的配置を初期構造にして、フラグメントの空間配置を乱雑化するための分子動力学計算あるいはモンテカルロ計算を行う。分子動力学計算あるいはモンテカルロ計算は、周知のものである。
・ステップＳ３−３では、分子動力学計算あるいはモンテカルロ計算で得られる構造から、ステップ４）以降の処理に使うものを記録する。すなわち、分子動力学計算あるいはモンテカルロ計算では構造の微少変化を逐次的に計算するので、定期的な間隔をおいたステップの構造（スナップショット）のみを記録して、類似構造の処理をしないようにする（空間配置サンプリング）。
こうして記録された構造がステップＳ１で与えられたフラグメントの乱雑空間配置である。
なお、分子動力学計算法あるいはモンテカルロ計算法は、巨大分子や分子集合体の構造を乱雑化するための手段としてしばしば使用される。このことは、例えば参考文献１に記載されている。
なお、分子動力学計算法あるいはモンテカルロ計算法は、巨大分子や分子集合体の構造を乱雑化するための手段としてしばしば使用される。このことは、例えば参考文献１に記載されている。
・ステップＳ４は、フラグメントの連結を決定する部分である。このステップが本発明の重要な部分である。
・ステップＳ４−１では、フラグメントを暫定的に連結する（図２−３））。
・ステップＳ４−１−１では、連結される可能性のある原子対の距離を計算しテーブル化する。図２の例においては、異なるフラグメント上にある不飽和のＳｉ原子とＯ原子の任意の対が連結される可能性がある。
・ステップＳ４−１では、フラグメントを暫定的に連結する（図２−３））。
・ステップＳ４−１−１では、連結される可能性のある原子対の距離を計算しテーブル化する。図２の例においては、異なるフラグメント上にある不飽和のＳｉ原子とＯ原子の任意の対が連結される可能性がある。
・ステップＳ４−１−２では、ステップＳ４−１−１で得られた距離テーブルを参照に、距離の短いものから順に不飽和原子対を連結する。この時、途中で原子価が飽和した原子には新たに連結を作成しないようにする。
・ステップＳ４−１−３では、ステップＳ−４−１−２での操作によって全ての不飽和原子価を飽和させられない場合を判断する。ステップＳ４−１−２のように連結を決定し不飽和原子価を飽和させられない場合として、たとえば最後に連結する原子対が同一フラグメント上にあるという場合がある。
・ステップＳ４−１−３では、ステップＳ−４−１−２での操作によって全ての不飽和原子価を飽和させられない場合を判断する。ステップＳ４−１−２のように連結を決定し不飽和原子価を飽和させられない場合として、たとえば最後に連結する原子対が同一フラグメント上にあるという場合がある。
・ステップＳ４−２では、後のステップの構造緩和の操作の後に適切な構造が得られるように、ステップＳ４−１で決定した連結のネットワークを改変する。すなわち、ステップＳ４−１で決定したネットワークは、往々にして結合が不自然に絡まった構造を与え易いのでこれを回避する。
このために、長距離の連結が出来るだけなくなる様に連結原子対の組み合わせを交換する（図２−４））。
このために、長距離の連結が出来るだけなくなる様に連結原子対の組み合わせを交換する（図２−４））。
・ステップＳ４−２−１では、最長の連結が短くなるように連結の対を交換する。具体的な交換方法は、最長の連結を含む２個の連結で連結原子対の交換を試行して、長い方の連結が交換の後に短くなるならこの交換を採用する。改めて連結の距離を調べ、最長の連結に対してこの操作を行い、最長の連結が変わらなくなるまで繰り返す。
・ステップＳ４−２−２では、最長の連結対は変更しないことにしておいて、２番目に長い連結が短くなるように、ステップＳ４−２−１の連結交換の操作を行う。同様にして１，２，・・・ｋ番目に長い連結対を変更しないことにしておき、ｋ＋１番目に長い連結対を最小化するための連結交換を行う。
ステップＳ４−２−１からステップＳ４−２−４の操作を５回繰り返す。
・ステップＳ４−２の上述の一連の処理は、図３に示したフローチャートのように整理して示すことができる。
ステップＳ４−２−１からステップＳ４−２−４の操作を５回繰り返す。
・ステップＳ４−２の上述の一連の処理は、図３に示したフローチャートのように整理して示すことができる。
・ステップＳ５では、構造緩和のための構造最適化と分子動力学計算を行う（図２−５））。
・Ｓ５−１では、次ぎに行う分子動力学計算に使う初期構造を得るための構造最適化を行う。
・Ｓ５−１では、次ぎに行う分子動力学計算に使う初期構造を得るための構造最適化を行う。
・ステップＳ５−２では、ステップ５−１で得られた最適化構造を初期構造にし、アニーリングのための分子動力学計算あるいはモンテカルロ計算を行う。
・ステップ６は、ステップ５）までに得られた分子構造から分子構造モデルとして適切なものを選択する部分である。
・ステップ６は、ステップ５）までに得られた分子構造から分子構造モデルとして適切なものを選択する部分である。
・ステップＳ６−１では、結合長、結合角などの微視的な構造パラメータが標準値付近の適当な範囲内に収まっているものを選択する。
・ステップＳ６−２では、実験的に物性値が測定されている場合には、この測定結果が計算によって良く再現されるような構造を選択する。
・ステップＳ６−２では、実験的に物性値が測定されている場合には、この測定結果が計算によって良く再現されるような構造を選択する。
図４は、以上の計算に必要なソフトウェア構成と、それぞれのソフトウェアでのデータ授受を可視化して描いた図である。
手作業以外の処理は、４種の独立したプログラム（分子構造最適化プログラム、分子動力学計算あるいはモンテカルロ計算のプログラム、ネットワーク生成プログラム、物性計算プログラム）で処理され、これらのプログラムはパーソナルコンピュータやワークステーション上で稼動する。ネットワーク生成プログラムは上記のステップ３）と４）の操作（フローチャートではＳ４の操作）を行う。物性計算プログラムが使用されない場合もある。また、分子動力学計算あるいはモンテカルロ計算が物性計算を兼ねる場合もある。
手作業以外の処理は、４種の独立したプログラム（分子構造最適化プログラム、分子動力学計算あるいはモンテカルロ計算のプログラム、ネットワーク生成プログラム、物性計算プログラム）で処理され、これらのプログラムはパーソナルコンピュータやワークステーション上で稼動する。ネットワーク生成プログラムは上記のステップ３）と４）の操作（フローチャートではＳ４の操作）を行う。物性計算プログラムが使用されない場合もある。また、分子動力学計算あるいはモンテカルロ計算が物性計算を兼ねる場合もある。
以下に、本件発明の予測法を用いて、具体的に有機シロキサン膜のモデリングを行った例を説明する。
この例のシロキサン膜モデルは、ユニットセル内に図５に示した５種のフラグメントを持ち、これらの結合価不飽和のＳｉ原子とＯ原子を連結することによって作成される。
この例のシロキサン膜モデルは、ユニットセル内に図５に示した５種のフラグメントを持ち、これらの結合価不飽和のＳｉ原子とＯ原子を連結することによって作成される。
こうして作成されるシロキサン膜モデルは、７個のＣＨ３ＳｉＯ３と３個のＨＳｉＯ３と１個のＳｉＯ４を持つ網状構造になり、原子組成としてはＳｉ１１Ｏ１７Ｃ７Ｈ２４（原子個数比はＳｉ：Ｏ：Ｃ：Ｈ＝１：１．５４：０．６７：２．１４）を持つことになる。
このモデリングの対象となっている実物の有機シロキサン膜は、テトラメチルシクロテトラシロキサン分子を原料としてプラズマＣＶＤ法によって生成されたもので、その物性値に関しては表１に示したような実験観測の結果が得られている。上述のフラグメントの種類と組成は、実験観測の結合種と原子組成が再現されるように選んだものである。
この例においては、フラグメントの乱雑構造生成（ステップ２））のために分子動力学法を使用した。１フェムト秒の時間刻みで５５０ピコ秒間の構造変化を求め、１ピコ秒毎の構造を記録する（空間配置サンプリング）ことにより５００個（初めの５０ピコ秒分を除く）のフラグメントの乱雑空間配置を得た。
この分子動力学計算は、温度が２９８．１５Ｋで圧力が１００００ａｔｍの定温定圧分子動力学計算（参考文献２）である。なお、ステップ３）とステップ４）を経てフラグメントを連結し、不適当な連結を持つものを除去した結果、得られた試行構造は３４３個となった。
構造緩和の分子動力学計算（ステップ５））では、温度が２９８．１５Ｋで圧力が１ａｔｍの定温定圧の条件で、１フェムト秒の時間刻みで１０ピコ秒の構造変化を求めた。
ステップ１）と５）の構造最適化計算およびステップ２）と５）の分子動力学計算にはＭＭ３分子力場（参考文献３）を使用した。また、構造最適化計算および分子動力学計算にはプログラムパッケージＴｉｎｋｅｒ（参考文献４）を使用した。
ステップ１）と５）の構造最適化計算およびステップ２）と５）の分子動力学計算にはＭＭ３分子力場（参考文献３）を使用した。また、構造最適化計算および分子動力学計算にはプログラムパッケージＴｉｎｋｅｒ（参考文献４）を使用した。
図６に本発明方法を使用して生成した３４３個の試行構造の密度の分布を表した。図６から実験値（１．３８ｇ／ｃｍ３）に近い密度を有する試行構造がいくつか得られていることがわかる。
表２に、上記３４３個の試行構造のなかで、結合長の偏差の小さい構造（Ｓｉ−Ｏ、Ｓｉ−Ｈ、Ｓｉ−Ｃ、Ｃ−Ｈの結合種それぞれに関して最大値と最小値との差が０．０５オングストローム以下の構造）を６種選択し、リストとした。実際には、Ｓｉ−Ｈ、Ｓｉ−Ｃ、Ｃ−Ｈは結合長偏差が小さく、表の値はＳｉ−Ｏ結合のものである。なお、シロキサン高分子では、結合がイオン結合的であるので、結合角、二面角がかなり大きな範囲に分布するので、分子構造モデル選定の参照にはしていない。表２から結合長が、標準値の近傍に分布する構造が得られたことがわかる。
表２に、上記３４３個の試行構造のなかで、結合長の偏差の小さい構造（Ｓｉ−Ｏ、Ｓｉ−Ｈ、Ｓｉ−Ｃ、Ｃ−Ｈの結合種それぞれに関して最大値と最小値との差が０．０５オングストローム以下の構造）を６種選択し、リストとした。実際には、Ｓｉ−Ｈ、Ｓｉ−Ｃ、Ｃ−Ｈは結合長偏差が小さく、表の値はＳｉ−Ｏ結合のものである。なお、シロキサン高分子では、結合がイオン結合的であるので、結合角、二面角がかなり大きな範囲に分布するので、分子構造モデル選定の参照にはしていない。表２から結合長が、標準値の近傍に分布する構造が得られたことがわかる。
図７には、これら多数の試行構造から得られた分子構造モデルの一つ（表２における構造２のもの）を示す。この構造は、実測の有機シロキサン膜と近い密度を示す構造である。
図８は、この構造に対して量子化学的手法を用いて計算された赤外吸収スペクトルｂ）と、実測の有機シロキサン膜の赤外吸収スペクトルａ）を示す。実測のスペクトルが良く再現されていることがわかる。このことは、この構造が適切なものであることを（原子や原子団の組成が実物とほぼ等しく、結合長、結合角、二面体角の構造的パラメータの分布が実物とほぼ等しいこと）を示している。
この例は、本発明の方法を使用して、化学結合のネットワークと立体的分子構造を生成した結果、網状高分子の構造を良く表現する分子構造モデルが得られたことを示している。
上述の３４３個の試行構造を得るのに費やした時間は２日の程度である。本発明の方法を使用しない場合、図７に示した構造のようにアモルファス網状高分子の構造を良く表現する分子構造モデルをこの時間の範囲で得ることは容易ではない。
上述の３４３個の試行構造を得るのに費やした時間は２日の程度である。本発明の方法を使用しない場合、図７に示した構造のようにアモルファス網状高分子の構造を良く表現する分子構造モデルをこの時間の範囲で得ることは容易ではない。
（参考文献）
[1] F.Nardi and R.C.Wade,in Molecular Dynamics:From Classical to Quantum Method，ed.P.B.Balbuena and J.M.Seminario,Elsevier,1999,pp.859-898
[2] H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. DiNola, and J. R. Haak, J. Chem. Phys., 81, 3684 (1984).
[3] N. L. Allinger, Y. H. Yuh, and J.-H. Lii, J. Am. Chem. Soc., 111, 8551 (1989)., J.-H. Lii and N. L. Allinger, J. Am. Chem. Soc., 111, 8566 (1989)., J.-H. Lii, and N. L. Allinger, J. Am. Chem. Soc., 111, 8576 (1989).
[4] Ponder, J. W. TINKER: Software Tools for Molecular Design, 3.9 ed.; Washington University School of Medicine: Saint Louis, MO, 2001.
[1] F.Nardi and R.C.Wade,in Molecular Dynamics:From Classical to Quantum Method，ed.P.B.Balbuena and J.M.Seminario,Elsevier,1999,pp.859-898
[2] H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. DiNola, and J. R. Haak, J. Chem. Phys., 81, 3684 (1984).
[3] N. L. Allinger, Y. H. Yuh, and J.-H. Lii, J. Am. Chem. Soc., 111, 8551 (1989)., J.-H. Lii and N. L. Allinger, J. Am. Chem. Soc., 111, 8566 (1989)., J.-H. Lii, and N. L. Allinger, J. Am. Chem. Soc., 111, 8576 (1989).
[4] Ponder, J. W. TINKER: Software Tools for Molecular Design, 3.9 ed.; Washington University School of Medicine: Saint Louis, MO, 2001.
Claims (1)
- 与えられた構成原子団に基づいて、アモルファス網状高分子の分子構造モデルを予測する方法であって、
１)アモルファス網状高分子を構成する構成原子団の個々の構造の最適化を行うステップと、
２)分子動力学法あるいはモンテカルロ法を使用して構成原子団の乱雑な空間的配置を決めるステップと、
３) 与えられた空間的配置を構成する構成原子団を暫定的に化学結合により連結するステップと、
４)構成原子団の連結結合のうちの２組の間で連結原子対を交換することによって、絡み合いを持たないように連結結合のネットワークを変更するステップと、
５) 出来上がったネットワークの分子構造に対して構造最適化計算と分子動力学計算およびモンテカルロ計算を行うステップと、
６）上記５つのステップによって得られるアモルファス網状高分子の試行構造から、結合長、結合角などの微視的な構造パラメータが適切な値を取っているものを分子構造モデルとして選択するステップと、
７）実験的に得られる物性値を再現することを条件に、アモルファス網状高分子の分子構造モデルとしてより適切な分子構造モデルを選択するステップを備えることを特徴とする分子構造予測方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005107361A JP2006282929A (ja) | 2005-04-04 | 2005-04-04 | 分子構造予測方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005107361A JP2006282929A (ja) | 2005-04-04 | 2005-04-04 | 分子構造予測方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2006282929A true JP2006282929A (ja) | 2006-10-19 |
Family
ID=37405172
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2005107361A Pending JP2006282929A (ja) | 2005-04-04 | 2005-04-04 | 分子構造予測方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2006282929A (ja) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007031541A (ja) * | 2005-07-26 | 2007-02-08 | Sumitomo Chemical Co Ltd | 高分子材料設計方法および装置 |
| JP2010256300A (ja) * | 2009-04-28 | 2010-11-11 | Sumitomo Rubber Ind Ltd | 高分子材料モデルの作成方法 |
| EP2645279A2 (en) | 2012-03-29 | 2013-10-02 | Fujitsu Limited | A program, apparatus, and method for determining bond types between atoms |
| CN104200071A (zh) * | 2014-08-15 | 2014-12-10 | 浙江师范大学 | 一种预测羟基取代的多溴代联苯醚甲状腺激素效应的方法及模型的建立方法 |
| CN106096258A (zh) * | 2016-06-07 | 2016-11-09 | 上海理工大学 | 一种含氟聚氨酯表面疏水性的分析方法 |
| JP2018070692A (ja) * | 2016-10-25 | 2018-05-10 | 東洋ゴム工業株式会社 | 架橋ゴムの分子モデルにおけるモデル間の再現性を向上させる分子モデル作成方法、及び同方法を実行する装置 |
| JP2020030796A (ja) * | 2018-08-23 | 2020-02-27 | タタ コンサルタンシー サービシズ リミテッドＴＡＴＡ Ｃｏｎｓｕｌｔａｎｃｙ Ｓｅｒｖｉｃｅｓ Ｌｉｍｉｔｅｄ | 原子及び原子合金材料の構造及び特性を予測するシステム及び方法 |
| CN111192628A (zh) * | 2020-02-27 | 2020-05-22 | 西安理工大学 | 一种聚二甲基硅氧烷的氧气扩散系数计算方法 |
| CN115831271A (zh) * | 2022-09-23 | 2023-03-21 | 哈尔滨工业大学 | 一种CuZrAl多晶模型的分子动力学模拟构建方法 |
Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1997006496A1 (fr) * | 1995-08-04 | 1997-02-20 | Hitachi, Ltd. | Systeme de conception de materiaux polymeres |
| JPH09127030A (ja) * | 1995-11-01 | 1997-05-16 | Sekisui Chem Co Ltd | 無定形高分子のガラス転移温度の予測方法 |
| JPH09223166A (ja) * | 1996-02-15 | 1997-08-26 | Toray Ind Inc | 分子または分子集合体の構造の解析方法および解析装置 |
| JP2002189978A (ja) * | 2000-12-19 | 2002-07-05 | Nippon Steel Corp | 分子材料の解析装置、分子材料の解析方法及び記憶媒体 |
| JP2003012566A (ja) * | 2001-07-02 | 2003-01-15 | Ricoh Co Ltd | 分子シミュレーション方法 |
| JP2003106989A (ja) * | 2001-09-27 | 2003-04-09 | Asahi Kasei Corp | 高分子材料の密着性を予測する方法、装置および粗視化モデルを記述するパラメーターを決定する方法、装置 |
| JP2003105090A (ja) * | 2001-09-28 | 2003-04-09 | Asahi Kasei Corp | 高分子材料の形態予測方法およびその装置 |
| JP2003213003A (ja) * | 2002-01-23 | 2003-07-30 | Sumitomo Chem Co Ltd | 重合反応のシミュレーション装置、重合反応のシミュレーションプログラム、記録媒体、および、シミュレーション方法 |
| JP2003528288A (ja) * | 1999-05-26 | 2003-09-24 | ザ・リージェンツ・オブ・ジ・ユニバーシティ・オブ・カリフォルニア | 高分子の三次元形状の決定方法 |
| JP2003272980A (ja) * | 2002-03-14 | 2003-09-26 | Ube Ind Ltd | 薄膜内部構造シュミレーション法及びその装置 |
| JP2003317118A (ja) * | 2002-04-22 | 2003-11-07 | Japan Science & Technology Corp | 分子間相互作用を構造に基づき判定するための支援システム |
| JP2004192275A (ja) * | 2002-12-10 | 2004-07-08 | Japan Science & Technology Agency | シミュレーション方法及びその装置、プログラム |
| JP2004527027A (ja) * | 2000-11-02 | 2004-09-02 | プロテイン メカニクス， インコーポレイテッド | 分子モデリングにおける解析的ヤコビアン計算のための方法 |
| JP2004287812A (ja) * | 2003-03-20 | 2004-10-14 | Japan Research Institute Ltd | 物性値算出方法、物性値算出装置、画像作成装置、コンピュータプログラム及び記録媒体 |
| JP2005032058A (ja) * | 2003-07-08 | 2005-02-03 | Nippon Zeon Co Ltd | 高分子材料の相分離構造を予測する方法及びその予測方法により高分子材料を製造する方法 |
| JP2005524129A (ja) * | 2002-03-26 | 2005-08-11 | カウンシル・オブ・サイエンティフィック・アンド・インダストリアル・リサーチ | ３次元分子構造の最適モデルを構築する方法及びシステム |
| JP2006028038A (ja) * | 2004-07-13 | 2006-02-02 | Ishikawajima Harima Heavy Ind Co Ltd | 高分子化合物の結合評価方法 |
| JP2006072578A (ja) * | 2004-08-31 | 2006-03-16 | Japan Science & Technology Agency | 実験装置と連動させる仮想実験インタフェース |
-
2005
- 2005-04-04 JP JP2005107361A patent/JP2006282929A/ja active Pending
Patent Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1997006496A1 (fr) * | 1995-08-04 | 1997-02-20 | Hitachi, Ltd. | Systeme de conception de materiaux polymeres |
| JPH09127030A (ja) * | 1995-11-01 | 1997-05-16 | Sekisui Chem Co Ltd | 無定形高分子のガラス転移温度の予測方法 |
| JPH09223166A (ja) * | 1996-02-15 | 1997-08-26 | Toray Ind Inc | 分子または分子集合体の構造の解析方法および解析装置 |
| JP2003528288A (ja) * | 1999-05-26 | 2003-09-24 | ザ・リージェンツ・オブ・ジ・ユニバーシティ・オブ・カリフォルニア | 高分子の三次元形状の決定方法 |
| JP2004527027A (ja) * | 2000-11-02 | 2004-09-02 | プロテイン メカニクス， インコーポレイテッド | 分子モデリングにおける解析的ヤコビアン計算のための方法 |
| JP2002189978A (ja) * | 2000-12-19 | 2002-07-05 | Nippon Steel Corp | 分子材料の解析装置、分子材料の解析方法及び記憶媒体 |
| JP2003012566A (ja) * | 2001-07-02 | 2003-01-15 | Ricoh Co Ltd | 分子シミュレーション方法 |
| JP2003106989A (ja) * | 2001-09-27 | 2003-04-09 | Asahi Kasei Corp | 高分子材料の密着性を予測する方法、装置および粗視化モデルを記述するパラメーターを決定する方法、装置 |
| JP2003105090A (ja) * | 2001-09-28 | 2003-04-09 | Asahi Kasei Corp | 高分子材料の形態予測方法およびその装置 |
| JP2003213003A (ja) * | 2002-01-23 | 2003-07-30 | Sumitomo Chem Co Ltd | 重合反応のシミュレーション装置、重合反応のシミュレーションプログラム、記録媒体、および、シミュレーション方法 |
| JP2003272980A (ja) * | 2002-03-14 | 2003-09-26 | Ube Ind Ltd | 薄膜内部構造シュミレーション法及びその装置 |
| JP2005524129A (ja) * | 2002-03-26 | 2005-08-11 | カウンシル・オブ・サイエンティフィック・アンド・インダストリアル・リサーチ | ３次元分子構造の最適モデルを構築する方法及びシステム |
| JP2003317118A (ja) * | 2002-04-22 | 2003-11-07 | Japan Science & Technology Corp | 分子間相互作用を構造に基づき判定するための支援システム |
| JP2004192275A (ja) * | 2002-12-10 | 2004-07-08 | Japan Science & Technology Agency | シミュレーション方法及びその装置、プログラム |
| JP2004287812A (ja) * | 2003-03-20 | 2004-10-14 | Japan Research Institute Ltd | 物性値算出方法、物性値算出装置、画像作成装置、コンピュータプログラム及び記録媒体 |
| JP2005032058A (ja) * | 2003-07-08 | 2005-02-03 | Nippon Zeon Co Ltd | 高分子材料の相分離構造を予測する方法及びその予測方法により高分子材料を製造する方法 |
| JP2006028038A (ja) * | 2004-07-13 | 2006-02-02 | Ishikawajima Harima Heavy Ind Co Ltd | 高分子化合物の結合評価方法 |
| JP2006072578A (ja) * | 2004-08-31 | 2006-03-16 | Japan Science & Technology Agency | 実験装置と連動させる仮想実験インタフェース |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007031541A (ja) * | 2005-07-26 | 2007-02-08 | Sumitomo Chemical Co Ltd | 高分子材料設計方法および装置 |
| JP2010256300A (ja) * | 2009-04-28 | 2010-11-11 | Sumitomo Rubber Ind Ltd | 高分子材料モデルの作成方法 |
| EP2645279A2 (en) | 2012-03-29 | 2013-10-02 | Fujitsu Limited | A program, apparatus, and method for determining bond types between atoms |
| CN104200071A (zh) * | 2014-08-15 | 2014-12-10 | 浙江师范大学 | 一种预测羟基取代的多溴代联苯醚甲状腺激素效应的方法及模型的建立方法 |
| CN106096258A (zh) * | 2016-06-07 | 2016-11-09 | 上海理工大学 | 一种含氟聚氨酯表面疏水性的分析方法 |
| JP2018070692A (ja) * | 2016-10-25 | 2018-05-10 | 東洋ゴム工業株式会社 | 架橋ゴムの分子モデルにおけるモデル間の再現性を向上させる分子モデル作成方法、及び同方法を実行する装置 |
| JP2020030796A (ja) * | 2018-08-23 | 2020-02-27 | タタ コンサルタンシー サービシズ リミテッドＴＡＴＡ Ｃｏｎｓｕｌｔａｎｃｙ Ｓｅｒｖｉｃｅｓ Ｌｉｍｉｔｅｄ | 原子及び原子合金材料の構造及び特性を予測するシステム及び方法 |
| CN111192628A (zh) * | 2020-02-27 | 2020-05-22 | 西安理工大学 | 一种聚二甲基硅氧烷的氧气扩散系数计算方法 |
| CN111192628B (zh) * | 2020-02-27 | 2023-03-24 | 西安理工大学 | 一种聚二甲基硅氧烷的氧气扩散系数计算方法 |
| CN115831271A (zh) * | 2022-09-23 | 2023-03-21 | 哈尔滨工业大学 | 一种CuZrAl多晶模型的分子动力学模拟构建方法 |
| CN115831271B (zh) * | 2022-09-23 | 2023-08-08 | 哈尔滨工业大学 | 一种CuZrAl多晶模型的分子动力学模拟构建方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2006282929A (ja) | 分子構造予測方法 | |
| Kobayashi et al. | Neural network potential for Al-Mg-Si alloys | |
| RU2727631C1 (ru) | Системы и способы прогнозирования структуры и свойств атомарных веществ и их сплавов | |
| Marti et al. | Complete-graph tensor network states: a new fermionic wave function ansatz for molecules | |
| Ghoniem et al. | Multiscale modelling of nanomechanics and micromechanics: an overview | |
| Sokolov et al. | Time-dependent N-electron valence perturbation theory with matrix product state reference wavefunctions for large active spaces and basis sets: Applications to the chromium dimer and all-trans polyenes | |
| Pedersen et al. | Simulations of hydrogen diffusion at grain boundaries in aluminum | |
| Mao et al. | Elasticity of a filamentous kagome lattice | |
| Javadi et al. | Design and optimization of microstructure of auxetic materials | |
| Zinchenko et al. | Nanoporous gold formation by dealloying: A Metropolis Monte Carlo study | |
| Park et al. | Dibaryons with two strange quarks and total spin zero in a constituent quark model | |
| Masur et al. | Fragment-based direct-local-ring-coupled-cluster doubles treatment embedded in the periodic Hartree–Fock solution | |
| Ozdemir et al. | A novel design framework for generating functionally graded multi-morphology lattices via hybrid optimization and blending methods | |
| Tsai et al. | Genetic algorithm for discovery of globally stable phases in block copolymers | |
| Shi et al. | Effect of single-site mutations on hydrophobic-polar lattice proteins | |
| Kazemi et al. | Stereolithography process optimization for tensile strength improvement of products | |
| Huang et al. | Cluster analysis of accelerated molecular dynamics simulations: A case study of the decahedron to icosahedron transition in Pt nanoparticles | |
| Biswas et al. | Materials modeling by design: applications to amorphous solids | |
| Kolluri et al. | Behavior of vacancies and interstitials at semicoherent interfaces | |
| Moat et al. | A class of aperiodic honeycombs with tuneable mechanical properties | |
| CN116913414A (zh) | 基于分子生成模型的融合亲脂性药物分子生成方法及系统 | |
| Soper et al. | Empirical potential structure refinement of semi-crystalline polymer systems: Polytetrafluoroethylene and polychlorotrifluoroethylene | |
| Cherman et al. | Alternate 1/N c expansions and SU (3) breaking from baryon lattice results | |
| Nakamura et al. | An efficient molecular orbital approach for self-consistent calculations of molecular junctions | |
| Albright et al. | Explicit SU (12) family and flavor unification model with natural fermion masses and mixings |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080312 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20080312 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100422 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100506 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20100914 |