2026年11月1日日曜日

X線結晶構造解析~その5~:構造因子について

構造因子(Structure factor)F(h)とは、結晶の電子密度(Electron density)ρ(r)のフーリエ変換(Fourier transform)として定義される関数です。これは逆格子空間(Reciprocal space)の関数であり、逆格子空間のベクトルbは逆の長さ(例えばÅ⁻¹)の次元を持ちます。

回折実験、特にX線回折実験は、物質の電子分布を表す電子密度ρ(r)に関する情報を得ることを目的としています。X線は物質の電子と相互作用して回折を起こします。

構造因子は回折理論における主要な関数の1つであり、実験で観測される回折強度I(h)は構造因子の大きさの二乗に比例します。 I(h) ~ |F(h)|²

この関係は回折理論の主要な結果ですが、同時にその主要な問題点でもあります。実験によって得られるのは回折強度Iから計算される構造因子の大きさ|F(h)|のみであり、構造因子という複素関数の位相(phase)φ(h)を得ることはできません。この位相問題(Phase problem)は、あらゆる回折実験における中心的な問題であり、何十年もの間、X線結晶学の進歩を妨げてきました。現代では重原子誘導体化法や異常分散法(SAD/MAD)などの実験的手法や、直接法・分子置換法などの計算手法によって解決策が確立されています。

構造因子は、結晶学における他の重要な関数である電子密度ρ(r)(実空間の関数)およびパターソン関数(Patterson function)P(u)(実空間の関数)と関連しています。パターソン関数P(u)は、回折強度I(h)のフーリエ変換として定義され、電子密度の畳み込み二乗(convolution square)です。情報間の関係を図示すると、電子密度ρ(r)から構造因子F(h)への変換(およびその逆)はフーリエ変換であり、回折強度I(h)(≈|F(h)|²)からパターソン関数P(u)への変換もフーリエ変換です。しかし、回折強度から電子密度への直接的な道筋は位相問題によって阻まれています。

数値計算においては、構造因子F(h)は通常、以下の原子座標に基づいた級数表現で表されます。

$$F(\mathbf{h}) = \sum_{j} f_j(\mathbf{h}) \exp(2\pi i \mathbf{h} \cdot \mathbf{r}_j)$$

ここで、hは逆格子ベクトル(格子面の法線ベクトルに相当)、\(r_j\)はj番目の原子の位置ベクトル、\(f_j(h)\)はその原子の原子散乱因子(Atomic scattering factor)です。原子散乱因子は、通常、hの大きさ、つまりsinθ/λの関数として扱われ、sinθ/λが大きくなるとゼロに近づきます。これにより、高角度の回折強度が弱くなります。

また、原子の熱運動(振動)や小さな変位(変位パラメータ displacement parameters)は、原子散乱因子をさらに減衰させます。この減衰はデバイ・ワラー因子(Debye–Waller factor)B(等方性の場合)や異方性変位パラメータ(anisotropic displacement parameters)Uijまたはβij(異方性の場合)によって記述され、例えば \( f_j(s) * exp(-B * (sinθ/λ)²) \) のように散乱因子に乗じられます。異方性変位パラメータは、原子の振動や変位の方向依存性を楕円体で表現します。

結晶の対称性も構造因子に影響を与えます。

  • 点対称中心を持つ構造(Centrosymmetric structures)では、構造因子は実数になります。この場合、位相問題は「符号問題」(プラスかマイナスか)に還元されます。また、異常分散がない通常の条件では、結晶に点対称中心があるか否かにかかわらず、フリーデルの法則(Friedel’s law)が成り立ち、回折強度I(h)とI(-h)は等しくなります(I(h) = I(-h))。
  • 異常分散(Anomalous dispersion)が存在する場合(特に重い原子や特定の波長のX線を用いる場合)、非中心対称構造においてフリーデルの法則I(h) = I(-h)は一般的に成立しなくなります。これにより生じるhと-h間の強度差(ビイフート差)を利用して、結晶の絶対配置(absolute configuration)を決定することが可能になります。
  • 結晶格子中の並進対称要素(すべり面 glide planes やらせん軸 screw axes)は、特定の指数を持つ回折強度がゼロになる系統的消滅(Systematic extinctions)を引き起こします。これは、構造因子が特定のhに対してゼロになることを意味します。系統的消滅は、同じ強度対称性を持つ異なる対称要素を区別するのに役立ちます。
  • 「全天球」測定の理由: 異常分散を利用する場合、本来等価であるはずの \( {hkl} \)  と \(  {\bar{h}\bar{k}\bar{l}} \) (フリーデル対)のわずかな強度差を正確に評価する必要があります。そのため、統計的な精度を高める意味でも、対称性で省略せずに全天球(またはそれに準ずる広範囲)を測定することが実務上の鉄則となっています。

