JP6534668B2 - 回折データの解析方法、コンピュータプログラム及び記録媒体 - Google Patents

Description

本発明は、単結晶構造解析法において、回折データを、より簡便かつ効率よく解析する方法、この方法をコンピュータに実行させるコンピュータプログラム、及び記録媒体に関する。

従来、化合物の分子構造を決定する方法として、単結晶Ｘ線構造解析法が知られている。単結晶Ｘ線構造解析法によれば、分子の三次元画像を原子レベルで得ることができるため、この方法は生理活性物質等の機能性物質の研究等において極めて有用である。

単結晶Ｘ線構造解析法においては、単結晶にＸ線を照射し、回折Ｘ線を検出して回折データを収集した後、回折データを解析し、分子の構造を決定する。
回折データを解析する際、通常、検出された回折Ｘ線は位相情報が失われているため、フーリエ合成によって、直接、分子の構造（電子密度分布）を決定することができない。このため、従来、結晶構造モデルを初めに構築し、これを精密化して、回折データに適合する結晶構造を得ることにより、分子の構造を決定するという手法が採られてきた。
例えば、図１に示すように、単結晶Ｘ線構造解析法の一般的な手順によれば、回折データを収集し、空間群を決定した後、初期位相を決定するステップにより結晶構造モデルを構築し、次いで、それを精密化することにより、測定された回折データに適合する結晶構造を得、分子の構造を決定することができる（非特許文献１）。

このように、従来においては、回折データを解析する際は、適切な空間群及び初期位相を決定する必要がある。しかしながら、これらのステップを適切に行うためには、結晶学に関する専門知識が必要になるため、これらのステップを行うことは、結晶学に馴染みがない研究者等にとって大きな負担になっていた。
このため、回折データを、より簡便かつ効率よく解析し得る方法が望まれていた。

化学者のための基礎講座１２ Ｘ線構造解析、編集：社団法人日本化学会、発行所：株式会社朝倉書店、１９９９年３月２０日（初版第１刷）

本発明は、かかる実情に鑑みてなされたものであり、単結晶構造解析法において、回折データをより簡便かつ効率よく解析し得る方法、この方法をコンピュータに実行させるコンピュータプログラム、及び、このコンピュータプログラムを記録した記録媒体を提供することを目的とする。

本発明者らは、上記課題を解決すべく鋭意検討した結果、三次元骨格と、該三次元骨格によって仕切られて形成された、三次元的に規則正しく整列した細孔及び／又は中空とを有し、前記三次元骨格が結晶構造解析法によって解明されている多孔性化合物の単結晶の細孔及び／又は中空内に、構造を決定する化合物の分子が規則的に配列されてなる結晶構造解析用試料を用いて得られた回折データを解析する際、用いた多孔性化合物の単結晶の空間群及び結晶構造に関する回折データを初期値として利用することにより、回折データを、より簡便かつ効率よく解析し得ることを見出し、本発明を完成するに到った。

かくして本発明によれば、下記（１）の回折データの解析方法、（２）〜（５）のコンピュータプログラム、及び（６）の記録媒体が提供される。
（１）三次元骨格と、該三次元骨格によって仕切られて形成された、三次元的に規則正しく整列した細孔及び／又は中空とを有し、前記三次元骨格が結晶構造解析法によって解明されている多孔性化合物の単結晶の細孔及び／又は中空内に、構造を決定する化合物の分子が規則的に配列されてなる結晶構造解析用試料を用いて得られた回折データの解析方法であって、
結晶構造解析用試料の空間群として、前記多孔性化合物の単結晶の空間群と同一の空間群、又は、前記空間群より対称性が低い空間群を選択するステップ（Ｉ）、
前記多孔性化合物の単結晶の結晶構造に関する回折データを初期値として用いて、前記結晶構造解析用試料についての初期構造を決定するステップ（ＩＩ）、及び、
ステップ（ＩＩ）で得られた初期構造を精密化するステップ（ＩＩＩ）、
を含むことを特徴とする回折データの解析方法。
（２）回折データの解析用プログラムであって、（１）に記載の回折データの解析方法をコンピュータに実行させるプログラム。
（３）前記多孔性化合物の単結晶の空間群を基に、より対称性が低い空間群を導出する処理（Ｉ）、
前記多孔性化合物の単結晶の空間群と同一の空間群、及び、処理（Ｉ）で導出された空間群からなる群から選ばれる１の空間群を、前記結晶構造解析用試料についての回折データの解析に用いる空間群として決定する処理（ＩＩ）、
処理（ＩＩ）で決定された空間群及び前記多孔性化合物の単結晶の結晶構造に関する回折データを初期値として用いて、前記結晶構造解析用試料についての回折データを解析する処理（ＩＩＩ）、及び、
処理（ＩＩＩ）で得られた結晶構造を精密化する処理（ＩＶ）、
をコンピュータに実行させる（２）に記載のプログラム。
（４）さらに、処理（Ｉ）で導出された空間群をユーザに提示する処理を含み、処理（ＩＩ）における空間群の決定が、ユーザの判断に基づいて行われるものである、（３）に記載のプログラム。
（５）処理（ＩＩ）における空間群の決定が、あらかじめ設定した規則に従ってコンピュータによって行われるものである、（３）に記載のプログラム。
（６）前記（２）〜（５）のいずれかに記載のコンピュータプログラムを記録したことを特徴とするコンピュータに読み取り可能な記録媒体。

本発明によれば、単結晶構造解析法において、回折データを、より簡便かつ効率よく解析し得る方法、この方法をコンピュータに実行させるコンピュータプログラム、及びこのコンピュータプログラムを記録したことを特徴とするコンピュータに読み取り可能な記録媒体が提供される。
本発明によれば、ゲスト分子を包接する前の多孔性化合物の単結晶の結晶構造に関する回折データ（結晶学データの座標値）を初期構造として解析することで、初期位相を従来法で決定することなく、初期位相を定めることができる。

従来一般的に行われてきたＸ線結晶構造解析の手順を表す図である。 本発明の方法の手順の一例を表す図である。 本発明のプログラムの処理手順の一例を表す図である。 単結晶の細孔が延在する方向を表す図である。 ｓｈｅｌｘｓ（初期構造を得るプログラム）を実行した場合に得られる、錯体部分の構造を示す図である。 本発明の方法により初期構造を与えることにより、グアイアズレンの構造の殆どが既に見えた状態を示す図（構造解析を開始する際に得られる状態図）である。

以下、本発明を、１）回折データの解析方法、２）コンピュータプログラム、及び、３）記録媒体に項分けして詳細に説明する。

１）回折データの解析方法
本発明の回折データの解析方法は、三次元骨格と、該三次元骨格によって仕切られて形成された、三次元的に規則正しく整列した細孔及び／又は中空とを有し、前記三次元骨格が結晶構造解析法によって解明されている多孔性化合物の単結晶の細孔及び／又は中空内に、構造を決定する化合物の分子が規則的に配列されてなる結晶構造解析用試料を用いて得られた回折データの解析方法であって、以下のステップ（Ｉ）〜（ＩＩＩ）を含むことを特徴とする。
（Ｉ）ステップ（Ｉ）
結晶構造解析用試料の空間群として、前記多孔性化合物の単結晶の空間群と同一の空間群、又は、前記空間群より対称性が低い空間群を選択するステップ。
（ＩＩ）ステップ（ＩＩ）
前記多孔性化合物の単結晶の結晶構造に関する回折データを初期値として用いて、前記結晶構造解析用試料についての初期構造を決定するステップ。
（ＩＩＩ）ステップ（ＩＩＩ）
ステップ（ＩＩ）で得られた初期構造を精密化するステップ。
本発明の方法は、例えば、図２に示す手順に従って行うことができる。

〔回折データの収集〕
本発明に用いる回折データは、三次元骨格と、該三次元骨格によって仕切られて形成された、三次元的に規則正しく整列した細孔及び／又は中空とを有し、前記三次元骨格が結晶構造解析法によって解明されている多孔性化合物の単結晶の細孔及び／又は中空内に、構造を決定する化合物（以下、「化合物（Ａ）」ということがある。）の分子が規則的に配列されてなる結晶構造解析用試料を用いて得られたものである。

