アッベの回折限界を導出する — 光学系はなぜλ/2より細かく見えないのか

顕微鏡の倍率を上げていくと、あるところで像が「大きくなるだけで細かくならない」瞬間が来ます。100倍で見えていた細胞内の粒が、1000倍にしてもぼやけた丸のまま大きくなるだけ。カメラのレンズでも同じで、高画素センサに替えても、絞りを絞りすぎると解像感が落ちます。倍率をいくら上げても、レンズが「見せられる細かさ」には物理的な上限がある — これが回折限界です。

この上限を最初に正しく定式化したのが、ドイツの物理学者エルンスト・アッベ（Ernst Abbe）でした。彼が1873年に示した式は驚くほど単純です。

$$ d_{\min} = \frac{\lambda}{2\,\mathrm{NA}} $$

波長 $\lambda$ と開口数 $\mathrm{NA}$ だけで、見分けられる最小の構造周期が決まってしまう。レンズの枚数を増やしても、収差をゼロにしても、この壁は動きません。

この限界を理解すると、驚くほど広い分野が同じ言葉でつながります。半導体リソグラフィでは、露光波長を248nm→193nm→13.5nm（EUV）と短くし、NAを液浸で1.35まで上げてきた歴史がまるごとこの式の分子と分母を殴り続けた歴史です。天文観測では、口径 $D$ の望遠鏡の角分解能 $1.22\lambda/D$ が同じ物理から出てきて、「なぜ大口径望遠鏡を作るのか」の答えになります。さらに超解像顕微鏡（2014年ノーベル化学賞）は、この式を「破った」のではなく、式が前提としている条件の外に出ることで回避しました。その回避のロジックを理解するには、まず限界そのものをきちんと導いておく必要があります。

この1枚が本記事の全体像です。左のように粗い縞なら1次回折光が浅い角度で飛ぶのでレンズの口径に入り、像に縞が再現されます。右のように縞を細かくすると回折角が大きくなり、1次光が口径から外れた瞬間に構造の情報が失われ、像は明るさだけの一様な面になります。「ぼやけて見える」のではなく「消える」という切れ方が、回折限界の性格をよく表しています。

本記事の内容

• 回折限界の直感 — なぜ「光を集める角度」が分解能を決めるのか
• 回折格子の式 $\sin\theta = \lambda/d$ からアッベ限界 $d=\lambda/(2\mathrm{NA})$ を導出（垂直照明・斜め照明の両方）
• レイリー基準 $0.61\lambda/\mathrm{NA}$ とスパロー限界 $0.47\lambda/\mathrm{NA}$ の導出と、係数が違う理由
• 光学伝達関数（OTF）のカットオフ $2\mathrm{NA}/\lambda$ による統一的な理解
• Pythonで MTF曲線・2点分解のプロファイル・構造化照明の周波数折り返しを可視化

前提知識

この記事を読む前に、以下の記事を読んでおくと理解がぐっと楽になります。特に単スリットのフラウンホーファー回折は、本記事の議論の土台そのものです。

• フラウンホーファー回折とは？単スリットの強度分布を導出する
• 波動光学入門 — 光を波として扱うとは
• 光学機器の仕組み — 顕微鏡・望遠鏡の結像

回折限界とは — 「レンズを大きくすれば無限に細かく見える」わけではない

まず、幾何光学の世界観を思い出してください。幾何光学では、物体の1点から出た光線はレンズで曲げられ、像面の1点に集まります。点は点に写る。この世界観では、収差さえゼロにすれば像は無限に鮮明になり、分解能に上限はありません。倍率を上げれば上げるだけ細かく見えるはずです。

ところが実際の顕微鏡では、そうなりません。理想的に作った対物レンズで点光源を見ても、像は点にならず、中心の明るい円盤とそれを囲む同心円状のリングになります。これがエアリーディスクです。原因は、光が波だからです。レンズの縁で光の束が「切り取られる」ことで回折が起き、どんなに完璧なレンズでも点を点に集めきれません。

左が幾何光学の世界観で、点光源はレンズを通って像面の1点に集まります。しかし実際の像は中央のように中心の明るい円盤と同心リングになり、右の断面を見ると強度がゼロになる暗環が周期的に現れています。第1暗環の位置 $v_1=3.8317$ が後で出てくるレイリー基準の係数 0.61 の正体で、この1つの数字が望遠鏡の $1.22\lambda/D$ まで支配します。

ここで大事なのは、回折の原因がレンズの品質ではなく、レンズの「有限の大きさ」そのものにあるという点です。無限に大きなレンズなら回折は起きませんが、そんなものは作れません。有限の口径で光を切り取る以上、回折は必ず起きる。だから回折限界は「技術が未熟だから届かない限界」ではなく、「波動光学が課す原理的な壁」なのです。

アッベがすごかったのは、この壁を「エアリーディスクの大きさ」ではなく、「物体のどこまで細かい周期構造が像に伝わるか」という観点で捉え直したことです。彼はイェーナのカール・ツァイス社で顕微鏡の設計に携わっていて、経験的な試行錯誤ではなく理論から設計したかった。そこで「顕微鏡は物体の何を運んでいるのか」を突き詰めた結果、次の描像にたどり着きました。

物体を照らすと、物体の細かい構造が光を回折させる。その回折光のうち、対物レンズの口径に入ったものだけが像を作る。細かい構造ほど回折角が大きく、口径から外れる。外れた瞬間、その構造の情報は像から完全に消える。

この「情報が入るか外れるか」というイエス・ノーの判定が、アッベの回折限界の本質です。エアリーディスクの重なり具合のような「程度問題」ではなく、口径に入らなければゼロ、というデジタルな切断になっている。ここが、後で出てくるレイリー基準との根本的な性格の違いになります。

では、その「回折角」を定量化していきましょう。鍵になるのは、物体を回折格子だと思うことです。

アッベの発想 — 物体を「回折格子の重ね合わせ」とみる

なぜいきなり回折格子が出てくるのか。それは、どんな物体も周期構造の足し算で書けるからです。これはフーリエ解析の考え方そのものです。物体の透過率分布 $t(x)$ を

$$ t(x) = \int_{-\infty}^{\infty} \hat{t}(f)\, e^{i 2\pi f x}\, df $$

と展開すれば、$t(x)$ は空間周波数 $f$ [本/mm] の正弦波格子を無数に重ね合わせたものとみなせます。細かい構造は高い $f$（＝短い周期 $d=1/f$）の成分に対応します。

顕微鏡の光学系が線形であれば、各成分は独立に伝わります。つまり、「周期 $d$ の格子が像に伝わるか？」という問いにだけ答えれば、任意の物体について何が見えるかが決まる。アッベの戦略はここにあります。複雑な物体を相手にせず、格子1枚に絞って考えるわけです。

さらに、周期 $d$ の格子から出る回折光は、離散的な方向にしか進みません。$0$次、$\pm1$次、$\pm2$次…と、はっきり分かれた角度に飛びます。だから「入るか外れるか」の判定が非常にクリアになります。ここが、アッベの議論が美しく決まる理由です。

