(57)【特許請求の範囲】
【請求項１】
長鎖アルキル鎖をもつ化合物の分子集合構造を当該化合物の分子モデルを用いて電子計算機によるシミュレーションを行って予測する分子シミュレーション方法であって、
分子モデルを作成するステップと、
前記分子モデルを正方形セルに配置するステップと、
前記セル内に存在する分子に対して働く力を定義するステップと、
前記セルを定義するステップと、
上記各ステップにおいて設定した条件に対応するデータを電子計算機に入力するステップと、
電子計算機が該入力値を基に分子動力学法によるプログラム（ソフトウエア）に基づいてシミュレーションを行うステップと、
を含むことを特徴とする分子シミュレーション方法。
【請求項２】
前記分子に対して働く力としては，非結合性相互作用力として１２−６  Ｌｅｎｎａｒｄ−Ｊｏｎｅｓポテンシャルを用い、分子末端間においてのみそれ以外の部位より３倍以上９倍以下の値とし、また３体力ポテンシャルとして結合の狭角度の動きを組み込み、ポテンシャルパラメーターとしてはＤｒｅｉｄｉｎｇのパラメーターを用い、６０以上の分子モデルを３次元空間において、密度を０．４以上１．２以下に設定し、原子間の結合長さをＮＩＭＭ法により固定し、かつ，計算方法にリンクセル法を適用したことを特徴とする請求項１に記載の分子シミュレーション方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、棒状分子の高次構造を予測するための分子シミュレーション方法、並びに分子に任意の引力を与えた際に、高次構造がどのように変化し、また、温度に対してどのような挙動をするかを予測する結晶構造予測分子シミュレーション方法に関する。
【０００２】
【従来技術】
最近の情報化社会の発展とともに、電子情報と人間をつなげる接点であるディスプレイやハードコピーの技術の重要性がますます高まっている。昔から情報媒体として使われている“紙”の有効性は、単に携帯性や保存性だけではなく、人間の情報認識方法や思考方法と関係している。そのため、ＣＲＴ やＬＣＤ など従来の表示媒体だけでは十分ではなく、メモリー性があり形態的にも画像品質の点でもペーパーライクな表示媒体の実現が待たれている。
【０００３】
このような表示媒体への要求に対し、記録材料側からのアプローチとしてさまざまなリライタブル記録媒体が検討されてきた。特に実用性の高さから熱エネルギーの制御で書き替え可能な記録媒体の開発が盛んに行われている。たとえば、高分子中に脂肪酸粒子を分散した複合膜の透明・白濁変化を利用したリライタブル記録媒体は、すでに磁気カード用の表示媒体として広く使われている。さらにハードコピーに近い白地に黒発色画像の形成が可能なリライタブル記録媒体を得るため、ロイコ染料の可逆発色現象の利用が考えられた。
【０００４】
代表的なロイコ染料は一般的な感熱記録紙に用いられているフルオラン化合物であり、顕色剤と呼ばれる酸性物質の作用によりラクトン環が開き発色する。感熱紙ではこの発色状態を安定にとり得るロイコ染料と顕色剤との組み合わせが用いられているが、基本的にはロイコ染料の発色反応は可逆的であり、たとえば塩基性物質の作用によって消色することがわかっている。
【０００５】
リライタブル記録媒体では、十分な発色濃度が得られ、画像が完全に消去できる発色・消色特性をもつことは前提であるが、発色画像の熱安定性と消色速度も実用的に重要な特性となる。特に、使用環境や書き替え処理時間から考えると、発色画像の耐熱性は高く、なるべく速く消去できることが要求される。このような特性には、長鎖型顕色剤の分子集合構造形成能がかかわっていると考えられる。
【０００６】
しかしながら、走査型トンネル顕微鏡（ＳＴＭ）のような分子の形を計測できる手法をもってしても、実験的にこのような特性が分子のどのような性質に影響されているかを明らかにすることが非常に困難であった。さらに透過型電子顕微鏡(ＴＥＭ)で高次構造の解析を行う場合、結晶性が高い材料においては高次構造がある程度判別可能であるが、結晶性が低い材料では高次構造を明らかにすることができなかった。
【０００７】
長鎖型分子のシミュレーションによる分子集合構造の解析は、液晶分子においてこれまで幾つかの検討が行われ、実験結果を再現するシミュレーション方法が提案されている。例えば、青木圭子、秋山鉄夫、（物性研究、vol.68,No.3,p.316-319,(1997)）らは、分子のどのような性質が液晶の秩序相を形成するかを明らかにすべく、圧力一定のもとでの分子動力学シミュレーションを行っている。ここでは、液晶相が出現する必要条件の一つである「形状の異方性をもつ」という特徴だけを考慮した単純な模型でも、昇温過程においてスメクチックおよびネマチック液晶を経て異方性液体に相転移することを示した。
この計算例では液晶分子を粗視化して計算を行っており、大局的に現象を理解することはできるもののより具体的な、分子それぞれの個性を考慮した計算はなされていない。
【０００８】
また、岡崎進、（日本原子力研究所JAERI-Conf，p.54-72,(1998)）は、脂質二重層膜の長時間分子動力学シミュレーションを行い、脂質分子の構造や運動、さらに膜を特徴づける揺らぎを、ナノ秒スケールの計算で調べている。ここでの計算はアトミックスケールの計算であり、前者と異なって分子それぞれの個性を考慮した計算が可能となっているが、計算時間が数ヶ月と非常に長期間かかるという問題がある。
【０００９】
また、たとえこれらの方法を用いたとしても、長鎖型顕色剤の分子挙動が実験結果と合致するような分子動力学シミュレーションはできない。
【００１０】
また特願２０００−２６６９１６号の発明では、計算は、Δｔ＝１ｆｓを用い、始め１０００Ｋにて４ｎｓｅｃ計算を行い、熱平衡状態を実現させた上で、次いで温度を１ｎｓｅｃの間に４００Ｋまで低下させた。その後、４ｎｓｅｃ計算を続けた。しかしながら、これでは計算時間が１週間以上必要となり、パラメーターを変えて計算をすることを考えると多大な計算時間がかかることで、実用上問題があった。
【００１１】
【発明が解決しようとする課題】
上述したように、実験的な方法では分子の高次構造の解析を容易に行うことができないという問題があった。また従来のシミュレーション方法では前述したように粗視化したモデルでは、比較的短時間でシミュレーションによる高次構造解析が可能であるが、分子設計に応用できるわけではなく、一方アトミックスケールのシミュレーションでは膨大な計算時間がかかり、現実的には不可能であるといった問題があった。また本発明に近い発明である特願２０００−２６６９１６号の発明でも計算時間が長く実際的には使用できないものであった。
【００１２】
本発明の目的は、短時間の計算時間で分子それぞれの個性を考慮できるシミュレーション方法を提供することにある。
また、本発明の別の目的は、分子組成が変化したときに、または温度条件が変化したときに、分子が形成する高次構造がどのような状態になるかを明らかにできるシミュレーション方法を提供することにある。
【００１３】
【課題を解決するための手段】
本発明は、上記目的を達成するために以下のような態様を有する。
（１）長鎖アルキル鎖をもつ化合物の分子集合構造を当該化合物の分子モデルを用いて電子計算機によるシミュレーションを行って予測する分子シミュレーション方法であって、
分子モデルを作成するステップと、
前記分子モデルを正方形セルに配置するステップと、
前記セル内に存在する分子に対して働く力を定義するステップと、
前記セルを定義するステップと、
上記各ステップにおいて設定した条件に対応するデータを電子計算機に入力するステップと、
電子計算機が該入力値を基に分子動力学法によるプログラム（ソフトウエア）に基づいてシミュレーションを行うステップと、
を含むことを特徴とする分子シミュレーション方法。
分子モデルを配置させる空間は、密度から決定される。密度は実際の実験環境を考慮して、シミュレーションにより明らかにしたい物性から判断して決定され、一義的に決まるものではない。分子モデルの必要数は多いほど計算精度を上げることができて好ましいものである。しかしながら、現状の計算機速度では分子モデル数には限界がある。したがって、モデル数は少なすぎても境界条件の影響を受けやすいため適当ではない。そのために、６０以上１０００以下が適当な数字になる。しかしながら、分子モデル数はこれに限定されるものではない。
（２）前記分子に対して働く力としては，非結合性相互作用力として１２−６  Ｌｅｎｎａｒｄ−Ｊｏｎｅｓポテンシャルを用い、分子末端間においてのみそれ以外の部位より３倍以上９倍以下の値とし、また３体力ポテンシャルとして結合の狭角度の動きを組み込み、ポテンシャルパラメーターとしてはＤｒｅｉｄｉｎｇのパラメーターを用い、６０以上の分子モデルを３次元空間において、密度を０．４以上１．２以下に設定し、原子間の結合長さをＮＩＭＭ法により固定し、かつ，計算方法にリンクセル法を適用したことを特徴とする請求項１に記載の分子シミュレーション方法。
【００１４】
分子動力学法（Molecular Dynamics:ＭＤ）の概要
ＭＤ法に基づく計算方法の概要を以下、説明する。
ＭＤ法は基本的に古典力学のニュートンの運動方程式による数値解析であり、分子集団系において時間に依存する物理量が求められる。入力パラメーターとして粒子間の力を用いるが、２体力を用いる場合が多い。力はポテンシャルから容易に求められる。その他に、粒子数、体積、粒子の質量、タイムステップ（時間刻み）や、分子の初期配置等を入力する。
【００１５】
【発明の実施の形態】
本発明の好ましい実施の形態を以下に示す。
計算の手法は、次のとおりである。
（１）分子構造を引力が大きい部分とそれ以外の部分のビーズモデルに置き換える。
（２）分子モデルを複数個用いて、２次元セル内に配置する。
（３）初期の構造を高温でシミュレーションし、熱平衡状態構造を発生させる。
（４）．温度を適当な時間の温度勾配で低下させ、安定構造を作成する。
【００１６】
本発明におけるプログラムは、フォートランで作成し、Ｕｎｉｘワークステーションマシンであるコンパック製のＤＥＣαにて計算を実行した。計算を実行するマシンとしてはコンパック製のマシンに限定される必要はない。他社の、例えばＩＢＭ製のＲＳ６０００／ＳＰやＳＧＩ製のＲ１２０００等でもかまわない。
【００１７】
例えば計算対象となる長鎖アルキル基を有する化合物（ＣＨ3−（ＣＨ2）17−Ｃ6Ｈ4−ＯＨ：パラオクタデシルフェノール）をビーズモデルに置き換えると、次のように表すことができる。
○-○-○-○-○-○-○-○-○-○-○-○-○-○-○-○-○-○-◎
（ここで、○は「ＣＨ2」を、末端の○は「ＣＨ3」を、◎は「Ｃ6Ｈ4−ＯＨ」を示す。また、ここでは「ＣＨ2」をモノマーと定義する。）
そして、○の部分は次の式（１）に従って粒子間の相互作用を行なう（ε＝ε0とする）。
【００１８】
【数１】
0004443797.tif 000002
【００１９】
一方、末端分子間の相互作用（◎の部分）については、ε=６×ε0とする。つまり大きな引力が働くことになる。モノマー間には狭角度の調和ポテンシャルが３体間ポテンシャルとして入っている。角度の初期値は１０９．４７１度である。より正確には図１の様なモデルとなる。
【００２０】
ＭＤ法の場合、この分子を等間隔に分子を１００本、２００Åの２次元セル内に設置し、ボンド長はＮＩＭＭ法により１．５３Åに固定して計算を実行した。計算速度の高速化させるためにリンクセル法を適用した。温度は束縛法により制御した。角度のポテンシャルは下記の式（２）のとおりである。
【００２１】
【数２】
0004443797.tif 000003
【００２２】
リンクセル法について
リンクセル法とは、系を正方形のセルで分割して、力の計算時間を短縮する方法のことである。２次元の系では、実際相互作用を計算するのは、ある粒子の所属するセルとその周りに隣接している８つのセル内だけでよくなる。通常、粒子数が２倍になると計算時間は４倍になるが、この方法によって計算時間は約２倍で済むことになる。
【００２３】
ＮＩＭＭ法について
ＮＩＭＭ法とは、ボンドを拘束する方法である。一般的なものにはGaussの方法とShakeの方法がある。それらの長所を生かし、短所を改善したものがＮＩＭＭ法である。Shakeの方法では拘束力を与える未定乗数λｎを決定するのに反復計算が必要になるが、ＮＩＭＭ法では、反復計算なしにλｎを決定できる。（λｎを決定することによってボンド長を固定できる。）
【００２４】
Ｓｈａｋｅの方法について
Shakeの方法とは、ボンド長にある許容誤差εを認めて
【００２５】
【数３】
0004443797.tif 000004
｛上式において、Rn（ｔ）は粒子ｎの時刻tでの位置、Dｎはボンド長（今回は１．５３Å）を表している。｝
【００２６】
となるまで、未定乗数λｎの決定のための繰り返し計算を行う。Shakeの方法の欠点はこの繰り返し計算にある。
【００２７】
計算は、Δｔ＝１ｆｓを用い、始め１０００Ｋにて１ｎｓｅｃ計算を行い熱平衡状態を実現させた上で、次いで温度を１ｎｓｅｃの間に４００Ｋまで低下させた。その後、４ｎｓｅｃ計算を続けた。得られた計算結果は、分子がラメラ状に並んでいる構造を形成した。
得られた構造の例を図２、３に示す。
図２は発色状態にある顕色剤の構造の計算結果を示す図であり、図３は消色状態にある顕色剤の構造についての計算結果を示す図である。
この構造は、透過型電子顕微鏡写真および走査型トンネル顕微鏡の実験から推定される構造と酷似した。
【図面の簡単な説明】
【図１】  パラオクタデシルフェノールのビーズモデルを示す図
【図２】  ＭＤによる計算結果を示す図（発色状態の顕色剤）
【図３】  ＭＤによる計算結果を示す図（消色状態の顕色剤）
【図１】
0004443797.tif 000005
【図２】
0004443797.tif 000006
【図３】
0004443797.tif 000007