静電場によって誘起される分子の電子密度変化を量子力学計算によって求めておき、分子の多中心双極子分極率および超分極率を最適化法によって決定する方法を開発した。この方法では、分子集合体内のイオンなどによって形成される不均一な静電場を、分子表面に置いた場の点電荷によって代表させ、分子の電子密度変化を分子表面における一電子ポテンシャルの変化として評価する。分子の多中心双極子分極率と超分極率には、原子上または結合軸に沿った誘起双極子を用い、分子表面の一電子ポテンシャルを再現するように最適化する。量子力学計算は、電子相関を考慮した非経験的分子軌道法を用い、基底関数には分子分極率の実験値を再現するセットを使用した。分子軌道計算から求めた一電子ポテンシャルを参照データとして、多中心双極子分極率および超分極率を決定する非線型最適化プログラムの開発を行なった。このプログラムの開発によって、分子一般について多中心分極率と超分極率を決定することが可能となった。水分子、飽和炭化水素、アルコールおよびベンゼンに適用し以下の成果を得た。 (1)原子上に静電場方向の誘起双極子を置き、更に、結合軸に沿った誘起双極子を置くことによって、簡単ではあるが一電子ポテンシャルを非常に良く再現することのできるモデルを見つけることができた。このモデルでは、参照の量子力学一電子ポテンシャルからの平均二乗誤差は、10%以下であった。(2)一連の飽和炭化水素とアルコールの計算結果から、-CH3と-CH2基のC、HとC-Hおよび水酸基のO、HとO-Hにおいて、分子をこえた多中心分極率と超分極率パラメータの移行性が観測された。比較的小さい分子で得られたパラメータを大型分子に用いることが可能となる。これらの成果は、分子集合体のシミュレーションで用いる誘起分極ポテンシャルを開発する上で極めて重要なものである。
|