そこで次に、その回折角を与える基本式 $\sin\theta = \lambda/d$ を、丁寧に導いておきましょう。

回折格子の式 sinθ = λ/d を導く

周期 $d$ の格子に、波長 $\lambda$ の平面波を垂直に当てます。格子の各周期が同じように光を散乱するので、隣り合う周期から出た光の光路差が波長の整数倍になる方向で、光は強め合います。

隣接する2つの散乱点は距離 $d$ だけ離れています。この2点から角度 $\theta$ の方向に進む光を考えると、下側の点から出る光は上側より

$$ \Delta = d \sin\theta $$

だけ余分に進みます。図で言えば、間隔 $d$ の底辺に対して、$\theta$ 方向への射影が $d\sin\theta$ になるという、直角三角形1個ぶんの幾何です。

破線は $\theta$ 方向へ進む光にとっての共通の波面で、右の散乱点から出た光はここに到達済みです。左の散乱点だけが橙色の線分ぶん余計に走らねばならず、その長さが直角三角形の斜辺 $d$ に対する射影 $d\sin\theta$ になります。この線分が波長の整数倍のとき、隣り合う周期からの光が位相を揃えて強め合います。

強め合う条件は、この光路差が波長の整数倍になることです。

$$ \begin{equation} d \sin\theta_m = m\lambda \quad (m = 0, \pm 1, \pm 2, \ldots) \end{equation} $$

$m$ を回折次数と呼びます。$m=0$ は $\theta_0 = 0$、つまり真っ直ぐ抜ける光（0次光）です。1次回折光は

$$ \sin\theta_1 = \frac{\lambda}{d} $$

の方向に出ます。この式が本記事の心臓部です。格子が細かい（$d$ が小さい）ほど、$\sin\theta_1$ が大きくなり、1次光は大きく外側へ飛ぶ。逆に粗い格子なら1次光は0次光のすぐ横に出ます。

ここで、$\sin\theta_1 \le 1$ でなければ解が存在しないことに注目してください。$d < \lambda$ になると $\lambda/d > 1$ となり、1次回折光そのものがこの世に存在できなくなります。正確には、そのとき1次光は伝搬波ではなくエバネッセント波（表面から指数関数的に減衰する波）になり、物体表面から波長程度離れると消えてしまいます。この事実は後で「なぜ $\lambda/2$ という壁が出てくるのか」の根っこになるので、頭の片隅に置いておいてください。

左の扇形を見ると、$d=3\lambda$ では1次光がわずか19度なのに、$d=1.05\lambda$ では72度まで倒れ、格子面すれすれに逃げていくのがわかります。右の曲線 $\sin\theta_1=\lambda/d$ は $d$ が小さいほど急激に立ち上がり、$d=\lambda$ で1に到達します。ここから左の赤い領域は $\sin\theta_1>1$ で解が存在せず、1次光は伝搬波になれません。NA=0.95のレンズなら $d=1.05\lambda$ が垂直照明での限界で、これは伝搬の限界 $d=\lambda$ のすぐ手前です。

なお、屈折率 $n$ の媒質中では波長が $\lambda/n$ に縮むので、格子方程式は

$$ n\, d \sin\theta_m = m\lambda $$

と書けます（$\lambda$ は真空中の波長）。油浸レンズが効く理由がここに隠れています。

さて、格子から出た回折光の角度がわかりました。次に問うべきは、「その光は対物レンズに入るのか？」です。

垂直照明のときの限界 d = λ/NA

対物レンズが物体側で受け入れられる最大の角度を $\alpha$ とします（レンズの縁を見込む半角）。物体が屈折率 $n$ の媒質中にあるとき、この受光能力を表す量が開口数（Numerical Aperture）です。

$$ \begin{equation} \mathrm{NA} = n \sin\alpha \end{equation} $$

$\mathrm{NA}$ が大きいほど、大きく傾いた光まで拾えます。空気中（$n=1$）では $\mathrm{NA}<1$ ですが、屈折率1.515の浸液を使えば $\mathrm{NA}=1.4$ 程度まで届きます。

ここでアッベの決定的な問いを立てます。像に格子の縞が現れるためには、何が必要か？

答えは「2つ以上の回折次数がレンズに入ること」です。理由は干渉にあります。レンズを通った光は像面で重なり合いますが、1本の平面波だけでは干渉縞は作れません。0次光だけが入った場合、像面には一様な明るさが広がるだけで、縞は一切現れません。0次光と1次光の2本が入って初めて、両者の干渉によって周期的な明暗（＝格子の像）が生まれます。

これは非常に強い主張です。像のコントラストが「弱くなる」のではなく、完全にゼロになる。細かすぎる格子は、ぼやけて見えるのではなく、のっぺりした一様な明るさとして見える。この「消え方」がアッベ限界の特徴です。

3枚とも同じ縦軸で描いてあります。上段の0次光だけが通った場合、強度は位置によらず一定でコントラストはちょうど 0.000、下のグレースケール帯も一様な灰色です。中段のように+1次光が加わるとコントラストは 0.800 まで立ち上がり、下段で±1次まで揃うと 1.000 に達します。「情報が届かないと像が暗くなる」のではなく「明るさは残るのに縞だけが消える」という点が重要です。

では条件を書き下しましょう。垂直照明（$\theta_i = 0$）では、0次光は $\theta=0$ なので必ず入ります。したがって条件は「1次光が入ること」です。

$$ n \sin\theta_1 \le n\sin\alpha = \mathrm{NA} $$

左辺に格子方程式 $n d\sin\theta_1 = \lambda$ から得られる $n\sin\theta_1 = \lambda/d$ を代入すると、

$$ \frac{\lambda}{d} \le \mathrm{NA} $$

$d$ について解くために両辺に $d$ を掛け、$\mathrm{NA}$ で割ると、

$$ \begin{equation} d \ge \frac{\lambda}{\mathrm{NA}} \end{equation} $$

これが垂直照明でのアッベ限界です。$\lambda=550\,\mathrm{nm}$、$\mathrm{NA}=1.4$ の油浸対物なら $d_{\min}=393\,\mathrm{nm}$ となります。

瞳を真上から見た図です。0次光は必ず中心にいて、±1次光は中心から $\lambda/d$ だけ離れた位置に現れます。格子が粗い左では余裕で瞳の中に収まり、中央でちょうど縁に接し、右では白抜きの丸が示すとおり瞳の外へ出てしまいます。1次光が使えるのは中心から縁までの半径 NA ぶんだけで、これが $d_{\min}=\lambda/\mathrm{NA}$ の由来です。

ところが、教科書によく載っているアッベの式は $\lambda/(2\mathrm{NA})$、つまりこの半分です。この2倍の差はどこから来るのでしょうか。答えは、照明の当て方にあります。

斜め照明で限界が半分になる — d = λ/(2NA)

垂直照明の議論には、実は無駄がありました。0次光が瞳（レンズの開口）のど真ん中を通っているのです。0次光は「格子を素通りした光」なので、それ自体は構造の情報を持っていません。情報を持っているのは1次光のほう。にもかかわらず、0次光が瞳の中央という一等地を占領し、1次光は瞳の端までの距離しか使えていません。