（ｉ）多孔性化合物の単結晶
多孔性化合物の単結晶は、内部に、三次元骨格と、該三次元骨格によって仕切られて形成された、三次元的に規則正しく整列した細孔及び／又は中空とを有し、前記三次元骨格が結晶構造解析法によって解明されているものである。本発明においては、結晶構造解析法によって解明された、前記単結晶の三次元骨格に関する回折データを初期値として利用する。

前記三次元骨格は、単結晶内部において、三次元的な広がりを有する骨格状の構造体である。三次元骨格は、１若しくは２以上の分子鎖、又は、１若しくは２以上の分子鎖及び骨格形成性化合物によって構成されたものである。
「分子鎖」とは、共有結合及び／又は配位結合によって組織化された組織体をいう。この分子鎖内には、分岐構造や環状構造があってよい。
１の分子鎖によって構成された三次元骨格としては、例えば、「ジャングルジム」状に組織化された骨格が挙げられる。
２以上の分子鎖によって構成された三次元骨格としては、２以上の分子鎖が、水素結合、π−πスタッキング相互作用、ファンデルワールス力等の相互作用により、全体として一つに組織化された骨格、例えば、２つの分子鎖が、「ちえのわ」状に絡みあってなる骨格が挙げられる。このような三次元骨格としては、後述する、多核金属錯体１、２の三次元骨格が挙げられる。

「骨格形成性化合物」とは、分子鎖の一部を構成するものではないが、水素結合、π−πスタッキング相互作用、ファンデルワールス力等の相互作用により、三次元骨格の一部を構成する化合物をいう。例えば、後述する多核金属錯体における骨格形成性芳香族化合物が挙げられる。
「三次元的に規則正しく整列した、細孔及び／又は中空」とは、結晶構造解析によって、細孔や中空を確認することができる程度に乱れなく、規則的に整列している細孔や中空をいう。
「細孔」、「中空」は単結晶内における内部空間を表す。筒状に伸びている内部空間を「細孔」といい、それ以外の内部空間を「中空」という。

細孔の大きさは、細孔が延在する方向に対して、最も垂直に近い結晶面と平行な面（以下、平行面ということがある。）における細孔の内接円（以下、単に「細孔の内接円」ということがある。）の直径と相関がある。内接円が大きければ、細孔も大きくなり、内接円が小さければ、細孔も小さくなる。

「細孔が延在する方向」は、以下の方法により決定することができる。
すなわち、まず、対象の細孔を横切る適当な方向の結晶面Ｘ（Ａ面、Ｂ面、Ｃ面かそれぞれの対角面など）を選ぶ。そして、結晶面Ｘ上に存在し、かつ、三次元骨格を構成する原子を、ファンデルワールス半径を用いて表すことで、結晶面Ｘを切断面とする細孔の断面図を描く。同様に、当該結晶面Ｘと一単位胞ずれた結晶面Ｙを切断面とする細孔の断面図を描く。次に、それぞれの結晶面における細孔の断面形状の中心間を、立体図において直線（一点鎖線）で結ぶ（図３参照）。このとき得られる直線の方向が、細孔が延在する方向である。

また、「細孔の内接円の直径」は、以下の方法により求めることができる。
すなわち、まず、上記と同様の方法により、前記平行面を切断面とする細孔の断面図を描く。次に、その断面図において細孔の内接円を描き、その直径を測定した後、得られた測定値を実際のスケールに換算することで、実際の細孔の内接円の直径を求めることができる。
さらに、前記平行面を、一単位胞分、徐々に平行移動させながら、各平行面における細孔の内接円の直径を測定することで、最も狭い部分の内接円の直径と、最も広い部分の内接円の直径が求められる。

単結晶の細孔の内接円の直径は、２〜３０Åが好ましく、３〜１０Åがより好ましい。

また、細孔の形状が真円とは大きく異なる場合、上記平行面における細孔の内接楕円の短径及び長径から、単結晶の包接能を予測することが好ましい。
単結晶の細孔の内接楕円の長径は、２〜３０Åが好ましく、３〜１０Åがより好ましい。また、単結晶の細孔の内接楕円の短径は、２〜３０Åが好ましく、３〜１０Åがより好ましい。

単結晶の細孔容積は、論文（Ａ）：Ａｃｔａ Ｃｒｙｓｔａｌｌｏｇｒ． Ａ ４６，１９４−２０１（１９９０）に記載の手法により求めることができる。すなわち、計算プログラム（ＰＬＡＴＯＮ ＳＱＵＥＥＺＥ ＰＲＯＧＲＡＭ）により算出したＳｏｌｖｅｎｔ Ａｃｃｅｓｓｉｂｌｅ Ｖｏｉｄ（単位格子内の空隙体積）をもとに「単結晶の体積×単位胞における空隙率」を用いて計算することができる。
単結晶の細孔容積（一粒の単結晶中のすべての細孔の容積）は、１×１０^−７〜０．１ｍｍ^３が好ましく、１×１０^−５〜１×１０^−３ｍｍ^３がより好ましい。

また、単結晶が中空を有する場合、その中空の大きさも、細孔容積と同様に、上記論文（Ａ）に記載の手法により求めることができる。

単結晶は、立方体または直方体形状を有するものが好ましい。その一辺は、好ましくは１０〜１０００μｍ、より好ましくは、６０〜２００μｍである。このような形状、大きさの単結晶を用いることで、良質の結晶構造解析用試料が得られ易くなる。

結晶構造解析用試料の調製に用いる多孔性化合物の単結晶は、その三次元骨格が結晶構造解析法によって解明されているものである。
「三次元骨格が結晶構造解析法によって解明されている」とは、結晶構造解析用試料の調製に用いる単結晶が、既に結晶構造解析が行われたものあることを意味するのではなく、結晶構造解析用試料の調製に用いる単結晶のものと同じ三次元骨格を有する多孔性化合物の単結晶を用いて、既に結晶構造解析が行われたことを意味する。

結晶構造解析用試料の調製に用いる多孔性化合物の単結晶は、三次元骨格（いわゆるホスト分子）のみからなるものであってもよいし、三次元骨格と、細孔及び／又は中空内に、交換可能な分子（いわゆるゲスト分子）とを有するものであってもよい。
また、三次元骨格を解明する際に用いられた多孔性化合物の単結晶と、結晶構造解析用試料の調製に用いる多孔性化合物の単結晶は、三次元骨格に変化がない限り、ゲスト分子の有無やゲスト分子の種類に関して相違していてもよい。

例えば、多孔性化合物の単結晶（１）の細孔にゲスト分子（１）が包接されたものの結晶構造が知られている場合、多孔性化合物の単結晶（１）の細孔にゲスト分子（１）が包接されたもの、多孔性化合物の単結晶（１）の細孔にゲスト分子（２）が包接されたもの、多孔性化合物の単結晶（１）の細孔にゲスト分子が包接されていないもの、のいずれも、結晶構造解析用試料の調製に使用することができる。そして、得られた結晶構造解析用試料を用いて回折データを収集した後、既知の結晶構造解析結果（多孔性化合物の単結晶（１）の細孔にゲスト分子（１）が包接された状態の結晶構造解析結果）を利用して、回折データを解析することができる。具体的には、多孔性化合物の単結晶（１）の細孔にゲスト分子（１）が包接された状態の結晶構造解析結果の原子座標から、多孔性化合物の単結晶（１）の原子座標のみを抽出したファイルを作成し、多孔性化合物の単結晶（１）の細孔にゲスト分子（２）が包接された回折データと共に使用することで、初期位相が予め決まった状態で結晶構造解析を行うことができる。

単結晶は、管電圧が２４ｋＶ、管電流が５０ｍＡで発生させたＭｏＫα線（波長：０．７１Å）を照射し、回折Ｘ線をＣＣＤ検出器で検出したときに、少なくとも１．５Åの分解能で分子構造を決定できるものが好ましい。かかる特性を有する単結晶を用いることで、良質の結晶構造解析用試料が得られ易くなる。

多孔性化合物の単結晶としては、前記多孔性化合物の単結晶の細孔及び／又は中空内に、化合物（Ａ）の分子を、結晶構造解析によって構造を決定することができる程度に乱れなく、前記多孔性化合物の単結晶の細孔及び中空内に規則正しく収容させることができるものであれば、特に限定されない。例えば、多核金属錯体の単結晶や、尿素結晶等が挙げられる。なかでも、細孔や中空の大きさや、細孔や中空内の環境（極性等）を制御し易いことから、多核金属錯体の結晶が好ましい。

