二光子吸収計算機 - TPA係数の計算

二光子吸収計算機

二光子吸収計算機は、3つのレーザー設定（波長、強度、パルス幅）を、βと表記されcm/GWで示される単一の例示的な数値に変換します。二光子吸収とは、分子が1個ではなく、ほぼ同時に2個の光子を吸収する過程であり、高い光強度でのみ起こります。このページが生成する数値は計算方法を学ぶための補助であり、実在する物質の測定値や予測係数ではありません。その理由を以下で説明します。

二光子吸収とは？

通常の単光子吸収では、分子は光子を一個吸収して、より高いエネルギー状態へ遷移します。二光子吸収（TPA）は異なり、分子は互いにフェムト秒以内の間隔で二個の光子を吸収し、それらの合計エネルギーを使って同じ遷移を起こします。フェムト秒は一秒の十億分の一の、さらに百万分の一です。

この二重吸収はまれです。2個の光子が同じ分子にほぼ同時に到達する必要があるため、意味を持つほど頻繁に起こるのは光が極めて強い場合だけです。通常の光は、太陽光でさえ、はるかに弱すぎます。エネルギーを短く強いパルスに集中させるレーザーによって、TPAの観測と利用が可能になります。

物理学者マリア・ゲッパート＝メイヤーは、博士論文の一部として1931年に二光子吸収を予測しました。当時の光源は、それを示すのに十分な強度を持っていませんでした。1960年に発明されたレーザーが状況を変えました。1961年、ベル研究所の研究者ウィルヘルム・カイザーとC. G. B. ギャレットは、ルビーレーザーとユウロピウムを添加した結晶を使い、初めて二光子吸収を観測しました。ゲッパート＝メイヤーは後に、核構造に関する別の研究により1963年のノーベル物理学賞を共同受賞しました。

二光子吸収係数の式

計算機は、次の固定された式を適用します。

$\beta = K \times \frac{I \times \tau}{\lambda^2}$

ここで、

• β — 二光子吸収係数（単位：cm/GW）
• K — 固定定数。この計算機では1.5に設定
• I — 光の強度。W/cm²からGW/cm²へ変換（10⁹で割る）
• τ — パルス幅（単位：フェムト秒（fs））
• λ — 光の波長（単位：ナノメートル（nm））

この式が物理法則ではなく仮の式である理由

計算機もこの記事も、その点を明示しています。実際の光学では、βは物質固有の性質です。これは通常のベール・ランベルト吸収に対応する二光子版の式 dI/dz = −βI² に現れ、物質のバンドギャップ、線形屈折率、カネエネルギー、光子エネルギーだけに依存します。ビーム強度とパルス幅はこの式に現れません。レーザー出力を2倍にしても、物質のβは2倍になりません。

上の式では強度とパルス幅が増えるとβも増えるため、物理的なこの量を表しているとはいえません。また、単位もcm/GWには換算されません。GW/cm² × fs ÷ nm²は、長さを電力で割った単位ではありません。したがって、この数値は計算方法の例示であり、計算機にはその旨の警告が常時表示されます。

実際のβを得るには、通常はZスキャン法による直接測定、またはSheik-Bahae、Hagan、Van Strylandによる二バンドスケーリング則のような物質モデルが必要です。β = K·√E_p·F₂(2ħω/E_g) / (n₀²·E_g³)。ここでE_gはバンドギャップ、n₀は線形屈折率、E_pはカネエネルギー、F₂は比2ħω/E_gの固定関数です。このモデルには、この計算機が入力として求めない物質パラメーターが必要です。

各変数の挙動

**波長（λ）**は分母で2乗されるため、波長が長くなると数値は逆二乗関係に従って低下します。強度とパルス幅を変えずに800 nmから1000 nmへ変えると、数値は (800/1000)² = 0.64 倍となり、半分を超える低下ではなく36%の低下になります。二光子実験の多くでは、近赤外光が生体組織によく浸透するため、およそ700–1000 nmの光が使われます。この理由から、800 nm付近で波長可変なチタンサファイア（Ti:Sapphire）レーザーが最も一般的な光源です。

**強度（I）**は、実験室で二光子吸収が顕著になるにはおよそ10¹⁰～10¹⁴ W/cm²に達する必要があり、計算機はこの範囲外の値に警告を出します。地表での太陽光はわずか0.1 W/cm²程度で、弱すぎます。この強度に達するには、エネルギーを連続的に広げるのではなく、短いバーストに集中させるパルスレーザーが必要です。

**パルス幅（τ）**は、TPA実験では通常10～1000フェムト秒の範囲です。短いパルスは、一定量のエネルギーをより短い時間に集中させるためピーク強度を高めますが、より複雑で高価なレーザー装置も必要になります。

計算例

代表的なチタンサファイアレーザーのパラメーターを使って式を計算すると、

• 波長：800 nm
• 強度：10¹² W/cm²
• パルス幅：100 fs
• K = 1.5

まず、強度をGW/cm²に変換します。10¹² ÷ 10⁹ = 1000 GW/cm²。

次に式を適用します。

β = 1.5 × (1000 × 100) ÷ 800² = 150,000 ÷ 640,000 = 0.234375。計算機では 0.2344 cm/GW と表示されます。

強度とパルス幅を同じに保ち、波長を1000 nmに変えると、

β = 1.5 × (1000 × 100) ÷ 1000² = 150,000 ÷ 1,000,000 = 0.15。0.1500 cm/GW と表示されます。

これは36%の低下であり、上で説明した逆二乗関係と一致します。

計算機に表示される内容

ページには、係数がcm/GW単位で1つの数値として表示されます。0.001～999の値は小数点以下4桁で表示され、それより小さい値や大きい値は、たとえば1.5000 × 10^-5のような科学的記数法に切り替わります。その下のパネルには、変換後の強度を含む計算過程が段階的に繰り返し表示され、手計算で確認できます。

