大学化学のための数学 Unit 09

フーリエ解析(NMR・FTIR・X線回折)

1. フーリエ級数

任意の周期関数は正弦波・余弦波の和として表現できます。振動分光・結晶のX線回折・NMRの解析の数学的基盤です。

フーリエ級数展開(周期 \(2L\))

\[f(x) = \frac{a_0}{2} + \sum_{n=1}^\infty \left(a_n \cos\frac{n\pi x}{L} + b_n \sin\frac{n\pi x}{L}\right)\]

フーリエ係数:

\[a_n = \frac{1}{L}\int_{-L}^{L} f(x)\cos\frac{n\pi x}{L}\,dx, \quad b_n = \frac{1}{L}\int_{-L}^{L} f(x)\sin\frac{n\pi x}{L}\,dx\]

複素フーリエ級数

\[f(x) = \sum_{n=-\infty}^{\infty} c_n e^{in\pi x/L}, \quad c_n = \frac{1}{2L}\int_{-L}^{L} f(x) e^{-in\pi x/L}\,dx\]
矩形波のフーリエ級数近似 -4 -2.67 -1.33 1.33 2.67 4 -1.50 -0.90 -0.30 0.30 0.90 1.50 x f(x) N=1 N=5 N=15
図1: 矩形波のフーリエ級数。項数を増やすほど原関数に近づくが、不連続点付近にギブス現象が残る。

パーセバルの等式(エネルギー保存)

\[\frac{1}{2L}\int_{-L}^{L} |f(x)|^2\,dx = \sum_{n=-\infty}^\infty |c_n|^2\]

2. フーリエ変換

非周期関数へのフーリエ解析の拡張です。時間領域と周波数領域を結びつけます。

フーリエ変換と逆変換

\[\tilde{f}(\nu) = \int_{-\infty}^{\infty} f(t) e^{-2\pi i \nu t}\,dt\] \[f(t) = \int_{-\infty}^{\infty} \tilde{f}(\nu) e^{2\pi i \nu t}\,d\nu\]

重要な変換ペア

時間領域 \(f(t)\)周波数領域 \(\tilde{f}(\nu)\)
\(e^{-at}\) (\(a>0, t\geq0\))\(\dfrac{1}{a + 2\pi i\nu}\)
\(e^{-at^2}\)\(\sqrt{\pi/a}\,e^{-\pi^2\nu^2/a}\)
\(\delta(t-t_0)\)\(e^{-2\pi i\nu t_0}\)
\(\cos(2\pi\nu_0 t)\)\(\frac{1}{2}[\delta(\nu-\nu_0)+\delta(\nu+\nu_0)]\)

ガウス関数の自己相似性

\[\mathcal{F}\{e^{-ax^2}\} = \sqrt{\frac{\pi}{a}}e^{-\pi^2\xi^2/a}\]

ガウス関数はフーリエ変換してもガウス関数になります。時間分解能が高い(幅が狭い)パルスほど、周波数帯域が広い(不確定性原理の時間-周波数版)。

3. NMR・分光への応用

NMR FID 信号のモデル 1.33 2.67 4 5.33 6.67 8 -1.10 -0.66 -0.22 0.22 0.66 1.10 t (ms) 信号強度 FID: e^{-t/T₂}·cos(ω₀t)
図2: NMR の FID 信号モデル \(e^{-t/T_2}\cos(\omega_0 t)\)。このフーリエ変換がスペクトルを与える。

FID からスペクトルへ

NMR の自由誘導減衰(FID)信号をフーリエ変換するとケミカルシフトのスペクトルが得られます:

\[\text{FID: } s(t) = e^{-t/T_2} e^{i\omega_0 t} \xrightarrow{\mathcal{F}} \text{スペクトル: }\tilde{s}(\omega) = \frac{T_2}{1+i(\omega-\omega_0)T_2}\]

実部がローレンツ型のスペクトル線型、線幅 \(\Delta\nu = 1/(\pi T_2)\) です。

IR・ラマン分光(FTIR)

フーリエ変換赤外分光計(FTIR)では干渉計で得られたインターフェログラムをフーリエ変換してスペクトルを得ます:

