結晶構造解析 第2回：原子の『箱』と『配置ルール』を特定する - Indexingと空間群推定
前回の「第1回：粉末パターンから何がわかる？EXPO2014構造解析の全体像」では、なぜ物質の結晶構造を知ることが重要なのか、そして粉末X線回折データからEXPO2014のようなソフトウェアを使ってどのように構造解析を進めていくのか、その全体的な流れとステップの概要をご紹介しました。
第2回となる今回は、いよいよ具体的な解析プロセスに入っていきます。解析全体の流れの中でも、特に最初の関門であり、構造解析の土台を築く上で非常に重要な2つのステップ、「Indexing（指数付け）」と「空間群決定」に焦点を当てて解説します。
これらのステップを一言でいうと、
Indexing: 実験データ（粉末パターン）から、原子が詰まっている「箱（ユニットセル）」の正確な形と大きさを特定する作業。
空間群決定: その「箱」の中で、原子たちがどのような「配置ルール」に従って配置されているかを見つけ出す作業。
にあたります。この「箱」と「ルール」が正確に決まらなければ、その中に原子を正しく配置していくことはできません。
この記事では、これらIndexingと空間群決定がどのような原理に基づいて行われるのか、そしてEXPO2014ソフトウェアがどのようにこれらのステップを実行し、解析者をサポートしてくれるのかを、順に見ていきます。
3. ステップ1：Indexing (格子定数の決定)
3-1. Indexingとは？ - 構造解析の最初の関門
EXPO2014を用いた結晶構造解析の最初の具体的なステップが「Indexing (インデキシング、指数付け)」です。これは、実験で得られた粉末X線回折パターンから、結晶構造の基本単位であるユニットセル（単位格子）の「形」と「大きさ」を決定するプロセスです。ユニットセルは、原子が詰まった「箱」のようなものだと考えると、Indexingはこの「箱」の寸法（辺の長さ a,b,c と、辺同士がなす角度 α,β,γ。これらを合わせて格子定数と呼びます）を特定する作業に相当します。
なぜこのステップが最初に来るのでしょうか？ それは、ユニットセルの形と大きさが分からないと、その中に原子がどのように配置されているかを議論することができないからです。Indexingは、いわば原子の世界の設計図を描くための「用紙サイズと罫線」を決める作業であり、構造解析全体の基礎となります。もしここで間違った格子定数を選んでしまうと、その後の解析ステップ（空間群決定や構造決定）がすべて正しく進められなくなってしまうため、非常に重要なステップです。
3-2. Indexingの原理 - 1次元データから3次元情報へ
粉末X線回折パターンは、横軸が回折角度 2θ、縦軸が回折X線強度という1次元的なグラフです。一方、求めたい格子定数は3次元のユニットセルの情報です。Indexingは、この1次元のパターン情報から、3次元のユニットセル情報を再構築するプロセスと言えます。
その原理を簡単に見てみましょう。
ピーク位置から格子面間隔dへ: 回折パターンに現れるピークの角度 2θ は、結晶中の特定の原子面（格子面）からのX線の反射（回折）に対応しています。有名なブラッグの法則(2d sinθ = n λ)を使うと、各ピークの角度 2θ (実験データ) とX線の波長 λ (既知) から、そのピークに対応する格子面の間隔 d を計算できます。
d 値と格子定数・ミラー指数 (hkl) の関係: 結晶学では、格子面間隔 d は、ユニットセルの格子定数 (a, b, c, α, β, γ) と、その格子面を表す3つの整数ミラー指数 (hkl) を使って、数学的に関係づけられます。この関係式は結晶の対称性（結晶系）によって異なりますが、例えば最も単純な立方晶（a=b=c, α=β=γ=90°）の場合は、以下のようになります。
$${1/d^2=(h^2+k^2+l^2)/a^2}$$
他の結晶系ではより複雑な式になりますが、基本的には「観測された d 値」と「未知の格子定数」、「未知のミラー指数 (hkl)」を結びつける式が存在します。
逆問題としてのIndexing: 粉末回折パターンからは、複数のピークに対応する d 値のリストが得られます。Indexingの仕事は、これらの複数の d 値すべてを、共通の格子定数 (a, b, c, α, β, γ) を用いて矛盾なく説明できるような、ミラー指数 (hkl) の組み合わせを見つけ出すことです。これは一種の逆問題であり、ピースに絵柄のないジグソーパズルを解く、または1枚の写真から 3Dモデルを生成するような作業になります。
3-3. EXPO2014での実行プロセス
EXPO2014では、このIndexingプロセスを効率的に進めるための機能が備わっています。
ピーク検出 (Peak Search): まず、解析に使用する回折ピークの位置 (2θ) を回折パターンから正確に読み取る必要があります。EXPO2014は、パターン中のバックグラウンドを除去し、ピーク形状を解析することで、ピーク位置を自動で検出する機能を持っています。ただし、ノイズが多い場合や弱いピーク、重なったピークなどでは誤検出もありうるため、自動検出されたピークリストが妥当かどうかをユーザーが確認し、必要に応じて手動でピークを追加・削除・修正することが重要です。ピーク位置の精度はIndexingの成功に直結します。
Indexingプログラムの実行: 精度の高いピーク位置（d値のリストに変換したもの）が得られたら、それを入力としてIndexingプログラムを実行します。EXPO2014には、実績のある複数のIndexingアルゴリズムが組み込まれており、ユーザーは選択して実行できます。
N-TREOR: EXPOシリーズで開発・改良されてきたプログラム。
DICVOL14: 広く使われている実績のあるプログラム。
McMaille: 別のアルゴリズムに基づくプログラム。
これらのプログラムは、入力されたd値リストを説明できる格子定数の候補を、様々な結晶系（立方晶から三斜晶まで）を仮定しながら探索します。
実験データからのノイズやバックグラウンドの除去で、たとえ回折パターンが既知である物質で構造解析の練習するときでも、その後の解析難易度が大きく変わります。実際の現場ではEXPO2014だけでなく、別のソフトウェアも使いながら慎重にデータの前処理をしています。
特に、未知物質のピークサーチは、単相かどうかもわからないので、どのピークが対象の結晶相のものか判定するのに試行錯誤の沼にはまります。解析だけでなく、合成条件をちょっと変えて実験をやり直してピークの強度の変化を細かく見ていったりもします。今後執筆するチュートリアル回で、実際のコツを紹介していく予定です。
3-4. 解の見極め方 - 正しい「箱」を選ぶ
Indexingプログラムを実行すると、多くの場合、一つの格子定数だけでなく、複数の候補解が出力されます。ここから最も確からしい、真の格子定数を選び出す必要があります。そのための判断材料として、以下の点を考慮します。
Figure of Merit (FOM): 各候補解の「良さ」を示す評価指標です。EXPO2014では、主に以下のFOMが表示されます。
M20 (de Wolff figure of merit): 最も古典的で広く使われる指標。一般的に大きいほどよく、M20が10以上であれば、その解は信頼性が高いとされますが、データの質によってはそれ以下でも正しい場合や、高くても間違っている場合もあります。
FN (Smith-Snyder figure of merit): 大きいほどよく、異なる結晶系同士でも比較しやすい工夫がなされています。
FOMnew (N-TREOR独自)、McM20 (McMaille独自)など、プログラム固有のFOMもあります。
未指数付けピーク数 (NIX: Number of Unindexed Lines): 提案された格子定数では説明できない（=ミラー指数 (hkl) が割り当てられない）ピークが、入力したピークリストの中にいくつ残っているかを示します。この数が少ないほど、良い解である可能性が高まります（理想はゼロ）。
パターンとの視覚的な比較: 候補の格子定数が正しいと仮定した場合に、どのような角度 (2θ) に回折ピークが現れるはずかを計算し、それを実測の回折パターンの下に並べて表示（ティックマーク）してみます。ティックマークで示された計算上のピーク位置が、実測パターンのピークとよく一致しているか、逆に、実測パターンにピークがあるのにティックマークがない（未指数付け）、あるいはティックマークがあるのに実測ピークがない（系統的消滅の可能性、または解が間違い）といった点を目で見て確認することが非常に重要です。
複数プログラムの結果比較: N-TREOR, DICVOL, McMailleなど、異なるアルゴリズムで得られた結果を比較し、共通して見つかる解があれば、その信頼性はより高いと考えられます。
これらの情報を総合的に判断し、最も確からしい格子定数（ユニットセル）を決定します。
M20は低角度側から20本のピークにおいてFoMを計算します。これは、低角度（＝大きなd値）ほど隣接反射との間隔が広く、重なりなどの危険が小さいため、「信頼できる情報」であるという設計思想に由来しています。
FNは低角度に限らず指数付けに成功した全N本を評価に使います。M20は立方晶などの高対称なものは高く、三斜晶のように低対称なものは低くなるという傾向があります。FNでは、このバイアスを「理論的に可能な反射を何%説明できたか」という無次元の網羅率を導入することで低減させています。
M20がざっとふるいにかける一次審査、FNが上位を並べ替える二次審査という使い分けをすることもあります。
3-5. Indexingを成功させるためのTips
Indexingは時として非常に難しいステップであり、自動で実行しても正しい解が得られない場合があります。成功率を上げるために、以下の点を意識すると良いでしょう。
ピーク選択の精度: 何よりも正確なピーク位置 (2θまたはd値) リストを入力することが重要です。特に低角側のピークは重なりが少なく指数付けに有効ですが、ピーク位置がブロードだったり非対称だったりする場合は、かえって解析を妨げる可能性もあります。
不純物ピーク: もし試料に不純物が含まれていると、その回折ピークがIndexingを妨害します。不純物が少量であれば、強度の弱いピークを除外したり、逆に強度の強いピークだけ（例えば上位20本程度）を選んで解析したりすることで、うまくいく場合があります。N-TREORなどは、少数の不純物ピークが含まれていても頑健に動作するように設計されています。
試行錯誤: 最初から完璧なピークリストを用意するのは困難です。うまくいかない場合は、疑わしいピーク（非常に弱い、ショルダーになっているなど）をリストから除外したり、逆に検出閾値を下げて弱いピークを追加したりしながら、何度か試行錯誤することが有効です。
3-6. まとめと次のステップ
Indexingは、粉末X線回折パターンという1次元の情報から、結晶構造の基本となる3次元のユニットセルの形と大きさを決定する、構造解析の最初の重要なステップです。EXPO2014はこのステップを強力に支援しますが、その原理を理解し、ピーク選択や解の評価を適切に行うことが成功の鍵となります。
無事に最も確からしい格子定数が決定できたら、次はその「箱」の中に原子がどのような対称性を持って配置されているのかを探る「ステップ2：空間群の決定」へと進みます
4. ステップ2：空間群の決定
4-1. 空間群とは？ - 原子配置の対称性ルール
前のステップ「Indexing」で、原子が詰まる基本の「箱」であるユニットセルの形と大きさが決まりました。次のステップは、その箱の中で原子がどのような「配置ルール」(＝対称性)で配置されているかを明らかにすることです。この原子配置の対称性を記述するのが「空間群 (Space Group)」です。
結晶の中では、原子はランダムに存在するのではなく、多くの場合、ある原子（または原子のグループ）を特定の操作（回転、鏡映、反転など）で動かすと、別の原子（または原子のグループ）とぴったり重なるような、規則正しい対称性を持っています。空間群は、これらの対称操作の組み合わせを数学的に表現したもので、全部で230種類に分類されています。
なぜ空間群を決める必要があるのでしょうか？ それは、空間群が分かると、ユニットセル全体の中で、原子の座標を独立に決定する必要がある最小限の範囲（非対称単位、Asymmetric Unit） がどこであるかが分かるからです。例えば、ある対称操作によってある原子(#1)と別の座標にある原子(#2)が互いに移り変わる関係にある場合、原子#1の座標さえ決まれば、原子#2の座標は自動的に決まります。このように、空間群を知ることで、決定すべき原子座標の数を大幅に減らすことができ、構造解析を効率的に進めることができます。また、後述する構造決定（特にDirect Methods）の手法自体も、空間群の情報に基づいて行われます。
4-2. 空間群決定の手がかり - 消滅則
空間群を特定するための最も重要な手がかりは、X線回折パターンに現れる「系統的消滅 (Systematic Absence)」、通称「消滅則」に起因する現象です。
結晶が特定の対称操作、具体的には並進成分を含む対称操作（らせん軸、映進面） を持っている場合、その対称性によって特定のミラー指数 (hkl) を持つ回折X線の強度が、完全にゼロになってしまうことがあります。例えば、「a軸に平行な2₁らせん軸（180°回転して軸方向にユニットセルの1/2だけ進む操作のことです）」が存在すると、(h00) という指数を持つ反射のうち、h が奇数（1, 3, 5, ...）のものは必ず強度がゼロになります。
このように、「どのような (hkl) の反射が系統的に消えているか」という消滅のパターンを調べることで、結晶がどのような並進成分を含む対称操作を持っているかを推定でき、それによって空間群（あるいは、いくつかの空間群をまとめたグループである「消滅則記号 (Extinction Symbol)」）を絞り込むことができます。
4-3. 粉末回折データ特有の難しさ
原理的には消滅則を調べれば空間群が分かりそうですが、粉末X線回折データの場合、事はそう単純ではありません。主な理由は以下の2点です。
ピークの重なり: 粉末パターンでは、異なるミラー指数 (hkl) を持つ複数の反射が、偶然ほぼ同じ回折角度 2θ に現れる「ピークの重なり（オーバーラップ）」が頻繁に起こります。もし、消滅するはずの反射が、消滅しない別の強い反射と完全に重なってしまうと、見かけ上、その反射が存在するように見えてしまい、消滅則を正しく判定できません。
強度の不確かさ: 弱い反射は、バックグラウンドノイズに埋もれてしまったり、強度推定の誤差が大きくなったりしがちです。観測されなかったピークが、本当に強度がゼロ（系統的消滅）なのか、それとも単に弱すぎて検出できなかっただけなのかを区別するのは困難な場合があります。
これらの理由から、粉末回折データから単純に「反射があるかないか」だけで消滅則を判断しようとすると、誤った結論に至る可能性があります。
4-4. EXPO2014のアプローチ：確率的に推定する
そこでEXPO2014では、より洗練されたアプローチを採用しています。単純な有無の判定ではなく、観測された強度データ全体を使って、統計的・確率的に最も確からしい消滅則（ひいては空間群）を推定しようと試みます。
その手順の概要は以下の通りです。
強度データの準備: まず、Le Bail法などを用いて、重なり合ったピークを分解し、個々の反射 (hkl) の積分強度を推定します。この際、強度の推定誤差も考慮されます。
確率計算: Indexingで決定された結晶系（単斜晶、直方晶など）において可能な全ての消滅則記号について、推定された反射強度データがそれぞれの消滅則記号が示す「消滅」および「非消滅」のパターンとどれくらい整合的であるかを、確率論に基づいて計算します。強度の推定が不確かな反射の影響を適切に重み付けすることで、より信頼性の高い確率評価を目指します。
Le Bail(ル・ベイユ)法は原子座標を一切仮定せずに、与えられたユニットセルと空間群の拘束だけを使って各 hkl 反射の積分強度を抽出する方法です。ここまでの内容から、「空間群はまだ決まってないのでは？」と思いましたね？それはその通りです。
EXPO2014では、まず「Indexingの結果選ばれた格子の結晶系が持つ空間群のうち、最も高いLaue対称性(二次元回折パターンの左右上下対称性)を持ち、かつ（格子心や面心による消滅以外に）系統的な消滅則を持たない空間群」を初期値として仮定して各ピークの強度から各hkl反射の積分強度を計算しています。これは、系統的な消滅則を持つ空間群を初期値に選ぶと、「この反射は消えるはず」と決め打ちすることになり、もし間違っていた場合に困ってしまいます。まずは回折があるものとして強度を計算して、そのあとその結晶系に属する「すべての可能な消滅パターン（消滅則記号）」について、「観測された強度データが、その消滅パターンとどれくらい整合的か」を確率として計算・評価します。最後に、この確率計算の結果に基づいて、確率が高い順に消滅則記号と、それに対応する空間群の候補をリストアップしてユーザーに提示します。
全パターンを網羅的に計算するのではなく、安全めな計算から絞り込むことで効率的に計算していく工夫がなされている、ということです。
EXPO2014では、Indexingと空間群推定の間に、ユーザーに化学組成の入力を求める段階があります。これはLe Bail法の後、Wilsonプロットなどで「規格化された構造因子振幅」を計算するときに使います。この規格化の時、実験による回折X線強度から原子の種類による散乱能力の違い（=原子散乱因子と言います。）や熱振動の影響を取り除く必要があります。そのためには、ユニットセル内にどの種類の原子がいくつ含まれているか、つまり化学組成の情報が不可欠なのです。組成が分からなければ、平均的な散乱強度や原子散乱因子を正しく計算できず、規格化計算そのものができません。つまり、組成情報は空間群の確率評価を行うための基礎となる規格化強度を計算するために必要となります
4-5. 結果の解釈と空間群の選択
EXPO2014は、上記の確率計算の結果を、確からしい順に並べた空間群（または消滅則記号）の候補リストとしてユーザーに提示します。
候補リスト: 各行には、候補となる空間群とその元になる消滅則記号、そして計算された確率 (FoM: Figure of Merit) が表示されます。通常、最も確率の高いものが一番上に表示されます。
複数の対応: 一つの消滅則記号に対して複数の空間群が対応することがよくありえます。例えば、「P 1 n 1」という消滅則記号は、P2/n と Pn という二つの空間群に対応します。これらは消滅則だけでは区別できません。
補助情報: 確率だけでなく、ユーザーが判断を下すための補助情報も提供されます。
Nabs: 系統的に消滅していると判断された反射のうち、他の反射と重なっていないものの数。
Nasym: その空間群の場合に、非対称単位に含まれる原子数の推定値（化学組成から期待される原子数と比較できる）。
CSD頻度: その空間群が、既知の結晶構造データベース (The Cambridge Structural Database, CSD) にどれくらいの頻度で登録されているかの情報（一般的か、珍しいか）。
キラル情報: その空間群がキラル（Sohnke群）かどうか。もし解析対象の分子が光学活性（キラル）であることが分かっていれば、キラルな空間群しか取り得ません。
視覚的確認: ユーザーはリスト中の候補を選択し、その消滅則に従って消えるはずの反射の位置を、実際の回折パターン上に表示させて確認することができます。消えるはずの位置に強いピークがあれば、その候補は疑わしいと考え、解析をやり直すこともあります。
これらの情報を総合的に検討し、最も妥当と思われる空間群をユーザーが選択します。
4-6. 注意点と次のステップへ
空間群の決定は、特に粉末データの場合、曖昧さが残ることがあります。
確率の解釈: 計算された確率が最も高い候補が常に正しいとは限りません。特に、複数の候補の確率が近い値を示す場合は、慎重な判断が必要です。
候補の絞り込み: 消滅則だけでは区別できない空間群が複数残る場合（さきほどのP2/n と Pnの例など）、最終的な判断は後の構造決定・精密化ステップの結果を見て行うか、あるいは両方の可能性を試す必要があります。通常は両方を試します。
データが理想的であれば、多くの場合、EXPO2014は非常に確からしい空間群を提案してくれます。適切な空間群が選択できれば、いよいよ結晶構造解析の核心部分である「ステップ3：構造決定」、すなわちユニットセル内の原子の具体的な座標を求める段階に進みます。
今回の記事はここまでです。