そこで、照明を斜めから当てて、0次光を瞳の片端に寄せてしまう。すると1次光は瞳の反対の端まで使えるので、許容できる回折角が2倍になります。これがアッベ限界が半分になるカラクリです。

斜め入射のときの格子方程式を導いておきましょう。入射角を $\theta_i$（法線から測る）とすると、格子に到達する時点ですでに隣接周期間に光路差 $d\sin\theta_i$ が生じています。出射側の光路差は先ほどと同じく $d\sin\theta_m$ なので、正味の光路差は差し引きになり、

$$ \begin{equation} n\,d\,(\sin\theta_m – \sin\theta_i) = m\lambda \end{equation} $$

となります。垂直照明 $\theta_i=0$ とすれば、先ほどの式に戻ることを確認してください。

さて、瞳が受け入れる方向は $|n\sin\theta| \le \mathrm{NA}$ の範囲です。この帯域を最大限使いたいので、0次光を片端ぎりぎり、1次光を反対の端ぎりぎりに配置します。0次光は入射方向そのままに進むので $n\sin\theta_0 = n\sin\theta_i$。これを瞳の片端に置くと

$$ n\sin\theta_i = -\mathrm{NA} $$

1次光を反対の端に置くと

$$ n\sin\theta_1 = +\mathrm{NA} $$

この2つを斜め入射の格子方程式（$m=1$）に代入します。

$$ d\,\bigl(n\sin\theta_1 – n\sin\theta_i\bigr) = \lambda $$

代入すると、括弧の中は $\mathrm{NA} – (-\mathrm{NA}) = 2\mathrm{NA}$ になるので、

$$ d \cdot 2\mathrm{NA} = \lambda $$

両辺を $2\mathrm{NA}$ で割れば、目的の式が得られます。

$$ \begin{equation} \boxed{\;d_{\min} = \frac{\lambda}{2\,\mathrm{NA}}\;} \end{equation} $$

これがアッベの回折限界です。$\lambda=550\,\mathrm{nm}$、$\mathrm{NA}=1.4$ なら $d_{\min}=196\,\mathrm{nm}$。垂直照明の393nmからちょうど半分になりました。

左右の違いは0次光をどこに置くかだけです。左の垂直照明では0次光が瞳の中心に居座るため、0次と1次の間隔として使えるのは半径ぶんの NA しかありません。右のように照明を傾けて0次光を左端まで寄せると、1次光を右端まで持っていけるので間隔は直径ぶんの 2NA になります。この紫の矢印の長さが格子方程式の括弧 $n\sin\theta_1-n\sin\theta_i$ そのもので、長さが2倍になれば解ける周期は半分になります。

実務的には、照明側にもコンデンサレンズがあり、その開口数を $\mathrm{NA}_{\mathrm{cond}}$ とすると、より一般的な式は

$$ d_{\min} = \frac{\lambda}{\mathrm{NA}_{\mathrm{obj}} + \mathrm{NA}_{\mathrm{cond}}} $$

と書けます。$\mathrm{NA}_{\mathrm{cond}}=0$（点光源による垂直照明・完全コヒーレント）なら $\lambda/\mathrm{NA}_{\mathrm{obj}}$、$\mathrm{NA}_{\mathrm{cond}}=\mathrm{NA}_{\mathrm{obj}}$（コンデンサを対物と同じまで開く）なら $\lambda/(2\mathrm{NA})$。顕微鏡でコンデンサの絞りを絞りすぎると解像が落ちるのは、まさにこの式のためです。「明るさを稼ぐために絞りを開ける」のではなく、分解能そのものが照明の開き角で決まっている。アッベの理論が顕微鏡の設計思想を変えた理由がここにあります。

そして、$\mathrm{NA}=n\sin\alpha \le n$ という上限があるので、原理的な最良値は

$$ d_{\min} \ge \frac{\lambda}{2n} $$

空気中なら $\lambda/2$、水中なら $\lambda/2.66$。タイトルの「$\lambda/2$ より細かく見えない」はこの意味です。$\lambda=550\,\mathrm{nm}$ の緑色光なら 275nm。これがレンズを使う限りの絶対的な壁になります。

ここまでで「周期構造がどこまで伝わるか」がわかりました。次は、この結論を空間周波数の言葉に翻訳しておきます。すると、レイリー基準との関係が見通しよくなります。

空間周波数で言い換える — カットオフ 2NA/λ

周期 $d$ の格子は、空間周波数 $f = 1/d$ の成分です。アッベ限界 $d \ge \lambda/(2\mathrm{NA})$ を $f$ について書き直すと、

$$ \begin{equation} f \le f_c = \frac{2\,\mathrm{NA}}{\lambda} \end{equation} $$

$f_c$ をカットオフ空間周波数と呼びます。光学系は $f_c$ 以下の空間周波数を通し、それ以上は完全に遮断するローパスフィルタとして振る舞う、というのがアッベ理論の最終形です。信号処理を知っている人にはおなじみの構図でしょう。実際、両者は数学的に同じもので、「レンズの瞳＝周波数領域の窓」という対応がそのまま成り立ちます。

この描像を定量化するのが光学伝達関数（OTF: Optical Transfer Function）です。インコヒーレント照明（蛍光顕微鏡など、各点が独立に光る場合）では、OTFは瞳関数の自己相関で与えられます。

$$ \mathrm{OTF}(f) = \frac{\displaystyle\int P(\bm{u})\,P^*(\bm{u}+\lambda f\,\hat{\bm{e}})\,d\bm{u}}{\displaystyle\int |P(\bm{u})|^2 d\bm{u}} $$

ここで $P$ は瞳関数（開口内で1、外で0）です。自己相関なので、2つの円をずらして重なり面積を測る操作に相当します。円の半径が $\mathrm{NA}/\lambda$ なら、ずらし量が直径 $2\mathrm{NA}/\lambda$ を超えた瞬間に重なりが消えます。これがカットオフが $2\mathrm{NA}/\lambda$ になる幾何学的な理由です。アッベが「0次光と1次光を瞳の両端に置く」と言ったのと、まったく同じことを言っています。

上段の紫色の領域が2つの瞳の重なりで、ずらし量が増えるほど痩せていきます。ずらし量0.5NAで重なり比は 0.685、1.0NAで 0.391、1.6NAで 0.104 と単調に減り、下段の曲線がそのままMTFになっています。半径 $\mathrm{NA}/\lambda$ の円を2枚ずらしているので、ずらし量が直径 $2\mathrm{NA}/\lambda$ に達したところで重なりが消える — カットオフの値が $2\mathrm{NA}/\lambda$ である理由は、この単純な作図に尽きます。

半径 $\rho_0=\mathrm{NA}/\lambda$ の2円の重なり面積を計算すると、正規化空間周波数 $u=f/f_c$ に対して

$$ \begin{equation} \mathrm{MTF}(u) = \frac{2}{\pi}\left(\arccos u – u\sqrt{1-u^2}\right), \quad 0 \le u \le 1 \end{equation} $$