\[I(\tilde{\nu}) = \int_{-\infty}^{\infty} B(\delta) e^{2\pi i \tilde{\nu} \delta}\,d\delta\]

\(B(\delta)\):干渉信号(光路差 \(\delta\) の関数)、\(\tilde{\nu}\):波数。

X線回折(構造因子)

電子密度分布 \(\rho(\mathbf{r})\) と回折強度の関係もフーリエ変換です:

\[F(\mathbf{Q}) = \int \rho(\mathbf{r}) e^{i\mathbf{Q}\cdot\mathbf{r}}\,d^3r\]

4. 畳み込みとパワースペクトル

畳み込み(Convolution)

\[(f * g)(t) = \int_{-\infty}^{\infty} f(\tau) g(t-\tau)\,d\tau\]

畳み込み定理:

\[\mathcal{F}\{f * g\} = \tilde{f} \cdot \tilde{g}\]

スペクトルの線型(分解能関数との畳み込み)はこの関係で記述されます。

パワースペクトル(ウィーナー-ヒンチンの定理)

\[S(\nu) = |\tilde{f}(\nu)|^2 = \mathcal{F}\{C(\tau)\}\]

\(C(\tau) = \int f(t)f(t+\tau)dt\) は自己相関関数。分子動力学シミュレーションで IR スペクトルを計算する際に使います。

🎯 練習問題

まず自分で解いてから、「答えを見る」を押して確かめましょう。

問1 ★ フーリエ変換

フーリエ変換は、時間領域の関数を何領域の関数に変換するか。

答えを見る

答え周波数領域

逆変換は、周波数領域から時間領域へ戻す。

問2 ★★ デルタ関数のフーリエ変換

\(\delta(t - t_0)\) のフーリエ変換は \(e^{-2\pi i\nu t_0}\) である。この絶対値はいくつか。

答えを見る

答え1

\(|e^{i\theta}| = \sqrt{\cos^2\theta + \sin^2\theta} = 1\)。

問3 ★★ ガウス関数のフーリエ変換

ガウス関数をフーリエ変換すると、どの形の関数になるか。

  • アローレンツ関数
  • イガウス関数
  • ウデルタ関数
答えを見る

答えイ

ガウス関数は、フーリエ変換してもガウス関数になる。

問4 ★★ 畳み込み定理

畳み込み定理によると、\(\mathcal{F}\{f * g\}\) はどれに等しいか。

  • ア\(\tilde{f} + \tilde{g}\)
  • イ\(\tilde{f} * \tilde{g}\)
  • ウ\(\tilde{f}\cdot\tilde{g}\)
答えを見る

答えウ

畳み込みのフーリエ変換は、それぞれのフーリエ変換の積になる。

問5 ★★ NMR の線幅

NMR のスペクトル線の線幅は \(\Delta\nu = 1/(\pi T_2)\) で表される。\(T_2 = 0.10\) s のとき、線幅は何 Hz か。\(\pi = 3.14\) とし、小数第2位まで答えよ。

答えを見る

答え3.18 Hz

\(1/(3.14 \times 0.10) = 3.184\ldots\)

問6 ★★★ 時間幅と周波数幅

パルスの時間幅を狭くすると、そのパルスにふくまれる周波数の広がりはどうなるか。

  • ア広くなる
  • イ変わらない
  • ウ狭くなる
答えを見る

答えア

時間幅と周波数幅の積は一定に保たれる(時間と周波数の不確定性関係)。

✅ この単元のまとめ

  • フーリエ級数は周期関数を三角関数の和に分解し、フーリエ係数が各周波数成分の振幅を与える。
  • フーリエ変換は時間領域と周波数領域を相互変換し、NMR(FID→スペクトル)、FTIR、X線回折に応用される。
  • ガウス関数のFT はガウス関数で、幅の積が一定(時間-周波数の不確定性関係)。
  • 畳み込み定理 \(\mathcal{F}\{f*g\} = \tilde{f}\cdot\tilde{g}\) はスペクトルの線型解析に使われる。