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法 をはじめ、
反応経路探索を自動化・高速化する優れた手法が次々と提案されています。
これらの手法をどのように実務に落とし込むか、
具体的な計算戦略とともに解説できればなと。

それでは、また。

Catch Key Points of a Paper ~0260~

論文のタイトル: Kinetic or Thermodynamic Product? Case Studies on the Formation of Regioisomers of Tetraphenyladamantanes

著者: Tim Berking, Wolfgang Frey, Clemens Richert*

雑誌名: Synthesis

巻: Volume 55, Issue 16, pp. 2473–2482

出版年: 2023

DOI: https://doi.org/10.1055/a-2097-0092

---

背景

1: アダマンタン誘導体と結晶化シャペロン

• アダマンタンは、ダイヤモンド格子を構成する最小単位に相当する、高度に対称的で剛直な炭化水素です。
• アダマンタンの橋頭位を置換した誘導体は、ナノ構造、多孔質材料、医薬品などのビルディングブロックとして重要です。
• 特に、アルコキシフェニル基を持つテトラアリールアダマンタン（TAA）は、液体化合物を結晶内に取り込む「結晶化シャペロン」として注目されています。
• これらは、有害物質や揮発性化合物の捕捉、あるいは液体の絶対立体配置を決定するためのツールとして利用されます。

2: 合成上の課題と研究のギャップ

• 代表的なシャペロンであるTBro（ブロモ基含有TAA）は、これまでFriedel–Craftsアルキレーションの熱力学的安定生成物と考えられてきました。
• しかし、従来の合成法では収率が20%以下と低く、効率的な製造が困難でした。
• 原料のブロモアニソールには複数の反応点（ortho/para位）があるため、配向異性体（レジオアイソマー）が生成する可能性が高いことが知られていました。
• これまで、どの異性体が真に最も安定な生成物（熱力学的支配生成物）であるかは不明確なままでした。

3: 本研究の目的と期待される成果

• 本研究の目的は、強ブレンステッド酸存在下での反応を検証し、真の熱力学的安定異性体を特定することです。
• ハロゲン（フッ素、塩素、臭素、ヨウ素）の種類が、生成物の配向性に与える微妙な影響を明らかにすることを目指します。
• 得られた知見を基に、新しい結晶化シャペロンを高収率かつ高純度で提供する手法を確立します。
• 熱力学的安定異性体の構造と特性（溶解性、結晶化能）を詳細に解析し、応用可能性を検討します。

---

方法

1: 研究デザインと合成スキーム

• 本研究は、アダマンタン-1,3,5,7-テトラオール（TOA）と各種3-ハロアニソールを用いたFriedel–Craftsアルキレーションを基礎としています。
• 反応条件（温度、酸触媒の種類と量、反応時間）を変化させ、生成する異性体の分布を調査しました。
• さらに、一度生成した異性体混合物を強酸で処理し、異性化が進行するかを確認する熱力学的平衡試験を実施しました。
• この手法により、速度論的に有利な生成物（初期にできるもの）と熱力学的に有利な生成物（最終的に安定なもの）を区別しました。

2: 使用した原料、試薬、選定基準

• 原料として、TOAと3-フルオロ、3-クロロ、3-ブロモ、3-ヨードアニソールの4種を使用しました。
• 酸触媒には、トシル酸（TosOH）およびより強力なトリフルオロメタンスルホン酸（TfOH）を選択しました。
• 反応溶媒として、過剰量のハロアニソール（80当量）自体を用いています。
• 酸の強度が異性化速度や副反応（エーテル結合の切断）に与える影響を評価基準としました。

3: 主要な評価項目と分析手法

• 生成物の同定と異性体比率の決定には、¹H NMR分光法を主に使用しました。
• 各異性体の化学シフトを帰属するため、個別のピークを詳細に解析しています。
• 単離した化合物の正確な三次元構造を確認するため、X線結晶構造解析を行いました。
• 高分解能質量分析（EIまたはESI）により、目的物の組成を確認しました。

4: 最適化プロトコルの詳細