多核金属錯体としては、配位性部位を２つ以上有する配位子の複数個、及び中心金属としての金属イオンの複数個を含むものが挙げられる。

配位性部位を２つ以上有する配位子（以下、「多座配位子」ということがある。）は、前記三次元骨格を形成し得るものである限り、特に限定されず、公知の多座配位子を利用することができる。
ここで、「配位性部位」とは、配位結合が可能な非共有電子対を有する、配位子中の原子又は原子団をいう。例えば、窒素原子、酸素原子、硫黄原子、リン原子等のヘテロ原子；ニトロ基、アミノ基、シアノ基、カルボキシル基等の原子団；等が挙げられる。なかでも、窒素原子又は窒素原子を含む原子団が好ましい。
なかでも、配位子の平面性が高く、強固な三次元骨格が形成され易いことから、多座配位子としては、芳香環を有するものが好ましい。
一般的に、配位子の中心から、配位性部位までの距離が長い多座配位子を用いると、相対的に細孔や中空が大きい多核金属錯体の単結晶が得られ、配位子の中心から、配位性部位までの距離が短い多座配位子を用いると、相対的に細孔や中空が小さい多核金属錯体の単結晶が得られる。

また、比較的大きな細孔や中空を有する単結晶を容易に得ることができることから、多座配位子としては、配位性部位を２つ以上有する多座配位子が好ましく、配位性部位を３つ有する配位子（以下、「三座配位子」ということがある。）がより好ましく、３つの配位性部位の非共有電子対（軌道）が擬同一平面上に存在し、かつ、３つの配位性部位が、三座配位子の中心部に対して等間隔放射状に配置されているものがより好ましい。

ここで、「擬同一平面上に存在する」とは、各非共有電子対が、同一平面上に存在する状態の他、若干ずれた平面、例えば、基準となる平面に対して、２０°以下で交差するような平面に存在する状態も含む意味である。
また、「３つの配位性部位が、三座配位子の中心部に対して等間隔放射状に配置されている」とは、配位子の中心部から等間隔で放射状に延びる線上に、３つの配位性部位が前記中心部から略等距離に配置されている状態をいう。

三座配位子としては、例えば、下記式（１）

（式中、Ａｒは、置換基を有していてもよい３価の芳香族基を表す。Ｘ^１〜Ｘ^３は、それぞれ独立に、２価の有機基、 又はＡｒとＹ^１〜Ｙ^３とを直接結ぶ単結合を表す。Ｙ^１〜Ｙ^３は、それぞれ独立に、配位性部位を有する１価の有機基を表す。）で示される配位子が挙げられる。

式（１）中、Ａｒは３価の芳香族基を表す。
Ａｒを構成する炭素原子の数は、通常３〜２２、好ましくは３〜１３、より好ましくは３〜６である。

Ａｒとしては、６員環の芳香環１つからなる単環構造を有する３価の芳香族基や、６員環の芳香環が３個縮合してなる縮合環構造を有する３価の芳香族基が挙げられる。

６員環の芳香環１つからなる単環構造を有する３価の芳香族基としては、下記式（２ａ）〜式（２ｄ）で示される基が挙げられる。また、６員環の芳香環が３個縮合してなる縮合環構造を有する３価の芳香族基としては、下記式（２ｅ）で示される基が挙げられる。なお、式（２ａ）〜式（２ｅ）において、「＊」は、それぞれ、Ｘ^１〜Ｘ^３との結合位置を表す。

Ａｒは、式（２ａ）、式（２ｃ）〜式（２ｅ）で示される芳香族基の任意の位置に置換基を有するものであってもよい。かかる置換基としては、メチル基、エチル基、イソプロピル基、ｎ−プロピル基、ｔ−ブチル基等のアルキル基；メトキシ基、エトキシ基、ｎ−プロポキシ基、ｎ−ブトキシ基等のアルコキシ基；フッ素原子、塩素原子、臭素原子等のハロゲン原子；等が挙げられる。これらの中でも、式（２ａ）又は（２ｂ）で示される芳香族基が好ましく、式（２ｂ）で示される芳香族基が特に好ましい。

Ｘ^１〜Ｘ^３は、それぞれ独立に、２価の有機基、又はＡｒとＹ^１〜Ｙ^３とを直接結ぶ単結合を表す。

２価の有機基としては、Ａｒとともに、π電子共役系を構成し得るものが好ましい。Ｘ^１〜Ｘ^３で表される２価の有機基がπ電子共役系を構成することで、式（１）で示される三座配位子の平面性が向上し、より強固な三次元ネットワーク構造が形成され易くなる。
２価の有機基を構成する炭素原子の数は、２〜１８が好ましく、２〜１２がより好ましく、２〜６がさらに好ましい。

２価の有機基としては、炭素数２〜１０の２価の不飽和脂肪族基、６員芳香環１つからなる単環構造を有する２価の有機基、６員芳香環が２〜４個縮合してなる縮合環構造を有する２価の有機基、アミド基〔−Ｃ（＝Ｏ）−ＮＨ−〕、エステル基〔−Ｃ（＝Ｏ）−Ｏ−〕、これらの２価の有機基の２種以上の組み合わせ等が挙げられる。

炭素数２〜１０の２価の不飽和脂肪族基としては、ビニレン基、アセチレン基（エチニレン基）等が挙げられる。
６員環の芳香環１つからなる単環構造を有する２価の有機基としては、１，４−フェニレン基等が挙げられる。
６員環の芳香環が２〜４個縮合してなる縮合環構造を有する２価の有機基としては、１，４−ナフチレン基、アントラセン−１，４−ジイル基等が挙げられる。
これらの２価の有機基の２種以上の組み合わせとしては、下記のものが挙げられる。

これらの芳香環は、環内に、窒素原子、酸素原子、硫黄原子等のヘテロ原子を含んでいてもよい。
また、２価の有機基は、置換基を有するものであってもよい。かかる置換基としては、Ａｒの置換基として先に示したものと同じものが挙げられる。
これらの中でも、Ｘ^１〜Ｘ^３で表される２価の有機基としては、下記のものが好ましい。

Ｙ^１〜Ｙ^３は、それぞれ独立に、配位性部位を有する１価の有機基を表す。
Ｙ^１〜Ｙ^３で表される有機基としては、Ａｒ、Ｘ^１〜Ｘ^３とともに、π電子共役系を構成し得るものが好ましい。
Ｙ^１〜Ｙ^３で表される有機基がπ電子共役系を構成することで、式（１）で示される三座配位子の平面性が向上し、強固な三次元骨格が形成され易くなる。
Ｙ^１〜Ｙ^３を構成する炭素原子の数は、５〜１１が好ましく、５〜７がより好ましい。

Ｙ^１〜Ｙ^３としては、下記式（３ａ）〜式（３ｆ）で示される有機基が挙げられる。なお、式（３ａ）〜式（３ｆ）において、「＊」は、Ｘ^１〜Ｘ^３との結合位置を表す。

Ｙ^１〜Ｙ^３は、式（３ａ）〜式（３ｆ）で示される有機基の任意の位置に、置換基を有するものであってもよい。かかる置換基としては、Ａｒの置換基として先に例示したものと同様のものが挙げられる。
これらの中でも、式（３ａ）で表される基が特に好ましい。

式（１）で示される三座配位子中の、Ａｒ、Ｘ^１〜Ｘ^３、Ｙ^１〜Ｙ^３を適宜選択することで、単結晶の細孔や中空の大きさを調節することができる。この方法を利用することで、目的の分子を包接し得る大きさの細孔や中空を有する単結晶を効率よく得ることができる。

式（１）で示される三座配位子としては、強固な三次元骨格が形成され易いことから、平面性及び対称性が高く、かつ、π共役系が配位子全体に広がっているものが好ましい。このような三座配位子としては、下記式（４ａ）〜式（４ｆ）で示される配位子が挙げられるが、これらに限定されるものではない。

これらの中でも、式（１）で示される三座配位子としては、上記式（４ａ）で示される２，４，６−トリス（４−ピリジル）−１，３，５−トリアジン（ＴＰＴ）が特に好ましい。