が得られます。MTF（Modulation Transfer Function）はOTFの絶対値で、「入力コントラストの何割が出力に残るか」を表します。$u=0$ で 1、$u=1$ で 0 になり、途中は単調減少です。数値で言うと $u=0.25$ で 0.685、$u=0.5$ で 0.391、$u=0.75$ で 0.144。カットオフ手前でも、細かい構造のコントラストはかなり削られていることがわかります。この「削られ具合」こそが、次に扱うレイリー基準の正体です。

なお、コヒーレント照明（レーザー照明など、物体面での光が完全に位相相関している場合）では話が変わります。コヒーレント系では振幅が線形に伝わるので、伝達関数は瞳関数そのものになり、カットオフは $\mathrm{NA}/\lambda$ です。これは先ほど導いた垂直照明のアッベ限界 $d=\lambda/\mathrm{NA}$ と一致します。つまり、

• コヒーレント（垂直照明・点光源）：$f_c^{\mathrm{coh}} = \mathrm{NA}/\lambda$、$d_{\min}=\lambda/\mathrm{NA}$
• インコヒーレント（全方位照明・$\sigma=1$）：$f_c^{\mathrm{inc}} = 2\mathrm{NA}/\lambda$、$d_{\min}=\lambda/(2\mathrm{NA})$

という対応になります。「垂直照明と斜め照明で2倍違う」という話と、「コヒーレントとインコヒーレントで2倍違う」という話は、同じ物理の別の言い方なのです。よく混乱するポイントなので、瞳の中で0次光がどこにいるかを思い浮かべて整理してください。

ここまでは「周期構造」の話でした。次に、天文学の伝統である「2つの点をどこまで見分けられるか」に移ります。同じ物理から、なぜ違う係数が出てくるのかを見ていきます。

レイリー基準 — エアリーディスクの重なりで決める

夜空の二重星を望遠鏡で見るとき、2つの星が分かれて見えるか、1つの楕円に見えるかは接眼倍率ではなく口径で決まります。この経験則を定式化したのがレイリー卿でした。彼の発想はアッベとまったく違います。周期構造ではなく、2つの点光源を主役に据えたのです。

まず、1つの点光源の像を求めましょう。円形開口を通った光の遠方回折像は、フラウンホーファー回折の結果として

$$ \begin{equation} I(v) = \left[\frac{2J_1(v)}{v}\right]^2, \qquad v = \frac{2\pi\,\mathrm{NA}}{\lambda}\,r \end{equation} $$

と書けます。$J_1$ は1次のベッセル関数、$r$ は物体面に換算した像面上の距離です。この $I(v)$ がエアリーパターンです。中心の明るい円盤（エアリーディスク）の周りに、暗環と明環が交互に並びます。

第1暗環は $J_1(v)=0$ の最初の非自明な根で決まります。数値的に

$$ v_1 = 3.8317 $$

です。$v$ の定義から $r$ に戻すと、

$$ r_1 = \frac{v_1 \lambda}{2\pi\,\mathrm{NA}} = \frac{3.8317}{6.2832}\cdot\frac{\lambda}{\mathrm{NA}} = 0.6098\,\frac{\lambda}{\mathrm{NA}} $$

丸めて $0.61\lambda/\mathrm{NA}$。これがエアリーディスクの半径です。

レイリー基準は、次のように定義されます。

一方の点像の主極大が、もう一方の点像の第1暗環に重なるとき、2点はぎりぎり分解されている。

つまり2点の間隔 $s$ がエアリーディスク半径に等しいとき、

$$ \begin{equation} \boxed{\;s_{\mathrm{Rayleigh}} = 0.61\,\frac{\lambda}{\mathrm{NA}}\;} \end{equation} $$

望遠鏡の文脈では、$\mathrm{NA}\approx \sin\alpha \approx (D/2)/f$（$D$ は口径、$f$ は焦点距離）を代入して角度に直します。

$$ \theta_{\min} = \frac{s}{f} = 0.61\frac{\lambda}{f}\cdot\frac{2f}{D} = 1.22\,\frac{\lambda}{D} $$

おなじみの $1.22\lambda/D$ が出てきました。係数 1.22 は $3.8317/\pi$ に他なりません。ベッセル関数の零点という、一見天下り的な数字が、天文学の教科書に載る係数の正体です。口径8mの望遠鏡なら $\theta_{\min}=8.39\times10^{-8}\,\mathrm{rad}=0.017$ 秒角、ハッブル（2.4m）なら 0.058 秒角。大口径を追い求める理由が数字で見えます。

さて、この基準はどれくらい「見える」状態なのでしょうか。2つの点が非コヒーレントに光っている場合、像は2つのエアリーパターンの強度の単純和になります。この和をレイリー間隔で計算すると、中心の谷の強度はピークの 73.5% になります。つまり26.5%のくぼみができている。「明確に2つに見える」というより「なんとか2つと言える」レベルで、実はかなり控えめな基準です。

望遠鏡の視野で二重星を見たときの見え方を計算したものです。黄色い線は水平断面の強度で、左端では谷が深く2つの星が明確に分かれています。レイリー基準（左から2番目）では確かに2つの山が判別できますが、中心はピークの 0.735 までしか下がっておらず、少しでも大気の揺らぎやノイズがあれば怪しくなる程度の分かれ方です。右2枚では谷が完全に消え、1つのぼやけた円盤にしか見えません。

そしてもう1つ重要なのは、レイリー基準は物理法則ではなく、人が決めた便宜的なしきい値だということです。「第1暗環に重なったとき」という条件に必然性はありません。ベッセル関数の零点というきりのいい場所を選んだだけです。だから係数0.61には、アッベの0.5のような導出上の必然性がない。ここが2つの基準の性格の違いです。

では、しきい値をもっと攻めるとどうなるか。次に見るスパロー限界がその答えです。

スパロー限界 — 谷が消える瞬間

レイリー基準では「谷が26.5%残っている」状態でした。もっと近づけると谷は浅くなり、やがて完全に平らになって消えます。この「谷が消える直前」を分解の限界とするのがスパロー（Sparrow）限界です。「2つに見えると言い張れる最後の瞬間」を採る、より攻めた基準です。

数学的な条件は明快です。合成強度プロファイルを

$$ I_{\mathrm{sum}}(x) = I\!\left(x-\tfrac{s}{2}\right) + I\!\left(x+\tfrac{s}{2}\right) $$

とすると、対称性から中心 $x=0$ は常に停留点（1階微分ゼロ）です。谷があるということは中心が極小、つまり2階微分が正。谷が消えるのは、2階微分がゼロになる瞬間です。

$$ \left.\frac{d^2 I_{\mathrm{sum}}}{dx^2}\right|_{x=0} = 0 $$

エアリーパターンについてこれを数値的に解くと、正規化変数で $v_s = 2.9756$ が得られます。$s = v_s\lambda/(2\pi\mathrm{NA})$ に戻すと、