• プロトコルA（長時間）: トシル酸1当量を用い、140°Cで16時間加熱する条件です。
• プロトコルB（短時間）: トシル酸10当量を用い、140°Cでわずか2時間反応させる条件です。
• プロトコルC（異性化）: TfOH 10当量を用い、25°Cという温和な条件で65時間攪拌するプロセスです。
• これらのプロトコルを比較することで、各ハロゲンにおける最適収率条件を導き出しました。

---

結果

1: 臭素およびヨウ素置換体の逆転現象

• 臭素（Br）およびヨウ素（I）の場合、para位置換体（iTBroおよびiTIM）が真の熱力学的安定生成物であることが判明しました。
• 初期に報告されていたortho置換体（TBro）は、実際には速度論的生成物に過ぎませんでした。
• 強酸を用いた異性化プロトコルにより、iTBroの収率は従来の20%から最大91%まで飛躍的に向上しました。
• ヨウ素置換体（iTIM）についても、短時間プロトコルで66%の良好な単離収率が得られました。

2: 塩素およびフッ素における配向性の違い

• 塩素（Cl）の場合、ortho位とpara位の安定性が近く、複雑な異性体混合物が生成し、純粋な単一生成物の入手は困難でした。
• 対照的に、フッ素（F）ではortho位置換体（TFM）が圧倒的に優先され、52%の収率で得られました。
• フッ素置換体では異性化条件下でもortho位（TFM）が優勢であり、他のハロゲンとは異なる挙動を示しました。
• これはフッ素の原子サイズが小さく、電子的な要因が立体障害を上回るためと考えられます。

3: 構造特性と結晶化能の評価

• 熱力学生成物であるiTBroは、従来のTBroと比較して一般的な有機溶媒への溶解性が著しく低いことがわかりました。
• X線構造解析の結果、iTBroはジクロロメタン（CH₂Cl₂）を安定に取り込んだ溶媒和結晶を形成しました。
• iTBroの結晶構造では、遠位のハロゲン原子間の強力なパッキング力が溶解性の低下に寄与しています。
• 一方で、フッ素置換体（TFM）は新たな結晶化シャペロンとしての可能性を示唆する結晶構造を示しました。

---

考察

1: 熱力学的安定性と立体障害の関係

• 臭素やヨウ素のような大きなハロゲン原子では、アダマンタン核に近いortho位よりも、遠いpara位の方が立体的反発が少ないため安定化します。
• これに対し、原子半径の小さいフッ素では、立体障害の影響が小さいためortho位が優先される結果となりました。
• この発見は、アダマンタンのFriedel–Crafts反応において、ハロゲンの種類が配向性を決定する「微妙なバランス」を支配していることを示しています。

2: 酸触媒の役割と反応メカニズム

• トシル酸は初期のアルキレーションに有効ですが、異性化の進行は極めて緩やかでした。
• TfOH（トリフル酸）は、室温で速やかに「レトロFriedel–Crafts反応（脱アリール化と再アリール化）」を誘起し、平衡を達成させます。
• 強酸を用いることで安定なWeiland中間体が生成しやすくなり、温和な条件下での異性化が可能になったと考察されます。

3: 先行研究との比較と学術的意義

• Stetterらによる古典的な研究（1960年代）を、現代的な視点と精密なNMR分析で補完・拡張しました。
• これまで「唯一の生成物」と思われていたものが、実は平衡状態の一部に過ぎないことを明らかにした点に意義があります。
• 本研究は、アダマンタン化学における速度論と熱力学の制御に関する重要な教訓を提供しています。

4: 合成応用への道筋

• 高純度・高収率で得られるiTIMやiTBroは、パラジウム触媒を用いたクロスカップリング反応の優れたビルディングブロックとなります。
• アリール腕をさらに伸長させることで、構造的に興味深い巨大で剛直なナノ分子の合成が可能になります。
• これは材料科学分野において、新たな三次元構造設計を容易にする成果です。

5: 研究の限界点

• 今回特定された熱力学安定生成物（iTBro, iTIM）は溶解性が非常に低いため、結晶化シャペロンとしての汎用性には課題が残ります。
• 良溶媒ではない分析対象（ゲスト分子）との共結晶化は困難である可能性が高いです。
• 塩素置換体のように、安定性の差がわずかな場合は、依然として単一異性体の分離が合成上のボトルネックとなります。

---

結論