まとめると、構造因子は結晶の電子密度分布を逆格子空間で記述するものであり、回折強度はその大きさの二乗に比例します。構造因子の位相が実験的に直接得られないという位相問題は、結晶構造解析の主要な課題です。構造因子の数学的な表現には、原子散乱因子、原子位置、熱運動(変位パラメータ)が含まれ、結晶対称性は構造因子の性質(実数であるか、特定の消滅則を持つかなど)に影響を与えます。

芳香属性再考~その2~:「二重結合2本分」はホンマに2倍なんか?

今回は、芳香族性の議論を、1,4-シクロヘキサジエンから問い直してみます。

「ベンゼンは芳香族だから安定である。」

化学を学ぶと、何度も出会う説明ですね。代表的な根拠として、よく登場するのが水素化熱の比較です。シクロヘキセンの二重結合を1本、水素化したときの発熱を測る。ベンゼンには2重結合が3本あるはずなので、水素化熱を3倍した予想と比べる。

すると、ベンゼンの水素化は予想よりも発熱が小さい。

「ほら、これが芳香族安定化です」

ホンマかいな。その「とりあえず3倍」、どこで審査通ったん?

そもそも、二重結合1本の水素化熱を、何本分でも単純に足してよいのだろうか。

ヘスの法則は、指定した反応のエンタルピーを組み合わせる関係を保証しています。

しかし、異なる分子にある二重結合が、いつも同じ水素化熱を持つことまでは保証していないのではないか?

また、反応熱の差を求めることと、その全量を特定の電子的相互作用に帰属させることも、別の作業として区別しなくてよいのか?

芳香族性に明確な定義を与えるためには、こうした参照系の振る舞いを曖昧にしたまま進むべきではないやろと考えます。そこで、今回はベンゼンを直接扱う前に、より単純な六員環の比較から始めることにしました。

まず考えたのは、次の反応です。

\[
2\,\mathrm{cyclohexene}
\longrightarrow
\mathrm{cyclohexadiene}+\mathrm{cyclohexane}
\]

「シクロヘキセン3分子」と「シクロヘキサジエン、シクロヘキセン、シクロヘキサンを各1分子」の比較から、両辺に共通するシクロヘキセンを消去した形です。

この反応は、炭素数、水素数、二重結合の総数を保存しています。もし二重結合の寄与を独立で一定のものとして扱えるなら、この比較は加成性を検査する一つの基準になると考えました。

ジエンからシクロヘキセンへの水素化熱を \(h_1\)、シクロヘキセンからシクロヘキサンへの水素化熱を \(h_2\) とすれば、上の反応のエンタルピーは、ヘスの法則により厳密に

\[
\Delta H=h_2-h_1
\]

となります。ここで注意していただきたいのは、検査するのは、その右辺がゼロになるかどうか。つまり、1本目と2本目で、水素化の発熱はホンマに同じなんかということです。ヘスの法則そのものを疑っているわけではありません。

ジエンには1,3体と1,4体がある。今回、重点を置いたのは1,4-シクロヘキサジエンです。先に1,3体を持ってくると、「そら共役してるからやろ」で話が終わりかねません。通常の直接的なπ共役鎖が、飽和炭素によって途切れている。ここでも非加成性が出るなら、「共役しているから」という一言で片づけることはできないと言えるでしょう。

もちろん、直接的なπ共役がないことは、分子内の相互作用が何もないことを意味しないです。σ骨格を介する相互作用、結合周辺の環境、環の形の影響が数字に出る可能性はある。この区別が、今回の出発点になります。

NIST Chemistry WebBookに掲載された気相の標準生成エンタルピーから、次の値を採用しました。[1–3]