また、多核金属錯体の多座配位子として、市販品を用いることもできる。例えば、２０１２年９月発行のシグマアルドリッチ社パンフレット（材料科学の基礎 第７号−多孔性配位高分子（ＰＣＰ）／金属有機構造体（ＭＯＦ）の基礎）には、ＰＣＰ／ＭＯＦ用配位子およびリンカー用化合物として、ピラジン、１，４−ジアザビシクロ［２．２．２］オクタン、１，２−ジ（４−ピリジル）エチレン、４，４’−ビピリジル、４，４’−ビフェニルジカルボン酸、ベンゼン−１，３−ジカルボン酸、ピラジン−２，３−ジカルボン酸、ピラジン−３，５−ジカルボン酸等が記載されており、これらを多核金属錯体の多座配位子として用いることができる。

多核金属錯体の中心金属としての金属イオンは、前記多座配位子と配位結合を形成して、三次元骨格を形成し得るものである限り特に限定されない。なかでも、鉄イオン、コバルトイオン、ニッケルイオン、銅イオン、亜鉛イオン、銀イオン、パラジウムイオン、ルテニウムイオン、ロジウムイオン、白金イオン等の周期表第８〜１２族の金属のイオンが好ましく、２価の、周期表第８〜１２族の金属イオンがより好ましい。なかでも、大きな細孔や中空を有する単結晶が得られ易いことから、亜鉛（ＩＩ）イオン、コバルト（ＩＩ）イオンが好ましい。

多核金属錯体の中心金属には、前記多座配位子の他に、単座配位子が配位していてもよい。かかる単座配位子としては、塩化物イオン（Ｃｌ⁻）、臭化物イオン（Ｂｒ⁻）、ヨウ化物イオン（Ｉ⁻）、チオシアン酸イオン（ＳＣＮ⁻）等の１価の陰イオン；アンモニア、モノアルキルアミン、ジアルキルアミン、トリアルキルアミン、エチレンジアミン等の電気的に中性の配位性化合物；等が挙げられる。

また、多核金属錯体は、反応溶媒（多核金属錯体の合成に用いた溶媒）、置換溶媒（反応溶媒と置き換えられた他の溶媒をいう。以下にて同じ。）、後述する骨格形成性芳香族化合物を含むものであってもよい。

「骨格形成性芳香族化合物」とは、三次元骨格を構成する分子鎖と相互作用（ただし、共有結合、配位結合を除く。）し、三次元骨格の一部を構成し得る芳香族化合物をいう。
多核金属錯体が骨格形成性芳香族化合物を含むことで、三次元骨格がより強固になり易く、化合物（Ａ）の分子を包接した後であっても、三次元骨格がより安定化する場合がある。

骨格形成性芳香族化合物としては、縮合多環芳香族化合物が挙げられる。例えば、下記式（５ａ）〜式（５ｉ）で示されるものが挙げられる。

多核金属錯体としては、例えば、以下の化合物が挙げられる。
（１）配位子及び金属イオンのみからなる化合物〔多核金属錯体（α）〕
（２）前記多核金属錯体（α）と、骨格形成性芳香族化合物とからなる化合物〔多核金属錯体（β）〕
（３）前記多核金属錯体（α）又は多核金属錯体（β）に、溶媒分子等のゲスト分子が包接されてなる化合物〔多核金属錯体（γ）〕

本発明に用いる多核金属錯体は、化合物（Ａ）の分子をその細孔や中空内に取り込んだ後においても結晶性を失わず、かつ、比較的大きな細孔や中空を有するものが好ましい。

このような特性を有する多核金属錯体は、前記式（１）で示される三座配位子を用いることで、簡便に得ることができる。
前記式（１）で示される三座配位子を用いることで得られる多核金属錯体としては、下記式（６ａ）〜（６ｃ）で示される多核金属錯体が挙げられる。

式（６ａ）〜式（６ｃ）中、Ｍは、２価の、周期表第８〜１２族の金属イオンを表し、Ｘは、１価の陰イオン性単座配位子を表し、Ｌは、前記式（１）で示される三座配位子を表し、「ｓｏｌｖ」は、合成時に用いた溶媒分子等のゲスト分子を表し、「ＳＡ」は、骨格形成性芳香族化合物を表し、ａ、ｂ、ｃは任意の自然数を表す。

式（６ａ）〜式（６ｃ）で示される多核金属錯体において、Ｌとして、前記式（４ａ）で示されるＴＰＴを用いた多核金属錯体は、これまでに、溶媒などのゲスト分子を取り込んだ形での分子構造が単結晶Ｘ線構造解析で決定されており、本発明に用いる多核金属錯体として特に適している。

このような多核金属錯体としては、下記式（７ａ）〜式（７ｄ）で示される多核金属錯体が挙げられる。

式（７ａ）〜（７ｄ）中、「ｓｏｌｖ」、「ＳＡ」、ａ、ｂ、ｃは、前記と同じ意味を表す。
式（７ａ）で示される多核金属錯体としては、特開２００８−２１４５８４号公報、Ｊ．Ａｍ．Ｃｈｅｍ．Ｓｏｃ．２００４，ｖ．１２６，ｐｐ１６２９２−１６２９３に記載の［（ＺｎＩ_２）_３（ＴＰＴ）_２（ＰｈＮＯ_２）_５．５］_ｎ（多核金属錯体１）や、多核金属錯体１中の反応溶媒分子の全部又は一部を置換溶媒に交換したものが挙げられる。
式（７ｂ）で示される多核金属錯体としては、特開２００８−２１４３１８号公報に記載の［（ＺｎＢｒ_２）_３（ＴＰＴ）_２（ＰｈＮＯ_２）_５（Ｈ_２Ｏ）］_ｎ（多核金属錯体２）や、多核金属錯体２中の反応溶媒分子の全部又は一部を置換溶媒に交換したものが挙げられる。
式（７ｃ）で示される多核金属錯体としては、特開２００６−１８８５６０号公報に記載の［（ＺｎＩ_２）_３（ＴＰＴ）_２（ＴＰＨ）（ＰｈＮＯ_２）_３．９（ＭｅＯＨ）_１．８］_ｎ（多核金属錯体３）や、［（ＺｎＩ_２）_３（ＴＰＴ）_２（ＰＥＲ）（ＰｈＮＯ_２）_４］_ｎ（多核金属錯体４）や、これらの多核金属錯体中の反応溶媒分子の全部又は一部を置換溶媒に交換したものが挙げられる。
式（７ｄ）で示される多核金属錯体としては、ＷＯ２０１１／０６２２６０号公報に記載の［（Ｃｏ（ＮＣＳ）_２）_３（ＴＰＴ）_４（ＤＣＢ）_２５（ＭｅＯＨ）_５］_ｎ（多核金属錯体５）や、多核金属錯体５中の反応溶媒分子の全部又は一部を置換溶媒に交換したものが挙げられる。

また、多核金属錯体としては、上記の式（６ａ）〜（６ｃ）で示されるものの他に、多孔性配位高分子（ＰＣＰ）や金属有機構造体（ＭＯＦ）と称される公知の多核金属錯体を用いることもできる。例えば、２０１２年９月発行のシグマアルドリッチ社パンフレット（材料科学の基礎 第７号−多孔性配位高分子（ＰＣＰ）／金属有機構造体（ＭＯＦ）の基礎）には、
［Ｃｕ_２（ｂｚｄｃ）_２（ｐｙｚ）］_ｎ
（「ｂｚｄｃ」は、２，３−ピラジンジカルボン酸を表し、「ｐｙｚ」は、ピラジンを表す。ｎは任意の数を表す。）、
［Ｚｎ_２（１４ｂｄｃ）_２（ｄａｂｃｏ）］_ｎ
（「１４ｂｄｃ」は、１，４−ベンゼンジカルボン酸を表し、「ｄａｂｃｏ」は、１，４−ジアザビシクロ［２．２．２］オクタンを表し、ｎは任意の数を表す。）、
［Ｃｕ（ｄｈｂｐｃ）_２（ｂｐｙ）］_ｎ
（「Ｈ_３ｄｈｂｐｃ」は、４，４’−ジヒドロキシビフェニル−３−カルボン酸を表し、「ｂｐｙ」は、４，４’−ビピリジルを表し、ｎは任意の数を表す。）、
［Ｃｒ（ｂｔｃ）_２」_ｎ
（「Ｈ_３ｂｔｃ」は、１，３，５−ベンゼントリカルボン酸を表し、ｎは任意の数を表す。）等の多核金属錯体が記載されており、本発明においては、これらの単結晶を用いることができる。