$$ \begin{equation} \boxed{\;s_{\mathrm{Sparrow}} = 0.4736\,\frac{\lambda}{\mathrm{NA}} \approx 0.47\,\frac{\lambda}{\mathrm{NA}}\;} \end{equation} $$

レイリーの0.61より2割ほど小さい値です。同じ光学系、同じ物理なのに、「分解できた」の定義を変えるだけで数字が2割動く。これが「分解能」という言葉を扱うときに一番気をつけるべき点です。カタログスペックを比べるときは、どの基準の数字なのかを必ず確認する必要があります。

3つの基準を並べてみましょう。

係数が 0.47〜0.61 の狭い範囲に収まっていることに注目してください。定義がバラバラなのに、桁も前後関係も変わらない。これは3つが同じ物理（瞳の有限性）を別々の切り口で測っているからです。実務では「$\lambda/(2\mathrm{NA})$ 程度」と覚えておけば十分で、係数の細かい違いは定義の問題にすぎません。

ここで自然な疑問が湧きます。スパロー限界 $0.474\lambda/\mathrm{NA}$ は、アッベのカットオフ周期 $0.5\lambda/\mathrm{NA}$ より小さい。カットオフを超えた細かさなのに、なぜ分解できるのでしょうか。次節でこの矛盾を解きます。

3つの基準はどう統一されるか — OTFで見る

矛盾に見えたものの正体は、「2点」と「格子」のスペクトルの違いです。

周期 $d$ の無限格子は、空間周波数 $f=1/d$ に線スペクトルを持ちます。エネルギーが1本の周波数に集中しているので、その周波数がカットオフを超えれば情報は完全にゼロ。これがアッベの「イエス・ノー」の判定です。

一方、2つの点光源はどうでしょう。$o(x) = \delta(x-s/2)+\delta(x+s/2)$ のフーリエ変換は

$$ \hat{o}(f) = e^{-i\pi f s} + e^{i\pi f s} = 2\cos(\pi f s) $$

すべての空間周波数に成分を持つ、広帯域のスペクトルです。だから、間隔 $s$ に対応する周波数 $1/s$ がカットオフを超えていても、通過帯域内に残った低周波成分だけで谷を再構成できる。像の強度分布は

$$ \hat{I}(f) = \hat{o}(f)\cdot\mathrm{OTF}(f) $$

を逆変換したもので、$|f|\le f_c$ の帯域を積分した結果として谷が現れます。谷の有無は「$1/s$ が通ったかどうか」ではなく、「通過帯域全体の重ね合わせがどんな形になるか」で決まる。だからスパロー限界がアッベのカットオフ周期をわずかに下回っても、何も矛盾しません。

逆に、レイリー間隔 $0.61\lambda/\mathrm{NA}$ に対応する周波数は $f = 1.64\,\mathrm{NA}/\lambda = 0.82 f_c$ です。この周波数でのMTFは 0.089 — わずか9%まで落ちています。レイリー基準はカットオフのかなり手前にあり、「MTFが1割弱でも人間の目には2つに見える」という経験的な線引きだったことがわかります。

まとめると、統一的な理解は次のようになります。

1. 光学系の本質は $f_c = 2\mathrm{NA}/\lambda$ のローパスフィルタである。これがアッベの結論であり、動かせない物理。
2. レイリー基準・スパロー限界は「フィルタを通した後の像を、人がどう判定するか」の規約である。物体の種類（点か格子か）と判定基準（谷が何%か）に依存するので、係数が変わる。
3. 3つとも $f_c$ を親として派生している。だから係数は 0.47〜0.61 という狭い範囲に収まり、$\lambda/\mathrm{NA}$ というスケール則は共通する。

この整理ができると、「うちの装置は $0.4\lambda/\mathrm{NA}$ を達成した」といった主張を正しく読めるようになります。$f_c$ を超えていなければ基準の取り方の問題、超えているなら何らかの追加情報を使っている、という切り分けができるからです。

では実際に、これらをPythonで描いて確かめましょう。

Pythonでの実装

3つの可視化を行います。(1) MTF曲線とカットオフ、(2) 2点光源の間隔を変えた像プロファイル、(3) 構造化照明による周波数の折り返し。まずは共通の準備からです。

import numpy as np
import matplotlib
import matplotlib.pyplot as plt
from scipy.special import j1, jn_zeros
from scipy.optimize import brentq

# 日本語フォント設定
for cand in ["Hiragino Sans", "Yu Gothic", "Noto Sans CJK JP", "IPAexGothic", "Meiryo"]:
    if any(cand == f.name for f in matplotlib.font_manager.fontManager.ttflist):
        plt.rcParams["font.family"] = cand
        break
plt.rcParams["axes.unicode_minus"] = False

# エアリーパターン I(v) = [2 J1(v)/v]^2  （v=0 は極限値1）
def airy(v):
    v = np.asarray(v, dtype=float)
    out = np.ones_like(v)
    m = np.abs(v) > 1e-12
    out[m] = (2.0 * j1(v[m]) / v[m]) ** 2
    return out

# 円形瞳のMTF（インコヒーレント）: u = f/fc
def mtf_circ(u):
    u = np.clip(np.asarray(u, dtype=float), 0.0, 1.0)
    return (2.0 / np.pi) * (np.arccos(u) - u * np.sqrt(1.0 - u ** 2))

v1 = jn_zeros(1, 1)[0]           # 第1暗環 3.8317
print(f"J1の第1零点 v1 = {v1:.4f},  v1/pi = {v1/np.pi:.4f}")
print(f"レイリー係数 = v1/(2pi) = {v1/(2*np.pi):.4f}")

出力は J1の第1零点 v1 = 3.8317, v1/pi = 1.2197、レイリー係数 = 0.6098 となります。手計算で追った $1.22\lambda/D$ の 1.22 と、$0.61\lambda/\mathrm{NA}$ の 0.61 が、どちらもベッセル関数の零点1個から出ていることが数値でも確認できました。係数を暗記する必要はなく、$v_1=3.8317$ さえ覚えておけば両方復元できます。

(1) MTF曲線とカットオフ

光学系がローパスフィルタであることを、MTF曲線で可視化します。3つの基準がこの曲線のどこに位置するかも重ねて示します。

u = np.linspace(0, 1, 400)
M = mtf_circ(u)

# 各基準に対応する正規化空間周波数 u = f/fc  （f = 1/s, fc = 2NA/λ）
u_abbe = 1.0                       # 0.5 λ/NA
u_ray = 1.0 / (0.6098 * 2)         # 0.61 λ/NA
u_spa = 1.0 / (0.4736 * 2)         # 0.47 λ/NA（>1: 帯域外）

fig, ax = plt.subplots(figsize=(9, 5.5))
ax.plot(u, M, lw=2.5, color="#1f77b4", label="MTF（円形瞳・インコヒーレント）")
ax.fill_between(u, 0, M, color="#1f77b4", alpha=0.12)
ax.axvline(u_abbe, color="crimson", ls="--", lw=2,
           label=f"アッベ限界 $d=0.5\\lambda$/NA（MTF=0）")