分子 気相の標準生成エンタルピー/kJ mol⁻¹ 不確かさ
シクロヘキサン −123.10 ±0.79
シクロヘキセン −4.32 ±0.98
1,4-シクロヘキサジエン +104.75 ±0.59

したがって、298.15 Kの気相反応について、

\[
\Delta H=104.75-123.10-2(-4.32)
=\boxed{-9.71\ \mathrm{kJ\,mol^{-1}}}
\]

となりました。ゼロちゃいますやん。

NISTには別の文献値も掲載されてるので、採用する文献値や不確かさの確認は必要です。

とはいえ、採用したデータの中心値はゼロから外れています。、「二重結合2本やから、はい2倍」で済ませてよいのか、見過ごせない数字が出てきました。


では、この非加成性を量子化学計算にも聞いてみましょう。

まずはPySCFでB3LYP/def2-SVPによる構造最適化と振動解析を行い、def2-TZVPで基底関数を大きくした単一点計算も試しました。さらに、同じ構造と熱補正でωB97Xにも聞いてみます。

PySCFでの比較 反応エンタルピー/kJ mol⁻¹
B3LYP/def2-SVPで最適化・振動解析 −7.67
B3LYP/def2-TZVP一点計算+上記の熱補正 −6.40
ωB97X/def2-TZVP一点計算+同じ熱補正 −6.58

お前らもゼロとは言わんのかい。

基底関数を大きくすると多少動く。でも、B3LYPからωB97Xへの変更は約0.18 kJ/mol。今回のずれが、B3LYP一人の独創的な見解という感じではありません。

たしかに、二重結合の単純加算からのずれは、計算でも負の値として現れました。

とはいえ、ここで納得すると、あとで面倒くさい人が出てきます。

「しょせんB3LYPやろ」

「構造を使い回したからちゃうん」

「熱補正が悪さしてるんちゃうん」

分かる。だいたい私みたいな人間が言うてます。

なので、念のため、同じ3分子を対象にGaussianでもやり直しました。

固定構造の単一点計算だけでなく、B3LYPとωB97Xそれぞれで構造最適化と振動解析。さらに、標準G4まで投入。

DFTではグリッドと収束条件を揃え、G4にはDFT用の基底関数を無理に当てはめず、G4による所定の複合手順を用いました。G4はCCSD(T)などを組み合わせる手法であり、単に別のDFT汎関数を試すこととは異なります。ただし、構造と振動の扱いまでDFTから完全に独立しているわけではないです。[4]

結果は次のとおりでした。

Gaussianでの計算 電子成分 ZPE 0→298 K補正 反応エンタルピー
B3LYP/def2-TZVP・固定構造単一点 −6.505 +0.323 −0.241 −6.422
ωB97X/def2-TZVP・固定構造単一点 −6.664 +0.323 −0.241 −6.582
B3LYP/def2-TZVP・最適化+振動解析 −6.435 −0.152 −0.155 −6.742
ωB97X/def2-TZVP・最適化+振動解析 −6.526 −0.197 −0.104 −6.826
標準G4 −5.784 −0.173 −0.194 −6.152

単位はすべてkJ mol⁻¹です。

G4でも、ゼロには戻りません。今回のGaussian計算は、どれも−6〜−7 kJ/mol付近です。しかも、中身を分けても、振動や熱補正が大きく変化して、発熱を小さくしたわけではなさそうです。主な差は、電子成分の段階ですでに出ています。ただし、この電子成分の差を純粋な「電子相関の効果」と呼ぶことはできないですが。。。

ここまでで、少なくとも「B3LYPだけの癖」「前回の構造や熱補正を固定したせい」「G4の未完了」「集計時のZPEの取り違え」という説明は、今回のデータからは支持されないことが確認できます。

では、実測由来の−9.71 kJ mol⁻¹と、G4の−6.15 kJ mol⁻¹の差は何なのか。

差は約3.56 kJ mol⁻¹です。

採用した生成エンタルピーの掲載不確かさを独立なものとして二乗和で伝播すると、

\[
\sqrt{0.59^2+0.79^2+(2\times0.98)^2}
\simeq2.19\ \mathrm{kJ\,mol^{-1}}
\]

となります。