多核金属錯体の合成方法は特に限定されず、公知の方法を利用することができる。
例えば、２０１２年９月発行のシグマアルドリッチ社パンフレット（材料科学の基礎 第７号−多孔性配位高分子（ＰＣＰ）／金属有機構造体（ＭＯＦ）の基礎）には、多座配位子等を含有する溶液と、金属イオン等を含有する溶液を混合する溶液法；耐圧容器内に、溶媒、多座配位子、金属イオン等を入れ、耐圧容器を密封した後、溶媒の沸点以上に加熱して水熱反応を行う水熱法；容器内に、溶媒、多座配位子、金属イオン等を入れ、マイクロ波を照射するマイクロ波法；容器内に、溶媒、多座配位子、金属イオン等を入れ、超音波を照射する超音波法；溶媒を用いることなく、多座配位子、金属イオン等を機械的に混合する固相合成法；等が記載されており、これらの方法を用いて、多核金属錯体の単結晶を得ることができる。

これらの中でも、特別の装置等を要しないことから、溶液法が好ましく用いられる。
溶液法としては、例えば、多座配位子の第１の溶媒の溶媒溶液に、金属イオン含有化合物の第２の溶媒の溶媒溶液を加え、このまま、０〜７０℃で、数時間から数日間、静置する方法が挙げられる。

金属イオン含有化合物は、特に制限されない。例えば、式：ＭＸ_ｎで示される化合物が挙げられる。ここで、Ｍは金属イオンを表し、Ｘは対イオンを表し、ｎはＭの価数を表す。

前記Ｘの具体例としては、Ｆ⁻、Ｃｌ⁻、Ｂｒ⁻、Ｉ⁻、ＳＣＮ⁻、ＮＯ_３ ⁻、ＣｌＯ_４ ⁻、ＢＦ_４ ⁻、ＳｂＦ_４ ⁻、ＰＦ_６ ⁻、ＡｓＦ_６ ⁻、ＣＨ_３ＣＯ_２ ⁻等が挙げられる。

用いる反応溶媒（第１の溶媒及び第２の溶媒）としては、ベンゼン、トルエン、キシレン、クロロベンゼン、１，２−ジクロロベンゼン、ニトロベンゼン等の芳香族炭化水素類；ｎ−ペンタン、ｎ−ヘキサン、ｎ−ヘプタン等の脂肪族炭化水素類；シクロペンタン、シクロヘキサン、シクロヘプタン等の脂環式炭化水素類；アセトニトリル、ベンゾニトリル等のニトリル類；ジメチルスルホキシド（ＤＭＳＯ）等のスルホキシド類；Ｎ，Ｎ−ジメチルホルムアミド、ｎ−メチルピロリドン等のアミド類；ジエチルエーテル、テトラヒドロフラン、１，２−ジメトキシエタン、１，４−ジオキサン等のエーテル類；メタノール、エタノール、イソプロピルアルコール等のアルコール類；アセトン、メチルエチルケトン、シクロヘキサノン等のケトン類；エチルセロソルブ等のセロソルブ類；ジクロロメタン、クロロホルム、四塩化炭素、１，２−ジクロロエタン等のハロゲン化炭化水素類；酢酸メチル、酢酸エチル、乳酸エチル、プロピオン酸エチル等のエステル類；水；等が挙げられる。これらの溶媒は一種単独で、あるいは二種以上を組み合わせて用いることができる。

比較的大きな多核金属錯体の単結晶を得たい場合には、前記第１の溶媒と第２の溶媒として、互いに相溶性を有さない（すなわち、２層分離する）ものを用いることが好ましい。例えば、第１の溶媒として、ニトロベンゼン、ジクロロベンゼン、ニトロベンゼンとメタノールの混合溶媒、ジクロロベンゼンとメタノールの混合溶媒を用い、第２の溶媒としてメタノールを用いる方法が挙げられる。
また、上記多核金属錯体１〜５については、それぞれ、上記文献に記載された方法にしたがって合成することができる。

（ｉｉ）結晶構造解析用試料
本発明に用いる結晶構造解析用試料は、前記多孔性化合物の単結晶の細孔及び／又は中空内に、化合物（Ａ）の分子が規則的に配列されてなるものである。
「化合物（Ａ）の分子が、規則的に配列される」とは、化合物（Ａ）の分子が、結晶構造解析によって構造を決定することができる程度に乱れなく、多孔性化合物の単結晶の細孔及び中空内に規則正しく収容されていることをいう。

結晶構造解析用試料は、管電圧が２４ｋＶ、管電流が５０ｍＡで発生させたＭｏＫα線（波長：０．７１Å）を照射し、回折Ｘ線をＣＣＤ検出器で検出したときに、少なくとも１．５Åの分解能で分子構造を決定できるものが好ましい。

結晶構造解析用試料は、化合物（Ａ）の構造を決定することができるものであれば、前記多孔性化合物の単結晶中のすべての細孔及び中空内に化合物（Ａ）の分子が取り込まれている必要はない。例えば、前記多孔性化合物の単結晶中の細孔及び中空内の一部に、化合物（Ａ）の溶媒溶液に用いた溶媒が取り込まれたものであっても良い。

結晶構造解析用試料は、化合物（Ａ）の分子の占有率が１０％以上のものであることが好ましい。
占有率は、結晶構造解析により得られる値であり、理想的な包接状態におけるゲスト分子〔化合物（Ａ）の分子〕の量を１００％としたときの、単結晶中に実際に存在するゲスト分子の量を表すものである。

結晶構造解析用試料は、前記多孔性化合物の単結晶を、化合物（Ａ）を含む溶媒溶液と接触させることにより得ることができる。
化合物（Ａ）の大きさは、化合物（Ａ）が単結晶の細孔及び／又は中空に入り得る大きさのものである限り、特に限定されない。化合物（Ａ）の分子量は、通常、２０〜３，０００、好ましくは１００〜２，０００である。
本発明においては、あらかじめ、核磁気共鳴分光法、質量分析法、元素分析等により、化合物（Ａ）の分子の大きさをある程度把握し、適当な細孔や中空を有する単結晶を適宜選択して用いることも好ましい。

化合物（Ａ）を含む溶媒溶液の溶媒としては、用いる単結晶を溶解せず、かつ、キラル化合物（Ａ）を溶解するものであれば、特に限定されない。

用いる溶媒の具体例としては、ベンゼン、トルエン、キシレン、クロロベンゼン、１，２−ジクロロベンゼン、ニトロベンゼン等の芳香族炭化水素類；ｎ−ブタン、ｎ−ペンタン、ｎ−ヘキサン、ｎ−ヘプタン等の脂肪族炭化水素類；シクロペンタン、シクロヘキサン、シクロヘプタン等の脂環式炭化水素類；アセトニトリル、ベンゾニトリル等のニトリル類；ジメチルスルホキシド（ＤＭＳＯ）等のスルホキシド類；Ｎ，Ｎ−ジメチルホルムアミド、ｎ−メチルピロリドン等のアミド類；ジエチルエーテル、テトラヒドロフラン、１，２−ジメトキシエタン、１，４−ジオキサン等のエーテル類；メタノール、エタノール、イソプロピルアルコール等のアルコール類；アセトン、メチルエチルケトン、シクロヘキサノン等のケトン類；エチルセロソルブ等のセロソルブ類；ジクロロメタン、クロロホルム、四塩化炭素、１，２−ジクロロエタン等のハロゲン化炭化水素類；酢酸メチル、酢酸エチル、乳酸エチル、プロピオン酸エチル等のエステル類；水；等が挙げられる。これらの溶媒は一種単独で、あるいは二種以上を組み合わせて用いることができる。

多孔性化合物の単結晶と、前記化合物（Ａ）を含む溶媒溶液を接触させる方法は特に限定されない。例えば、前記単結晶を、化合物（Ａ）を含む溶媒溶液に浸漬させる方法、前記単結晶をキャピラリーの中に詰めた後、化合物（Ａ）を含む溶媒溶液を、そのキャピラリー内を通過させる方法等が挙げられる。