• TOAとハロアニソールのFriedel–Crafts反応における、真の熱力学的安定生成物を初めて網羅的に特定しました。
• 臭素・ヨウ素ではpara置換体、フッ素ではortho置換体が熱力学的に最安定であり、塩素は混合物を与えます。
• 強酸（TfOH）による異性化プロトコルの確立により、これまで低収率だった化合物の高効率合成が可能となりました。

将来の展望

• 実践への提言: 溶解性の高いTFM（フッ素体）は新たな結晶化シャペロンとして有望であり、今後の応用研究が期待されます。
• 将来の方向性: 低溶解性の問題を克服するため、側鎖を導入した誘導体の開発が望まれます。

---

用語集

• アダマンタン (Adamantane): 籠状の構造を持つ飽和炭化水素。非常に剛直で、化学的に安定。
• Friedel–Crafts アルキレーション: 酸触媒を用いて、芳香環にアルキル基（本研究ではアダマンタン基）を導入する有機反応。
• 配向異性体 (Regioisomer): 分子式は同じだが、置換基が付く位置（ortho位、para位など）が異なる分子。
• 熱力学支配生成物: 平衡状態において、自由エネルギーが最も低く、最も安定な生成物。
• 速度論支配生成物: 反応のエネルギー障壁が最も低く、最も速く生成する生成物。必ずしも一番安定とは限らない。
• 結晶化シャペロン (Crystallization chaperone): 結晶化しにくい液体分子などを、自らの結晶格子の中に閉じ込めて一緒に結晶化させる「介添え」分子。
• ブレンステッド酸: プロトン（H⁺）を供与する物質。トシル酸やトリフル酸など。

---

TAKE HOME QUIZ

1：研究の背景と目的

問 1: これまでの研究で、結晶化シャペロンとして知られる TBro（ブロモ基を持つ誘導体）は、反応においてどのような性質の生成物（速度論的または熱力学的）であると考えられていましたか。また、本研究ではそれがどのように修正されましたか。

問 2: テトラアリールアダマンタン（TAA）誘導体が、化学において重要視されている主な理由は何ですか。1つ挙げてください。

セクション2：ハロゲンの種類と配向性

問 3: 3-ハロアニソールを原料とした場合、ハロゲンの種類（フッ素、臭素、ヨウ素）によって、最終的に得られる主生成物の置換位置（ortho位またはpara位）はどのように異なりますか。具体的に説明してください。

問 4: 塩素（Cl）置換体において、純粋な単一異性体を得ることが困難であった理由は何ですか。

セクション3：合成手法と反応条件

問 5: 本研究で「異性化（Isomerization）」を促進するために使用された強力なブレンステッド酸は何ですか。また、従来のトシル酸（TosOH）と比較してどのような利点がありましたか。

問 6: 熱力学的安定異性体である iTBro は、従来の TBro と比較して、物理的特性（特に溶解性）においてどのような違いがありますか。また、その違いが結晶化シャペロンとしての利用にどのような影響を与えますか。

---

解答と解説

問 1: 以前は TBro が熱力学的に安定な生成物であると考えられていました。しかし、本研究により、実際には iTBro（para-位置換体）が真の熱力学的安定生成物であり、TBro は速度論的生成物であったことが明らかになりました。

問 2: 常温で液体の化合物を結晶内に取り込み、その絶対立体配置を決定するための「結晶化シャペロン」として機能するためです。

問 3: 臭素（Br）およびヨウ素（I）の場合、アダマンタン核から遠い para位置換体（iTBro、iTIM）が熱力学的に安定な主生成物となります。一方、フッ素（F）の場合は、ortho位置換体（TFM）が優先的に生成されます。

問 4: 塩素の場合、ortho位置換体とpara位置換体の安定性の差が非常に小さく、複雑な異性体混合物が形成されるためです。

問 5: トリフルオロメタンスルホン酸（TfOH / トリフル酸）が使用されました。トシル酸では異性化に140℃で1週間以上かかりましたが、TfOHを使用すると25℃という温和な条件で、65時間以内に効率よく異性化を完了させることができました。

問 6: iTBro は TBro よりも有機溶媒への溶解性が著しく低いという特徴があります。このため、それ自体が溶媒としての性質を持たないゲスト分子との共結晶化（熱的結晶化）に使用することが難しくなるという制約が生じます。