だからといって、「だいたい誤差や」で解散するつもりもありません。文献間の共通の基準や、不確かさの意味も確認せずに、都合よく誤差棒を振り回したらあかん。

計算の身元確認をしたら、次は実験値の身元確認も要る。忙しいな。

通常の直接的なπ共役鎖を持たない1,4-シクロヘキサジエンの比較でも、水素化熱の加成性は自明ではなく、標準G4でも約−6 kJ/molの反応差が残る。 これは、ベンゼンの議論に進む前に押さえておくべき事実です。

すなわち、ここまでで言えることは、案外シンプルです。

通常の直接的なπ共役鎖がない1,4-シクロヘキサジエンでも、水素化熱の「一本いくら」は、そのままでは通用しない。

少なくとも今回の計算では、標準G4まで使っても約−6 kJ/molの差が残りました。

では、その差は何なのか。

環の形が変わるからか。ひずみの増減か。σ骨格を介した相互作用か。それらが混ざった結果か。

ここを調べずに、ベンゼンとの反応熱の差へ進んで、その全額に「芳香族安定化エネルギー」と名前をつけてよいのか。

ベンゼンの安定性を否定する話ではありません。

まず前提がちゃんと合うてるか見ようや、という話です。

次は実験値の出どころをたどりつつ、環を大きくしたら、この「2本なのに2倍にならない問題」がどう変わるのかを調べます。

六員環が窮屈なだけなら、部屋を広くしたら落ち着くかもしれない。広くしても残るなら、別の事情がある。

現在、8・10・12・14・16員環まで計算を延長中です。

ベンゼンの話をするつもりが、まずシクロアルケンの家の間取りを調べている。

芳香族性、思ってたより玄関までが遠いみたいですわ。


参照資料:

1. NIST Chemistry WebBook:Cyclohexane(気相熱化学)。採用値:Prosen, Johnson et al.(1946)。

2. NIST Chemistry WebBook:Cyclohexene(気相熱化学)。採用値:Steele, Chirico et al.(1996)。

3. NIST Chemistry WebBook:1,4-Cyclohexadiene(気相熱化学)。採用値:Steele, Chirico et al.(1989)。

4. L. A. Curtiss, P. C. Redfern, K. Raghavachari, Gaussian-4 theory, J. Chem. Phys. 126, 084108(2007)。


To be continued

2026年10月2日金曜日

有機合成化学は、現代のリベラルアーツになりつつあるのか

今回は長々と最近なんとなーく思っていることを書き出して、ChatGPTと議論した内容をまとめてみた。

まず、かつて有機合成化学は「ものを作る技術」であった。

それがいまや「世界を考える技芸」へと昇華しつつあるということだ。


有機合成化学は、何のためにあるのか。


医薬品をつくるため。

農薬をつくるため。

色素や香料、高分子、機能性材料をつくるため。


化学に素養のある諸君には想像に難くないだろう。


近代以降の有機化学は、これほどまでに「役に立つ学問」であり続けたことも珍しいと思えるほど、産業と深く結びついてきた。


19世紀に始まった有機合成は、染料、医薬、香料、ビタミン、農薬、プラスチック、電子材料へとその領域を広げ、文字通り人間の周囲に存在する物質世界を書き換えてきた。あのK. C. Nicolaouは有機合成を、 “the art and science” of constructing organic substances と表現している。有機合成は、科学の一分野であると同時に、「art」でもある。


私は最近、この「art」という言葉を以前より重く感じる。

近年は特に、有機合成化学は、産業のための高度な道具だけではなくなってきている。

それはひょっとすると、世界の構造を考えるための一種の教養――現代のリベラルアーツ――になりつつあるのではないか。

ふと、そんな考えに至ったのだ。


リベラルアーツ(Liberal Arts)の「Arts」は、もともと「芸術」ではなかった

ここでいうArtsは、美術館に展示される絵画や、ホールで演奏される音楽だけを意味しない。

Arts の源流にあるラテン語 ars は、技芸、術、体系化された実践知といった意味を持つ。

古代から中世にかけて形成された「自由七科」とも呼ばれるリベラルアーツは、文法・修辞学・弁証法からなるTriviumと、算術・幾何学・音楽・天文学からなるQuadriviumによって構成されていた。つまり、そこには現在なら日本式に言う「理系」と分類される数学や天文学も当然のように含まれていた。

