X線結晶構造解析～その４～：X線回折における「位相問題」

位相問題（Phase problem）とは、X線回折実験において、回折強度 I(b) から結晶の電子密度 ρ(r) を直接決定することができないという、回折理論の主要な問題であり、かつあらゆる回折実験における中心的な問題です。

これは、回折実験によって測定されるのは回折強度 I(b)のみであり、この強度I(b)は構造因子（Structure factor）F(b)の大きさの二乗に比例する（I(b) ~ |F(b)|²）という事実に基づいています。構造因子F(b)は結晶の電子密度ρ(r)のフーリエ変換として定義される複素関数です。複素関数である構造因子F(b)は、その大きさ |F(b)| と位相 φ(b) の両方によって完全に記述されます。しかし、実験から得られるのは強度の大きさ|F(b)|だけであり、構造因子の位相φ(b)は実験的に直接得ることができません。

電子密度ρ(r)は、構造因子F(b)を逆フーリエ変換することによって計算できます。

$$\rho(\mathbf{r}) = \frac{1}{V} \sum_{\mathbf{h}} |F(\mathbf{h})| \exp(i\phi(\mathbf{h})) \exp(-2\pi i \mathbf{h} \cdot \mathbf{r})$$

つまり、ρ(r)を求めるためには、構造因子F(b)の大きさだけでなく、その位相φ(b)も必要になります。位相情報が失われているため、実験結果である強度I(b)から、目的とする電子密度ρ(r)への直接的な道筋が存在しないのです。情報間の関係を図示すると、電子密度ρ(r)から構造因子F(b)への変換とその逆はフーリエ変換ですが、回折強度I(b)（≈|F(b)|²）から電子密度ρ(r)へのパスは位相問題によって阻まれています。

この位相問題は、何十年もの間、X線結晶学の進歩を妨げてきました。回折実験で得られる回折点が逆格子ベクトルhに対応することから、すべての単結晶回折実験は逆格子によって記述されます。構造因子F(h)は通常、原子位置r_jと原子散乱因子f_j(h)を用いた級数で表されます。

位相問題へのアプローチや関連する概念がソースにいくつか示されています。

• 直接法は、「電子密度は常に正である（負の領域はない）」や「原子は離散的に存在する」といった物理的制約から、位相間の確率的な関係（三位相和など）を導き出して位相を推定する手法です。H. Hauptman と J. Karle はこの業績で1985年にノーベル化学賞を受賞しました。
• パターソン関数 P(u) は、回折強度I(b)のフーリエ変換として定義され、電子密度の畳み込み二乗にあたります。パターソン関数は、この位相問題に対する最初のアプローチとしてパターソンによって導入されました。
• 対称心を持つ構造（Centrosymmetric structures）の場合、構造因子は実数となり、位相問題は構造因子の符号（sign）を決定する「符号問題」に還元されます。このような構造では、異常分散がない限りフリーデルの法則 I(h) = I(-h) が成り立ちます。
• 異常分散（Anomalous dispersion）が存在する場合、フリーデルの法則は一般的に成り立たなくなります（I(h) ≠ I(-h)）。この現象を利用して、絶対配置（absolute configuration）を決定することが可能です。異常分散は、特に重い原子や特定の波長のX線で顕著になりますが、酸素のような比較的軽い原子でも正確な測定を行えば絶対配置決定に利用できる場合があります。このためには、フリーデル対応の反射（hと-h）を測定する必要があり、これは通常、半天球ではなく全天球のデータを測定することを意味し、データ量が2倍に増加します。ソースでは、異常分散を利用した多波長回折実験がマクロ分子の位相決定（phasing）に利用されることが示唆されています。
• 系統的消滅（Systematic extinctions）は、結晶の並進対称要素（すべり面やらせん軸）によって特定の逆格子ベクトルhに対する構造因子F(h)がゼロになる現象であり、強度測定において特定の反射が観測されない形で現れます。系統的消滅は、位相問題を直接解決するものではありませんが、空間群決定に不可欠な情報であり、空間群情報は構造解析（位相決定を含む）の重要な前提となります。
• 回折データセットの品質も位相問題の解決に影響します。例えば、観測されない反射（"less thans"）があまりに多い場合（例えば50%を超える場合）、位相問題を解決することが困難になる危険性があり、構造を解くことができない可能性があります。これは「完全性 (Completeness)」という指標に関連します。高角側のデータ（高解像度データ）が不足すると、電子密度マップの解像度が落ち、原子の「形」がぼやけてしまい、位相決定のアルゴリズム（特に直接法）が収束しにくくなります。