ax.axvline(u_ray, color="darkorange", ls="-.", lw=2,
           label=f"レイリー間隔 $s=0.61\\lambda$/NA（MTF={mtf_circ(u_ray):.3f}）")
ax.axvline(u_spa, color="green", ls=":", lw=2,
           label="スパロー間隔 $s=0.47\\lambda$/NA（帯域外）")
ax.set_xlim(0, 1.15); ax.set_ylim(0, 1.02)
ax.set_xlabel("正規化空間周波数 $u = f / f_c$（$f_c = 2$NA$/\\lambda$）")
ax.set_ylabel("コントラスト伝達率 MTF")
ax.set_title("光学系は $f_c=2$NA$/\\lambda$ のローパスフィルタ")
ax.legend(loc="upper right", fontsize=9)
ax.grid(alpha=0.3)
plt.tight_layout()
plt.show()

このグラフから3つのことが読み取れます。第一に、MTFは $u=1$（$f_c=2\mathrm{NA}/\lambda$）できれいにゼロに落ちる。ここがアッベの壁で、これより細かい周期は像に一切現れません。第二に、カットオフ直前ではMTFが急速に落ちるので、「限界ぎりぎりの構造は見えるが極端にコントラストが低い」という実感と一致します。第三に、レイリー間隔に対応する $u=0.82$ でのMTFは 0.089 しかなく、レイリー基準が「9%のコントラストでも2点と判定する」かなり緩い規約であることがわかります。スパロー間隔は $u=1.06$ と帯域外に出てしまいますが、これは前節で述べたとおり2点像が広帯域スペクトルを持つためで、矛盾ではありません。

(2) 2点光源の像プロファイル — レイリー・スパロー・未分解

次に、実際に2点の像がどう見えるかを描きます。間隔を変えて谷の深さがどう変化するかを追いましょう。

def two_points(x, vs):
    """間隔 vs（正規化）で置いた2つの非コヒーレント点像の強度和"""
    return airy(x - vs / 2.0) + airy(x + vs / 2.0)

# スパロー限界を数値的に求める（中心の2階微分がゼロになる間隔）
def curvature(vs, h=1e-4):
    f0 = two_points(np.array([0.0]), vs)[0]
    fp = two_points(np.array([h]), vs)[0]
    fm = two_points(np.array([-h]), vs)[0]
    return (fp - 2 * f0 + fm) / h ** 2

vs_sparrow = brentq(curvature, 1.0, 3.8)
print(f"スパロー間隔 v_s = {vs_sparrow:.4f} -> 係数 {vs_sparrow/(2*np.pi):.4f} λ/NA")
print(f"レイリー間隔 v_s = {v1:.4f} -> 係数 {v1/(2*np.pi):.4f} λ/NA")

出力は スパロー間隔 v_s = 2.9756 -> 係数 0.4736 λ/NA、レイリー間隔 v_s = 3.8317 -> 係数 0.6098 λ/NA です。「谷の2階微分がゼロ」という条件だけから、教科書の 0.47 という係数が数値的に再現されました。導出が正しかったことの確認になっています。

続いて、間隔を4通り変えたプロファイルを描きます。

x = np.linspace(-9, 9, 2000)
cases = [
    (v1 * 1.6,   "十分離れている（$1.0\\lambda$/NA）", "#2ca02c"),
    (v1,         "レイリー基準（$0.61\\lambda$/NA）",  "#ff7f0e"),
    (vs_sparrow, "スパロー限界（$0.47\\lambda$/NA）",  "#1f77b4"),
    (vs_sparrow * 0.7, "未分解（$0.33\\lambda$/NA）",  "#d62728"),
]

fig, axes = plt.subplots(2, 2, figsize=(11, 7), sharex=True)
for ax, (vs, label, c) in zip(axes.ravel(), cases):
    prof = two_points(x, vs)
    ax.plot(x, prof / prof.max(), color=c, lw=2.2, label="2点の合成像")
    ax.plot(x, airy(x - vs / 2) / prof.max(), color="gray", lw=1, ls="--")
    ax.plot(x, airy(x + vs / 2) / prof.max(), color="gray", lw=1, ls="--",
            label="個々の点像")
    dip = two_points(np.array([0.0]), vs)[0] / prof.max()
    ax.set_title(f"{label}\n中心/ピーク = {dip:.3f}", fontsize=10)
    ax.axvline(vs / 2, color="k", lw=0.6, alpha=0.4)
    ax.axvline(-vs / 2, color="k", lw=0.6, alpha=0.4)
    ax.grid(alpha=0.3); ax.legend(fontsize=8)
    ax.set_xlabel("像面位置 $v$（正規化）"); ax.set_ylabel("規格化強度")
plt.tight_layout()
plt.show()

4枚のプロファイルを見比べると、分解能の「連続的な劣化」がはっきり見えます。十分離れているとき（左上）は谷が深く、2つの点は疑いようもなく分かれています。レイリー基準（右上）では中心の強度がピークの 0.735、つまり26.5%のくぼみが残っており、「注意して見れば2つ」という程度です。スパロー限界（左下）では中心/ピークが 1.000 になり、プロファイルの頂上が完全に平らな台地になります。ここが「2つ」と主張できる最後の瞬間で、これより近づけると（右下）中心が最大値になり、単一のぼやけた点にしか見えません。実験で「分解できた」と言うときに、自分がどのレベルの谷を根拠にしているかを意識すべき理由が、この4枚に凝縮されています。

さらに、谷の深さを間隔の関数として連続的に描いておくと、3つの基準の位置関係が一目でわかります。

vs_list = np.linspace(1.0, 7.0, 400)
dips = []
for vs in vs_list:
    p = two_points(x, vs)
    dips.append(1.0 - two_points(np.array([0.0]), vs)[0] / p.max())
dips = np.array(dips)
coef = vs_list / (2 * np.pi)   # λ/NA 単位の間隔

fig, ax = plt.subplots(figsize=(9, 5))
ax.plot(coef, dips * 100, lw=2.5, color="#1f77b4")
ax.axvline(0.4736, color="green", ls=":", lw=2, label="スパロー限界 0.47")
ax.axvline(0.6098, color="darkorange", ls="-.", lw=2, label="レイリー基準 0.61")
ax.axvline(0.5, color="crimson", ls="--", lw=2, label="アッベ限界 0.50")
ax.axhline(26.5, color="gray", lw=1, ls="--")
ax.text(0.75, 28, "レイリー基準の谷の深さ 26.5%", fontsize=9, color="gray")
ax.set_xlabel("2点の間隔 $s$ [$\\lambda$/NA]")
ax.set_ylabel("谷の深さ（ピークに対する%）")
ax.set_title("2点像の谷の深さと3つの分解能基準")
ax.legend(); ax.grid(alpha=0.3)
plt.tight_layout()
plt.show()

この曲線から、分解能の基準が「どこで線を引くか」の問題であることが決定的にわかります。谷の深さは間隔とともに滑らかに増えていくだけで、どこにも不連続なジャンプはありません。スパロー限界（0.47）は谷の深さが0になる点、レイリー基準（0.61）は26.5%になる点、アッベ限界（0.50）は谷の深さとは無関係に周期構造の観点から引かれた線。3本の縦線が近接していることが、係数の細かい違いに本質的な意味がないことを物語っています。

