(57)【要約】【課題】軟Ｘ線集光光学系の高精度で効率的なアライメント及び形状評価が可能な高精度波面計測手法を提供せんとする。【解決手段】高ＮＡ集光素子の焦点面上においてピンホールを走査し、ピンホールの光源と反対側である下流側に設置された撮像手段により，前記ピンホールの走査で変化する回析光の強度分布を計測し、前記計測される回析光の強度分布および各走査位置の情報に基づき、前記高ＮＡ集光素子の焦点面と前記撮像手段の観察面との間で、タイコグラフィ位相回復法による光の伝播・逆伝播計算の繰り返しにより焦点面上の光の波動場分布を求め、該焦点面上の光の波動場に基づき、回析積分の逆計算で前記高ＮＡ集光素子の出口波面を求めることを特徴とする。【選択図】図６
2017037002.tif 000002
【特許請求の範囲】
【請求項１】
高ＮＡ集光素子の出口波面計測方法であって、
高ＮＡ集光素子の焦点面上においてピンホールを走査し、
ピンホールの光源と反対側である下流側に設置された撮像手段により，前記ピンホールの走査で変化する回析光の強度分布を計測し、
前記計測される回析光の強度分布および各走査位置の情報に基づき、前記高ＮＡ集光素子の焦点面と前記撮像手段の観察面との間で、タイコグラフィ位相回復法による光の伝播・逆伝播計算の繰り返しにより焦点面上の光の波動場分布を求め、
該焦点面上の光の波動場に基づき、回析積分の逆計算で前記高ＮＡ集光素子の出口波面を求めることを特徴とする高ＮＡ集光素子の出口波面計測方法。
【請求項２】
前記タイコグラフィ位相回復法による光の伝播・逆伝播計算において、観察面上の設定している測定強度固定の領域の外側に位相とともに強度も変更可能な領域を設定するとともに、前記集光素子の出口形状寸法を拘束条件として加える請求項１記載の高ＮＡ集光素子の出口波面計測方法。
【請求項３】
前記タイコグラフィ位相回復法による光の伝播・逆伝播計算において、
前記計測される回折光の強度分布を、下記［１］〜［７］に示す高ＮＡ近似計算法の（ξ、η）座標系に補間した形で出力し、
該出力結果を用いて、下記式（８）の（ξ、η）座標系でタイコグラフィ位相回復法による光の伝播・逆伝播計算を繰り返し行う請求項１又は２記載の高ＮＡ集光素子の出口波面計測方法。
［１］  観察面Ｄ上の点（ｕ，ｖ，ｚ）における複素波動場を、Reighley-Sommerfeld第Ｉ種回折公式の式（１）のＵ（ｕ，ｖ，ｚ）と定義する。Ｕ（ｘ，ｙ，０）は焦点面Ｓ上の波動場とする。
【数１】
2017037002.tif 000024
焦点面に（ｘ，ｙ）座標系を設定
観察面に（ｕ，ｖ）座標系を設定
各座標系の原点を貫いて光が伝播する方向をＺ軸の正
λは光の波長
［２］  （ｘ，ｙ，０）から点（ｕ，ｖ，ｚ）までの距離は式（２）で表わされる。
【数２】
2017037002.tif 000025
［３］  新たに焦点面のピンホール中心（０，０，０）から観察面（ｕ，ｖ，ｚ）までの距離ｒ２を式（３）で定義する。
【数３】
2017037002.tif 000026
［４］  振幅係数について、
【数４】
2017037002.tif 000027
を用いて
【数５】
2017037002.tif 000028
と近似する。
［５］  位相係数について、式（４）及びｒ２＞ｚから求まる近似式（５）を設定するとともに焦点面上の２次項を無視し、これを近似条件とする。
【数６】
2017037002.tif 000029
［６］  以上から式（６）が得られる。ここで式（７）とおいて、式（８）のフーリエ変換の形の式が得られる。
【数７】
2017037002.tif 000030
［７］  （ｕ、ｖ）座標系と（ξ、η）座標系の補間の式は、下記式（９）、式（１０）となる。
【数８】
2017037002.tif 000031
【請求項４】
高ＮＡ集光素子の出口波面計測システムであって、
ピンホールを有する走査物体と、
前記走査物体を、前記高ＮＡ集光素子の焦点面上において走査する走査駆動装置と、
前記走査物体の光源と反対側である下流側に設置され、前記ピンホールの走査で変化する回析光の強度分布を計測する撮像装置と、
前記撮像手段で計測される回析光の強度分布および各走査位置の情報に基づき、前記高ＮＡ集光素子の焦点面と前記撮像手段の観察面との間で、タイコグラフィ位相回復法による光の伝播・逆伝播計算の繰り返しにより焦点面上の光の波動場分布を求め、該焦点面上の光の波動場に基づき、回析積分の逆計算で前記高ＮＡ集光素子の出口波面を求める演算装置と、
を備えてなることを特徴とする高ＮＡ集光素子の出口波面計測システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、軟Ｘ線集光素子のアライメントや形状誤差測定に好適に用いることのできる高ＮＡ集光素子の出口波面計測方法、及び出口波面計測システムに関する。
【背景技術】
【０００２】
従来、Ｘ線を用いて物質や現象を分析するにあたっては、その強度密度ならびに集光サイズが非常に重要であり、集光光学系の高性能化が常に必要とされている。
【０００３】
例えば楕円関数を回転した形状の回転楕円ミラーは、ＮＡが大きく軟Ｘ線を高効率かつシングルナノメートルに集光することができる優れたＸ線集光素子である（非特許文献１）。しかしながら、いまだ理想的な性能を持つ回転楕円ミラーは実用化されていない。理想的な回転楕円ミラーの実現には、反射面である内面の３次元非球面形状を高精度に計測する手法の確立が求められている。
【０００４】
また、回転楕円ミラーの発展形として、上流のミラーにより軟Ｘ線ビームをリング状に拡大し、拡大した軟Ｘ線を下流の回転体ミラーにより集光する２段集光システムが提案されている（特許文献１参照）。このシステムによれば、理想的な回折限界性能１０ｎｍでロスのない高効率な集光ビームの形成が可能となる。しかし、多数のミラーを反射させるこのシステムは、光学設計の自由度が高い反面、アライメントが複雑で且つ時間がかかる。アライメントの為のＸ線の利用には時間的に限界があり、可視光レーザを用いたオフ2ライン時での高精度のアライメント及び光学系評価が求められている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１５−１７９５７号公報
【非特許文献】
【０００６】
【非特許文献１】Takahiro Saito, Yoshinori Takei, Hidekazu Mimura, Development of Surface Profile Measurement Method for Ellisoidal X-ray Mirror using Phase retrieval, Proc. SPIE, 8501, 850103, 2012
【発明の概要】
【発明が解決しようとする課題】
【０００７】
そこで、本発明が前述の状況に鑑み、解決しようとするところは、軟Ｘ線集光光学系の高精度で効率的なアライメント及び形状評価が可能な高精度波面計測手法を提供する点にある。
【課題を解決するための手段】
【０００８】
本発明者は、従来顕微鏡の手法として知られているタイコグラフィ（Ptychography）による位相回復法（H. M. L. Faulkner and J. M. Rodenburg, Movable aperture lensless transmission microscopy: a novel phase retrieval   algorithm, Phys. Rev. Lett. 93, 023903, 2004.）に着目し、集光素子の焦点面上にピンホールを走査させることを着想した。そして、このようにタイコグラフィ法を用いることにより，強度が直接測定できない波動場上（ここではＡ（ｘ，ｙ））もしくは近傍において，光の強度もしくは位相を変調させる透過関数(Ｏ（ｘ，ｙ）)をもつ強度遮蔽物体や位相物体を走査し、走査毎に変化するもう片方の波動場の（ここではＵ（ｕ，ｖ））の強度分布を撮像手段で測定することで，情報量を増やし，波動場Ａ（ｘ，ｙ）と走査させた物体の透過関数を求め、さらに回折積分の逆計算をすることで集光素子の出口波面を求めることができることを見出し、本発明を完成するに至った。
【０００９】
すなわち本発明は、高ＮＡ集光素子の出口波面計測方法であって、高ＮＡ集光素子の焦点面上においてピンホールを走査し、ピンホールの光源と反対側である下流側に設置された撮像手段により、前記ピンホールの走査で変化する回析光の強度分布を計測し、前記計測される回析光の強度分布および各走査位置の情報に基づき、前記高ＮＡ集光素子の焦点面と前記撮像手段の観察面との間で、タイコグラフィ位相回復法による光の伝播・逆伝播計算の繰り返しにより焦点面上の光の波動場分布を求め、該焦点面上の光の波動場に基づき、回析積分の逆計算で前記高ＮＡ集光素子の出口波面を求めることを特徴とする高ＮＡ集光素子の出口波面計測方法を提供する。
【００１０】
ここに、焦点面とは集光点又はその近傍で光軸に対して垂直な面をいい、撮像手段が計測する光の強度分布は、前記ピンホールから広がる回析光を撮像した回析像の強度プロファイルをいう。本明細書では、光源側を上流、撮像手段側を下流とする。また、計算される焦点面上の「波動場分布」とは、回復される光の複素波動場であり、光波動場の回析積分の逆計算は、具体的には、Reighley-Sommerfeld第Ｉ種回析公式に波数ｋ＝−２π／λを代入することで求められる。出口波面は光波動場（強度分布及び波面）として求められる。
【００１１】
本発明におけるタイコグラフィ位相回復法による光の伝播・逆伝播計算の概要について説明する。前記焦点面上のピンホールを通過する光の透過関数は、一般的に複素数であらわされ，通過した時の光の強度の減衰、位相変化量を与える。透過関数は、次の式（１８）、式（１９）のように強度変調項と位相変調項に分離される。
【００１２】
【数１】
2017037002.tif 000003
【００１３】
ピンホールの中心位置を（ｔ，ｓ）とすると、位置（ｔ，ｓ）での透過関数はＯ（ｘ−ｔ，ｙ−ｓ）で表現される。波動場分布Ａ（ｘ，ｙ）のもつビームが透過関数Ｏ（ｘ−ｔ，ｙ−ｓ）を通過した時のビームの波動場Ｂ（ｘ，ｙ，ｔ，ｓ）は式（２０）と表現される。
【００１４】
【数２】
2017037002.tif 000004
【００１５】
透過関数Ｏ（ｘ−ｔ，ｙ−ｓ）が既知であると、ピンホール位置（ｔ，ｓ）の時に繰り返し計算で用いるＵ面（観察面）の波動場、強度、位相をそれぞれＵ（ｘ，ｙ，ｔ，ｓ）、Ｉ（ｘ，ｙ，ｔ，ｓ）、φ（ｘ，ｙ，ｔ，ｓ）とおき、測定される強度をＩexp（ｘ，ｙ，ｔ，ｓ）とおく。このときタイコグラフィの繰り返し計算方法を図１に示す。図ではピンホール中心場所をｐ１（ｔ１，ｓ１）とｐ２（ｔ２，ｓ２）としている。測定されるＩexp（ｘ，ｙ，ｔ，ｓ）の強度データを用いて、ピンホールの各場所において二つの波動場Ｂ（ｘ，ｙ，ｔ，ｓ）とＵ（ｘ，ｙ，ｔ，ｓ）の間で繰り返し計算を行う。
【００１６】
すなわち、まず走査物体の場所がｐ１のとき、図のようにＩexp（ｘ，ｙ，ｔ１，ｓ１）のデータを用いてＢ（ｘ，ｙ，ｔ１，ｓ１）とＵ（ｘ，ｙ，ｔ１，ｓ１）の間で繰り返し計算を行う。算出される強度分布を既知強度データに置き換える拘束をかけ、繰り返し計算を実施する。Ｂ面は強度データが未知でありＵ面の強度データのみの拘束であるが、ある一定値に収束する。そして、収束したＢ（ｘ，ｙ）の値を用いて式（２０）からＯ（ｘ−ｔ，ｙ−ｓ）を利用してＡ（ｘ，ｙ）を求める。Ａ（ｘ，ｙ）が求まると既知であるＯ（ｘ−ｔ２，ｙ−ｓ２）を用いて、走査物体の場所がｐ２のときのＢ（ｘ，ｙ，ｔ２，ｓ２）の初期値を仮に設定する。
【００１７】
次に、Ｉexp（ｘ，ｙ，ｔ２，ｓ２）を用いてＢ（ｘ，ｙ，ｔ２，ｓ２）とＵ（ｘ，ｙ，ｔ２，ｓ２）の間で繰り返し計算して収束させる。この結果、再び求めたＡ（ｘ，ｙ）から走査物体の場所がｐ１の時のＢ（ｘ，ｙ，ｔ１，ｓ１）を求め、繰り返し計算の初期値と利用し，Ｂ（ｘ，ｙ，ｔ１，ｓ１）とＵ（ｘ，ｙ，ｔ１，ｓ１）の間で繰り返し計算を行う。このサイクルを全体が収束するまで実施すると最終的に求まった波動場分布Ａ（ｘ、ｙ）は測定された二つの強度データＩexp（ｘ，ｙ，ｔ１，ｓ１）とＩexp（ｘ，ｙ，ｔ２，ｓ２）を満たす波動場分布となる。
【００１８】
Ａ（ｘ，ｙ）を既知とし，透過関数であるＯ（ｘ，ｙ）を未知にしたときも同じ考え方でＯ（ｘ，ｙ）を収束させることができる。このように本発明では、タイゴグラフィー位相回復法により走査物体の各場所で繰り返し計算を実施し，その収束した値を別の場所における繰り返し計算の初期値に利用することで測定されたすべての強度情報を満たすように収束させることができる。ピンホールは２次元的に走査し、走査回数を多くすることで大量の強度データから波動場分布を求めることができる。図２に走査回数Ｍ回の時のタイコグラフィ位相回復アルゴリズムを示す。データが多くなると、はじめにＡ（ｘ，ｙ）を既知、Ｏ（ｘ，ｙ）を未知として収束計算を実施し、Ｏ（ｘ，ｙ）を求め、求めたＯ（ｘ，ｙ）を固定してＡ（ｘ，ｙ）を求めることができ、これを繰り返すことでＡ（ｘ，ｙ），Ｏ（ｘ，ｙ）の両方を求めることができる。
【００１９】
焦点面上で計算される波動場の分解能は、撮像手段の観察面上での領域の広さに相当する。焦点面上での分解能を向上するためには、フーリエ変換の性質を考えると観察面を大きくする必要があるが、撮像手段の限界がある。そこで本発明者は、図３に示すように、観察面上の設定している測定強度固定の領域の外側に、位相とともに強度も変更可能な領域を設定するとともに、集光光学素子の周りの強度がゼロであることを利用し、模式図４に示すように、タイコグラフィ位相回復法の反復計算中において、前記集光素子の出口形状寸法の波動場を計算し、出口周りを強度ゼロに置き換えて焦点面に戻すループを導入した。
【００２０】
集光素子の出口周りの強度がゼロである設定は、焦点面上において外側の強度を低減させるように作用する。その結果、撮像手段の観察面上において測定強度固定外の領域の強度が大きくなることを防ぐことができる。これにより、図５に示すように測定強度よりも外側の領域を含めた波動場が回復し、焦点面上での回復される波動場の空間分解能を向上させることができた。同時に、ピンホールを表す透過関数の分解能も向上した。
【００２１】
また、前記タイコグラフィ位相回復法による光の伝播・逆伝播計算において、前記計測される回折光の強度分布を、下記［１］〜［７］に示す高ＮＡ近似計算法の（ξ、η）座標系に補間した形で出力し、該出力結果を用いて、下記式（８）の（ξ、η）座標系でタイコグラフィ位相回復法による光の伝播・逆伝播計算を繰り返し行うことが好ましい。これによりタイコグラフィ位相回復法の計算を汎用の高速フーリエ変換ソフトで効率よく行うことが可能であり、且つ強度の回復精度を向上させることができる。
【００２２】
［１］  観察面Ｄ上の点（ｕ，ｖ，ｚ）における複素波動場を、Reighley-Sommerfeld第Ｉ種回折公式の式（１）のＵ（ｕ，ｖ，ｚ）と定義する。Ｕ（ｘ，ｙ，０）は焦点面Ｓ上の波動場とする。
【００２３】
【数３】
2017037002.tif 000005
焦点面に（ｘ，ｙ）座標系を設定
観察面に（ｕ，ｖ）座標系を設定
各座標系の原点を貫いて光が伝播する方向をＺ軸の正
λは光の波長
【００２４】
［２］  （ｘ，ｙ，０）から点（ｕ，ｖ，ｚ）までの距離は式（２）で表される。
【００２５】
【数４】
2017037002.tif 000006
【００２６】
［３］  新たに焦点面のピンホール中心（０，０，０）から観察面（ｕ，ｖ，ｚ）までの距離ｒ２を式（３）で定義する。
【００２７】
【数５】
2017037002.tif 000007
【００２８】
［４］  振幅係数について，
【００２９】
【数６】
2017037002.tif 000008
【００３０】
を用いて
【００３１】
【数７】
2017037002.tif 000009
【００３２】
と近似する。
【００３３】
［５］  位相係数について、式（４）及びｒ２＞ｚから求まる近似式（５）を設定するとともに焦点面上の２次項を無視し、これを近似条件とする。
【００３４】
【数８】
2017037002.tif 000010
【００３５】
［６］  以上から式（６）が得られる。ここで式（７）とおいて、式（８）のフーリエ変換の形の式が得られる。
【００３６】
【数９】
2017037002.tif 000011
【００３７】
［７］  （ｕ、ｖ）座標系と（ξ、η）座標系の補間の式は、下記式（９）、式（１０）となる。
【００３８】
【数１０】
2017037002.tif 000012
【００３９】
また本発明は、高ＮＡ集光素子の出口波面計測システムであって、ピンホールを有する走査物体と、前記走査物体を、前記高ＮＡ集光素子の焦点面上において走査する走査駆動装置と、前記走査物体の光源と反対側である下流側に設置され、前記ピンホールの走査で変化する回析光の強度分布を計測する撮像装置と、前記撮像手段で計測される回析光の強度分布および各走査位置の情報に基づき、前記高ＮＡ集光素子の焦点面と前記撮像手段の観察面との間で、タイコグラフィ位相回復法による光の伝播・逆伝播計算の繰り返しにより焦点面上の光の波動場分布を求め、該焦点面上の光の波動場に基づき、回析積分の逆計算で前記高ＮＡ集光素子の出口波面を求める演算装置とを備えてなることを特徴とする高ＮＡ集光素子の出口波面計測システムをも提供する。
【発明の効果】
【００４０】
以上にしてなる本願発明に係る高ＮＡ集光素子の出口波面計測方法は、軟Ｘ線集光光学系の高精度で効率的なアライメント及び形状評価を可能とする。すなわち、可視光レーザで効率のよい高精度なアライメントを行ったうえで、波面誤差をもとに複雑な非球面３次元形状を持つ回転楕円ミラーの内面全面の形状誤差を計測・評価することが可能である。この形状計測結果はマンドレル作成や電鋳法といったミラー作成プロセスの改良に利用できるだけでなく、回転楕円ミラーの形状修正加工にも活用できる。また、軟Ｘ線集光素子のアライメント、形状評価手法としてだけでなく，その計測結果を活用することによって従来にない理想的な軟Ｘ線集光光学系を提供することにつながり、Ｘ線光学分野に大きく貢献するものである。
【図面の簡単な説明】
【００４１】
【図１】走査物体の位置が２か所のときのタイコグラフィ計算の概要を示す説明図。
【図２】走査物体が位置がＭか所のときのタイコグラフィ位相回復法のアルゴリズムを示す説明図。
【図３】位相回復計算における焦点面での分解能の向上を示す図。
【図４】位相回復計算における拘束条件の追加を示す説明図。
【図５】ＣＣＤ面上における振幅分布の回復を示す図（領域サイズ：３６．８６４×３６．８６４ｍｍ）。
【図６】本発明の代表的実施形態に係るタイコグラフィ位相回復による波面計測システムの概略構成を示す説明図。
【図７】タイコグラフィ位相回復計算における座標系を示す模式図。
【図８】焦点面上の波動場からの逆伝播計算により求めた出口波面分布を示す図。
【図９】タイコグラフィ位相回復計算における強度情報更新順序の改善を示す説明図。
【図１０】想定する回転楕円ミラー下流開口での振幅位相分布を示す図（外径：５．２２９ｍｍ，内径：２．８５９ｍｍ）
【図１１】計算により求めたＣＣＤ面上の仮想強度（一部を抜粋）を示す図（領域サイズ：１８．４３２×１８．４３２ｍｍ）。
【図１２】ループ回数と入力強度データと回復強度の差の２乗和の関係を示すグラフ。
【図１３】回復された焦点面上での強度と位相（ＣＣＤ面のデータはピンホール中心が焦点位置にあるとき）を示す図。
【図１４】使用した回転楕円ミラーの写真。
【図１５】アライメント調整による回転楕円ミラーの集光波面誤差の改善を示す図。
【図１６】アライメント回数と波面誤差のコマ収差成分のＲＭＳ値の関係を示すグラフ。
【図１７】同一条件で測定された回転楕円ミラーの波面誤差プロファイルを示す図。
【図１８】回転楕円ミラーを回転させる前後の波面計測結果を示す図（外径４．６ｍｍ，内径２．９ｍｍ）。
【図１９】波面誤差プロファイルのゼルニケ解析結果を示すグラフ。
【図２０】下流開口波面と回転楕円ミラーの対応を示す図。
【図２１】波面誤差から決定した回転楕円ミラーの形状誤差を示す図であり、（ａ）はミラー内面全面の形状誤差を示す図、（ｂ）はミラー上流端から４ｍｍ断面における真円度を示す図、（ｃ）は長手方向の形状誤差プロファイルを示すグラフをそれぞれ示している。
【図２２】真円度測定装置により計測された回転楕円ミラーの真円度を示す図。
【発明を実施するための形態】
【００４２】
次に、本発明の実施形態を添付図面に基づき詳細に説明する。
【００４３】
本発明に係るタイコグラフィ位相回復法を用いた集光素子の波面計測法は、集光素子に対して理想的なコヒーレント球面波を入射させ、その集光点のさらに後方に撮像手段（例えばＣＣＤカメラ）を設置し、集光点付近で光軸に対して垂直な面（焦点面）内でピンホールを走査し、投影される回折像の強度プロファイル変化を観察し、タイコグラフィ位相回復計算によって、焦点面上の光の複素波動場を回復させ、さらに光波動場の回折積分を逆計算することで、集光素子の下流開口における光波動場（強度分布及び波面）を求めるものである。
【００４４】
図６に本発明に係る高ＮＡ集光素子の出口波面計測システムの一例を示す。本例は、回転楕円ミラーに対してタイコグラフィ位相回復を行うために構築された光学系であり、全体は暗室内部の防振台の上に固定されている。光は図中左上のレーザから出射され、最終的に図中右下の撮像手段であるＣＣＤカメラに入射する。光源として使用したのは直線偏光を持つ波長６３２．８ｎｍのＨｅ−Ｎｅレーザである。
【００４５】
レーザから発せられるガウシアンビームは、まずＮ／Ｄフィルタによって出力を調整され、その後、直径５０μｍのピンホールに入射される。集光素子の出口波面誤差と形状誤差を対応付けるためには、入射するコヒーレント光が点光源から発せられた理想的な球面波であることが好ましく、ピンホールは光が球面波に近づく位置に設けられる。集光素子に入射した光は、その下流に存在する焦点に集光される。
【００４６】
タイコグラフィ位相回復計算では、上記した高ＮＡ近似計算法の（ξ、η）座標系で光の伝播・逆伝播計算を繰り返し行うことが好ましく、このために焦点近傍に挿入し走査される走査物体は透過面積が小さいことが好ましく、そのピンホールは直径２０μｍに設定されている。走査物体は、走査駆動装置である並進３軸の駆動ステージに固定され、焦点面上において走査される。
【００４７】
集光素子によって焦点付近に集められた光は、走査物体のピンホール内側のみを通過し，その回折パターンが下流側のＣＣＤカメラによって観察される。ＣＣＤカメラは、線形性が保証される冷却ＣＣＤカメラが好適であり、本例ではBitran社（http://www.bitran.co.jp/）製のＢＵ−５３ＬＮが使用されている。ＢＵ−５３ＬＮの仕様を表１に示す。
【００４８】
【表１】
2017037002.tif 000013
【００４９】
駆動ステージに固定した走査物体（ピンホール）を焦点面上で走査しながら、順次ＣＣＤカメラによって強度プロファイルを取得する。光は焦点面上全体に広がっているが，その強度は強度中心に著しく集中している。したがって、ピンホールの内側に強度中心を含まない場合はＣＣＤカメラに入射する光の強度が非常に弱くなる。したがって、走査位置に応じて露光時間を１００ｍｓから３００００ｍｓまで計測可能な最大強度に対して飽和しない範囲で変化させている。ＣＣＤカメラの計測値は露光時間に対して線形性を持つため、位相回復計算において取得強度を露光時間によって割ることで、同一条件で撮像を続けた場合と同様の結果として用いることができる。
【００５０】
演算装置によるタイコグラフィ位相回復法を用いた光の伝播・逆伝播計算は、本例では、次のように行われる。図７に示すように、焦点面上に原点０ｆの直交座標系（ｘ，ｙ）、ＣＣＤカメラのピクセルアレイ上に原点０ｃの直交座標系（ｕ，ｖ）をとり、両座標系の原点はｚ軸（光軸）が垂直に貫き、その距離をｆ２とする。焦点面上にピンホールを中心が０ｐ（ｘ＝Ｘ、ｙ＝Ｙ）に位置するよう設け、その上で焦点面上にｘ２ｙ２直交座標系（ｘ２＝ｘ−Ｘ、ｙ２＝ｙ−Ｙ）、ＣＣＤ面上にｕ２ｖ２直交座標系（ｕ２＝ｕ−Ｕ，ｖ２＝ｖ−Ｖ）を原点がピンホールの中心を通るようにそれぞれ再設定する。
【００５１】
ｘｙ平面の波動場をＡ（ｘ２，ｙ２）、ｕｖ平面の波動場をＵ（ｘ２，ｙ２）とし、焦点面に挿入した走査物体が波動場に与える影響を複素関数Ｏ（ｘ２，ｙ２）とする。ピンホールの場合、この関数はピンホールの内側においては「１」、外側においては「０」の値をとる実数関数となる。焦点面においてピンホール内部を透過した光の波動場は，Ｂ（ｘ２，ｙ２）を用いて式（１１）で表わされる。これが出口波動場となる。
【００５２】
【数１１】
2017037002.tif 000014
【００５３】
Ｕ（ｕ２，ｖ２）の波動場はＢ（ｘ２，ｙ２）からのみ影響を受け、Ｂ（ｘ２，ｙ２）のうち光が実際に存在する領域はピンホールの内側のみであるから、高ＮＡ近似回折計算を適用し、ＣＣＤ面上の波動場は以下の式（１２）ように表される。
【００５４】
【数１２】
2017037002.tif 000015
【００５５】
ｘ２ｙ２平面の逆空間であるξη平面を以下の式（１３）に設定することで、式（１２）は出口波動場Ｂ（ｘ２，ｙ２）のフーリエ変換Ｆ［Ｂ（ｘ２，ｙ２）］（ξ,η）を用いて簡略化される。
【００５６】
【数１３】
2017037002.tif 000016
【００５７】
式（１３）において両辺の強度のみに注目すると、式（１４）、式（１５）となる。
【００５８】
【数１４】
2017037002.tif 000017
【００５９】
式（１５）からわかるように、ピンホールの中心を原点とする焦点面上波動場Ａ（ｘ２，ｙ２）にピンホールの透過関数Ｏ（ｘ２，ｙ２）を掛け合わせた出口波動場Ｂ（ｘ２，ｙ２）のフーリエ変換の強度はＣＣＤカメラによる撮像結果によって再生される。できることがわかる。このＣＣＤ撮像結果を元にはじめに式（１６）を計算し，これを拘束条件として反復計算を行う。
【００６０】
【数１５】
2017037002.tif 000018
【００６１】
この際、ＣＣＤ面上の撮像範囲の２倍のサイズを現す行列を用意し、撮像範囲内のピクセルにのみ拘束条件をかけて計算を行う。焦点面におけるピクセルサイズはＣＣＤ面での行列の大きさに反比例するため、このことは焦点面上の分解能を１／２にする効果がある。これにより，より緻密な焦点面上の波動場が得られる。一方で，撮像範囲外のピクセルは拘束条件が存在しないため、実験上発生し得ない強度が回復されるが、それを防ぐために本手法では，集光素子の下流開口直径を拘束条件として加えている。
【００６２】
このようにして求められた焦点面上の光の波動場から、さらに集光素子の下流開口における波動場を取得するために、焦点面上の波動場に対して回折積分の逆計算を行う。具体的には、式（１）、式（２）のReighley-Sommerfeld第Ｉ種回折公式を用いて波数ｋ＝−２π／λを代入することによって得られる。回復される波動場はミラーに反射した光のみであり、回転楕円ミラーは筒状の形状をしており，下流開口における波面は図８のような輪帯状となる。
【００６３】
焦点面に挿入される走査物体がピンホールという非常に透過面積が小さい物体であるため、ミラーに反射せず内部を通過した球面波の一部はほとんどが遮られ、ＣＣＤカメラに届かない。したがって、回復した波面は反射光のみの情報を持っていることになる。図８の輪帯状の下流開口波動場は、外周側が回転楕円ミラーのうち下流側で反射した光を表わし，内周側がミラーの上流側で反射した光を表している。
【００６４】
走査物体の走査手順は、特に限定されないが、ピンホールの位置が焦点面上を渦巻き方式で中心部から外周部に向かって走査し、この順で得られるＣＣＤ強度情報を元に反復計算を行うことが好ましい。例えばピンホールを図９（ａ）に示すようにｘｙ方向にラスタースキャンしながら測定し、その順序でＣＣＤ強度情報の更新を行った場合、ある確率で図９（ｂ）のように焦点面上波動場に位相飛び（光渦）が発生する。これはピンホール位置を強度の低い外周部から集光点に向けて移動しながら計算して行くことで前列までの位相回復結果と矛盾が蓄積することが原因と考えられる。位相回復過程を見ると光渦の場所は外側に移動していき最終的に消滅するが、計算の収束に時間を要することとなる。また光渦の領域は強度がゼロになるためにＣＣＤ面上においても局所的に正確に位相が回復されない問題がある。
【００６５】
これに対し、図９（ｃ）に示すように渦巻き方式で強度情報の更新を行うと、タイコグラフィ計算では焦点面で強度が高い中心付近にピンホールがあるときほどＣＣＤ強度情報の信頼性が高いことから、信頼性が高い情報をもとに先に収束させながら信頼性が低い情報を採用して行くことになり，回復波動場に矛盾が生じにくく、前述のように矛盾が一点に収束してしまった場合でも，光渦領域を外側に向けて迅速に押し出して行くことができ、光渦発生による位相回復の停滞を防止し、効率よく演算を行うことができる。
【００６６】
次に、本例の高ＮＡ集光素子の出口波面計測システムを用いてタイコグラフィ位相回復シミュレーションを行った結果について説明する。
【００６７】
本シミュレーションは、回転楕円ミラーの出口波面を仮想的に与えてＣＣＤカメラで測定される強度データを理論的に求めたうえ、ＣＣＤ面と焦点面の間の反復計算により位相回復を行い、焦点面波動場から出口波面を求める。そして、仮想的に与えた上記出口波面と位相回復法により求められた出口波面とを比較する。想定した回転楕円ミラーの出口波面を図１０に示す。また、表２にシミュレーションにおいて想定した実験条件を示す。ピンホールの各位置において、波動場を順方向の回折積分によってＣＣＤ面まで伝播させ、図１１に示す仮想のＣＣＤカメラでの測定データを作成する。この仮想測定データはピンホールのスキャン回数の数だけ強度データが存在する。これら仮想測定データを用いて、上記した式（１１）〜のタイコグラフィ位相回復を行い、焦点面の波動場、出口波面を求める。
【００６８】
【表２】
2017037002.tif 000019
【００６９】
図１２に全体のループ回数と，式（１７）で計算される入力強度データと反復計算の中での回復強度プロファイルの一致度を示す。このように約２０回のループ回数を行うことで計算がほぼ収束していることがわかる。
【００７０】
【数１６】
2017037002.tif 000020
【００７１】
また、図１３に位相回復計算において出力された焦点面の位相と振幅、ＣＣＤ観察面での振幅プロファイルを示す。図８に位相回復計算の最終結果を示す。入力した波面誤差と比較した結果、差異はＲＭＳで０．００２１λであり、高精度に波面決定できることがわかる。
【００７２】
次に、本例の高ＮＡ集光素子の出口波面計測システムを用いて軟Ｘ線集光用回転楕円ミラーの形状・アライメント誤差の計測実験を行った結果について説明する。
【００７３】
表３および図１４に使用した軟Ｘ線集光用回転楕円ミラーの仕様と写真を示す。本ミラーは，形状誤差約５０ｎｍ（Ｐ−Ｖ）の精度で作製されたマンドレルを，久米らが開発した常温の電析出条件を用いた高精度電鋳法により形状転写を行ったものであり、内面の真円度としてミラー中央部で約１００ｎｍ（Ｐ−Ｖ）である。
【００７４】
【表３】
2017037002.tif 000021
【００７５】
ミラーの設置は、シャックハルトマンセンサーを用いた波面計測により、入射ビームの曲率半径を見ながら回転楕円ミラーの設計値であるＲ=６４４０ｍｍの位置にミラー上流端を設置する。ＣＣＤカメラに映る回転楕円ミラーの反射光とミラー中心を通る球面波の中心を一致させ、ミラーの姿勢をおおよそビームに対して平行にする。また、ＣＣＤカメラに映る反射光が対称性の良いリング状となるように姿勢を調整する。そして、演算装置による波面回復結果から、ゼルニケ解析によりコマ収差成分を算出し、入射角度誤差を予測してコマ収差成分が低減するようにミラーの姿勢をアライメント調整する。この演算、アライメントを数回繰り返す。
【００７６】
表４にその他の実験条件を示す。回転楕円ミラーの形状計測の長手方向の分解能５ｍｍを達成するためには出口波面において分解能０．１４８ｍｍを実現する必要がある。これは焦点面上において直径７３μｍの範囲でスキャンを行う必要があることを示している。そこで本実験ではピンホールのスキャン範囲を９３μｍとしている。
【００７７】
【表４】
2017037002.tif 000022
【００７８】
（アライメント調整結果）
図１５にアライメントを行った結果を示す。このように各アライメントを実施することにより、傾きに起因するコマ収差成分が低減していることがわかる。アライメント回数と波面誤差のコマ収差成分の関係を図１６に示す。アライメント調整４回でコマ収差成分がＲＭＳで０．００２０λ（１．２８ｎｍ）に低減している。このことはミスアライメントに起因する波面誤差を想定する軟Ｘ線の波長４．２ｎｍ以下に抑えられたことを示している。アライメント精度としてミラーの傾きに換算すると１．３９μｒａｄに相当し、可視光とピンホールのみで構成される簡便な光学系によってこの精度が得られたことは非常に有用である。アライメントが終了した段階で得られた波面誤差はＲＭＳで０．０２９７λ（λ＝６３２．８ｎｍ），Ｐ−Ｖで０．２７３λであった。
【００７９】
（波面誤差プロファイルの再現性、確からしさの評価）
最終的に得られた波面誤差プロファイルの再現性を評価するために、連続２回測定を行った。図１７に同一条件において行った２回の波面誤差プロファイルを示す。２回連続測定における波面誤差プロファイルの差の２乗平均平方根は０．００８２１λ（λ=６３２．８ｎｍ）であった。また、評価対象とするゼルニケ解析における８１次項以下の成分でみると０．００４９１λであり、高い再現性を示している。
【００８０】
図１８に回転楕円ミラーを光軸周りに９０度回転させた時に得られた波面誤差プロファイルを示す。この計測では，回転楕円ミラーの焦点位置、ピッチング、ヨーイングの角度調整を再度行っており，完全に独立した計測となっている。回転角度の決定は波面誤差プロファイルと同時に回復される下流開口上の強度プロファイルの形状をもとに行う。輪帯状の強度プロファイルの外周円は回転楕円ミラーの下流開口のエッジ形状を反映しているため、凹凸形状のフィッティングを行うことにより１°レベルでの回転角度決定が可能である。
【００８１】
図１８に示すように、ミラーの光軸中心の回転に対して連動して回転する波面誤差成分と回転しない成分が存在する。回転する成分は、ミラーの形状誤差に起因する成分である。回転しない成分は、ミラーに入射する前の波面における波面誤差を含むシステムエラー成分であり、非点収差成分がシステムエラーとして多く含まれていると考えられる。図１８（ｃ），（ｄ）はミラー回転前後の非点収差を表しており，これらはほぼ変化していない。非点収差は９０度の回転によって反転するという特徴があるため、この成分は明らかにミラーに起因するものでないことがわかる。
【００８２】
非点収差を含んだままの波面誤差同士の一方を回転させ、ミラーの向きを揃えたうえで比較すると、ＲＭＳで０．０４２０λ（λ＝６３２．８ｎｍ）もの差異があった。他方、双方の非点収差の平均分布をシステムエラーとして同定し、回転前後の波面誤差分布から除去した結果を図１８（ｅ），（ｆ）に示す。システムエラーを除けば互いの波面誤差プロファイルは似通っていることがわかる。
【００８３】
非点収差成分を除去した回転前後それぞれの波面誤差プロファイルのゼルニケ解析結果を図１９に示す。球面収差を示すｚ９が良く一致している。光軸周りに９０°回転した前後の波面誤差プロファイルの差分の２乗平均平方根は０．０１２８λ（λ=６３２．８ｎｍ）であった。また、評価対象とするゼルニケ解析における８１次項以下の成分でみると０．００７８４λ（λ=６３２．８ｎｍ）であった。
【００８４】
表５に，再現性／確からしさ評価をまとめた結果を示す。ミラーを入射波面に対して回転させてシステムエラーを分離することにより、ミラーの形状に起因する波面誤差を計測できた。以上の結果から、回転楕円ミラーによる波面計測における確からしさはλ／１００レベルであるといえる。
【００８５】
【表５】
2017037002.tif 000023
【００８６】
（回転楕円ミラーの３次元形状誤差プロファイルの算出）
計測された波面誤差プロファイルから回転楕円ミラーの３次元形状誤差プロファイルを算出する。光源からミラー反射面，下流開口断面へと進む光線は、回転楕円ミラーの形状誤差が小さい場合、重なることはない。そのためミラー反射面と出口波面は１対１で対応すると近似することができる。図２０に示すように，下流開口における波面誤差の半径方向成分と周方向成分は、それぞれ回転楕円ミラーの長手方向形状誤差と周方向形状誤差に相当する。図２１（ａ）にミラー内面全面の形状誤差をプロットした。この図よりミラー上流端１０ｍｍ及び下流端５ｍｍの領域にそれぞれ形状誤差が集中していることがわかる。
【００８７】
また、図２１（ｂ）に回転楕円ミラー周方向のプロファイルを示す。これは、上流端から４ｍｍの距離の断面における形状プロファイルの真円からの差を表示している。この図より、上流端付近で真円度は１μｍを超えており、楕円状の形状を持つことがわかる。他方、中央部における断面をとると１００ｎｍ以下の真円度が得られる。また、図２１（ｃ）に回転楕円ミラー長手方向のプロファイル断面を示す。ミラー中央部３０ｍｍに限定すれば、長手方向においてもＰ−Ｖ１００ｎｍ前後の形状誤差プロファイルが得られる。
【００８８】
この形状誤差プロファイルの妥当性を評価するために、真円度測定器（小坂研究所製：ＥＣ１５５０Ｈ）によりミラー内面の真円度を計測した。結果を図２２に示す。プロファイルを比較すると、扁平方向及びくびれ位置においてよく似た傾向がみられ、真円度の値もオーダーが一致していることがわかる。また、ミラー中央部における真円度測定器による計測結果も１００ｎｍ以下であり、波面計測法に基づく値とオーダーが一致した。これらより本波面計測法がミラー形状誤差に基づく情報を回復できたと確認できた。
【００８９】
以上、本発明の実施形態について説明したが、本発明はこうした実施例に何ら限定されるものではなく、軟Ｘ線集光光学系以外の高ＮＡ集光素子の出口波面計測、アライメント調整、形状誤差測定にも広く適用でき、本発明の要旨を逸脱しない範囲において種々なる形態で実施し得ることは勿論である。
【図１】
2017037002.tif 000032
【図２】
2017037002.tif 000033
【図３】
2017037002.tif 000034
【図４】
2017037002.tif 000035
【図５】
2017037002.tif 000036
【図６】
2017037002.tif 000037
【図７】
2017037002.tif 000038
【図８】
2017037002.tif 000039
【図９】
2017037002.tif 000040
【図１０】
2017037002.tif 000041
【図１１】
2017037002.tif 000042
【図１２】
2017037002.tif 000043
【図１３】
2017037002.tif 000044
【図１４】
2017037002.tif 000045
【図１５】
2017037002.tif 000046
【図１６】
2017037002.tif 000047
【図１７】
2017037002.tif 000048
【図１８】
2017037002.tif 000049
【図１９】
2017037002.tif 000050
【図２０】
2017037002.tif 000051
【図２１】
2017037002.tif 000052
【図２２】
2017037002.tif 000053