まとめると、位相問題は、X線回折実験で測定される強度が構造因子の大きさの二乗であるために、構造因子の位相情報が失われるという基本的な問題です。電子密度を得るには構造因子の位相が必要であり、この問題は結晶構造解析における中心的な課題となっています。パターソン関数、対称性の利用（符号問題）、そして異常分散（特にマクロ分子の位相決定や絶対配置の決定）などが、この位相問題に取り組むための手法として関連付けられています。

X線結晶構造解析～その３～：逆格子 (Reciprocal Lattice)の概念と役割

逆格子は、結晶構造解析において非常に重要な概念であり、結晶の直接格子 (direct lattice) とは別の、関連する空間（逆空間 (reciprocal space) と呼ばれます）における格子です。具体的には、結晶の構造やX線回折現象を記述するために導入されます。

1. 定義と関係:

   • 逆格子は、基底ベクトル a*、b*、c* によって定義される空間です。この空間はフーリエ空間 (Fourier space) とも呼ばれます。
   • 逆空間のベクトルは、逆長さの次元を持ちます (例: Å⁻¹)。
   • 逆格子ベクトル h (h = ha* + kb* + lc*) は、直接格子の面と密接な関係があります。具体的には、ベクトル h は直接格子の特定の面（ベクトル1/ha, 1/kb, 1/lcで単位胞ベクトルを横断する面）の法線の方向を持ち、その大きさ ⎪h⎪ はその面から原点までの距離 d の逆数に等しくなります (d = 1/⎪h⎪)。この性質が、算術的な計算において平面の記述よりもベクトルの使用を有利にし、逆格子を導入する理由となります。
   • 逆格子ベクトルの定義式（数学的厳密性）
     
     $$\mathbf{a^*} = \frac{\mathbf{b} \times \mathbf{c}}{V}, \quad \mathbf{b^*} = \frac{\mathbf{c} \times \mathbf{a}}{V}, \quad \mathbf{c^*} = \frac{\mathbf{a} \times \mathbf{b}}{V}$$

     ここで $${V = \mathbf{a} \cdot (\mathbf{b} \times \mathbf{c})}$$

2. 回折における役割:

   • X線回折実験は、逆格子を基に記述することができます。
   • 回折は、エワルド球 (Ewald sphere) 上に、結晶の格子面に対応する逆格子ベクトル h が乗るときに起こります。
   • すべての回折斑点（diffraction spots）は、逆空間における格子ベクトル h に割り当てることができます。
   • 逆格子ベクトル h は、結晶中の回折を記述する上で支配的な役割を果たします。

3. 構造因子との関連:

   • 構造因子 (Structure Factor) F(b) は、逆空間の関数です。その引数である b (または結晶学では通常 h と表記される) は、逆格子ベクトルを表します。
   • 実験では回折X線の強度 I が測定されますが、I は構造因子の絶対値（構造振幅 ⎪F(b)⎪）の二乗に比例します (I(b) ∼ ⎪F(b)⎪²)。
   • 回折実験から得られるのは構造振幅 ⎪F(b)⎪ のみであり、構造因子の位相 ϕ(b) を直接得ることはできません。これが結晶構造解析の中心的な課題である位相問題 (phase problem) です。

4. 実験的な観測:

   • フィルム法 (回転法、ワイゼンベルグ法、プリセッション法) は、逆格子の層や断面を視覚化するのに用いられました。回転写真では層線が逆格子の層に対応し、プリセッション写真では歪みのない逆格子の像が得られます。
   • フィルム上の回折斑点の位置や距離から、逆格子定数 (a*, b*, c*, α*, β*, γ*) を推定し、それから直接格子の格子定数を計算することができました。
   • 現代の自動回折計（点検出器や面積検出器）は、特定の逆格子点 (hkl) における回折強度を自動的に測定します。

5. 制限球と解像度:

   • 特定の波長 λ のX線では、逆格子ベクトル h の大きさが 2/λ より小さい、すなわち半径 2/λ の球の内部にある逆格子点のみが回折を起こし観測可能です。この球は制限球 (limiting sphere) と呼ばれ、逆空間に存在します。
   • 制限球の半径 (2/λ) は、回折実験の解像度 (resolution) の限界を定めます。高い解像度（小さいd値）は、逆空間でより原点から離れた（⎪h⎪が大きい）回折点を観測することに対応します。波長の短いMoKα線はCuKα線よりも制限球が大きくなり、より高解像度のデータを収集できます。最小面間隔（解像度） \(  d_{\min} \) は、 \( \theta = 90^\circ \)  のとき \( \lambda/2 \) となります。