重要なのは、それらが単なる職業訓練ではなかったことである。

何か特定の製品を作るためでも、特定の職に就くためだけでもない。

世界を数えること。

形を理解すること。

言葉を論理的に扱うこと。

天体の運動を読み解くこと。

そうした知的技芸を身につけることで、人間が物事を自分自身で考えられるようになる。

つまりLiberal Artsとは、本来かなり文字通りに、人を「自由」にするためのArtsだった。

だとすれば、成熟しつつある現代の有機合成化学をそこに加えることは、それほど突飛な話ではないのではないか。


有機合成を学んだ人なら、おそらく一度は奇妙な感覚を経験したことがある。


分子をつくることは、世界について考えることである


紙の上に複雑な天然物の構造式が描かれている。

環がいくつも縮環している。

不斉中心が並んでいる。

官能基が互いに邪魔をするような位置に配置されている。


最初は、ほとんど「こんなものをどうやって作るのか」としか思えないこともあるだろう。

しかし、しばらく眺めていると、ある結合が目に入る。

ここを切ればいいのではないか。

その一本の結合を逆向きに切ると、突然巨大な構造が二つの意味のある断片に分かれる。

さらに切る。

もう一度切る。

すると、理不尽なほどに構築不能に見えた分子が、知っている分子の集合体、既知反応の連鎖へと変わっていく。


これは単なる化学反応の暗記ではない。

複雑なものを、どこから分解すれば理解可能になるのか。

という思考そのものである。


E. J. Coreyが体系化したretrosynthetic analysisの核心もそこにあるのではないだろうか。CoreyとChengの著書が『The Logic of Chemical Synthesis』と題されているのは象徴的だ。複雑な分子を単純な前駆体へと論理的に分解し、到達可能な経路を生成する。その方法論は、化学であると同時に、問題解決についての一種の形式知でもある。


すなわち、有機合成を学ぶとは、

「この反応では何ができるか」

を覚えることだけではない。


むしろ、

「複雑なものを、どの単位で見るか」

「どこを切断すれば本質が現れるか」

「どの順序なら矛盾なく組み上がるか」

という思考様式を身につけることでもあるわけだ。


数学者が世界を数や構造で見るように。

物理学者が対称性や保存則で見るように。

合成化学者は世界を、結合可能性と変換可能性によって見る。


炭素―炭素結合を見ると、「そこに存在する線」ではなく、「そこをどう形成したのか」が見える。

不斉炭素を見ると、単なる四面体ではなく、「どの面から試薬を近づければ、この立体配置だけを得られるのか」という空間的な物語が見える。

分子が名詞ではなく、動詞になるのである。


時に、合成経路には「美しい」としか呼べないものがある

そして有機合成には、どうしても「正しい」「効率的」という言葉だけでは説明できない価値がある。

それが、美しさである。


例えば、十数工程必要だと思われた骨格が、一つのcascade reactionによって一挙に立ち上がる。

離れた場所に存在していた原子が、反応の瞬間に連続して結合し、複数の環と複数の不斉中心が一挙に形成される。

ほとんど同じエネルギーを持つ二つの遷移状態のうち、一方だけを触媒が選び取り、99:1という比率で一つのエナンチオマーを作る。

最終生成物を見ただけでは思いつかないような一時的な対称性を導入し、途中まで左右同じものとして扱った後、最後の数工程でその対称性を壊す。


そこには明らかに、「うまくいった」という以上の何かがある。

迂回せず、力ずくでもなく、分子自身の性質を利用して、最小限の操作で構造が立ち上がっていく。

ある種の良い数学の証明を読んだときに感じる感覚に近い。

あるいは、人によっては優れた建築物を見たときの感覚に近い。

ほかにも解き方はあったはずなのに、これ以外ではありえなかったように見える。

それを私たちは時に「美しい」といい、時に「エレガント」と呼ぶ。


実際、全合成をArtとして語ることは単なる比喩ではない。Nicolaouは、天然物全合成を知的挑戦そのものとして行う理由や、それが若い研究者に与えるproblem-solving trainingについて明確に論じている。Organic synthesis: the art and science of replicating the molecules of living nature and creating others like them in the laboratory

