光学と量子エレクトロニクス 第8章 §8.4 ★★★ 線形のままなら照明の工夫で2倍まで(構造化照明)。2準位の蛍光分子を誘導放出で飽和させると応答が非線形になって上限が消え(STED)、1度に1分子だけ光らせれば位置は光子数の平方根で細かく決まる(位置決め)
蛍光の顕微鏡では、物体は蛍光分子の数の密度 c ( x ) c(\bm{x}) c ( x ) で、分子は照明の光を吸って励起状態に上がり、少し長い波長の光を出して下りる。分子どうしは互いに無関係に光るので、像はインコヒーレントな像(式 (8.4) )である。この節では、分子を基底状態 ∣ g ⟩ | g \rangle ∣ g ⟩ と励起状態 ∣ e ⟩ | e \rangle ∣ e ⟩ の2準位系として扱う(実際の分子では振動の準位が細かく付いているが、それは吸収と放出の波長をずらす役だけを持つとみなす)。
構造化照明:線形のままの2倍
照明を一様にせず、波数 K i l l \bm{K}_{\mathrm{ill}} K ill の縞 1 + cos ( K i l l ⋅ x ) 1+\cos(\bm{K}_{\mathrm{ill}}\cdot\bm{x}) 1 + cos ( K ill ⋅ x ) にする。分子が出す光の強さは、弱い照明では照明の強度に比例するので、蛍光の分布は c ( x ) ( 1 + cos ( K i l l ⋅ x ) ) c(\bm{x})\bigl(1+\cos(\bm{K}_{\mathrm{ill}}\cdot\bm{x})\bigr) c ( x ) ( 1 + cos ( K ill ⋅ x ) ) で、そのスペクトルは c ~ ( q ) \tilde c(\bm{q}) c ~ ( q ) に 1 2 c ~ ( q ∓ K i l l ) \frac12\tilde c(\bm{q}\mp\bm{K}_{\mathrm{ill}}) 2 1 c ~ ( q ∓ K ill ) を足したものになる。像では ∣ q ∣ ≤ 2 k 0 N A \lvert \bm{q}\rvert \le2k_0\mathrm{NA} ∣ q ∣ ≤ 2 k 0 NA しか見えないが、ずらされた項のおかげで、物体のスペクトルの∣ q ∓ K i l l ∣ ≤ 2 k 0 N A \lvert \bm{q}\mp\bm{K}_{\mathrm{ill}}\rvert \le2k_0\mathrm{NA} ∣ q ∓ K ill ∣ ≤ 2 k 0 NA の部分も像に入る。縞の向きと位相を変えて何枚か撮り、ずらされた項を分けてもとの位置へ戻すと、物体のスペクトルが ∣ q ∣ ≤ 2 k 0 N A + ∣ K i l l ∣ \lvert \bm{q}\rvert \le2k_0\mathrm{NA}+\lvert \bm{K}_{\mathrm{ill}}\rvert ∣ q ∣ ≤ 2 k 0 NA + ∣ K ill ∣ までわかる。照明の縞も同じレンズで作るので ∣ K i l l ∣ ≤ 2 k 0 N A \lvert \bm{K}_{\mathrm{ill}}\rvert \le2k_0\mathrm{NA} ∣ K ill ∣ ≤ 2 k 0 NA で、見える範囲は最大で 4 k 0 N A 4k_0\mathrm{NA} 4 k 0 NA 、限界の2倍になる(M. G. L. Gustafsson (2000), J. Microsc. 198 , 82)。2倍で止まるのは、蛍光の分布が照明について線形で、照明そのものが定理の円板の中にあるからである。
STED:2準位の飽和で上限を外す
励起した分子が光を出す前に、蛍光より長い波長の強い光(STED 光)を当てると、分子は誘導放出で基底状態へ戻り、蛍光の波長の光を出さない。STED 光の強度が 0 0 0 になる点を中心に置いたドーナツ形の分布にすれば、中心のごく近くの分子だけが蛍光を出す。この方法は Hell と Wichmann が提案した(S. W. Hell and J. Wichmann (1994), Opt. Lett. 19 , 780)。
蛍光の寿命 τ f l \tau_{\mathrm{fl}} τ fl より十分短い STED 光のパルスを当てる場合を考える。誘導放出の断面積を σ \sigma σ 、光子の流れの密度を I ( t ) / ( ℏ ω ) I(t)/(\hbar\omega) I ( t ) / ( ℏ ω ) とすると、パルスの間の励起状態の占有確率 P e P_e P e は
d P e d t = − σ I ( t ) ℏ ω P e ⟹ η = P e ( 後 ) P e ( 前 ) = exp ( − F F s ) , F = ∫ I d t , F s = ℏ ω σ (8.9) \frac{\mathrm{d} P_e}{\mathrm{d} t}=-\frac{\sigma I(t)}{\hbar\omega}P_e
\quad\Longrightarrow\quad
\eta=\frac{P_e(\text{後})}{P_e(\text{前})}=\exp\Bigl(-\frac{F}{F_{\mathrm s}}\Bigr),\qquad F=\int I\,\mathrm{d} t,\quad F_{\mathrm s}=\frac{\hbar\omega}{\sigma}\tag{8.9} d t d P e = − ℏ ω σ I ( t ) P e ⟹ η = P e ( 前 ) P e ( 後 ) = exp ( − F s F ) , F = ∫ I d t , F s = σ ℏ ω ( 8.9 )
で減る(パルスの間の自然放出と、誘導放出で基底状態に戻った分子が STED 光で再び励起されることは無視した)。F F F はパルスのエネルギーの面密度(フルエンス)で、F s F_{\mathrm s} F s を飽和フルエンス と呼ぶ。ドーナツの中心の近くで F ( x ) = F max sin 2 ( π N A ∣ x ∣ / λ 0 ) F(\bm{x})=F_{\max}\sin^2(\pi\mathrm{NA}\,\lvert \bm{x}\rvert /\lambda_0) F ( x ) = F m a x sin 2 ( π NA ∣ x ∣ / λ 0 ) とし(中心で 0 0 0 、半径 λ 0 / ( 2 N A ) \lambda_0/(2\mathrm{NA}) λ 0 / ( 2 NA ) で最大。1次元の模型)、ζ = F max / F s \zeta=F_{\max}/F_{\mathrm s} ζ = F m a x / F s と置く。励起の点の像を半値全幅 d 0 = λ 0 / ( 2 N A ) d_0=\lambda_0/(2\mathrm{NA}) d 0 = λ 0 / ( 2 NA ) の Gauss 関数 exp ( − 4 ln 2 x 2 / d 0 2 ) \exp(-4\ln2\,x^2/d_0^2) exp ( − 4 ln 2 x 2 / d 0 2 ) で近似すると、蛍光を出す分子の分布(実効の点の像)は
h e f f ( x ) = exp ( − 4 ln 2 x 2 d 0 2 ) exp ( − ζ sin 2 π N A x λ 0 ) ≈ exp ( − 4 ln 2 x 2 d 2 ) , d = λ 0 2 N A 1 + ξ ζ , ξ = π 2 16 ln 2 = 0.890 (8.10) \begin{gathered}
h_{\mathrm{eff}}(x)=\exp\Bigl(-\frac{4\ln2\,x^2}{d_0^2}\Bigr)\exp\Bigl(-\zeta\sin^2\frac{\pi\mathrm{NA}\,x}{\lambda_0}\Bigr)
\approx\exp\Bigl(-\frac{4\ln2\,x^2}{d^2}\Bigr),\\
d=\frac{\lambda_0}{2\mathrm{NA}\sqrt{1+\xi\,\zeta}},\qquad \xi=\frac{\pi^2}{16\ln2}=0.890
\end{gathered}\tag{8.10} h eff ( x ) = exp ( − d 0 2 4 ln 2 x 2 ) exp ( − ζ sin 2 λ 0 π NA x ) ≈ exp ( − d 2 4 ln 2 x 2 ) , d = 2 NA 1 + ξ ζ λ 0 , ξ = 16 ln 2 π 2 = 0.890 ( 8.10 )
である。近似は sin 2 u ≈ u 2 \sin^2u\approx u^2 sin 2 u ≈ u 2 で、2つの Gauss 関数の積の幅を足し合わせた(1 / d 2 = 1 / d 0 2 + ζ π 2 N A 2 / ( 4 ln 2 λ 0 2 ) 1/d^2=1/d_0^2+\zeta\pi^2\mathrm{NA}^2/(4\ln2\,\lambda_0^2) 1/ d 2 = 1/ d 0 2 + ζ π 2 NA 2 / ( 4 ln 2 λ 0 2 ) )。ζ \zeta ζ を上げれば幅は 1 / ζ 1/\sqrt\zeta 1/ ζ で際限なく細くなる(図 8.3 )。係数 ξ \xi ξ は励起と STED の光の分布の形で変わり、定数を F s F_{\mathrm s} F s の定義に含めたd ≈ λ 0 / ( 2 N A 1 + F max / F s ) d\approx\lambda_0/(2\mathrm{NA}\sqrt{1+F_{\max}/F_{\mathrm s}}) d ≈ λ 0 / ( 2 NA 1 + F m a x / F s ) の形で書かれることが多い(V. Westphal and S. W. Hell (2005), Phys. Rev. Lett. 94 , 143903)。
STED では第二の前提が破れているので、定理 8.1 とは矛盾しない。STED 光の分布 F ( x ) F(\bm{x}) F ( x ) そのものは定理の円板の中の波数しか持たないが、η = e − ζ sin 2 ( ⋅ ) \eta=\mathrm{e}^{-\zeta\sin^2(\cdot)} η = e − ζ s i n 2 ( ⋅ ) は F F F の非線形の関数で、e − ζ sin 2 u \mathrm{e}^{-\zeta\sin^2u} e − ζ s i n 2 u を Fourier 級数に展開すると cos 2 m u \cos2mu cos 2 m u の項がすべての m m m について現れ、ζ \zeta ζ が大きいほど高い m m m の係数が大きくなる。実効の点の像は円板の外の波数を持ち、その大きさは2準位系の飽和の深さ ζ \zeta ζ で決まる。多重露光のしきい値(8.2 節)と同じく、応答が非線形なので円板の外の波数が現れる。
図 8.3 STED の実効の点の像(1次元の模型)。点線は式 (8.10) の左辺(sin 2 \sin^2 sin 2 を残した厳密な積)、細い実線は右辺の Gauss 関数の近似である。λ 0 = 600 n m \lambda_0=600\,\mathrm{nm} λ 0 = 600 nm , N A = 1.4 \mathrm{NA}=1.4 NA = 1.4 で、半値全幅は ζ = 0 \zeta=0 ζ = 0 で 214 n m 214\,\mathrm{nm} 214 nm 、ζ = 10 \zeta=10 ζ = 10 で 68.8 n m 68.8\,\mathrm{nm} 68.8 nm (近似 68.1 68.1 68.1 )、ζ = 100 \zeta=100 ζ = 100 で 22.6 n m 22.6\,\mathrm{nm} 22.6 nm (近似 22.6 22.6 22.6 )。
1分子の位置決め:事前の情報を使う
第三の前提は「物体について何も知らない」ことだった。視野の中で光っている分子が1つだけだと知っていれば、像の明るさの分布は点の像 ∣ h ∣ 2 \lvert h\rvert ^2 ∣ h ∣ 2 をその分子の位置 x 0 x_0 x 0 だけずらしたもので、見るべきものはスペクトルでなく、1つの数 x 0 x_0 x 0 である。x 0 x_0 x 0 は点の像の幅よりはるかに細かく決めることができる。点の像を標準偏差 s s s の Gauss 関数で近似する。Airy の点の像の半値全幅は 0.514 λ 0 / N A 0.514\lambda_0/\mathrm{NA} 0.514 λ 0 / NA (第6章の Airy の半値の点 k 0 N A r = 1.616 k_0\mathrm{NA}\,r=1.616 k 0 NA r = 1.616 から)で、同じ半値全幅の Gauss 関数はs = 0.514 λ 0 / ( 2.355 N A ) = 0.218 λ 0 / N A s=0.514\lambda_0/(2.355\,\mathrm{NA})=0.218\lambda_0/\mathrm{NA} s = 0.514 λ 0 / ( 2.355 NA ) = 0.218 λ 0 / NA である。検出した N p h N_{\mathrm{ph}} N ph 個の光子の位置は、それぞれ独立にこの分布に従う。1個の光子の位置の分布 p ( x ; x 0 ) ∝ exp ( − ( x − x 0 ) 2 / ( 2 s 2 ) ) p(x;x_0)\propto\exp(-(x-x_0)^2/(2s^2)) p ( x ; x 0 ) ∝ exp ( − ( x − x 0 ) 2 / ( 2 s 2 )) の Fisher 情報は J ( x 0 ) = 1 / s 2 J(x_0)=1/s^2 J ( x 0 ) = 1/ s 2 で、N p h N_{\mathrm{ph}} N ph 個では N p h / s 2 N_{\mathrm{ph}}/s^2 N ph / s 2 になるので、Cramér–Rao の不等式(相対エントロピーから情報と物理を組み立てる 第1部第7章 )から
V a r ( x ^ 0 ) ≥ s 2 N p h , V a r ( x ^ 0 ) ≥ s N p h (8.11) \mathrm{Var}(\hat x_0)\ge\frac{s^2}{N_{\mathrm{ph}}},\qquad \sqrt{\mathrm{Var}(\hat x_0)}\ge\frac{s}{\sqrt{N_{\mathrm{ph}}}}\tag{8.11} Var ( x ^ 0 ) ≥ N ph s 2 , Var ( x ^ 0 ) ≥ N ph s ( 8.11 )
で、光子の位置の平均を推定量にとると等号になる。背景の光や画素の大きさがあれば精度はこれより下がる。
近くにある多数の分子を分けるには、光で状態を切り替えられる蛍光分子を使い、1度に視野の中のまばらな一部だけを光らせて位置を決め、消してから別の一部を光らせることを繰り返す。集めた位置を1枚の図に描けば、分子の間隔が s s s より小さくても分子の並びを描くことができる。この方法は 2006 年に別々のグループから報告された(E. Betzig et al. (2006), Science 313 , 1642;M. J. Rust, M. Bates and X. Zhuang (2006), Nat. Methods 3 , 793)。STED と位置決めの開発に対して、2014 年の Nobel 化学賞が Betzig・Hell・Moerner に贈られた。
例題 8.5 STED の飽和の深さと位置決めの光子数 ★★
オリジナル λ 0 = 600 n m \lambda_0=600\,\mathrm{nm} λ 0 = 600 nm , N A = 1.4 \mathrm{NA}=1.4 NA = 1.4 とする。
(1) 励起の点の像の半値全幅 d 0 d_0 d 0 を求め、ζ = 100 \zeta=100 ζ = 100 の STED の実効の幅を式 (8.10) で求めよ。
(2) 同じ N A \mathrm{NA} NA の構造化照明で見える最小の周期を求め、(1) と比べよ。
(3) 1分子の位置決めで、N p h = 1000 N_{\mathrm{ph}}=1000 N ph = 1000 個の光子を集めたときの位置の精度の下限を式 (8.11) で求めよ。
解答 (1) d 0 = λ 0 / ( 2 N A ) = 214 n m d_0=\lambda_0/(2\mathrm{NA})=214\,\mathrm{nm} d 0 = λ 0 / ( 2 NA ) = 214 nm 。ζ = 100 \zeta=100 ζ = 100 では 1 + 0.890 × 100 = 9.49 \sqrt{1+0.890\times100}=9.49 1 + 0.890 × 100 = 9.49 で d = 22.6 n m d=22.6\,\mathrm{nm} d = 22.6 nm 、励起だけの像の約 9.5 9.5 9.5 分の 1 1 1 になる。sin 2 \sin^2 sin 2 を残した積の半値全幅を数値で求めても 22.6 n m 22.6\,\mathrm{nm} 22.6 nm で、近似の誤差は図 8.3 で見えないほど小さい。
(2) 見える波数の上限は 4 k 0 N A 4k_0\mathrm{NA} 4 k 0 NA なので、最小の周期は λ 0 / ( 4 N A ) = 107 n m \lambda_0/(4\mathrm{NA})=107\,\mathrm{nm} λ 0 / ( 4 NA ) = 107 nm 。STED の 22.6 n m 22.6\,\mathrm{nm} 22.6 nm よりはるかに粗い。構造化照明は線形のまま2倍にとどまり、STED は非線形の深さ ζ \zeta ζ で幅が決まる。
(3) s = 0.218 × 600 / 1.4 = 93.6 n m s=0.218\times600/1.4=93.6\,\mathrm{nm} s = 0.218 × 600/1.4 = 93.6 nm で、s / 1000 = 2.96 n m s/\sqrt{1000}=2.96\,\mathrm{nm} s / 1000 = 2.96 nm 。点の像の幅の 30 30 30 分の 1 1 1 程度の精度で位置が決まる。光子の数を 4 4 4 倍にしても精度は2倍にしかならない。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com