(3) 構造化照明が周波数空間を折り返す

最後に、超解像の代表格である構造化照明顕微鏡（SIM）が、なぜアッベ限界の2倍の分解能を出せるのかを周波数空間で描きます。

fc = 1.0                       # 正規化カットオフ
f_pat = 0.95 * fc              # 照明パターンの空間周波数（帯域内が条件）
f = np.linspace(-2.5, 2.5, 1200)

def band(f, center):
    """中心 center に移動した通過帯域（MTF形状）"""
    u = np.abs(f - center) / fc
    out = np.zeros_like(f)
    m = u <= 1.0
    out[m] = mtf_circ(u[m])
    return out

fig, axes = plt.subplots(2, 1, figsize=(10, 7.5), sharex=True)

ax = axes[0]
ax.plot(f, band(f, 0.0), lw=2.5, color="#1f77b4", label="通常照明の通過帯域")
ax.axvspan(-fc, fc, color="#1f77b4", alpha=0.10)
ax.axvline(fc, color="crimson", ls="--", lw=2, label="アッベのカットオフ $f_c$")
ax.axvline(-fc, color="crimson", ls="--", lw=2)
ax.set_ylabel("伝達率"); ax.set_ylim(0, 1.05)
ax.set_title("通常照明：$|f|\\leq f_c$ しか観測できない")
ax.legend(fontsize=9); ax.grid(alpha=0.3)

ax = axes[1]
ax.plot(f, band(f, 0.0), lw=2.0, color="#1f77b4", label="0次成分（従来と同じ）")
ax.plot(f, band(f, +f_pat), lw=2.0, color="#2ca02c", ls="--",
        label="+1次成分（$+f_p$ だけ折り返し）")
ax.plot(f, band(f, -f_pat), lw=2.0, color="#ff7f0e", ls="--",
        label="-1次成分（$-f_p$ だけ折り返し）")
env = np.maximum.reduce([band(f, 0.0), band(f, f_pat), band(f, -f_pat)])
ax.plot(f, env, lw=1.2, color="k", alpha=0.5)
ax.axvline(2 * fc, color="purple", ls="--", lw=2, label="SIMの実効カットオフ $2f_c$")
ax.axvline(-2 * fc, color="purple", ls="--", lw=2)
ax.set_xlabel("空間周波数 $f$ [$f_c$ 単位]"); ax.set_ylabel("伝達率")
ax.set_ylim(0, 1.05)
ax.set_title("構造化照明：帯域が $\\pm f_p$ に複製され、実効カットオフが約2倍に")
ax.legend(fontsize=9); ax.grid(alpha=0.3)
plt.tight_layout()
plt.show()

上下2枚の比較が、SIMの原理そのものです。上段の通常照明では、観測できる情報は $|f|\le f_c$ の1つの帯域に限られます。下段では、周波数 $f_p$ の縞模様で照明することにより、物体の周波数成分が $f \pm f_p$ にシフトして帯域内に入り込みます。数学的には、照明パターンとの掛け算が周波数領域での畳み込みになるためで、AM変調とまったく同じ仕組みです。照明パターン自身も光学系を通らねばならないので $f_p \le f_c$ が必要であり、その最大値 $f_p=f_c$ を採ると、実効カットオフは $f_c + f_p = 2f_c = 4\mathrm{NA}/\lambda$、分解能にして $\lambda/(4\mathrm{NA})$ まで届きます。ちょうど2倍です。

ただし、生の観測データでは3つの成分が重なって混ざっているので、そのままでは分離できません。実際のSIMでは照明パターンの位相を3段階（0, 2π/3, 4π/3）ずらして3枚撮影し、3元1次方程式を解いて0次・±1次成分を分離します。さらに縞の向きを3方向変えることで、2次元の周波数空間を等方的に埋めます。合計9枚の撮影が必要になるのは、この「情報を余分に取ってから解く」構造のためです。

具体例 — 数字で実感する

導いた式に実際の数値を入れて、身の回りの光学機器がどのあたりにいるかを確認しましょう。すべて $d_{\min}=\lambda/(2\mathrm{NA})$ で計算します。

対数目盛で並べると、実体顕微鏡の 1100 nm からEUV露光の 20.5 nm まで、同じ $\lambda/(2\mathrm{NA})$ という1本の式が50倍以上の範囲をカバーしているのがわかります。青とオレンジの棒の比はどの装置でも一定（0.61/0.5＝1.22倍）で、基準の取り方の差が装置によらない定数倍でしかないことも読み取れます。棒の長さを決めているのは波長とNAだけで、倍率はどこにも効いていません。

この表と図からいくつも学びがあります。まず、実体顕微鏡（10倍・NA0.25）は生物顕微鏡の5倍以上粗い。倍率は10倍と100倍で10倍差ですが、分解能を決めるのはNAです。「倍率を上げても分解能は上がらない」というのは、単なる標語ではなく、式に倍率が一切現れないという事実の言い換えです。

次に、Blu-rayのスポット径 291nm に対して、規格上のトラックピッチが 320nm というのは示唆的です。記録密度は回折限界ぎりぎりまで攻めて設計されている。CD（$\lambda$=780nm, NA=0.45）からDVD（650nm, 0.60）、BD（405nm, 0.85）への進化は、まさに $\lambda$ を下げ NA を上げる戦いでした。

そしてリソグラフィ。ArF液浸の 71.5nm に対し、実際の半導体は10nm以下のパターンを作っています。これは矛盾ではなく、マルチパターニング（1つのパターンを複数回の露光に分けて描く）や位相シフトマスクなどの技法で、1回の露光の限界を回避しているからです。業界では $\mathrm{CD}=k_1\lambda/\mathrm{NA}$ という式を使い、$k_1=0.5$ がアッベ限界に相当します。$k_1$ を 0.25 まで下げる努力が、そのまま微細化の歴史です。

最後に、原理的な壁も確認しておきましょう。$\mathrm{NA}\le n$ なので、空気中では $d_{\min}\ge\lambda/2 = 275\,\mathrm{nm}$（緑色光）、水中でも $207\,\mathrm{nm}$。細胞小器官（ミトコンドリアの内膜構造は数十nm、シナプス小胞は約40nm）は、通常の光学顕微鏡では原理的に見えない。これが超解像顕微鏡が生命科学に決定的なインパクトを与えた理由です。

では、その超解像は何をしているのでしょうか。最後に、限界の「回避」の論理を整理します。

超解像は回折限界をどう回避するのか

まず強調したいのは、超解像技術はマクスウェル方程式を破っていないという点です。アッベの導出を振り返ると、暗黙の前提が3つ入っています。

1. 線形結像である（物体の明るさと像の明るさが比例する）
2. 伝搬する光波だけを使う（エバネッセント波は遠方に届かない）
3. 物体について事前知識がない（何が写っているかわからない）

超解像は、このいずれかの前提を外すことで限界を回避します。分類すると次のようになります。