2018年にはさらに踏み込み、有機合成を “an advancing Art and Science for Its Own Sake” として維持する重要性を論じている。つまり、有機合成には応用先とは独立して発展させる価値がある、という主張である。


そして、我々研究者が資本主義経済社会の中で研究活動をする上で、避けて通れない問いがある。

「何に使えるのですか」

こう問われたとき、

「この分子を作ることで社会課題を解決できるようになります」

とだけ答える必要はないというイメージはあるが、そこに毎回の正しい正解はない。


しかし、新しい反応は、新しい思考方法でもあるように、

新しい全合成は、人間が物質をどこまで理解し、どこまで制御できるかという、一つの知性の到達点でもあるのではないだろうか。

例えば、現代においては創薬のモダリティが増えたからこそ、逆に見えてきたものもある。


この問題を考える背景には、製薬産業そのものの変化もある。

二十世紀の創薬では、合成低分子化合物が圧倒的な存在感を持っていた。

ところが現在、抗体、核酸医薬、遺伝子治療、細胞治療、ADC、二重・多重特異性分子、標的タンパク質分解など、治療モダリティは急速に多様化している。

2026年の Nature Reviews Drug Discovery のレビューも、過去25年間の創薬を、遺伝学・ゲノミクス・プロテオミクスの進展とともに、従来の合成低分子や天然物誘導体の外へ治療モダリティが拡張した時代として整理している。

ただし、ここから「低分子創薬は終わった」と結論を出すのは早計である。

FDAのCDERだけでも2024年には50のnovel drugが承認されているし、現在の製薬企業も低分子を主要モダリティの一つとして維持している。例えば武田薬品のR&D戦略も、一時は十数種類のモダリティへ展開した後、2025年時点ではsmall molecules、biologics、ADCs、allogeneic cell therapiesという四つへ重点化している。


したがって、ここで起きているのは「有機合成の没落」ではない。

一時的にそのような空気感や風潮があったにせよ、むしろ、有機合成だけが医薬品を作る唯一の中心言語だった時代から、複数の言語が並存する時代への移行と考えた方がよいだろう。

そして私は、この相対化こそが重要だと思う。

産業的有用性があまりにも巨大だった時代、有機合成には常に、

「何トン作れるか」

「何工程短縮できるか」

「原薬になるか」

「特許になるか」

という実践的な問いがついて回った。

もちろん、それらは現在も重要である。

しかし、現代は有機合成の価値がそれだけではなくなったことが、かえって見えやすくなった。


分子をどのように切断するか。

どの結合を最初につくるか。

どの立体情報をいつ導入するか。

分子がもともと持つ反応性をどこまで利用するか。

保護基を使うのか、それとも使わずに化学選択性で押し切るのか。


そこには、同じ目的物へ向かう無数の可能な道がある。

だから合成には「解答」だけでなく、文学のような作風が存在する。

論文のSchemeを見れば、ときに研究者の顔を知らなくても、その思想が見えることがある。


とすると、全合成はただ単に「分子を作ること」ではなく、「選択すること」でもある。

考えてみれば、全合成で研究者が実際に行っていることの大部分は、反応そのものではない。

選択である。


どの結合を切るか。

どの官能基を残すか。

何を後回しにするか。

どの不斉を利用し、どの不斉をいつ新しく作るか。

自然物の構造のどこを「本質」と見なし、どこを「偶然」と見なすか。


つまり全合成とは、物質を作る行為であると同時に、

複雑さに対して、何を重要と判断するかという営み

でもある。


だから、一つの標的分子に対して複数の全合成経路が存在する意味がある。

最初の全合成が達成された瞬間に、その分子についての問題が終了するわけではない。


別の研究者が別の場所で切る。

別の反応を鍵反応にする。

二十工程を十二工程にする。

十二工程を七工程にする。


あるいは工程数を減らすのではなく、自然物の生合成を思わせる経路を選ぶ。

結果として得られる分子は同じである。

しかし、そこに至るまでに提示された理解の仕方は異なる。

この点で全合成は、山頂に登ることよりも、登山道を設計することに近い。

頂上に到達したという事実だけなら一度でよい。

それでも人が別のルートを登るのは、道そのものに意味があるからだ。