（ｉｉｉ）回折データの収集
本発明の方法において、収集する回折データとしては、多孔質化合物の単結晶及び結晶構造解析用試料についての、Ｘ線回折データや中性子線回折データ等が挙げられる。
回折データを収集する際は、従来の単結晶構造解析における回折データの収集方法を利用することができる。
具体的には、図１に示す従来の手順において、単結晶の代わりに上記の結晶構造解析用試料をマウントする点を除き、その手順に従って回折データの測定を行い、データを収集することができる。図１中の各ステップの詳細は、例えば、前記非特許文献１に記載のとおりである。
なお、近年の結晶構造解析装置においては、回折データの収集における多くのステップがコンピュータにより自動化されている。本発明の方法においても、自動的に収集された回折データを利用することができる。

〔回折データの解析〕
結晶構造解析用試料の回折データの解析は、図２に示すように、結晶構造解析用試料の空間群として、前記多孔性化合物の単結晶の空間群と同一の空間群、又は、前記多孔性化合物の単結晶の空間群より対称性が低い空間群を選択するステップ（Ｉ）、前記多孔性化合物の単結晶の結晶構造に関する回折データ（結晶学データの座標値）を初期値として用いることにより、初期位相を従来法で決定することなく、結晶構造解析用試料の初期構造を決定するステップ（ＩＩ）、及び、ステップ（ＩＩ）で得られた初期構造を精密化するステップ（ＩＩＩ）を経て行われる。
結晶構造解析用試料の回折データの解析を行う前提として、前記多孔性化合物の単結晶の結晶解析データ（空間群、位相情報、結晶構造情報等）を有していることが必要である。

（ｉ）ステップ（Ｉ）
ステップ（Ｉ）は、結晶構造解析用試料の空間群として、前記多孔性化合物の単結晶の空間群と同一の空間群、又は、記多孔性化合物の単結晶の空間群より対称性が低い空間群を選択するステップである。

このステップは、具体的には次のようにして行うことができる。
まず、反射点の指数付けを行う。これにより、大雑把な結晶格子の格子定数とブラベ格子が決定される。
格子定数とは、結晶軸の長さや軸間角度のことをいい、単位格子の各稜間の角度 α，β，γ と、各軸の長さａ，ｂ，ｃを表す６個の定数である。例えば、［（ＺｎＩ_２）_３（ＴＰＴ）_２］_ｎで表される金属錯体の場合には、ａ＝３５、ｂ＝１５、ｃ＝３１、ｍｏｎｏｃｌｉｎｉｃＣである。
ブラベ格子とは、格子点の配列の対称性（晶系）と格子の型の組み合わせにより分類される結晶格子をいい、１４種のブラベ格子が存在する。
この解析は、Ｘ線解析装置メーカー独自のプログラム（例えば、Ｂｒｕｋｅｒであれば、ＡＰＥＸ等）や、汎用プログラム（ＨＫＬ２０００等）を用いて、コンピュータにより行うことができる。
この解析により、いくつかの候補が挙げられた場合には、親化合物が持つ格子定数（例えば、親化合物が、[（ＺｎＩ_２）_３（ＴＰＴ）_２]_ｎで表される金属錯体である場合には、ａ＝３５、ｂ＝１５、ｃ＝３１、ｍｏｎｏｃｌｉｎｉｃＣである。）に最も近いもの（最も近いもので対称性が低いもの）を選択することにより、格子定数及びブラベ格子を決定することができる。
また、コンピュータによる計算の結果、格子定数の見当がつけられない場合には、親化合物がもつ格子定数（例えば、親化合物が、［（ＺｎＩ_２）_３（ＴＰＴ）_２］_ｎで表される金属錯体である場合には、ａ＝３５、ｂ＝１５、ｃ＝３１、ｍｏｎｏｃｌｉｎｉｃＣである。）を直接入力し、このデータを基に実測データから、ゲスト化合物が包接したことによる誤差分をｒｅｆｉｎｅｍｅｎｔすることで、格子定数及びブラベ格子を決定することができる。

次に、空間群の決定を行う。空間群の決定は、上記で決定した結晶格子の格子定数とブラベ格子の情報を用いて、汎用プログラム（例えば、ＰＬＡＴＯＮやＢｒｕｋｅｒ社のＷＰＲＥＰ）を用いて、コンピュータにより行うことができる。

空間群とは、結晶構造における対称要素の集合によって作られる群をいう。
対称要素とは、対称中心、鏡映面、映進面、回転軸及びらせん軸をいい、これらの対称要素に基づく操作を対称操作という。ある対象物に対称操作を施したときに、その前後を区別することができないという性質を対称性という。

結晶の空間的対称性は、結晶で許される対称要素とブラベ格子の組み合わせにより、２３０種の空間群に分類される。

通常、誤った空間群に基づいて回折データの解析を行うと、構造解析に失敗したり、分子の構造が異常に歪んだりするため、回折データを解析する際は、２３０種の空間群の中から、真の空間群を決定することが求められる。

従来の空間群の決定方法においては、回折Ｘ線の方向や強度を基に計算を行い、さらに消滅則等の結晶学における専門知識が必要であった。
一方、本発明の方法においては、結晶構造解析用試料〔以下、「単結晶（α）」ということがある。〕の三次元骨格が、その調製に用いた多孔性化合物の単結晶〔以下、「単結晶（β）」ということがある。〕の三次元骨格と同じであると考えられるため、単結晶（α）の空間群を決定する際に、単結晶（β）の空間群を利用する。
具体的には、上記のように、単結晶（α）の空間群として、単結晶（β）の空間群と同一の空間群、又は、前記多孔性化合物の単結晶の空間群より対称性が低い空間群〔単結晶（β）の空間群から、任意の対称要素を除くことで得られる空間群〕を選択し、回折データの解析を行う。

例えば、前記多核金属錯体１｛［（ＺｎＩ_２）_３（ＴＰＴ）_２（ＰｈＮＯ_２）_５．５］_ｎ｝の空間群は、Ｃ２／ｃである。したがって、多核金属錯体１を用いて結晶構造解析用試料を作製した場合、基本的には、候補となる空間群としてＣ２／ｃを入力し、解析を行う。もし、ゲスト化合物が対称面上に存在して、偽対象データとなっている場合には、親空間群Ｃ２／ｃの部分群（Ｃｃ、Ｃ２、Ｐ２１、Ｐ−１、Ｐ１）を入力して計算を行うことで、真の空間群に容易にたどり着くことができる。
このように本発明の方法においては、候補となる空間群の数が限定されるため、効率よく真の空間群を決定することができる。
真の空間群であるか否かは、従来と同様、得られた構造に問題が無いか等により判断することができる。
空間群を記述する方法には、ヘルマン・モーガン記号（Ｈｅｒｍａｎｎ−Ｍａｕｇｕｉｎ）とシェーンフリース記号（Ｓｃｈｏｅｎｆｌｉｅｓ）の２つがあるが、どちらであってもよい。

（ｉｉ）ステップ（ＩＩ）
ステップ（ＩＩ）は、単結晶（β）の結晶構造に関する回折データ（結晶学データの座標値）を初期値として用いて、単結晶（α）の回折データを解析して、単結晶（α）の初期構造を決定するステップである。

本発明の方法においては、単結晶（α）と単結晶（β）が、共通の三次元骨格を有するため、単結晶（β）の結晶構造をモデルとして利用することができる。
単結晶Ｘ線結晶構造解析において、結晶構造を得ることは、結晶中に存在する原子周りの電子の密度を記述できる関数（構造因子Ｆ）を求めることと同義である。しかしながら、実測可能な回折点のデータからは複素関数であるＦの大きさ部分しか決定することはできない。完全なＦを記述するには、位相にあたる部分を求める必要があるが、これは、観測データに対し適当な近似値を与え、それから計算的に予測される回折点と実測データがどれくらい合致するかで判断しなくてはならない。

従来の手法では、この適当な近似値（初期位相）が得られない場合、構造解析に至らないという問題（いわゆる初期位相問題）があった。
本発明では、空間群や構造が既知の単結晶（α）（結晶スポンジ）に解析対象化合物を導入する（包接させる）ことで得られた結晶構造解析試料を用いる。そのため、上述した空間群や初期位相に関しては、”適当な近似値”が既知である。なぜなら、単結晶（α）と単結晶（β）は骨格構造がほぼ同じと考えることができるからである。本発明では、その”適当な近似値”を用いて解析を行うことで、初期位相問題が発生することなく、解析を行うことで、初期構造を決定することができる。