(a) 帯域を拡げる — 構造化照明（SIM） 前節で見たとおり、既知の縞模様で照明し、周波数を折り返して帯域内に取り込みます。前提3を外す方法で、「照明パターンを知っている」という事前情報が効いています。線形なので改善は2倍が上限（$\lambda/(4\mathrm{NA})$）ですが、生細胞をそのまま高速に撮れるという実用上の強みがあります。

(b) 非線形性を持ち込む — STED、非線形SIM STED（誘導放出抑制）顕微鏡は、励起ビームの周りにドーナツ型の消光ビームを重ね、周辺の蛍光分子を強制的に消してしまいます。消光の効率が強度に対して飽和する非線形性を持つため、消光ビームを強くするほど「光り残る」中心部を任意に細くできます。実効的な分解能は $d \approx \lambda/(2\mathrm{NA}\sqrt{1+I/I_{\mathrm{sat}}})$ で、原理上の下限がありません。前提1（線形性）を外した例です。

(c) 時間で分ける — PALM / STORM（単一分子局在化） 1つの分子だけがぽつんと光っていれば、その像がぼやけたエアリーディスクでも、重心を精密に求めれば位置は高精度でわかります。局在化精度は光子数 $N$ に対しておよそ $\sigma/\sqrt{N}$ で改善します。PSFの幅を 240nm（$\sigma\approx102\,\mathrm{nm}$）とすると、10nm の精度に必要な光子数はわずか $N=(102/10)^2\approx104$ 個。実際には背景光や画素化で悪化するため数百〜数千光子を集めますが、いずれにせよ回折限界の200nm級から1桁以上細かい位置決めが可能です。ランダムに少数の分子だけを光らせては消す、を数万回繰り返し、位置を積み上げて画像を再構成します。これは「分解」ではなく「局在化」で、前提3を外した（＝点光源であるという事前知識を使う）方法です。

(d) 近接場を使う — NSOM / 固浸レンズ そもそもエバネッセント波が届かないことが $\lambda/2$ の根源でした。ならば、消える前に拾えばいい。物体から波長以下の距離に微小開口や探針を近づければ、エバネッセント成分を検出できます。前提2を外す方法で、原理的な分解能は探針サイズで決まり波長に依存しません。ただし表面しか見られず、走査に時間がかかります。

(e) 開口を合成する — フーリエタイコグラフィ、干渉計 複数の照明角度で撮った低NA画像を、周波数空間で貼り合わせて広い帯域を合成します。天文学の開口合成（VLBI）と同じ思想で、大口径レンズを作る代わりに情報を後から統合します。

3つの戦略は、外している前提がそれぞれ違います。(a)のSIMは帯域そのものを $2f_c$ まで拡げますが、照明パターンも光学系を通る以上 $f_p\le f_c$ なので改善は2倍で頭打ちです。(b)のSTEDは消光の飽和という非線形性を使うため、消光強度を上げるほど中心のピークが際限なく細くなり、原理的な下限がありません。(c)のPALM/STORMはそもそも「分解」をあきらめ、単一分子の重心位置を $\sigma/\sqrt{N}$ の精度で求める問題に読み替えており、数百光子で10nm級に届きます。

こうして並べると、回折限界は「情報の壁」であって「材料の壁」ではないことがよくわかります。瞳を通らなかった空間周波数の情報は、そのままでは永遠に失われる。だから超解像は必ず、どこかから追加の情報（照明パターン、非線形応答、疎性、近接場、複数視点）を持ち込んでいます。何も足さずに解像度だけ上がる手法があったら、それは疑ってかかるべきです。

なお、単純なデコンボリューション（PSFの逆フィルタ）は、$f_c$ 以下でMTFに削られたコントラストを回復するだけで、$f_c$ を超えた成分は復元できません（そこではMTFがゼロなので、割り算が定義できない）。SNRが良ければ「解析接続で少し外挿できる」という議論はありますが、実用上は帯域内のコントラスト回復と考えるのが安全です。ここを混同すると、「ソフトで超解像した」という主張の妥当性を見誤ります。

まとめ

本記事では、アッベの回折限界を回折格子の議論から導出し、レイリー基準・スパロー限界との関係をOTFの観点から統一的に整理しました。

• アッベの発想は「物体＝周期構造の重ね合わせ」「像＝2つ以上の回折次数の干渉」。細かい構造ほど回折角が大きく、瞳から外れた瞬間に情報がゼロになる
• 格子方程式 $n d(\sin\theta_m – \sin\theta_i)=m\lambda$ から、垂直照明では $d_{\min}=\lambda/\mathrm{NA}$、斜め照明では $d_{\min}=\lambda/(2\mathrm{NA})$ が導かれる。差の2倍は「0次光を瞳の端に寄せられるか」の違い
• レイリー基準 $0.61\lambda/\mathrm{NA}$ はエアリーディスク半径（$J_1$ の第1零点 3.8317 由来）で、谷の深さ26.5%という慣習的なしきい値。望遠鏡では $1.22\lambda/D$ と同じもの
• スパロー限界 $0.474\lambda/\mathrm{NA}$ は谷の2階微分がゼロになる条件で、最も攻めた基準
• 統一的な見方は OTFのカットオフ $f_c = 2\mathrm{NA}/\lambda$。光学系はローパスフィルタであり、3つの基準はいずれもこの $f_c$ から派生した「判定規約」の違いにすぎない
• $\lambda/2$ の壁は $\mathrm{NA}\le n$ から来る。それより細かい構造の情報はエバネッセント波が担い、遠方には届かない
• 超解像は限界を破るのではなく、線形性・伝搬波のみ・事前知識なし、という前提のどれかを外して回避している

回折限界は「レンズの性能の話」ではなく「情報がどこまで運ばれるか」の話です。この視点を持つと、光学だけでなく、サンプリング定理・アンテナの指向性・レーダーの角度分解能まで、同じ数学が姿を変えて現れていることに気づけるはずです。

次のステップとして、以下の記事も参考にしてください。

• フラウンホーファー回折とは？単スリットの強度分布を導出する — 本記事のエアリーパターンの土台
• 回折格子の分光分解能を導出する — 同じ格子方程式を「波長を分ける」側から使う
• 光学収差とは？ザイデルの5収差を理解する — 回折限界に到達する前に効いてくる、もう一つの像劣化要因
• 波動光学入門 — 光を波として扱うとは — 回折の出発点であるホイヘンス・フレネルの原理

参考文献

• E. Abbe, “Beiträge zur Theorie des Mikroskops und der mikroskopischen Wahrnehmung”, Archiv für Mikroskopische Anatomie, 9, 413–468 (1873)
• M. Born and E. Wolf, Principles of Optics, 7th ed., Cambridge University Press — 8章（回折理論）と9章（部分コヒーレンス）
• J. W. Goodman, Introduction to Fourier Optics, 4th ed., W. H. Freeman — 6章（周波数解析による結像系の解析）
• M. G. L. Gustafsson, “Surpassing the lateral resolution limit by a factor of two using structured illumination microscopy”, Journal of Microscopy, 198, 82–87 (2000)