6. 対称性と消滅則:

   • 逆格子には、結晶の対称性に対応した強度対称性 (intensity symmetry) が現れます。フリーデルの法則 (Friedel's law) は、異常分散など特別な場合を除き、I(h) = I(-h) であることを示しており、X線回折が常に中心対称であることを意味します。これにより、観測すべき独立な反射は制限球の半分に減ります。
   • 独立な反射の集合（非対称単位 (asymmetric unit)）は、逆空間内の特定の領域に存在し、結晶のラウエ群（Laue group）によって定義されます。
   • 結晶の並進対称要素（すべり面やらせん軸）は、逆空間の特定の hkl において、回折強度が系統的にゼロになる現象を引き起こします。これは系統的消滅 (systematic extinctions) と呼ばれ、フィルム上では特定の回折斑点が消えたり、ディフラクトメーターデータでは非常に弱い強度として観測されたりします。系統的消滅は、すべり面と鏡面、回転軸とらせん軸を区別するのに役立ち、空間群 (space group) の決定に不可欠な情報を提供します。KAMTRAの例では、00l反射のl=2n+1が消滅していることから、c方向の2₁らせん軸の存在が示され、空間群の決定に利用されました。

7. 単位胞の内容と密度:

   • 逆格子の単位胞の体積 V* は、直接格子の単位胞の体積 V の逆数 (V* = 1/V) に等しいです。
   • 格子定数から計算される単位胞の体積 V は、単位胞あたりの分子数 Z やX線密度 ρx の計算に用いられます。

8. 原子変位パラメーター:

   • 原子の熱運動や位置の乱れを示す変位パラメーター（以前は温度因子と呼ばれていました）は、逆格子ベクトル h の成分 (h, k, l) と逆格子定数 (a*, b*, c*, α*, β*, γ*) を含む式で表現されます。これは、回折強度に対するこれらの効果を補正するために用いられます。

要するに、逆格子は、X線回折実験で直接観測される回折斑点と、結晶の原子配置という直接空間の情報との橋渡しをする概念です。回折パターンを逆空間で解釈することで、結晶の基本的な情報（対称性、単位胞のサイズと形状）が得られ、最終的に原子の電子密度分布を導出するための構造因子が得られます（ただし位相問題の解決が必要）。

X線結晶構造解析～その２～：X線回折について

X線回折とは、X線が物質と相互作用する際に起こる回折現象を利用した実験手法です。特に、X線は物質中の電子と相互作用するため、X線回折実験から得られる結果は、その物質の電子分布（電子密度）に関する情報に起因します。したがって、X線回折実験の主な目的は、サンプルの電子密度 ρ(r) に関する情報を得ることです。

結晶のような周期的な物質における回折理論では、前述のように電子密度 ρ(r)のフーリエ変換として構造因子 F(b) が定義されます。構造因子 F(b)は逆格子空間（フーリエ空間）のベクトル b の関数であり、一般的には複素関数です。つまり、F(b)は大きさ（振幅）|F(b)| と位相 φ(b) の両方を持っています。

X線回折実験において測定されるのは、回折されたX線の強度 I(b) です。回折理論の重要な結果として、この測定される回折強度 I(b)は、構造因子 F(b)の大きさの二乗に比例します。I(b) ≈ |F(b)|²。

繰り返しになりますが、これがX線回折における位相問題 (Phase problem) の根源です。実験によって測定されるのは強度 I(b)のみであるため、これから直接得られるのは構造因子の大きさ |F(b)| だけです。しかし、目的とする物質の電子密度 ρ(r)を計算するには、構造因子 F(b)を逆フーリエ変換する必要がありますが、そのためには大きさ |F(b)|だけでなく、位相 φ(b)の情報も不可欠です。実験では位相φ(b)を直接得ることができないため、回折強度 I(b)から電子密度 ρ(r)への直接的な道筋が存在しないのです。この位相問題は、何十年もの間、X線結晶学の進歩を妨げてきた中心的な課題です。

X線回折実験と構造解析のプロセスは多岐にわたります。

• サンプル: 主に単結晶が用いられます。
• X線源: 従来のX線管（封入管、回転対陰極、マイクロフォーカス） や、より強力で波長選択が可能なシンクロトロン放射光 が使用されます。実験目的やサンプルに応じて、CuKα線（比較的長波長）やMoKα線（短波長）が一般的です。吸収が大きいサンプルや高分解能データ取得には、MoKα線やシンクロトロン放射光が有利です。単色X線が必要な場合は、フィルターやモノクロメーターが使用されます。
• 検出器: 歴史的にはX線フィルム（回転法、ワイスナーグ法、プリセッション法など）が使われましたが、近現代では回折計を用いた自動測定が主流であり、特に面積検出器（イメージングプレート、CCD）が広く使われています。現在の最先端は ハイブリッドピクセル検出器（CPAD/HPC） です。CCDに比べてノイズが極めて少なく、読み出し速度が圧倒的に速いため、現代のシンクロトロン施設や最新の回折計の主流となっています。
• データ収集: 結晶を回折計にマウントし、格子定数や方位行列を決定した後、回折強度が測定されます。測定は自動化されており、スキャン方法や2θ範囲などが設定されます。データセットの質は、「分解能（d値）」 や「測定された反射の数」、「観測されない反射」の割合によって評価されます。データの質を高めるため、低温測定が行われることもあります。
• 構造解析: 測定された回折強度データから構造を決定・精密化します。
  • 空間群決定: 測定された強度の対称性（強度対称性）と、特定の反射が観測されない系統的消滅のパターンから、結晶の空間群を推定します。強度対称性はラウエ群に対応し、通常はフリーデルの法則 I(h) = I(-h) により中心対称的な情報が得られますが、特定の波長で異常分散が存在する場合、フリーデルの法則が破れることを利用して絶対配置を決定できます。系統的消滅は空間群の並進対称要素（らせん軸や映進面）を示します。
  • 構造決定（位相決定）: 位相問題を解決するための最も重要なステップです。初期構造モデルを得るために、パターソン法（回折強度のフーリエ変換であるパターソン関数 P(u)を利用。これは電子密度の畳み込み二乗に等しく、位相問題への最初のアプローチでした。重原子法などの基盤となります）や直接法 といった手法が用いられます。異常分散を利用した多波長回折実験は、マクロ分子などの位相決定に特に有効です。SHELXS, SIR, SHELXD, DIRDIF, PATSEEなどのソフトウェアが構造決定に用いられます。
  • 構造精密化: 初期モデルを改善し、測定強度に最もよく合うように調整します。一般的に最小二乗法が用いられ、原子位置、変位パラメータ（かつて温度因子と呼ばれていました）などが最適化されます。変位パラメータは原子の熱運動や位置の乱れを反映し、ORTEPなどの図で楕円体として表現されます。SHELXL, OLEX2, CRYSTALS/CAMERONといったプログラムが使用されます。
  • 結果の解析と表示: 得られた構造の幾何学的情報（距離や角度）の計算、ORTEPなどを用いた分子構造のグラフィカルな表示、そして結果を国際的な標準形式であるCIF (Crystallographic Information File) として出力し検証する作業が行われます。

X線回折実験は、逆格子空間の概念を用いて記述されることが多く、ラウエ条件やブラッグの法則、エワルド球 といった幾何学的な条件が回折現象を理解する上で重要です。観測可能な反射の範囲は、入射X線の波長によって決まる制限球 (limiting sphere) 内に限定されます。波長  \( \lambda \) が短くなるほどエワルド球の半径 ( \( 1/\lambda \) ) が大きくなり、より高分解能（逆格子の遠くまで）のデータが収集できるという点も、「MoKα線が有利」という理由を裏付けます。また、結晶の格子体積 V、化学式量 Mr、単位格子あたりの分子数 Z、X線密度 ρx の間には関係があり、空間群決定や組成の検証に利用されます。

初期のXRAYシステム から、現代のSHELX suite やWINGX といったソフトウェアシステムに至るまで、多くの計算機プログラムがX線回折データの処理と構造解析に不可欠な役割を果たしています。

X線結晶構造解析～その１～：歴史的背景と技術の進化

X線結晶構造解析は、X線と物質中の電子との相互作用（回折現象）を利用して、物質の原子配置（電子密度分布）を決定する手法です。X線結晶構造解析は、100年以上の歴史の中で、光源、検出器、そして計算機の劇的な進化とともに、様々な理論、手法、機器、ソフトウェアの発展の歴史があります。まずは、それらを概説します。

歴史的背景と技術の進化

• 初期の実験と基本理論: X線回折実験の目的は、試料の電子密度 ρ(r) の情報を得ることにあります。X線回折は、X線が物質の電子と相互作用することによって発生し、これは量子化学におけるいわゆる「物理的に観測可能な量」である電子密度と関連付けられます。回折理論の基本結果として、回折強度 I は構造因子 F(b) の絶対値の2乗に比例することが示されています。ここで F(b) は電子密度 ρ(r) のフーリエ変換です。しかし、実験から得られるのは F(b) の絶対値 ⎪F(b)⎪ のみであり、その位相 ϕ(b) は得られません。この位相問題は、すべての回折実験における中心的な問題であり、数十年にわたりX線結晶構造解析の進歩を妨げてきました。

• 構造因子と電子密度の関係:

  電子密度 \( \rho(\mathbf{r}) \) と構造因子 \( F(\mathbf{h}) \) はフーリエ変換のペアであり、定式化すると以下の関係にあります。

  $$F(\mathbf{h}) = \int_{V} \rho(\mathbf{r}) \exp(2\pi i \mathbf{h} \cdot \mathbf{r}) dV$$

  
  すなわち、実験で得られる強度は \( I(\mathbf{h}) \propto |F(\mathbf{h})|^2 \) であるため、複素数である \( F(\mathbf{h}) \) の位相成分が失われるという「位相問題」が不可避であることが分かります。なお、F(b) とF(h)の使い分けは、主に「連続的な空間」か「離散的な格子点」かという数学的・結晶学的な文脈の違いにあります。
  F(b) ：逆空間（フーリエ空間）全体の任意の点を示すベクトル。
  F(h) ：逆格子上の「格子点」のみを示すベクトル（ミラー指数 \( h, k, l \) に対応）。

• 位相問題への初期の取り組み: この位相問題へのアプローチとして、パターソン関数 P(u) が重要です。これはパターソンによって1934/35年に初めて導入された関数であり、電子密度のコンボリューション平方に相当します。重原子法やパターソン法と呼ばれる手法は、この関数に基づいて位相問題を解決するための重要な方法です。

• フィルム法と初期の機器: 構造解析の最初の段階では、結晶の対称性や格子定数に関する情報を得るためにフィルム法が用いられました。フィルム法では、検出器としてX線感応フィルムを使用し、回折強度は点検出器が普及する1960年代中頃までフィルム技術で測定されていました。初期のフィルム法には、回転法、ラウエ法、ワイスナーグ法（ノーマルビーム法、等傾斜法）、M.J. Buergerのプリセッション法などがあります。フィルム法は現在では主に空間群の決定や結晶品質の検査に用いられ、過去の層記録能力を持つフィルム法とは異なり、数十年間にわたり点検出器（シンチレーションカウンター）が使用され、一度に一つの反射しか記録できませんでした。フィルム法から自動回折計への移行により、強度測定はより高精度かつ高速になりました。結晶の取り付けには、ゴニオメーターヘッドが使用され、結晶をアライメントするためのゴニオメーターの調整（セッティング）法やプリセッション法なども用いられましたが、これはC.W. Bunnの手法やドラッグスドルフ、M.J. Buergerによる改良に基づいています。

• 検出器の進化: 1960年代の点検出器（シンチレーションカウンター）から現代のシンクロトロン放射光の間には、「エリア検出器（イメージングプレート、CCD、現在のCMOS/ハイブリッドピクセル検出器）」の登場という大きな飛躍があります。これにより、一度に一つの反射ではなく、数千の反射を同時に記録できるようになったことで、解析スピードは劇的に向上しました。

• ソフトウェアの進化: 初期には、直接法による結晶学的な位相問題の解決のための完全なプログラムスイートとして、1971年に初めて公開された有名なプログラム MULTAN があります。これは歴史的に言及されるべきプログラムの一つです。その後の「直接法」の覇者は、G. Sheldrickによる SHELXS/SHELXD です。結晶計算の初期のアイデアを統合したプログラムシステムとして、L. FarrugiaによるWINGX が挙げられます。歴史的な主流としては、WinGX以外にも CCP4 や Phenix といった大規模な統合プラットフォームが構造生物学分野を支えてきました。これは、構造決定の開始から出版段階までの計算やグラフィックスをカバーし、重要な既存の結晶解析プログラムを統合しています。最近では、Android搭載のモバイルデバイス上で精製や表示を可能にするSHELXLのバージョンSHELXle@DROID が開発されており、これは「オフラボ結晶解析計算」の時代に入りつつあることを示しています。

• 原子モデルと温度因子: 構造解析では、結晶中の原子位置を決定しますが、原子散乱因子 fj(s) は静止した原子モデルに基づいており、ハートリーやフォックなど様々な著者による計算が国際結晶学表Vol. C にまとめられています。結晶中の原子は常に振動しているため、構造因子には補正が必要であり、これは温度因子または変位パラメーターとして考慮されます。これは1914年にデバイによって導入されたデバイ-ワラー因子 B として知られ、原子の等方的な振動を仮定した等方性温度因子と、より一般的な原子の異方的な振動を記述する異方性温度因子があります。振動または平衡位置からの小さな変位の混合効果を表すため、「熱パラメーター」という呼称が「変位パラメーター」に置き換えられることが決定されました。原子の振動挙動の図示には、1960年代にC.K. Johnsonによって導入され、現在でも広く使われているORTEP表示 が用いられます。低温実験を行うと、原子の熱運動が減少し、高角反射の強度が強まるため、データセットの分解能が高くなるという利点があり、これは近年の一般的な手法となっています。

• 対称性と空間群決定: X線回折パターンには強度対称性が現れます。これはX線回折が常に中心対称であるというフリーデルの法則 (I(h) = I(-h)) として知られています。結晶の対称要素に並進要素（並進格子、映進面、らせん軸）が加わることで空間群が生まれます（全230種類）。特に有機構造で非常に頻繁に現れる空間群としてP21/cが挙げられます。並進要素は特定の反射系列を消失させる系統的消滅を引き起こし、これが空間群の決定に大いに役立ち、映進面と鏡映面、あるいはらせん軸と回転軸の区別を可能にします。強度対称性が同じである結晶クラスはラウエ群に分類されますが、これは強度データだけでは空間群を uniquely に決定できないという問題も生じさせます。空間群決定の手順は確立されており、強度対称性からラウエ群を決定し、系統的消滅を調べ、X線密度（単位胞あたりの分子数 Z の計算）などを考慮して、国際結晶学表を参照しながら行われます。このプロセスは現代の回折計ソフトウェアに統合されていますが、ユーザー自身の結晶学的知識が重要であると強調されています。

• 特定の構造解析例: ビタミンB12は結晶学の歴史において重要な役割を果たしました。150以上の原子を持つこの構造は、1950年代中頃にイギリスのドロシー・ホジキンらのグループによって重原子法を用いて解かれました。氏の業績は、単に「大きな分子を解いた」だけでなく、コンピューターが未発達な時代にパターソン写像を駆使して解いたという執念の歴史でもあります。これは当時としては巨大な問題であり、ホジキンは他の大型生体分子に関する研究と共に1964年にノーベル化学賞を受賞しました。NaCl構造も、初期のF心格子（面心立方格子）の顕著な例として挙げられ、史上初めて報告された結晶構造でした。

• 異常分散と絶対配置: 特定の条件下ではフリーデルの法則は成り立たず、これは異常分散と呼ばれる現象によって引き起こされます。異常分散は重い元素でより顕著ですが、この効果を利用して分子の絶対配置を決定することができます。Johannes Bijvoetとその共同研究者ら は、酒石酸ナトリウムルビジウムに関する研究で、絶対配置決定への異常分散の最初の成功例を示しました。それまで、右手系か左手系かは「便宜上の仮定（フィッシャー投影式）」に過ぎませんでしたが、彼が異常分散を用いて「フィッシャーの仮定が偶然にも正しかった」ことを証明したエピソードは、化学史上極めて重要です。この研究により、酒石酸の絶対配置が実験的に初めて導かれました。絶対配置の決定には、通常よりも精密な強度測定が必要であり、フリーデル対（hと-h）の反射を測定する必要があります。

• 現代光源の登場: X線発生において熱は深刻な問題であり、冷却システムが必要です。これまでの密封管の進歩は遅かったですが、回転対陰極やマイクロフォーカスX線源が導入されました。さらに、リキッドメタルジェットX線源 のような新しい技術は、溶融金属合金の連続的な流れを標的とすることで、この物理的限界を克服し、大幅に高い輝度を達成しています。X線源の輝度は長年にわたり発展しており、X線管から異なる世代のシンクロトロンへと大きく向上しています。シンクロトロン放射光は非常に輝度が高く、シャープで高モノクロなビームを提供するため、大量の反射を短時間で高分解能かつ高完全度で測定することが可能になりました。これは特に大きな分子構造や結晶品質の低い試料の解析に有利であり、近年ではタンパク質やウイルスのような大型構造のデータ収集の大多数がシンクロトロンで行われています。ただし、シンクロトロンは大型施設であり利用に制限があるという歴史的な文脈も存在します。

これらの要素は、X線結晶構造解析がどのように発展し、現在の状態に至ったかという歴史的な流れを示しています。フィルム法から点検出器、そしてエリア検出器とシンクロトロンへの機器の進化、位相問題へのアプローチとしてのパターソン法や直接法、結晶対称性とその系統的消滅の理解による空間群決定法の確立、原子散乱因子や温度因子による原子モデルの改良、そして特定の構造解析（ビタミンB12、酒石酸など）における画期的な進歩 など、様々な側面が互いに影響し合いながらこの分野の歴史を形作ってきました。

余談

結晶学において回折斑点を議論する場合、 \( \mathbf{h} \) が使われます。

• 意味: 逆格子ベクトルを表します。具体的には \( \mathbf{h} = h\mathbf{a^*} + k\mathbf{b^*} + l\mathbf{c^*} \) であり、 \( h, k, l \) は整数（ミラー指数）です。

• 対象: 無限に続く周期的な結晶格子が前提です。

• 性質: 結晶の周期性により、回折強度は逆格子点上でのみ「鋭いピーク」として現れます。つまり、連続的な \( F(\mathbf{b}) \) を逆格子点 \( \mathbf{h} \) でサンプリング（間引き）したものが \( F(\mathbf{h}) \) です。

X線回折が起こる条件（ラウエ条件）は、「逆空間のベクトル \( \mathbf{b} \) が、ちょうど逆格子ベクトル \( \mathbf{h} \) と一致したとき（ \( \mathbf{b} = \mathbf{h} \) ）」です。

• 理論の説明時 ( \( F(\mathbf{b}) \) ): 「回折強度は構造因子の2乗に比例する」といった、物理現象の原理を語るときは、空間全体を指す  \( \mathbf{b} \) を使うのが一般的です。

• 実務・解析時 ( \( F(\mathbf{h}) \) ): 「どの斑点を測るか」「ミラー指数  \( 002 \) の強度はいくつか」という具体的なデータ（反射）を扱うときは、格子点を示す  \( \mathbf{h} \) を使います。

【遷移状態の探し方】～その１～

皆さんは、分子が踊っているところを見たことがありますか？

有機反応を扱う研究者として、
X線結晶構造解析で分子構造を「静止画」として眺めるだけでは、
どうしても物足りなさを感じてしまうときがあります。

反応とは、分子が動き、形を変え、
エネルギー障壁を越えて別の構造へと変換される“過程”そのものです。
その過程の中で、反応の行方を決定づけているのが、
一瞬だけ存在する 遷移状態（Transition State, TS） です。

このTSは寿命が極めて短く、
実験的に直接観測することはほぼ不可能です。

それにもかかわらず、
反応速度、収率、立体選択性、さらには副反応の有無に至るまで、
そのすべてがTS構造によって支配されていると言っても過言ではありません。

では、私たち研究者はどのようにして、
この「見えない構造」を覗き込むことができるのでしょうか。

その答えの一つが、計算化学による遷移状態探索です。
量子化学計算を用いれば、反応座標上を分子がどのように移動し、
どのような構造を経由して生成物へ至るのかを、
原子レベルで追跡することが可能になります。

とはいえ、実務的には、TS探索は決して簡単ではありません。
初期構造の与え方を少し誤るだけで、
計算は容易に別の極小構造へと落ち込み、
本来求めたいTSには一向に収束しません。
「TS初期構造をどう作るか」という問題は、
計算化学を扱う多くの研究者が一度は直面する壁でしょう。

今回は、その壁を越えるための実践的手法として、
xTBに実装されたメタダイナミクス（MTD）法を用い、
遷移状態に近い構造を高速にスクリーニングする方法を紹介します。

本手法の狙いは、
厳密なTSをいきなり最適化することではありません。
まず「TSっぽい構造」を大量に、短時間で集め、
それを出発点としてTS探索を行う――
このアプローチにより、試行錯誤と計算時間を大幅に削減しつつ、
TS探索の成功率を飛躍的に高めることができます。

有機合成研究者の方が、
自らの反応機構を自分の手で可視化するための第一歩として、
本稿が役立てば幸いです。

それでは、長い前振りもここまでで、初回は「単分子」における
xTBを用いたMTDシミュレーションによる遷移状態候補構造の高速探索法について、
具体的な操作手順とパラメータ設計の考え方を交えながら解説していきます。

1. 計算戦略の全体像

本手法の基本方針は以下である：

高温・強バイアス・短時間のMTDで構造遷移を強制的に起こし、
遷移状態付近の構造を大量に回収 → それをTS最適化の初期構造に使う

このアプローチにより：

• 初期構造の試行錯誤を激減

• TS探索成功率の劇的向上

• 計算時間の大幅短縮

が実現できる。

2. 初期構造の準備

まず、初期構造を xyz形式 で準備する。

Avogadro上で作成できるのであれば、可視化もAvogadroで行うので、そちらがベターです。

3. 計算環境へのログインと作業ディレクトリの作成

xTBのインストールに関しては、こちらなどを参照してください。

作業ディレクトリの作成

mkdir mtd_test
cd mtd_test

ここに：

• 初期構造：xx.xyz

• 条件ファイル：metadyn.inp

を配置する。

もし計算サーバを利用する場合はファイル転送は：

• sftp

• FileZilla

等を利用しましょう。

---

4. metadyn.inp の基本構成と意味

公式マニュアル：
https://xtb-docs.readthedocs.io/en/latest/mtd.html

基本構造は以下：

$metadyn
   save=10
   kpush=0.1
   alp=0.2
$end

この3つのパラメータが MTD挙動のほぼ全てを決定する。

---

5. MTDパラメータの物理的意味と設計指針

5.1 save：ヒル追加頻度

save = 10

→ 10ステップごと（dt=1 fsなら10 fsごと）にヒルを追加

save	意味
5	超高速・荒い
10	高速探索用（推奨）
50	精密探索
100以上	低速・高精度

TSスクリーニング目的では：
save = 10 が最適。

5.2 kpush：バイアス強度（最重要）

kpush = 0.1

→ 構造をどれだけ強く押し出すか

kpush	挙動
0.001	非常に弱い
0.01	穏やか
0.05–0.1	高速スクリーニング
0.2以上	地形破壊

TSスクリーニング目的では：

kpush = 0.05 – 0.1

が最も効率が良いが、分子サイズなどにもよりけり。

最初は少し走らせて止めて、適宜チューニングする。

5.3 alp：ヒルの幅（探索解像度）

alp = 0.2

alp	特性
0.1	超精密
0.2–0.3	実用的
0.4以上	雑・荒い

高速探索では：

alp = 0.2 – 0.3

が最適。

---

6. TSスクリーニング向け・推奨初期MTD設定

遷移状態っぽい構造を高速収集するための実用設定：

metadyne.inp

$metadyn
   save = 10
   kpush = 0.05
   alp   = 0.3
$end

追加計算条件の目安

項目	推奨値
温度	330–350 K
タイムステップ	1 fs
計算時間	20–50 ps

パラメータ設定を追加する場合は以下のようになる。

metadyne.inp

$md

   time=50
   step=1
   temp=333
$end
$metadyn 
   save=10 
   kpush=0.005 
   alp=0.3 
$end 

7. MTDシミュレーションの実行

xtb --md --input metadyn.inp xx.xyz

正常終了すると：

xtb.trj

が生成される。

8. エネルギー解析と構造の可視化

8.1 エネルギー変化の抽出

grep energy xtb.trj > mtd.txt

→ エネルギーのみを抽出
→ Excel等でエネルギープロファイル作成可

これにより：

• 反転タイミング

• エネルギー極大点

を特定できる。

8.2 トラジェクトリの可視化

mv xtb.trj xtb.xyz

として：

• Avogadro

で読み込む。

• File → Open → xtb.xyz

• Extensions → Animation → 再生

fpsを下げるとゆっくり再生でき、
構造遷移の瞬間が直感的に把握できる。

9. xtb.xyz ファイル形式の理解

構造は以下の繰り返しで記述される：

原子数
エネルギー行
座標

例：

186
energy: -265.141180200059 ...
B   0.000000   0.000000  -2.326118
N  -0.585299  -2.333618  -0.817113
...

これが1ステップ分の構造を表す。

10. TS候補構造の抽出戦略（最重要）

よくある間違い：

❌ 最大エネルギー構造をそのままTSと見なす

→ これは高確率で失敗する。

正しい構造抽出法：

1. 構造遷移が起きる瞬間（研究者としての直感で）を特定

2. 反応座標（距離・角度・二面角）が：

   • 極値を通過

   • 速度が反転
     する瞬間の構造を複数抽出

3. 5構造程度をTS候補として保存

11. 抽出構造からTS最適化への接続

TS候補構造を保存後：

• 原子間距離

• 二面角

• 結合角

の 1〜3個を固定 してGaussianなどで厳密な構造最適化を行い、
そこから通常のTS探索手順に移行する。

この方法により：

• TS収束率：劇的に向上

• 試行回数：1/5〜1/10

となる。

12. パラメータ調整指針（実務者向け）

現象	対処
動かない	kpush ↑
変形しすぎ	kpush ↓
変化が粗すぎ	alp ↓
遷移が遅い	save ↓ + kpush ↑

まとめ

xTBのMTDは、

「単分子のTS探索を人力から半自動化するための最強ツール」

であると個人的には思っています。

特に：

• 構造変化が大きい反応（ヘリセンの反転障壁など）

• 反転・環反転・異性化

• 配座変換

において、劇的な効率改善をもたらすでしょう。

以上、【遷移状態の探し方】として、xTBを用いたメタダイナミクスによるTS候補構造の高速探索法を紹介しました。

本手法は、厳密なTS最適化そのものではなく、
TSに近い初期構造を効率よく準備することに主眼を置いています。
初期構造の質がTS探索の成否を大きく左右することを考えると、
本手法は、計算化学に不慣れな有機合成研究者にとっても、
非常に強力な武器になるはずです。

その２では、2分子間反応におけるTS探索の高速化法について紹介する予定です。
近年では、北海道大学の前田先生らによって開発された AFIR法 をはじめ、
反応経路探索を自動化・高速化する優れた手法が次々と提案されています。
これらの手法をどのように実務に落とし込むか、
具体的な計算戦略とともに解説できればなと。

それでは、また。