単結晶（α）の回折データを解析して、初期構造を決定する方法としては、直説法、重原子法、分子置換法等が挙げられる。これらの方法は、プログラムを用いて実行することができる。
用いるプログラムとしては、直説法、重原子法、分子置換法等により、単結晶（α）の回折データを解析して、初期構造を決定するものであれば、特に限定されない。例えば、ＳＨＥＬＸ、ＳＩＲ、ｓｕｐｅｒｆｌｉｐ、Ｘ−ＰＬＯＲ（モレキュラーシミュレーション社）や、ＡＭＯＲＥ（ＣＣＰ４（Ｃｏｌｌａｂｏｒａｔｉｖｅ Ｃｏｍｐｕｔａｔｉｏｎａｌ Ｐｒｏｊｅｃｔ，Ｎｕｍｂｅｒ４．Ａｃｔａ Ｃｒｙｓｔａｌｌｏｇｒ． Ｄ５０， ６７０−６７３（１９９４））のプログラム群の１つ）等の公知のプログラムが挙げられる。

（ｉｉｉ）ステップ（ＩＩＩ）
ステップ（ＩＩＩ）は、ステップ（ＩＩ）で得られた結晶構造を精密化するステップである。
ステップ（ＩＩＩ）は、従来の、構造の精密化ステップと同じものであり、最小二乗法やフーリエ合成（差フーリエ合成）等を繰り返すことで、測定された回折データに適合する結晶構造を得、分子の構造を決定する（非特許文献１）。
これらの方法は、プログラムを用いて実行することができる。用いるプログラムとしては、最小二乗法やフーリエ合成（差フーリエ合成）等を繰り返すことで、測定された回折データに適合する結晶構造を得、分子の構造を決定することができるものであれば、特に限定されない。ＳＨＥＬＸＬ、ＲＥＦＭＡＣ、Ｘｔａｌ等の公知のプログラムが挙げられる。

従来、回折データを解析する際は、適切な空間群及び初期位相を決定する必要があり、これらのステップを適切に行うためには、結晶学に関する専門知識が必要であった。
しかしながら、結晶構造が既知の多孔性化合物の単結晶を利用する本発明の方法によれば、結晶学に馴染みがない研究者等であっても、回折データを簡便かつ効率よく解析することができる。

上述したステップ（Ｉ）〜（ＩＩＩ）は、コンピュータにより、ステップごとのプログラムを連続的に実行させることにより、実施することができる。

２）コンピュータプログラム
本発明のコンピュータプログラムは、本発明の回折データの解析方法をコンピュータに実行させるものである。
本発明のコンピュータプログラムは、コンピュータにインストールされて、図３に示すように、回折データの収集、候補となる空間群の表示〔処理（Ｉ）〕、データの整理、及び回折データの解析に用いるための空間群の決定〔処理（ＩＩ）〕、単結晶（β）の結晶構造に関する回折データ（結晶学データ等）を初期値として用いることにより、単結晶（α）の初期構造の決定〔処理（ＩＩＩ）〕、初期構造の精密化〔処理（ＩＶ）〕を、連続的に実行するものである。
それぞれの処理ごとに、それぞれの処理を行うプラグラムが実行される。
本発明のコンピュータプログラムは、これらのプログラムの集合体であって、処理（Ｉ）〜（ＩＶ）を連続的に行うものである。

本発明のコンピュータプログラムは、主制御装置（ＣＰＵ）、入出力装置及び記憶装置からなる処理装置にインストールされ、実行されるプログラムである。
主制御装置（ＣＰＵ）は、プログラムを実行して演算処理を行う装置である。
入出力装置は、プログラムが格納された記録媒体の読み取り装置、インターネットとの通信手段、ユーザとのインターフェース（表示画面、キーボード等）を備える。
記憶装置は、展開されたプログラムのデータ、及び、実行中のデータを保管する装置である。

プログラムのデータには、前記処理（Ｉ）〜（ＩＶ）のそれぞれを実行するプログラムデータのほかに、後述するように、親化合物（ゲスト化合物を包接する前の多孔性化合物の単結晶）の結晶構造回折データ（親化合物の空間群、結晶構造）、結晶構造解析の対象となる化合物の回折データ、一般的な空間群に関するデータ等が挙げられる。
本発明のコンピュータプログラムは、例えば、記録媒体（ＣＤＲＯＭ）やインターネットから取得できるものである。

＜処理（Ｉ）＞
処理（Ｉ）は、本発明の方法のステップ（Ｉ）のうち、候補となる空間群の導出を行うものである。すなわち、結晶構造解析用試料の空間群として、前記多孔性化合物の単結晶の空間群と同一の空間群、又は、記多孔性化合物の単結晶の空間群より対称性が低い空間群を選択し、表示画面に候補となる空間群を表示する処理である。

処理（Ｉ）では、使用者が親化合物（ゲスト化合物を包接する前の多孔性化合物の単結晶）の空間群をコンピュータに入力することで、候補となる空間群を表示させる。
例えば、前出の［（ＺｎＩ_２）_３（ＴＰＴ）_２（ＰｈＮＯ_２）_５．５］_ｎ｝の場合には、単結晶の空間群はＣ２／ｃである。使用者が、空間群Ｃ２／ｃを入力すると、親空間群Ｃ２／ｃとともに、部分群（Ｃｃ、Ｃ２、Ｐ２１、Ｐ−１、Ｐ１）も表示される。
この解析は、Ｘ線解析装置メーカー独自のプログラム（例えば、Ｂｒｕｋｅｒであれば、ＡＰＥＸ等）や、汎用プログラム（ＨＫＬ２０００等）を用いて、コンピュータにより行うことができる。
コンピュータには、空間群に関するデータが記憶されたメモリー領域が設けられており、親空間群を入力すると、親空間群とともに親空間群の部分群も表示する機能を付与することができる。
本発明のコンピュータプログラムは、処理（Ｉ）の後、その処理結果をユーザに提示する処理をコンピュータに実行させるものであってもよいし、コンピュータに処理結果をユーザに提示する処理をさせることなく、処理（ＩＩ）を実行させるものであってもよい。

＜処理（ＩＩ）＞
処理（ＩＩ）は、本発明の方法におけるステップ（Ｉ）のうち、空間群を決定するため演算処理を実行するものである。すなわち、前記多孔性化合物の単結晶の空間群と同一の空間群、及び、処理（Ｉ）で導出された空間群からなる群から選ばれる１の空間群を回折データの解析用に決定する処理（ＩＩ）である。
空間群の決定は、上記で決定した結晶格子の格子定数とブラベ格子の情報を用いて、汎用プログラム（例えば、ＰＬＡＴＯＮやＢｒｕｋｅｒ社のＷＰＲＥＰ）を用いて、コンピュータにより行うことができる。
処理（ＩＩ）は、前記多孔性化合物の単結晶の空間群と同一の空間群、及び、処理（Ｉ）で導出された空間群からなる群から選ばれる１の空間群を回折データの解析用に決定するものである。
この決定は、ユーザの判断に基づいて行われるものであってもよいし、あらかじめ設定した規則に従ってコンピュータが行うものであってもよい

ユーザの判断に基づいて空間群を決定する場合、プログラムに、コンピュータが処理（Ｉ）を終えた後、その処理結果をユーザに提示する機能を付加することで、ユーザはその中から１の空間群を決定するようにしてもよい。
コンピュータが空間群を決定する場合、対称性が高い空間群を選択する等の規則をあらかじめ設けておくことで、コンピュータが空間群を決定することができる。例えば、回折データから、ゲスト化合物が対称面上に存在して、偽対象データとなっている可能性がある場合には、親空間群Ｃ２／ｃの部分群（Ｃｃ、Ｃ２、Ｐ２１、Ｐ−１、Ｐ１）を用いて、計算を行うようにプログラムされていてもよい。

＜処理（ＩＩＩ）＞
処理（ＩＩＩ）は、本発明の方法におけるステップ（ＩＩ）を実行するものである。すなわち、処理（ＩＩ）で決定された空間群及び前記多孔性化合物の単結晶の結晶構造に関する回折データを初期値として用いて、結晶構造解析用試料の初期構造を決定するものである。処理（ＩＩＩ）によって、構造の精密化のための初期構造が決定される。