入力欄の上には常に警告ボックスがあり、この数値は検証済みの物理係数ではないと示します。各入力は個別に検査されます。0以下の値はエラーとして拒否され、正の値が通常の実験範囲外にある場合は範囲警告が表示されますが、結果は生成されます。結果の横には、ビームが試料に入射して出射する様子を示す図があり、現在の波長、強度、係数が表示されます。

二光子吸収の応用

**二光子顕微鏡法。**吸収はレーザービームが強く集束された場所でのみ起こるため、この技術では試料内部の深い単一点にある蛍光分子だけを励起できます。焦点外の領域を損傷または退色させにくいため、生きた脳組織の画像化に広く使われています。ウィンフリート・デンク、ジェームズ・ストリックラー、ワット・ウェッブは、1990年にコーネル大学で初の二光子顕微鏡を製作しました。

**光線力学療法。**近赤外光は可視光よりも組織の深部に到達します。二光子励起を使うと、光感受性薬剤をレーザーの焦点でのみ活性化できるため、周囲の正常組織を温存しながら腫瘍を標的にできます。

**3D光学データ記録。**物質のブロック内部でレーザーを異なる深さに集束させると、その深さだけで二光子吸収によってデータを記録できます。これにより、DVDのような表面型フォーマットより高い記録密度が得られます。

**二光子微細加工。**二光子重合とも呼ばれるこの技術では、レーザーの焦点でのみ感光性樹脂を硬化させ、3Dプリント構造を、100ナノメートル未満の特徴寸法で作製します。

**光リミッティング。**βが高い物質は通常の光強度では透明なままですが、強いレーザーパルスの下では強く吸収します。この性質は、レーザーの偶発的な曝露から保護する保護眼鏡やセンサーフィルターに利用されています。

吸収係数を手計算する方法

1. 光の強度を10⁹で割り、W/cm²からGW/cm²に変換します。
2. その結果に、フェムト秒単位のパルス幅を掛けます。
3. 定数K（この簡略モデルでは1.5）を掛けます。
4. 波長をナノメートル単位で二乗した値で割ります。

計算機はこの同じ4段階を実行し、標準的な実験範囲外の入力に警告を出します。範囲は、波長が400–1200 nm、強度が10¹⁰–10¹⁴ W/cm²、パルス幅が10–1000 fsです。これらの範囲外の値でも数値は生成されます。

よくある質問

二光子吸収とは何ですか。 分子が1個ではなく、ほぼ同時に2個の光子を吸収して、より高いエネルギー状態に達する過程です。非常に高い光強度でのみ顕著になります。

係数βは何を測定しますか。 光学では、βは物質がどの程度強く二光子吸収を起こすかをcm/GW単位で表します。βが高いほど、一定の強度で物質はより強く吸収します。公表値には幅があり、溶融シリカは800 nm付近で約10⁻⁶ cm/GW、ヒ化ガリウムは1064 nmでおよそ26 cm/GWです。このページが出力する数値は、これらの測定値のいずれでもありません。

なぜ二光子吸収にはこれほど強い光が必要なのですか。 2個の光子が互いにフェムト秒以内に分子へ当たる必要があり、通常の光強度ではその可能性が非常に低いからです。パルスレーザーはエネルギーを短いバーストに集中させ、観測と利用が可能な程度までその可能性を高めます。

この計算機では波長が結果にどう影響しますか。 波長が長くなると、結果は波長の二乗に反比例して低下します。他の設定を変えずに800ナノメートルから1000ナノメートルへ移すと、上の計算例では0.234375から0.15へ、36%低下します。

この式は三光子吸収に使えますか。 いいえ。三光子吸収には固有の係数があり、吸収されるパワーは強度の二乗ではなく三乗に比例して増加します。別のモデルが必要です。

この簡略式はどの程度正確ですか。 正確ではなく、正確であることを意図した式でもありません。物質の性質を含まず、単位もcm/GWには換算されないため、出力を物質の係数として引用してはいけません。実際の用途には、測定値またはSheik-Bahaeの二バンド式のような物質モデルを使用してください。

レーザーの出力を上げると、実際の係数は変化しますか。 いいえ。βは物質と波長の性質であるため、強度やパルス幅が変化しても同じです。変化するのは物質が吸収する光の量であり、二光子吸収によって強度の2乗に比例する形でエネルギーが失われるためです。

参考文献

1. ゲッパート＝メイヤー, M.（1931）「Über Elementarakte mit zwei Quantensprüngen」Annalen der Physik, 401（3）, 273–294.
2. カイザー, W., & ギャレット, C. G. B.（1961）「CaF₂: Eu²⁺における二光子励起」Physical Review Letters, 7（6）, 229–231.
3. デンク, W., ストリックラー, J. H., & ウェッブ, W. W.（1990）「二光子レーザー走査蛍光顕微鏡法」Science, 248（4951）, 73–76.
4. ルーミ, M., & ペリー, J. W.（2010）「二光子吸収：測定と原理の概説」Advances in Optics and Photonics, 2（4）, 451–518.
5. シェイク＝バハエ, M., ヘイガン, D. J., & ヴァン・ストリランド, E. W.（1990）「二光子吸収に関連する固体の電子カー効果の分散とバンドギャップスケーリング」Physical Review Letters, 65（1）, 96–99.
6. シェイク＝バハエ, M., ハッチングス, D. C., ヘイガン, D. J., & ヴァン・ストリランド, E. W.（1991）「固体における束縛電子非線形屈折の分散」IEEE Journal of Quantum Electronics, 27（6）, 1296–1309.