(57)【要約】【課題】分子結晶や高分子結晶を対象とした分光スペクトルの帰属同定において、結晶相に特徴的なコンフォメーションに起因するピークを同定するための手段を提供する。【解決手段】分子の構造最適化を実行する際、分子内コンフォメーション等の一部の構造パラメータを保存して構造パラメータの一部のみ最適化することを特徴とする、分子軌道法による分光スペクトルピーク位置計算システム。【選択図】　図２
2004012317.tif 000002
【特許請求の範囲】
【請求項１】
分子構造を入力する入力手段と入力された構造を基に分光スペクトルピーク位置計算を行う演算手段と計算結果を出力する出力手段とを具備する分光スペクトルピーク位置計算装置において、前記分光スペクトルピーク位置計算を行う演算手段は、少なくとも分子構造の一部のみを構造最適化する演算ステップとその後に分子軌道法によって分光スペクトルピーク位置計算を行う演算ステップとを実行することを特徴とする分光スペクトルピーク位置計算装置。
【請求項２】
前記分子構造の一部のみを構造最適化する演算ステップは、結合角および二面角を固定し、結合長のみを最適化するものである請求項１記載の分光スペクトルピーク位置計算装置。
【請求項３】
分光スペクトルピーク位置を計算する対象が高分子であって、該高分子若しくはその部分構造が入力手段によって入力され、これについての分子構造の一部のみを構造最適化する演算ステップは、主鎖の結合角および二面角を固定し、主鎖の結合長のみを最適化するものである請求項１記載の分光スペクトルピーク位置計算システム。
【請求項４】
分光スペクトルピーク位置を計算する対象が高分子であって、該高分子若しくはその部分構造が入力手段によって入力され、これについての分子構造の一部のみを構造最適化する演算ステップは、主鎖の結合長及び側鎖を定義する結合長、結合角、二面角の全部または一部を最適化するものであることを特徴とする請求項１記載の分光スペクトルピーク位置計算装置。
【請求項５】
分子構造の一部のみを構造最適化する演算ステップに用いる計算手法が非経験的分子軌道法または密度汎関数法であることを特徴とする請求項１から４のいずれかに記載の分光スペクトルピーク位置計算装置。
【請求項６】
分光スペクトルピーク位置計算を行う演算ステップに用いる計算手法が非経験的分子軌道法または密度汎関数法を用いてされることを特徴とする請求項１から５のいずれかに記載の分光スペクトルピーク位置計算装置。
【請求項７】
核磁気共鳴化学シフトを求めるものである請求項１から６のいずれかに記載の分光スペクトルピーク位置計算装置。
【請求項８】
分子構造の入力を実行するための手順と入力された分子構造の一部のみを構造最適化するための手順とその後に分子軌道法若しくは密度汎関数法によって分光スペクトルピーク位置計算を行う手順と計算結果を出力する手順とを実行させるためのプログラム。
【請求項９】
請求項１記載の装置を用いて分子の分光スペクトルピーク位置を計算することを特徴とする分光スペクトルピーク位置の計算方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、分光スペクトルのピーク位置計算技術に関する。特に、分子結晶や高分子結晶に適して分光分析装置等による測定で得られた分光スペクトルデータを基にスペクトル成分を同定するための分光スペクトルピーク位置計算装置に関するものである。
【０００２】
【従来の技術】
通常、材料計測で得られる分光スペクトルは、特定の成分や構造に特有の幾つかの孤立ピーク波形の重なりとして現れる。そこで、試料に関する分光スペクトルを測定し、分光スペクトル中の各スペクトル成分の波長および強度を基に、試料に含まれる物質や構造を同定し、それらの量を測定することが広く行われている。広く研究された既知の物質や構造であれば、それらに起因する波長や強度は既に知られており、同定に利用される。
【０００３】
しかしながら、物質が既知でない場合や、物質そのものは既知であるがその微細構造について詳細な情報が得られない場合、物質の研究が過去に十分行われておらず分光スペクトルに関して信頼できる情報が得られない場合など、既存の情報を利用した同定が不可能な場合がある。このような場合、分子軌道法の理論計算による同定が有力な手段として用いられる。例えば、物質の分子構造を基に分子軌道法により遷移エネルギー、振動子強度、外場応答量などを計算し、測定された分光スペクトルと比較することで、スペクトル同定を行うことができる。近年の理論化学の進歩により、分子軌道法によるスペクトル計算は、妥当な分子座標さえ与えれば極めて高い精度でスペクトル同定に必要なデータを計算することが可能である。このような技術については、分子軌道法計算プログラム・パッケージＧａｕｓｓｉａｎの開発元であるＧａｕｓｓｉａｎ社のホームページで詳細な情報を容易に入手できる。例えば、核磁気共鳴スペクトルピーク位置計算については、例えば、ｈｔｔｐ：／／ｗｗｗ．ｇａｕｓｓｉａｎ．ｃｏｍ／ｎｍｒｃｏｍｐ．ｈｔｍの文書およびその参考文献が参考になる。
【０００４】
分子軌道法によるスペクトルピーク位置計算に必要な分子座標は、通常、Ｘ線回折や中性子線回折などの実験によるものか、あるいは理論計算による構造最適化計算の結果が用いられる。
【０００５】
しかし、結晶中の原子座標を得るための手段として広く用いられるＸ線回折では、観測手法の特性のため、原子座標にわずかな誤差があり、特に重原子と水素原子との間の結合距離が短く見積もられている場合がほとんどである。ここで、重原子とは水素原子・ヘリウム原子以外の炭素・窒素・酸素等の原子を指す。ところが、分子軌道法によって計算される分光スペクトルピーク位置は、計算に用いた分子座標、特に結合距離の微小な差違に極めて敏感に影響されるため、Ｘ線回折により求めた原子座標に基づいて分子軌道法によって分光スペクトルピーク位置計算を行っても、実測結果と照合するのに十分と言える精度の計算結果を得るのは難しい。そのため、精度の良い計算結果を得るには、理論計算による構造最適化が必要な場合がほとんどである。
【０００６】
一方、理論計算による構造最適化は孤立分子を扱うものであり、分子結晶中の状態を再現するものとは言い難い。特に、通常の構造最適化ではコンフォメーションに関する情報が失われるため、分子結晶や高分子結晶における特定の結晶相に固有なコンフォメーションに起因するスペクトル成分を同定することができず、分光スペクトルの帰属同定に十分活用することが困難であった。
【０００７】
これまで、分子軌道法を用いて分光スペクトルの同定を行うための技術として、例えば特開２００１−１５９６０３号公報が知られている。これは、分子軌道法を用いて計算したピーク位置およびピーク強度を用いてガウス型ないしローレンツ型のモデル関数を作成し、実測で得たスペクトルとの残差平方和を最小化するようにパラメータを決めることによりスペクトルの同定を行う技術である。しかし、この技術は構造情報から分子結晶や高分子結晶の分光スペクトルピーク位置を計算するために用いるには、分子軌道法を用いてピーク位置を計算する過程における計算精度やコンフォメーションについて十分に考慮されていないために分子結晶や高分子結晶の分光スペクトルの帰属同定に用いることは困難である。
【０００８】
【発明が解決しようとする課題】
本発明は、上記の問題を解決するためになされたものであり、分子結晶や高分子結晶について特定の結晶相に固有なコンフォメーションに起因するスペクトル成分を帰属同定でき、しかも核磁気共鳴分光スペクトルや可視分光スペクトル、紫外分光スペクトル、真空紫外分光スペクトルに限らず赤外分光スペクトルやラマン分光スペクトルなど広い分野の分光スペクトルに適用できるような分子軌道法による分光スペクトルピーク位置計算システムを提供することを目的とする。
【０００９】
【課題を解決するための手段】
前記の課題を解決するための本発明は以下のような構成を本旨とする。すなわち、
（１）分子構造を入力する入力手段と入力された構造を基に分光スペクトルピーク位置計算を行う演算手段と計算結果を出力する出力手段とを具備する分光スペクトルピーク位置計算装置において、前記分光スペクトルピーク位置計算を行う演算手段は、少なくとも分子構造の一部のみを構造最適化する演算ステップとその後に分子軌道法によって分光スペクトルピーク位置計算を行う演算ステップとを実行することを特徴とする分光スペクトルピーク位置計算装置、
（２）前記分子構造の一部のみを構造最適化する演算ステップは、結合角および二面角を固定し、結合長のみを最適化するものである前記（１）記載の分光スペクトルピーク位置計算装置、
（３）分光スペクトルピーク位置を計算する対象が高分子であって、該高分子若しくはその部分構造が入力手段によって入力され、これについての分子構造の一部のみを構造最適化する演算ステップは、主鎖の結合角および二面角を固定し、主鎖の結合長のみを最適化するものである前記（１）記載の分光スペクトルピーク位置計算システム、
（４）分光スペクトルピーク位置を計算する対象が高分子であって、該高分子若しくはその部分構造が入力手段によって入力され、これについての分子構造の一部のみを構造最適化する演算ステップは、主鎖の結合長及び側鎖を定義する結合長、結合角、二面角の全部または一部を最適化するものであることを特徴とする前記（１）記載の分光スペクトルピーク位置計算装置、
（５）分子構造の一部のみを構造最適化する演算ステップに用いる計算手法が非経験的分子軌道法または密度汎関数法であることを特徴とする前記（１）〜（４）のいずれかに記載の分光スペクトルピーク位置計算装置、
（６）分光スペクトルピーク位置計算を行う演算ステップに用いる計算手法が非経験的分子軌道法または密度汎関数法を用いてされることを特徴とする前記（１）から（５）のいずれかに記載の分光スペクトルピーク位置計算装置、
（７）核磁気共鳴化学シフトを求めるものである前記（１）から（６）のいずれかに記載の分光スペクトルピーク位置計算装置、
（８）分子構造の入力を実行するための手順と入力された分子構造の一部のみを構造最適化するための手順とその後に分子軌道法若しくは密度汎関数法によって分光スペクトルピーク位置計算を行う手順と計算結果を出力する手順とを実行させるためのプログラム、
（９）請求項１記載の装置を用いて分子の分光スペクトルピーク位置を計算することを特徴とする分光スペクトルピーク位置の計算方法、である。
【００１０】
【発明の実施の形態】
本発明の特徴は、分子軌道法を用いて分光スペクトルピーク位置の計算が実行される前に構造最適化を行い、しかも構造最適化ステップが分子の構造の一部のみを構造最適化する部分構造最適化を行う（以下単に部分構造最適化とも言う）ことに特徴がある。これにより、分子結晶や高分子結晶について特定の結晶相に固有なコンフォメーションに起因するスペクトルピーク位置を精度良く計算でき、実測によるスペクトル成分の帰属同定が容易となるという効果を有する。かかる部分構造最適化を実行する処理によって、特定の結晶相に固有なコンフォメーションを保存しつつ、分子構造をスペクトル計算に必要な精度にまで最適化することが可能となるためである。
【００１１】
本発明は、核磁気共鳴分光スペクトルや可視分光スペクトル、紫外分光スペクトル、真空紫外分光スペクトルに限らず赤外分光スペクトルやラマン分光スペクトルなど広い分野の分光スペクトルに広汎に適用できる。なぜなら、一定の構造データからこれらスペクトル計算を行うための分子軌道計算手順は、有償または無償の分子軌道計算プログラム・パッケージにて開示されているためである。
【００１２】
以下、本発明の一実施形態について詳細に説明する。
【００１３】
初めに、本発明に用いる部分構造最適化について説明する。
【００１４】
分子軌道法における分子座標の指定には、一般にＺ−ｍａｔｒｉｘと呼ばれる方法が好適に用いられる。この方法では、まず第１の原子の位置を座標の原点とする。次に第２の原子の位置は第１の原子に対する結合長で定義され、同時にこれら２つの原子のなすベクトルが分子座標系のｚ軸の方向として定義される。第３の原子の位置は第２の原子に対する距離および第３、第２、第１の原子とのなす結合角を用いて定義され、同時にこれら３つの原子のなす平面が分子座標系のｚｘ面として定義される。第４の原子の位置は、第３の原子との距離、第４、第３、第２の原子とのなす角度、第４、第３、第２、第１の原子とがなす二面角という３つのパラメータによって定義される。ここで二面角とは、第４、第３、第２の原子のなす平面と第３、第２、第１の原子のなす平面との角度として定義されるパラメータである。第５の原子以降の原子の位置は、第４の原子と同様に結合長・結合角・二面角の３つのパラメータでもって定義される。
【００１５】
この他、デカルト座標を用いて分子座標を指定する方法も考えられるが、分子のコンフォメーションを規定する際にはＺ−ｍａｔｒｉｘを用いるほうが容易であり、より好ましく用いられる。
【００１６】
このような構造の入力手段については、従来公知の入力手段を採用しえ、通常コンピュータに入力するためのハードウェアおよびソフトウェアを用いうる。
【００１７】
次に、本発明の計算装置は、入力された構造に基づいて分光スペクトルピーク位置計算を用いる演算手段を具備する。この演算手段に係るプログラムには、もちろん他のステップが含まれていても構わないことは言うまでもないが、大きく２つのステップを実行する。すなわち、分子構造の一部のみを構造最適化するステップと該ステップよりも後に分子軌道法によって分光スペクトルピーク位置計算を行うステップであり、例えば、かかるステップを行う手順を実行するためのプログラムが記憶された演算手段として構成されている。
【００１８】
構造最適化は、まず与えられた分子座標に対して分子軌道法により全エネルギーを計算し、求めたエネルギーをこれらの座標パラメータで微分することにより傾きベクトルを計算し、さらに、傾きベクトルの方向に座標パラメータを増減させ、次の分子軌道計算に用いる新たな分子座標を定義する。このサイクルを、傾きベクトルが十分小さな値に収束するまで繰り返すことにより、座標パラメータは最適な値に収束させることで最適構造が求められる。
【００１９】
通常の構造最適化では、全ての座標パラメータを最適化の対象とするが、本発明においては、座標パラメータの一部のみを最適化の対象とするものである。この手順を部分構造最適化と言う。かかる手法を用いることにより、コンフォメーションを固定して結合長のみを最適化する、あるいは反対に結合長を固定してコンフォメーションのみを最適化する、さらにはポリマー構造における主鎖を固定して側鎖のみ構造最適化するといった処理が可能である。
【００２０】
構造最適化に用いる理論計算手法は、分子力学法、半経験的分子軌道法、非経験的分子軌道法、密度汎関数法などが考えられるが、精度の高い分子構造を得られる非経験的分子軌道法あるいは密度汎関数法を用いるのが望ましい。これらの手法を用いて構造最適化を行う公知のプログラムにおいて、最適化を行う際に一部の座標パラメータを固定する、すなわち計算の対象としない、手順を実行することで部分構造最適化を行うことが可能である。
【００２１】
部分構造最適化の演算ステップにおいては、結合角および二面角を固定し、結合長のみを最適化する方法が好ましく用いられる。その理由は、既に述べたように、Ｘ線回折によって決定した結晶構造では重原子と水素原子との間の結合距離が短く見積もられている場合がほとんどであるため、結合長を最適化することにより、原子座標の精度を向上させることが可能になるからである。また、結合角および二面角を固定することにより、Ｘ線回折から得られたコンフォメーション情報を保存することが可能となるからである。
【００２２】
分光スペクトルピーク位置を計算する対象が高分子である場合、主鎖のコンフォメーションの影響について情報を得るのが重要であることから、主鎖について結合角、二面角を固定し結合長のみを最適化する方法が好ましく用いられる。さらに加えて側鎖の位置を最適化することにより、分光スペクトルピーク位置の計算精度を向上させ、側鎖の位置の影響について情報を得ることが可能である。
【００２３】
次に、本発明の分光スペクトルピーク位置計算について説明する。分子軌道法による理論計算からは、核磁気共鳴分光スペクトルや可視分光スペクトル、紫外分光スペクトル、真空紫外分光スペクトル、赤外分光スペクトル、ラマン分光スペクトルなど広い分野の分光スペクトルに関する情報を得ることができるが、これらは概ね次のような原理に基づいている。
【００２４】
分子軌道法によると、外部電磁場が分子の波動関数に与える摂動を考慮して分子の外場応答量を計算することができる。ここから、核磁気共鳴分光スペクトルのピーク位置を計算できる他、赤外分光スペクトル、ラマン分光スペクトルのピーク強度などを計算できる。
【００２５】
また、分子軌道法によると、分子の励起状態の波動関数を計算して遷移エネルギーや遷移モーメントを計算することができる。ここから、可視分光スペクトル、紫外分光スペクトル、真空紫外分光スペクトルなどのピーク位置やピーク強度などを計算できる。
【００２６】
また、構造最適化の手法を応用し、全エネルギーの座標パラメータによる２階微分の行列を求めることができる。この行列を用いると基準振動解析と呼ばれる解析が可能になり、ここから赤外分光スペクトル、ラマン分光スペクトルなどのピーク位置を計算できる。
【００２７】
本発明の分光スペクトルピーク位置計算システムの効果を十分に発揮させるためには、分光スペクトルピーク位置計算の手法として、精度の高い非経験的分子軌道法あるいは密度汎関数法を用いるのが好ましい。また、本発明の分光スペクトルピーク位置計算システムは、これらの分光スペクトルのうち、特に分子のコンフォメーションの情報を鋭敏に反映する核磁気共鳴分光スペクトルに対して好ましく用いることができる。
【００２８】
以上に述べた部分構造最適化や分光スペクトルピーク位置・ピーク強度計算は、有償あるいは無償で入手可能な分子軌道計算プログラム・パッケージに組み込まれた機能を利用することで実行可能である。プログラム・パッケージとしては、Ｇａｕｓｓｉａｎ社製のＧａｕｓｓｉａｎ９８、アイオワ州立大学製のＧＡＭＥＳＳなどを用いることができる。
【００２９】
分子結晶や高分子結晶の分光スペクトルピーク位置を計算する場合、初期構造の入力方法として次のようなものが考えられる。
【００３０】
第一に、Ｘ線回折や中性子線回折などによる結晶構造の測定結果が存在する場合、その構造から適当な一部を抜き出し、部分構造最適化のための初期構造とするのが良い。
【００３１】
第二に、実験による結晶構造が利用できず構造式しか存在しない場合、まず通常の構造最適化を行ってエネルギー的に最も安定な構造を求めた後、二面角を変化させて可能なコンフォメーションをいくつか発生させるという方法を取ることにより、コンフォメーション変化と分光スペクトルとの関係について考察することが可能となる。
【００３２】
出力手段は、従来公知の手段が採用しえ、特に限定されるものではない。
【００３３】
【実施例】
以下、実施例により本発明をさらに具体的に説明する。なお、本発明はこの実施例により限定されるものではない。
【００３４】
まず、結晶多形を示す分子“バルビタール”について、ある結晶相の炭素１３核核磁気共鳴分光スペクトルを測定し、ピーク位置の標準データとした。固体核磁気共鳴スペクトルの測定には、測定装置としてＣｈｅｍａｇｎｅｔｉｃｓ社製のＣＭＸ−３００を用い、測定手法としてＣＰＭＡＳ法を用いた。
【００３５】
実施例１
次の手順で“バルビタール”結晶の核磁気共鳴分光スペクトルのピーク位置を計算した。
（１）Ｘ線回折により得られた上記“バルビタール”の結晶構造座標データをＡｃｃｅｌｒｙｓ社製の分子設計システムＣｅｒｉｕｓ２に入力し、結晶構造の可視化を行った。この構造から１つの分子を抜き出し、構造最適化のための分子モデルを作成した。
（２）上記の分子モデルについて、以下に述べる方法で部分構造最適化を行った。プログラム・パッケージとしてＧａｕｓｓｉａｎ社製のＧａｕｓｓｉａｎ９８を用いた。波動関数のキーワードにはＢ３ＬＹＰ／６−３１Ｇ（ｄ）を用い、構造最適化計算のキーワードにはＯＰＴ＝Ｚ−ｍａｔｒｉｘを用いた。ＳＣＦの収束条件や構造最適化計算の力の収束条件など、計算時に必要となるしきい値等のパラメーターは、ソフトウエアの標準値を用いた。分子座標の定義にはＺ−ｍａｔｒｉｘを使用し、可変パラメータとして結合長を、固定パラメータとして結合角および二面角を指定した。Ｚ−ｍａｔｒｉｘの定義には、Ｃｅｒｉｕｓ２のＧａｕｓｓｉａｎ　Ｉ／Ｆモジュールを用いた。
（３）部分構造最適化が終了した分子モデルについて、以下に述べる方法で核磁気共鳴分光スペクトルのピーク位置計算を行った。プログラム・パッケージとしてＧａｕｓｓｉａｎ９８を用いた。波動関数のキーワードにはＨＦ／６−３１＋Ｇ（ｄ，ｐ）を用い、核磁気共鳴分光スペクトルのピーク位置計算のキーワードにはＮＭＲを用いた。ＳＣＦの収束条件など計算時に必要なパラメーターは、ソフトウエアの標準値を用いた。
【００３６】
以上のようにして計算したピーク位置と実測との比較を表１に示す。表１におけるＣ１，　Ｃ２等の記号は、図１の構造式における記号と対応している。表１に示すように、実施例における平均誤差は−３．０ｐｐｍ、最大誤差は−６．６ｐｐｍ、誤差二乗平均は１３．９であった。また、実測で観測されるＣ１とＣ３、Ｃ５とＣ７、Ｃ６とＣ８のピークの分離は計算でも再現されていた。図２に、実測のスペクトルに分子軌道法で計算したピーク位置を縦棒で重ね書きしたものを示す。図２において実測ピークの帰属を示すＣ１，　Ｃ２，　．．．等の記号は、図１の構造式における記号と対応している。実測ピークと計算ピークを結ぶ細い実線は、両者の対応関係を表す。
【００３７】
比較例１
構造最適化を実行せず、その他は実施例１と同様にして核磁気共鳴分光スペクトルのピーク位置計算を行った。結果を表１および図３に示す。
【００３８】
計算されたピーク位置と実測値とを比較したところ、平均誤差は−１４．６ｐｐｍ、最大誤差は−２６．７ｐｐｍ、誤差二乗平均は３０７．７であり、実測値との乖離が極めて大きい。特にアルキル鎖の炭素核（Ｃ５、Ｃ６、Ｃ７、Ｃ８）で誤差が大きい。すなわち、実用的な精度を有しているとは言えないものであった。
【００３９】
比較例２
プログラム・パッケージとしてＧａｕｓｓｉａｎ９８を、構造最適化計算のキーワードにＦＯＰＴを用いて、全ての座標パラメータの最適化を行った。その他は実施例１と同様にして核磁気共鳴分光スペクトルのピーク位置計算を行った。結果を表１および図４に示す。
【００４０】
計算したピーク位置と実測値とを比較したところ、平均誤差は−２．３ｐｐｍ、最大誤差は−６．７ｐｐｍ、誤差二乗平均は１２．４であり、乖離幅の点では比較例１よりも良好ではあるが、実測で観測されるＣ１とＣ３、Ｃ５とＣ７、Ｃ６とＣ８のピークの分離を認めることはできず、構造情報として不十分なものとしてしか得られなかった。
【００４１】
【表１】
2004012317.tif 000003
【００４２】
【発明の効果】
本発明は、次のような効果を有する。
【００４３】
第一に、分子結晶や高分子結晶について特定の結晶相に固有なコンフォメーションに起因するスペクトル成分を精度良く且つピーク情報が豊富である結果として得ることが可能である。
【００４４】
第二に、計算の対象となる分光スペクトルとして、核磁気共鳴分光スペクトルや可視分光スペクトル、紫外分光スペクトル、真空紫外分光スペクトルに限らず赤外分光スペクトルやラマン分光スペクトルなど広い分野に適用可能である。
【図面の簡単な説明】
【図１】実施例の項で、ピーク位置の帰属を説明する“バルビタール”の構造式である。
【図２】実施例１における実測スペクトルと計算ピーク位置との対応を示すスペクトルチャートである。
【図３】比較例１における実測スペクトルと計算ピーク位置との対応を示すスペクトルチャートである。
【図４】比較例２における実測スペクトルと計算ピーク位置との対応を示すスペクトルチャートである。
【図１】
2004012317.tif 000004
【図２】
2004012317.tif 000005
【図３】
2004012317.tif 000006
【図４】
2004012317.tif 000007