一方で、分子構造そのものにも造形美がある。

さらに言えば、有機化学には合成以前の美しさもある。

分子構造そのものの美しさである。


cubaneの異様な正六面体。

heliceneの螺旋。

fullereneの閉じた曲面。

taxolの密集した立体中心。

strychnineの、ほとんど悪意すら感じるほど入り組んだ骨格。


人間が「美しい形を作ろう」と設計したわけでもないのに、原子価、軌道、立体反発、電子状態といった自然法則の結果として、そこに奇妙な造形が現れる。


しかも化学者は、その形を眺めるだけではない。

その形を組み立て直す。

彫刻家が大理石から形を取り出すのだとすれば、合成化学者は形そのものを結合から作る。

建築家が重力という制約の中で空間を設計するのなら、合成化学者は原子価と熱力学と速度論という制約の中で分子を設計する。


自由ではない。

だが、制約があるからこそ創造性が生まれる。

これはまさに技芸――ars――と呼ぶにふさわしい。


そして、先の問い「役に立つか」から少しだけ離れてみる

科学には、ときどき奇妙な逆説がある。

役に立つことだけを要求され続けると、その学問が本来持っていた知的な広がりが見えなくなる。


幾何学は土地を測るためだけに存在するのではない。

天文学は暦を作るためだけに存在するのではない。

音楽理論は儀式の音楽を書くためだけに存在するのではない。


同じように、有機合成も医薬品や材料を製造するためだけに存在するわけではないのかもしれない。

有機合成には、人間が物質世界を理解するために発明してきた、ひとつの独特な思考法が蓄積されている。


結合を切って考える。

対称性を利用する。

反応性を予測する。

三次元空間を想像する。

複数の経路を比較する。

望む結果だけでなく、副反応まで先回りして想像する。


そして、複雑さを単純さへ回帰し、再び単純さから複雑さを構築する。

これはかなり高度な「思考訓練」である。


ここで、冒頭で触れた「自由七科」に戻ろう。

二十一世紀のQuadriviumに、(有機合成)化学を置くということだ。

中世のQuadriviumには、

算術

幾何学

音楽

天文学

があった。

すなわち、

数

空間

比例

宇宙

これらは世界を理解するための異なるレンズだった。


もし二十一世紀のLiberal Artsを新しく設計するとしたら、そこに私は化学を


とりわけ有機合成化学を


置いてみたい。


なぜなら有機合成は、

「世界は何でできているか」ではなく、「世界はどのように組み替えられるか」を教える学問だからである。


分析する科学は、存在するものを理解する。

合成する科学は、存在しなかったものまで思考する。


ここには決定的な違いがある。


有機合成を学ぶと、分子を見る目が変わる。

やがて、物事を見る目そのものが少し変わる。

完成された複雑なものを見ても、その姿を絶対視しなくなる。


「これは、何からできているのだろう」

「どこで分ければ理解できるだろう」

「別の順序で組み立てることはできないだろうか」


と考えるようになる。

それは分子についての問いである。

しかし同時に、かなり普遍的な問いでもある。


複雑な問題を分解し、可能な経路を考え、制約の中から最もエレガントな解を探す。


その訓練によって人間の思考が少し自由になるのであれば、

それは古い意味での Liberal Art と呼んでもよいのではないだろうか。


有機合成化学は、産業の道具であった。

これからも、おそらくそうあり続けるだろう。


だが同時に、

一本の結合をどこに引くかを何時間も考え、

一つの立体中心の向きを制御するために触媒を設計し、

数十個の原子からなる構造を、可能な無数の経路の中からただ一つの思想によって組み上げていく。


その営みを眺めていると、これはもう単なる「製造技術」ではないと思う。

有機合成とは、物質を使って行う論理学であり、

三次元空間の中で行う修辞学であり、

自然法則という制約の上に成立する建築である。


そして、そのどれもが美しい。

だから私は、有機合成化学を「役に立つ科学」とだけ呼ぶことに、少し物足りなさを感じる。

それはものを作る術であり、

と同時に、

ものの見方を作る術でもある。

その意味で、有機合成化学は今、

現代のリベラルアーツとして読み直すことのできる学問なのだと思う。


2026年6月17日水曜日

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

2026年5月18日月曜日

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