本発明では、初期位相が決定された場合に出現すると想定されるホスト分子（多孔性化合物の単結晶）の骨格構造が予めわかっているので、初期位相を決定するためのプログラムを実行することなく、先の求めた格子定数と空間群の情報を用いて、ダイレクトに結晶構造解析用試料の初期構造を決定することができる。

処理（ＩＩＩ）は、処理（ＩＩ）で決定された空間群及び前記多孔性化合物の単結晶の結晶構造を基に、プログラムを用いて実行することができる。用いるプログラムとしては、処理（ＩＩＩ）を実行できるものであれば、特に限定されない。例えば、ＳＨＥＬＸ、ＳＩＲ、ｓｕｐｅｒｆｌｉｐ、Ｘ−ＰＬＯＲ（モレキュラーシミュレーション社）や、ＡＭＯＲＥ（ＣＣＰ４（Ｃｏｌｌａｂｏｒａｔｉｖｅ Ｃｏｍｐｕｔａｔｉｏｎａｌ Ｐｒｏｊｅｃｔ，Ｎｕｍｂｅｒ４．Ａｃｔａ Ｃｒｙｓｔａｌｌｏｇｒ． Ｄ５０， ６７０−６７３（１９９４））のプログラム群の１つ）等の公知のプログラムが挙げられる。

＜処理（ＩＶ）＞
処理（ＩＶ）は、本発明の方法におけるステップ（ＩＩＩ）を実行するものである。すなわち、処理（ＩＩＩ）で得られた初期構造を精密化するものである。

初期構造の精密化処理では、反射データ（ｈｋｌファイル）と、構造の精密化前の初期構造のデータファイル（ｉｎｓファイル)を用いる。
具体的には、次のようにして実行される。
まず、ホスト化合物（多孔性化合物の単結晶）の骨格にあたる原子団の座標値を、構造の精密化前の初期構造の座標値にあてはめる。実際の作業としては、ホスト化合物（多孔性化合物の単結晶）の骨格にあたる原子団の座標値を、前記初期構造のデータファイル（ｉｎｓファイル)にコピーすればよい。
次に、この座標を初期値として構造の精密化を行うことができる。構造の精密化法としては、フーリエ法、最小二乗法、最尤法等が挙げられる。
この処理は、従来の精密化処理と同様のものであり、プログラムを用いて行うことができる。用いるプログラムとしては、精密化処理を行うことができるものであれば、特に限定されない。例えば、ＳＨＥＬＸＬ、ＲＥＦＭＡＣ、Ｘｔａｌ等の公知のプログラムが挙げられる。

本発明のコンピュータプログラムは、以上のようにして構造精密化して得られた分子構造について、分子全体の投影図、原子間距離、結合角等を表示する機能をさらに有していてもよい。

本発明のコンピュータプログラムによれば、本発明の方法を効率よく実行することができる。したがって、本発明のプログラムを利用することで、結晶学に馴染みがない研究者等であっても、回折データを簡便かつ効率よく解析することができる。

３）記録媒体
本発明の記録媒体は、本発明のコンピュータプログラムを記録したことを特徴とするコンピュータに読み取り可能なものである。
記録媒体としては、フレキシブルディスク（ＦＤ）、ＭＯディスク、ＣＤＲ、ＣＤＲＷ、ＤＶＤ−ＲＯＭ、ＤＶＤ−ＲＡＭ、外付けＨＤＤ、メモリカード、ＵＳＢメモリ、シリコンディスク、ＨＤＤ互換シリコンディスク等が挙げられる。また、本発明のコンピュータプログラムが、複数の記録媒体に分けて記録されたものであってもよい。

以下に、実施例を挙げて、本発明についてより具体的に説明する。ただし、本発明はこれらの実施例によりなんら限定されるものではない。

フムレン（２，６，６，９−テトラメチル−１，４−８−シクロウンデカトリエン）の立体構造は、結晶中で銀イオンと共存させて重原子法によって位相を決定する過程を経て、Ｊ． Ｃｈｅｍ． Ｓｏｃ．Ｂ，１１２−１２０ （１９６６）にて報告されている。

フムレンを、分子構造が既知の多孔性化合物Ａ（例えば、［（ＺｎＩ_２）_３（ＴＰＴ）_２（ＰｈＮＯ_２）_５．５］_ｎが使用できる。）の細孔内に包接させた単結晶を作製し、これを結晶構造解析用試料とした。この結晶構造解析用試料から収集した反射データファイルと、多孔性化合物Ａのみの原子座標を記述したデータファイルを使用し、公知のプログラム（例えば、ＳＨＥＬＸＬが使用できる。）を用いて構造解析を開始すると、新規に位相を決定すること無く、ホスト化合物の構造はフムレンの包接に最適化された原子位置に速やかに変化し、フムレン中のいくつかの原子が初期構造として観察された。
得られた初期構造に関するデータを使用し、公知のプログラムを用いて、構造精密化を行うことで、フムレンの分子構造を容易に決定することができる。
多孔性化合物Ａのみの原子座標を記述したデータファイルは、分子構造が既知の多孔性化合物Ａ、例えば、多孔性化合物Ａが、［（ＺｎＩ_２）_３（ＴＰＴ）_２（ＰｈＮＯ_２）_５．５］_ｎである場合、このものの結構構造解析データから、溶媒である（ＰｈＮＯ_２）に対応するデータを削除する操作を行うことにより、入手することができる。
細孔性錯体にグアイアズレン分子を包摂させた結晶のデータでは、ｓｈｅｌｘｓ（初期構造を得るプログラム）を実行した場合は、錯体部分の構造しか得られない（図５）。一方、本発明の方法により初期構造を与えることで、グアイアズレンの構造の殆どが既に見えた状態で、構造解析を開始することができる（図６）。

１：結晶面Ｘ
２：結晶面Ｙ
３：細孔
４：細孔が延在する方向

Claims (6)

1. 三次元骨格と、該三次元骨格によって仕切られて形成された、三次元的に規則正しく整列した細孔及び／又は中空とを有し、前記三次元骨格が結晶構造解析法によって解明されている多孔性化合物の単結晶の細孔及び／又は中空内に、構造を決定する化合物の分子が規則的に配列されてなる結晶構造解析用試料を用いて得られた回折データの解析方法であって、
   結晶構造解析用試料の空間群として、前記多孔性化合物の単結晶の空間群と同一の空間群、又は、前記空間群より対称性が低い空間群を選択するステップ（Ｉ）、
   前記多孔性化合物の単結晶の結晶構造に関する回折データを初期値として用いて、前記結晶構造解析用試料についての初期構造を決定するステップ（ＩＩ）、及び、
   ステップ（ＩＩ）で得られた初期構造を精密化するステップ（ＩＩＩ）、
   を含むことを特徴とする回折データの解析方法。
2. 回折データの解析用プログラムであって、請求項１に記載の回折データの解析方法をコンピュータに実行させるコンピュータプログラム。
3. 前記多孔性化合物の単結晶の空間群を基に、より対称性が低い空間群を導出する処理（Ｉ）、
   前記多孔性化合物の単結晶の空間群と同一の空間群、及び、処理（Ｉ）で導出された空間群から選ばれる１の空間群を、結晶構造解析用試料についての回折データの解析に用いる空間群として決定する処理（ＩＩ）、
   処理（ＩＩ）で決定された空間群及び前記多孔性化合物の単結晶の結晶構造に関する回折データを初期値として用いて、前記結晶構造解析用試料についての初期構造を決定する処理（ＩＩＩ）、及び、
   処理（ＩＩＩ）で得られた初期構造を精密化する処理（ＩＶ）、
   をコンピュータに実行させる、請求項２に記載のコンピュータプログラム。
4. さらに、処理（Ｉ）で導出された空間群をユーザに提示する処理を含み、処理（ＩＩ）における空間群の決定が、ユーザの判断に基づいて行われるものである、請求項３に記載のコンピュータプログラム。
5. 処理（ＩＩ）における空間群の決定が、あらかじめ設定した規則に従ってコンピュータによって行われるものである、請求項３に記載のコンピュータプログラム。
6. 請求項２〜５のいずれかに記載のコンピュータプログラムを記録したことを特徴とするコンピュータに読み取り可能な記録媒体。

Images