(57)【要約】  本発明は、結像装置、結像方法および結像システムを開示し、サンプル画像データの取得及び結像技術の分野に属する。結像装置は、荷電粒子源、収束システム、走査制御システム、検知モジュール、及び検知モジュールの下方に配置されているスペクトル分析モジュールを含む。検知モジュールは、複数の画素化探知器ユニットを含み、穴が開設されている。本発明に係る結像装置は、穴が開設された検知モジュール及びスペクトル分析モジュールを含み、穴が開設された検知モジュールにより回折パターンを取得し、スペクトル分析モジュールは、穴を通る荷電粒子線に対してスペクトル分析を行い、一回の走査で回折パターン及びスペクトル信号を同時に取得することを可能にする。本発明に係る結像方法は、中空積層結像方法に基づいたのであり、穴が開設された検知モジュールで得られた回折パターンを結像させることを実現でき、結像の効果が良好である。【選択図】図１
2020537123.tif 000002
【特許請求の範囲】
【請求項１】
荷電粒子を放出するための荷電粒子源（１）と、
荷電粒子線を制限し収束するための収束システム（２）と、
荷電粒子線のサンプルでの走査を制御するための走査制御システム（３）と、
サンプル（４）と、
荷電粒子を取り込み、荷電粒子の信号強度を検出するための検知モジュール（５）と、
検知モジュール（５）の下方に配置され、荷電粒子の分光学特徴を分析するためのスペクトル分析モジュール（６）と、を含み、
前記検知モジュール（５）は、複数の画素化探知器ユニット（７）を含み、当該検知モジュール（５）には、穴（８）が開設されている
ことを特徴とする結像装置。
【請求項２】
前記穴（８）は、円形状、方形状又は環形状を含むが、これらに限定されないことを特徴する請求項１に記載の結像装置。
【請求項３】
前記検知モジュール（５）は、方形状、円形状および島状を含むが、これらに限定されないことを特徴とする請求項１又は２に記載の結像装置。
【請求項４】
回折パターンを取得する工程と、
中空積層結像方法により、取得した回折パターンに対して積層結像を行う工程と、
を含むことを特徴とする結像方法。
【請求項５】
前記回折パターンを取得する工程は、結像装置に基づくものであり、前記結像装置は、荷電粒子源（１）、収束システム（２）、走査制御システム（３）、サンプル（４）、および検知モジュール（５）を含み、検知モジュール（５）は、複数の画素化探知器ユニット（７）を含み、穴（８）が開設されており、
前記結像方法は、
Ａ：収束システム（２）によって、荷電粒子源（１）により放出された荷電粒子線をサンプル（４）の表面に収束する工程と、
Ｂ：走査制御システム（３）によって荷電粒子線を制御してサンプル（４）を走査させる工程と、
Ｃ：荷電粒子線がサンプル（４）を透過して検知モジュール（５）に到達し、検知モジュール（５）における画素化探知器ユニット（７）は、対応する走査位置における荷電粒子の信号強度を検出し、対応する走査位置の回折パターンを取得する工程と、
Ｄ：荷電粒子線がサンプル（４）を走査する時、走査ビームスポットに対応する走査領域Ｒｉと他の走査領域Ｒｊには、非空集合Ｒｏｖｅｒｌａｐが存在し、
【数１】
2020537123.tif 000023
であり、Ｎが荷電粒子線によるサンプル（４）の走査領域の総数である工程と、
を含むことを特徴とする請求項４に記載の結像方法。
【請求項６】
中空積層結像方法により、取得した回折パターンに対して積層結像を行う工程は、具体的には、以下の工程を含むことを特徴とする請求項４又は５に記載の結像方法。
ａ：Ｐ（ｒ）をプローブ関数とし、Ｏ（ｒ）を物体の複素振幅分布関数とし、複数回の反復計算によってＯ（ｒ）を再構成し、最終的な反復計算で再構成するＯ（ｒ）を、物体の最終的な複素振幅分布関数とし、積層結像は、物体の最終的な複素振幅分布関数に基づいて画像を再構成する
ｂ：Ψｎ，ｍを物体を透過する出射波関数とし、Ψｎ，ｍを、プローブ関数Ｐ（ｒ）と物体の複素振幅分布関数Ｏ（ｒ）との乗積に定義すると、
【数２】
2020537123.tif 000024
式（１）において、ｎはＯ（ｒ）のｎ回目の反復を表し、ｍは荷電粒子線のサンプルでのｍ番目の走査位置を表し、Ｒｍは、荷電粒子線のサンプルでの、１番目の走査位置に対するｍ番目の走査位置の相対的座標ベクトルを表し、ｒは空間座標である
ｃ：物体を透過する出射波関数Ψｎ，ｍに対してフーリエ変換を行うことにより、当該出射波関数Ψｎ，ｍのファーフィールドでの振幅および位相分布を得る
【数３】
2020537123.tif 000025
式（２）において、
【数４】
2020537123.tif 000026
は出射波関数Ψｎ，ｍのファーフィールドでの振幅を表し、Φｎ，ｍは出射波関数Ψｎ，ｍのファーフィールドでの位相を表す
ｄ：実験装置によってサンプルを透過する出射波のファーフィールド光強度を採集し、Ｉｍとし、Ｉｍは、荷電粒子線がサンプルにおけるｍ番目の走査位置を走査した後にサンプルを透過する出射波のファーフィールド光強度を表す
ｅ：制約関数Ｍを設定する
ｆ：出射波関数Ψｎ，ｍのファーフィールドでの振幅
【数５】
2020537123.tif 000027
に代えて、
【数６】
2020537123.tif 000028
を制約関数Ｍに代入し、式（３）を得る
【数７】
2020537123.tif 000029
ｇ：Ψｎ，ｍ，ｎｅｗ’（ｒ）に対して逆フーリエ変換を行い、式（４）で示される新たな出射波関数Ψｎ，ｍ，ｎｅｗ（ｒ）を得る
【数８】
2020537123.tif 000030
ｈ：新たな出射波関数Ψｎ，ｍ，ｎｅｗ（ｒ）により計算し、式（５）で示される新たな物体の複素振幅分布関数を得る
【数９】
2020537123.tif 000031
式（５）において、αおよびβは調節可能なパラメータであり、αは、分母が０ではないことを確保するために用いられ、βは、帰還の強度を制御するために用いられる
ｉ：新たな物体の複素振幅分布関数Ｏｎ＋１（ｒ＋Ｒｍ）を式（１）に代入して再度の反復計算を開始する
ｊ：最終的な反復計算で得られる物体の複素振幅分布関数Ｏｚ（ｒ＋Ｒｍ）を画像再構成の基礎とし、Ｏｎ＋１（ｒ＋Ｒｍ）の振幅および位相を描画して再構成画像を得る。
【請求項７】
前記制約関数Ｍは、検知モジュール（５）の構造に関連することを特徴とする請求項６に記載の結像方法。
【請求項８】
Ｏ（ｒ）の１回目の反復において、物体の複素振幅分布関数Ｏ０（ｒ＋Ｒｍ）をランダム分布関数にすることを特徴とする請求項７に記載の結像方法。
【請求項９】
請求項１〜３のいずれか一項に記載の結像装置および請求項４〜８のいずれか一項に記載の結像方法を含むことを特徴とする結像システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、サンプル画像データの取得および結像技術の分野に属し、具体的には、結像装置、結像方法および結像システムに関する。
【背景技術】
【０００２】
従来の透過型電子顕微鏡（ＣＴＥＭ）では、収束された電子線によりサンプルの表面に入射するのではなく、平行電子線によりサンプルに入射し、通常、ＣＣＤを探知器として使用し、電子のＣＣＤ全体範囲内における各画素点での信号を一括収集することができ、例えば中国特許ＣＮ１０５５７５７４９Ａには、ＣＣＤによる結像で透過型電子回折パターンを取得する。走査型透過電子顕微鏡（ＳＴＥＭ）構造は、所定のエネルギーにより加速した電子線を収束し、その収束点をサンプルの表面に位置させた後、電子線を制御してサンプルを走査させ、さらに探針によってサンプルを透過する電子線信号を収集し、結像を行う装置であり、ＳＴＥＭは、ＣＴＥＭに比べて、その解像度が高く、他の付属設備を同時に取り付けることができるので、物理、化学、材料科学、生物学などの学科に広く用いられている。現代の走査型透過電子顕微鏡、例えば米国特許ＵＳ４０９９０５５Ａに記載されたシステムでは、電子線が収束システムにより制限されて集束し、走査システムの制御に応じてサンプルを走査し、探知器を用いて明視野又は暗視野のサンプルを通る信号を取り込み、所定の顕微鏡構造設計によって、明視野、暗視野などの使わない信号を同時に取得し、一回の走査で複数の画像を得ることができる。
【０００３】
従来の商業用走査型透過電子顕微鏡の構造は、高度に統合化した複雑なシステムであり、通常、他の付属設備、例えば様々なスペクトル分析装置を多数含み、典型的には電子エネルギー損失分光があり、例えば米国特許ＵＳ４７４３７５６ＡおよびＵＳ８３３４５１２において、通常、各走査点のエネルギー損失が異なる電子を区別することにより、サンプルの成分、化学状態などの情報を測定する。したがって、スペクトル分析を行うことは、重要な意味がある。
【０００４】
従来の商業用走査型透過電子顕微鏡の構造は、通常、一組の環状探針を用いて異なる集光角の散乱電子を検知することにより、コントラストの異なるパターンを形成するものであり、電子線は、結像に用いられない場合に、特性スペクトルの分析、例えば電子エネルギー損失分光分析などに用いられることが多い。ＨＡＡＤＦ、ＡＤＦ、ＡＢＦ、ＢＦなどは、環状探針を用いる公知の技術であり、これらの方法の差異は集光角が異なることにあり、異なる集光角によって、画像は異なるコントラストを有することができ、例えば重元素には、ＨＡＤＤＦ結像が適用される。環状探針を用いて結像する方法の利点は、結像の速度が速いことであるが、特性スペクトルの分析を結像とともに行うことができず、かつ結像過程において、透過型電子顕微鏡の構造では複数種類の電子信号を同時に取得することができない。その後、分割された環状探針が出現し、例えば米国特許ＵＳ２０１６２５４１１８Ａ１に提案される分割探針は、環状探針を複数の領域に分割し、各領域で独立した信号を取得できることを特徴とし、ＤＰＣおよび派生したｉＤＰＣ技術は、通常、当該分割された環状探針を採用し、非対称な信号でサンプル電位または電気量を測定するものである。透過型電子顕微鏡の構造で複数種類の電子信号を同時に得る効果を達成することができるが、最終的には画素化の効果を達成していないので、画素化探知器アレイではない。
【０００５】
画素化した探針アレイは、画素ドットモデルの信号を採集し、当該画素化信号に対して後続の処理、例えばＤＰＣ（微分位相コントラスト）、ｉＤＰＣ（積分微分位相コントラスト）、ＨＡＡＤＦ（高角環状暗視野像）、ＡＤＦ（環状暗視野像）、ＡＢＦ（環状明視野像）、ＢＦ（明視野像）などの結像方法を行うものである。
【０００６】
上記をまとめて、一回の走査で回折パターンおよびスペクトル信号を同時に取得し、かつ回折パターンを取得する過程中に画素化検知モジュールの効果を同時に達成することができることを如何に実現する、という研究は、非常に重要な意味がある。
【０００７】
積層結像方法は、コヒーレンス回折の走査に基づく結像方法であり、位相回復アルゴリズムの結像手段を組み合わせて、従来の透過型走査結像方法に比べて、当該方法は、従来の対物レンズによる結像を必要とせず、位相差などの欠陥を回避することができる。かつ、可視光域、ｘ射線、電子顕微鏡などの異なる波長条件下で、結像の品質および結像の解像度の点でいずれも明らかな優勢があることが既に確証された。積層結像は、画像取得装置を用い、目標画像を再構成するための画像を取得し、過程全体において、主に位相を再構成する。積層結像反復エンジンは、プローブ関数を移動することにより再構成画像を取得するアルゴリズムを採用する。補強型積層結像反復エンジンは、反復過程中において未知の初期プローブ関数を得ることができ、積層結像と共通の技術基礎を含む。よって、積層結像は、回折パターンの取得を基礎として行わなければならず、走査過程中の回折パターンに対して数学処理を行うことにより積層結像の結果を得る。
【０００８】
以上をまとめて、従来の走査型透過電子顕微鏡では、電子収束回折パターンを取得するとともにスペクトル分析を行うことができる装置がなく、当該機能を実現する装置を研究して得ても、従来技術における積層結像方法および派生した積層結像反復エンジン方法や補強型積層結像反復エンジン方法などにとって、如何に当該装置で取得した回折パターンに対して積層結像を行うかは、やはり難題である。
【０００９】
したがって、従来の積層結像方法を如何に改良してそれに様々な走査型透過電子顕微鏡装置に適応させるかは、非常に重要である。
【発明の概要】
【発明が解決しようとする課題】
【００１０】
従来の走査型透過電子顕微鏡では回折パターンおよびスペクトルデータを同時に取得することができないこと、及び従来技術における積層結像などの方法では当該装置に対して結像することができないこと、という従来技術の欠点を克服するために、本発明は、一回の走査で回折パターンおよびスペクトル信号を同時に取得することを実現し、データの取得効率を向上させる結像装置、及び、特殊な結像装置に対して積層結像することができ、結像解像度が高く、画像の品質が良い結像方法を提供する。
【課題を解決するための手段】
【００１１】
上記目的を実現するために、本発明に係る結像装置は、荷電粒子を放出するための荷電粒子源と、
荷電粒子線を制限し収束するための収束システムと、
荷電粒子線のサンプルでの走査を制御するための走査制御システムと、
サンプルと、
荷電粒子を取り込み、荷電粒子の信号強度を検出するための検知モジュールと、
検知モジュールの下方に配置され、荷電粒子の分光学特徴を分析するためのスペクトル分析モジュールと、を含み、
前記検知モジュールは、複数の画素化探知器ユニットを含み、当該検知モジュールには、穴が開設されている。
【００１２】
当該結像装置の結像過程において、まず、荷電粒子源により荷電粒子線を放出し、荷電粒子線は収束システムに進入して収束され、収束後の荷電粒子線はサンプル上を走査し、走査制御システムにより荷電粒子線を制御してサンプルを走査させ、サンプルの所定領域が照射されて走査され、走査過程において、各走査点での荷電粒子線の全部又は一部は、サンプルを透過し、サンプルを透過した荷電粒子線は検知モジュールに進入し、荷電粒子線の一部が検知モジュールにより取り込まれ、荷電粒子線の他の一部が検知モジュールの穴を通ってスペクトル分析モジュールに進入し、検知モジュールにより取り込まれた荷電粒子線は、検知モジュールにより検出されることができ、サンプルの走査範囲における複数の空間位置での信号強度を一括して与え、スペクトル分析モジュールに進入した荷電粒子線は、スペクトル分析モジュールでスペクトルデータ又はエネルギー損失分光などを取得する。穴が開設された検知モジュールを用いて回折パターンを取得し、穴を通る荷電粒子線を用いてスペクトル分析を行い、一回の走査で回折（散乱）パターンおよびスペクトル信号を同時に取得し、データの取得効率を高めることができる。
【００１３】
さらには、前記穴は、円形状、方形状又は環形状を含むが、これらに限定されない。
【００１４】
さらには、前記検知モジュールは、方形状、円形状および島状を含むが、これらに限定されない。
【００１５】
本発明に係る結像方法は、
回折パターンを取得する工程と、
中空積層結像方法により、取得した回折パターンに対して積層結像を行う工程と、を含む。まず、穴が設けられた結像装置により回折パターンを取得し、そして中空積層結像方法を利用して当該回折パターンを処理して積層結像を得る。
【００１６】
さらには、前記回折パターンを取得する工程は、結像装置に基づくものであり、前記結像装置は、荷電粒子源、収束システム、走査制御システム、サンプルおよび検知モジュールを含み、検知モジュールは複数の画素化探知器ユニットを含み、当該検知モジュールには、穴が開設されている。
【００１７】
当該方法は、以下の工程を含む。
【００１８】
Ａ：収束システムによって、荷電粒子源により放出された荷電粒子線をサンプルに収束する。
【００１９】
Ｂ：走査制御システムによって荷電粒子線を制御してサンプルを走査させる。
【００２０】
Ｃ：荷電粒子線がサンプルを透過して検知モジュールに到達し、検知モジュールにおける画素化探知器ユニットは、対応する走査位置における荷電粒子の信号強度を検出し、対応する走査位置の回折パターンを取得する。
【００２１】
Ｄ：荷電粒子線がサンプル（４）を走査し、荷電粒子線プローブの走査位置の間に、一部の重合が存在する。
【００２２】
走査ビームスポットに対応する走査領域Ｒｉと他の走査領域Ｒｊには、非空集合Ｒｏｖｅｒｌａｐが存在し、
【数１】
2020537123.tif 000003
であり、Ｎは荷電粒子線によるサンプルの走査領域での総数である。即ち、回折パターンを生成する走査領域のそれぞれが回折パターンを生成する１つ又は複数の他の走査領域と重なることが求められ、最後に得られる回折パターンに欠損がないことをさらに確保する。
【００２３】
さらには、中空積層結像方法により、取得した回折パターンに対して積層結像を行う工程は、具体的には、以下の工程を含む。
【００２４】
ａ：Ｐ（ｒ）をプローブ関数とし、Ｏ（ｒ）を物体の複素振幅分布関数とし、複数回の反復計算によってＯ（ｒ）を再構成し、最終的な反復計算で再構成するＯ（ｒ）を、物体の最終的な複素振幅分布関数とし、積層結像は、物体の最終的な複素振幅分布関数に基づいて画像を再構成する。
【００２５】
ｂ：Ψｎ，ｍを物体を透過する出射波関数とし、Ψｎ，ｍを、プローブ関数Ｐ（ｒ）と物体の複素振幅分布関数Ｏ（ｒ）との乗積に定義すると、
【数２】
2020537123.tif 000004
【００２６】
式（１）において、ｎはＯ（ｒ）のｎ回目の反復を表し、ｍは荷電粒子線のサンプルでのｍ番目の走査位置を表し、Ｒｍは、荷電粒子線のサンプルでの、１番目の走査位置に対するｍ番目の走査位置の相対的座標ベクトルを表し、ｒは空間座標である。
【００２７】
ｃ：物体を透過する出射波関数Ψｎ，ｍに対してフーリエ変換を行うことにより、当該出射波関数Ψｎ，ｍのファーフィールドでの振幅および位相分布を得る。
【数３】
2020537123.tif 000005
【００２８】
式（２）において、
【数４】
2020537123.tif 000006
は出射波関数Ψｎ，ｍのファーフィールドでの振幅を表し、Φｎ，ｍは出射波関数Ψｎ，ｍのファーフィールドでの位相を表す。
【００２９】
ｄ：実験装置によってサンプルを透過する出射波のファーフィールド光強度を採集し、Ｉｍとし、Ｉｍは、荷電粒子線がサンプルにおけるｍ番目の走査位置を走査した後にサンプルを透過する出射波のファーフィールド光強度を表す。
【００３０】
ｅ：制約関数Ｍを設定する。
【００３１】
ｆ：出射波関数Ψｎ，ｍのファーフィールドでの振幅
【数５】
2020537123.tif 000007
に代えて、
【数６】
2020537123.tif 000008
を制約関数Ｍに代入し、式（３）を得る。
【数７】
2020537123.tif 000009
【００３２】
ｇ：Ψｎ，ｍ，ｎｅｗ’（ｒ）に対して逆フーリエ変換を行い、式（４）で示される新たな出射波関数Ψｎ，ｍ，ｎｅｗ（ｒ）を得る。
【数８】
2020537123.tif 000010
【００３３】
ｈ：新たな出射波関数Ψｎ，ｍ，ｎｅｗ（ｒ）により計算し、式（５）で示される新たな物体の複素振幅分布関数を得る。
【数９】
2020537123.tif 000011
【００３４】
式（５）において、αおよびβは調節可能なパラメータであり、αは、分母が０ではないことを確保するために用いられ、βは、帰還の強度を制御するために用いられる。
【００３５】
ｉ：新たな物体の複素振幅分布関数Ｏｎ＋１（ｒ＋Ｒｍ）を式（１）に代入して再度の反復計算を開始する。
【００３６】
ｊ：最終的な反復計算で得られる物体の複素振幅分布関数Ｏｚ（ｒ＋Ｒｍ）を画像再構成の基礎とし、再構成画像はＯｎ＋１（ｒ＋Ｒｍ）の振幅および位相に対して描画して得られるものである。
【００３７】
式（１）〜式（５）により計算して新たな物体の複素振幅分布関数が得られた後、新たな物体の複素振幅分布関数Ｏｎ＋１（ｒ＋Ｒｍ）を採用して改めて式（１）に代入して再度の反復計算を開始し、最終的な反復計算で得られた物体の複素振幅分布関数Ｏｚ（ｒ＋Ｒｍ）を画像再構成の基礎とし、Ｏｎ＋１（ｒ＋Ｒｍ）の振幅および位相を描画して再構成画像を得る。
【００３８】
さらには、前記制約関数Ｍは、検知モジュールの構造に関連する。離散化の条件で、マトリクスで当該制約関数Ｍを表し、当該マトリクスにおいて、各要素値は従属変数に対応し、従属変数の関数値は１又は０であり（１又は０に限定されない）、独立変数は、対応する二次元位置座標である。
【００３９】
さらには、Ｏ（ｒ）の１回目の反復において、物体の複素振幅分布関数Ｏ０（ｒ＋Ｒｍ）をランダム分布関数にする。
【００４０】
本発明に係る結像システムは、結像装置および結像方法を含む。
【発明の効果】
【００４１】
本発明は従来技術に比べて、以下の利点を有する。
【００４２】
本発明に係る結像装置は、穴が開設された検知モジュールおよびスペクトル分析モジュールを含み、穴が開設された検知モジュールにより回折パターンを取得し、スペクトル分析モジュールは、穴を通った荷電粒子線に対してスペクトル分析を行い、一回の走査で回折パターンおよびスペクトル信号を同時に取得することを実現する。
【００４３】
本発明に係る結像装置では、検知モジュールは複数の画素化探知器ユニットを配列して組み合わせたものであり、画素化探知器ユニットは、特定の位置での荷電粒子の信号強度を取得することができ、さらに回折パターンを取得し、取得した回折パターンの精度が高い。
【００４４】
本発明に係る結像方法は、中空積層結像方法に基づいたものであり、穴が開設された検知モジュールで得られた回折パターンを結像させることを実現でき、結像の効果が良い。
【００４５】
本発明に係る結像方法は、穴が開設された結像装置に適用し、穴が開設された検知モジュールの穴開け形状は円形状、方形状、環形状などに限定されるものではなく、結像効果は、従来の結像方法を適用した従来の結像装置により得られる結像に比べて、画像解像度が高く、画像の品質が高い。
【００４６】
本発明に係る結像方法は、穴が開設された結像装置に適用し、Ｍ関数を採用して当該穴が開設された検知モジュールに対応させ、全ての構造および形状の検知モジュールで得られる回折パターンを結像させることができる。
【００４７】
本発明は、最終的な結像結果として位相コントラスト結像を採用し、位相コントラストは、軽元素結像の分野においてより優勢を有し、結像の効果がより良い。
【図面の簡単な説明】
【００４８】
【図１】本発明に係る結像装置の構造模式図である。
【図２】本発明に係る検知モジュールの構造模式図である。
【図３】本発明に係る結像方法−中空積層結像アルゴリズムの工程のフローチャットである。
【図４】実験の原始写真の情報である。
【図５】実験の原始写真における再構成に用いられる領域の情報である。
【図６】実施例９における初期関数Ｐ（ｒ）の振幅画像である。
【図７】実施例９における初期関数Ｐ（ｒ）の位相画像である。
【図８】実施例９におけるＭ関数画像である。
【図９】実施例９における再構成画像の振幅画像である。
【図１０】実施例９における再構成画像の位相画像である。
【図１１】実施例１０における再構成画像の振幅画像である。
【図１２】実施例１０における再構成画像の位相画像である。
【図１３】実施例１１におけるＭ関数画像である。
【図１４】実施例１１における再構成画像の振幅画像である。
【図１５】実施例１１における再構成画像の位相画像である。
【図１６】実施例１２におけるＭ関数画像である。
【図１７】実施例１２における再構成画像の振幅画像である。
【図１８】実施例１２における再構成画像の位相画像である。
【発明を実施するための形態】
【００４９】
以下に、図面および実施例を組み合わせて本発明をさらに説明する。
【００５０】
実施例１：
本実施例に係る結像装置は、図１を参照し、荷電粒子源１、収束システム２、走査制御システム３、サンプル４、検知モジュール５および検知モジュール５の下方に配置されているスペクトル分析モジュール６を含む。
【００５１】
荷電粒子源１は、荷電粒子を放出するために用いられる。収束システム２は、荷電粒子線を制限し収束するために用いられる。走査制御システム３は、荷電粒子線のサンプルでの走査を制御するために用いられ、走査方式として、荷電粒子線を移動させてもよいし、サンプルを移動させてもよいが、これらに限定されていない。検知モジュール５は、荷電粒子を取り込んで荷電粒子の信号強度を検出するために用いられる。検知モジュール５の下方に配置されているスペクトル分析モジュール６は、荷電粒子の分光学特徴を分析するために用いられる。図２を参照し、本実施例における検知モジュール５は、複数の画素化探知器ユニット７を含み、検知モジュール５は複数の画素化探知器ユニット７を配列して組み合わせてなるものであり、複数の画素化探知器ユニット７の配列方式は、方形状、円形状又は島状分布を含むがこれらに限定されるものではなく、各画素化探知器ユニット７は、いずれも特定の空間位置での荷電粒子の信号強度を与えることができる。検知モジュール５には、穴８が開設され、穴８は、検知モジュール５に到達した荷電粒子線を透過させて下方のスペクトル分析モジュール６に進入させるために用いられ、荷電粒子線が検知モジュール５に照射する時、荷電粒子線の一部が穴８の位置を通過し、下方のスペクトル分析モジュール６に到達することを確保しなければならない。
【００５２】
当該結像装置の結像過程において、まず、荷電粒子源１により荷電粒子線を放出し、荷電粒子線は収束システム２に進入して収束され、収束後の荷電粒子線はサンプル４を走査し、走査制御システム３により荷電粒子線を制御してサンプル４を走査させ、サンプル４の所定領域が照射されて走査され、走査過程において、各走査点での荷電粒子線の全部又は一部は、サンプル４を透過し、サンプル４を透過した荷電粒子線は検知モジュール５に進入し、荷電粒子線の一部が検知モジュール５により取り込まれ、荷電粒子線の他の一部が検知モジュール５の穴８を通ってスペクトル分析モジュール６に進入し、検知モジュール５により取り込まれた荷電粒子線は、検知モジュールにより検出されることができ、サンプルの走査範囲における複数の空間位置での信号強度を一括して与え、スペクトル分析モジュール６に進入した荷電粒子線は、スペクトル分析モジュールでスペクトルデータ又はエネルギー損失分光などを取得する。
【００５３】
本実施例に係る装置は、穴が開設された検知モジュールを用いて回折パターンを取得し、穴を透過した荷電粒子線を用いてスペクトル分析を行い、一回の走査で回折（散乱）パターンおよびスペクトル信号を同時に取得し、データの取得効率を高めることができる。
【００５４】
実施例２：
本実施例に係る結像装置は、実施例１に基づき、図２を参照し、穴８は、円形状、方形状又は環形状を含むがこれらに限定されていない。収束システム２および走査制御システム３は、荷電粒子線がサンプルを照射する領域の範囲、領域範囲の大きさおよび相対位置を正確に確認することができ、荷電粒子線がサンプルを走査する度に、検知モジュール５で得られた回折パターンおよびスペクトル分析モジュール６で得られたスペクトルデータは、荷電粒子線がサンプルのどの位置および領域を透過して生じるものであるかを特定することができ、荷電粒子線の照射範囲をＲに定義し（原点に対する照射領域における全ての点の座標の集合）、得られる回折パターンをＤに定義し（Ｄは各探知器ユニットで得られる信号の大きさ（光強度）の集合である）、得られた分光学データをＳに定義し（スペクトル分析モジュール６によって取得される）、結像装置は、走査過程において、走査するごとに（Ｒ、Ｄ、Ｓ）のデータ情報を同時に取得する。
【００５５】
実施例３：
本実施例に係る結像方法は、回折パターンを取得する方法および中空積層結像方法により取得した回折パターンに対して積層結像を行う方法を含む。
【００５６】
そのうち、回折パターンを取得する方法は、結像装置に基づくものであり、当該結像装置は、図１および図２を参照し、荷電粒子源１、収束システム２、走査制御システム３、サンプル４および検知モジュール５を含み、検知モジュール５は、複数の画素化探知器ユニット７を含んで構成され、当該検知モジュール５には、穴８が開設されている。そのうち、荷電粒子源１は、荷電粒子を放出するために用いられる。収束システム２は、荷電粒子線を制限し収束するために用いられる。走査制御システム３は、荷電粒子線のサンプルでの走査を制御するために用いられ、走査方式は、荷電粒子線を移動させる方式やサンプルを移動させる方式を含むがこれらに限定されていない。検知モジュール５は、荷電粒子を取り込んで荷電粒子の信号強度を検出するために用いられる。本実施例における検知モジュール５は、複数の画素化探知器ユニット７を含み、検知モジュール５は複数の画素化探知器ユニット７を配列して組み合わせてなるものであり、配列方式は、方形状、円形状又は島状分布を含むがこれらに限定されるものではなく、各画素化探知器ユニット７は、いずれも特定の空間位置での荷電粒子の信号強度を与えることができる。
【００５７】
当該方法は次の工程を含む。収束システム２によって、荷電粒子源１により放出された荷電粒子線をサンプル４に収束する。走査制御システム３によって荷電粒子線を制御してサンプル４を走査させる。荷電粒子線は、サンプル４を透過して検知モジュール５に到達し、検知モジュール５における画素化探知器ユニット７は、対応する走査位置における荷電粒子の信号強度を検出し、穴８が開設された検知モジュール５は、対応する走査位置の回折パターンを取得する。走査制御システム３によって荷電粒子線を制御してサンプル４を走査させる。収束システム２および走査制御システム３は、荷電粒子線がサンプルに照射する領域の範囲、領域範囲の大きさおよび相対位置を正確に確認することができ、荷電粒子線がサンプルを走査する度に、検知モジュール５で得られる回折パターンは、荷電粒子線がサンプルのどの位置および領域を透過して生じるものであるかを特定することができ、荷電粒子線の照射範囲をＲに定義し（原点に対する照射領域における全ての点の座標の集合）、得られる回折パターンをＤに定義し（Ｄは各探知器ユニットで得られる信号の大きさ（光強度）の集合である）、結像装置は、走査過程において、走査するごとに（Ｒ、Ｄ）のデータ情報を同時に取得する。回折パターンは、検知モジュール５における探知器ユニットで得られた荷電粒子の強度に基づく。
【００５８】
荷電粒子線がサンプル４を走査して一組のサンプルの回折パターン情報を取得し、当該情報は、走査位置情報および回折パターン情報、例えば｛（Ｒ１、Ｄ１）、（Ｒ２、Ｄ２）．．．（ＲＮ、ＤＮ）｝を含み、荷電粒子線によるサンプル４での走査領域Ｒｉと他の１つ又は複数の走査領域Ｒｊには非空集合Ｒｏｖｅｒｌａｐが存在し、
【数１０】
2020537123.tif 000012
であり、Ｎは荷電粒子線がサンプル（４）を走査する走査領域の総数であり、即ち、回折パターンを生成する走査領域のそれぞれは、回折パターンを生成する１つ又は複数の他の走査領域と重なることが求められる。
【００５９】
実施例４：
本実施例に係る結像方法は、実施例３に基づき、中空積層結像アルゴリズム（Ｈｏｌｌｏｗ−ＰｔｙＣｈｏｇｒＡｐｈｙ）に基づいて検知モジュール５で得られた回折パターンを処理して画像を取得するために用いられ、当該結像方法は、荷電粒子の回折パターンによって位相コントラストのサンプル画像を取得することができ、回折パターンの初期データおよび出射波の初期の推測の反復過程に基づき、出射波の振幅および位相情報を再構成し、
再構成する位相、振幅と結像との関係は以下のとおりである。
【００６０】
再構成の最終的な結果は、Ｏ（ｒ）関数であり、その独立変数は位置座標であり、値は、荷電粒子線が特定の位置座標でサンプルを透過する振幅および位相の変化量であると理解することができ、荷電粒子線の振幅および位相の変化量の相対値を異なるグレースケール値に換算することが約束され、その空間座標が異なる場合に振幅および位相変化の大きさが異なることによりグレースケール画像を取得することを、画像の再構成と理解してもよい。Ｏ（ｒ）における振幅変化の大きさによる結像は幅コントラスト像と呼ばれ、位相変化の大きさによる結像は位相コントラスト像と呼ばれる。
【００６１】
当該方法は図３を参照し、以下の工程を含む。
【００６２】
Ｐ（ｒ）をプローブ関数（Ｐｒｏｂｅｆｕｎｃｔｉｏｎ）とし、プローブ関数は、荷電粒子線がサンプルに到達する前の入射波関数であり、ｒは空間座標である。
【００６３】
Ｏ（ｒ）を物体の複素振幅分布関数とし、物体の複素振幅分布関数は、物体の構造および特性を十分に反映することができ、その値は、サンプルの特徴位置の、通過した荷電粒子線への作用（振幅および位相の変化として具現される）を表し、Ｏ（ｒ）は目標関数であり、複素関数であり、本実施例の方法は、複数回の反復計算によりＯ（ｒ）を再構成し、最後の反復計算によるＯ（ｒ）を物体の複素振幅分布関数とし、ｒは空間座標であり、Ｏｎ（ｒ）は、ｎ回目の反復における目標関数を表し、反復過程の全体において、Ｏｎ（ｒ）の値もそれに伴って変化し、最終的な再構成画像は、Ｏｎ（ｒ）の振幅および位相を描画して得られるものである。
【００６４】
Ψｎ，ｍを物体を透過する出射波関数とし、物体を透過する出射波関数Ψｎ，ｍをプローブ関数Ｐ（ｒ）と物体の複素振幅分布関数Ｏ（ｒ）との乗積に定義すると、
【数１１】
2020537123.tif 000013
【００６５】
式（１）において、ｎはＯ（ｒ）のｎ回目の反復を表し、ｍは荷電粒子線のサンプルでのｍ番目の走査位置を表し、Ｒｍは、荷電粒子線のサンプルでの、１番目の走査位置に対するｍ番目の走査位置の相対的座標ベクトルを表す。
【００６６】
物体を透過する出射波関数Ψｎ，ｍに対してフーリエ変換を行うことにより、当該出射波関数Ψｎ，ｍのファーフィールドでの振幅および位相分布を得る。
【数１２】
2020537123.tif 000014
【００６７】
式（２）において、
【数１３】
2020537123.tif 000015
は、ｎ回目の反復における出射波関数Ψｎ，ｍのファーフィールドでの振幅を表し、Φｎ，ｍは、ｎ回目の反復における出射波関数Ψｎ，ｍのファーフィールドでの位相を表し、ｉは虚数単位である。
【００６８】
出射波関数は、通常、複数形式で表し、荷電粒子線の性質を波関数で説明することができ、実部および虚部を含み、波関数の位相および振幅で表すと、
【数１４】
2020537123.tif 000016
は、空間におけるある位置（ｒ）での荷電粒子線の強度に対応し、実際に測定することができ、実験装置により荷電粒子線のサンプル上のｍ番目の走査位置でのファーフィールド光強度を採集し、Ｉｍとする。
【００６９】
制約関数Ｍを設定し、関数Ｍは、上記結像システムにおける検知モジュールの形状に関連し、制約関数Ｍは実数関数であり、回折パターン（Ｄ）の空間の大きさと同様に、空間座標の関数であり、その値は、回折パターン（Ｄ）における異なる空間座標の信号の反復過程での重みを制御し、実験中に具現すると、Ｍ関数の値は、検知モジュールの具体的な構造および形状に関連する。
【００７０】
出射波関数Ψｎ，ｍのファーフィールドでの振幅
【数１５】
2020537123.tif 000017
に代えて、
【数１６】
2020537123.tif 000018
を制約関数Ｍに代入し、式（３）を得る。
【数１７】
2020537123.tif 000019
【００７１】
Ψｎ，ｍ，ｎｅｗ’（ｒ）に対して逆フーリエ変換を行い、式（４）で示される新たな出射波関数Ψｎ，ｍ，ｎｅｗ（ｒ）を得る。
【数１８】
2020537123.tif 000020
【００７２】
新たな出射波関数Ψｎ，ｍ，ｎｅｗ（ｒ）により計算し、式（５）で示される新たな物体の複素振幅分布関数を得る。
【数１９】
2020537123.tif 000021
【００７３】
式（５）において、αおよびβは調節可能なパラメータであり、αは、分母が０ではないことを確保するために用いられ、βは、帰還の強度を制御するために用いられる。
【００７４】
式（１）〜式（５）により計算して新たな物体の複素振幅分布関数を得ると、新たな物体の複素振幅分布関数Ｏｎ＋１（ｒ＋Ｒｍ）を改めて式（１）に代入して再度の反復計算を開始し、最終的な反復計算で得られた物体の複素振幅分布関数Ｏｚ（ｒ＋Ｒｍ）を画像再構成の基礎とし、Ｏｎ＋１（ｒ＋Ｒｍ）の振幅および位相を描画して再構成画像を得る。
【００７５】
実施例５：
本実施例に係る結像方法は、実施例４に基づき、制約関数Ｍは検知モジュール５の構造に関連する。
【００７６】
離散化の条件下で、マトリクス関数で当該制約関数Ｍを表し、当該マトリクス関数における各要素値は従属変数に対応し、従属変数の関数値は１又は０であり、独立変数は対応する二次元位置座標である。
【００７７】
また、シミュレーションプログラムの計算において、直接にマトリクス関数で当該制約関数Ｍを表し、マトリクス関数は検知モジュール５の構造に関連し、マトリクス関数の位置座標は検知モジュール５の位置座標に対応し、マトリクス関数における従属変数の関数値が１である位置座標は、検知モジュール５における穴８の位置座標に対応する。
【００７８】
シミュレーションプログラムにおいて、マトリクス関数を１つの画像に対応させ、マトリクス関数における従属変数の関数値が１である位置座標を、画像中において全て白色で表示し、マトリクス関数における従属変数の関数値が０である位置座標を、画像中において全て黒色で表示すると、マトリクス関数Ｍに対応する画像中において、白色で表示する領域は検知モジュール５における穴の位置に対応する。
【００７９】
実施例６：
本実施例に係る結像方法は、実施例５に基づき、Ｏ（ｒ）の１回目の反復において、物体の初期複素振幅分布関数Ｏ０（ｒ＋Ｒｍ）はランダムに分布する。物体の初期複素振幅分布関数Ｏ０（ｒ＋Ｒｍ）とプローブ関数Ｐ（ｒ）との乗積を、物体を透過する出射波関数Ψｎ，ｍとすると、
【数２０】
2020537123.tif 000022
である。
【００８０】
実施例７：
本実施例に係る結像方法は、実施例５に基づき、Ｏ（ｒ）の１回目の反復において、物体の初期複素振幅分布関数Ｏ０（ｒ＋Ｒｍ）は、他の先験的条件がある場合、他の方式によって、サンプルの低解像度画像を取得することができ、この低解像度画像を基礎として高解像度画像を取得することが望まれ、このとき、低解像度画像をＯ（ｒ）の初期複素振幅分布関数Ｏ０（ｒ＋Ｒｍ）として直接使用する。
【００８１】
実施例８：
本実施例に係る結像方法は、実施例５に基づき、荷電粒子線はサンプルに収束され、サンプル走査を行った後、大部分の荷電粒子線は検知モジュールに到達するが、やはり検知モジュールの穴８を通る荷電粒子線が存在し、検知モジュールの穴８を貫通した荷電粒子線の信号強度が検知モジュールにより検知されないので、対応する回折パターンが欠損してしまう。補償として、検知モジュール５に開設される穴８の弧度、即ち回折パターンを取得する集光角の範囲は５ｍｒａｄ〜２２ｍｒａｄであり、この場合に、荷電粒子線はサンプルを絶えず走査し、各走査領域は、いずれも他の複数の領域と重なり、次の走査位置で欠損した回折パターンを改めて取得する。結像の品質を低下させることなく、穴が開設された検知モジュールにより完全な再構成結果を得るためには、穴８の中央が十分に小さく、エッジが十分に大きいという条件を満たさなければならない。
【００８２】
実施例９：
本実施例に係る結像方法は、実施例８に基づき、図４を、検知モジュール５で取得した回折パターンとしてアルゴリズムの実行可能性を検証し、当該回折パターンには、元の位相情報および振幅情報を含み、写真解像度は２５６＊２５６であり、中空積層結像方法によりその一部の領域（図５を参照）を再構成し、再構成領域における走査点の数は８１００個（９０＊９０）であり、その回折パターンの集光角は２２ｍｒａｄである。
【００８３】
コンピュータにより模擬される探知器が各走査点で得られた回折パターン解像度は５１２＊５１２であり、調節可能なパラメータαおよびβは、反復回数に対する目標関数の変化幅を制御し、自由に設定可能なパラメータであり、典型値は０．１〜０．０１であり、当該実施例において、パラメータはα＝０．０１、β＝０．０１である。
【００８４】
Ｐ（ｒ）関数は、平面での位置を変数とする二次元複素関数であり、その形式は唯一ではなく、積層結像アルゴリズムにおいて、Ｐ（ｒ）関数は反復において変化せず、当該実施例は、Ｐ（ｒ）が変化しない条件下でコンピュータで模擬するものである。当該実施例で用いるＰ（ｒ）関数画像は図６および図７に示すように、図６はＰ（ｒ）関数の振幅画像であり、図７はＰ（ｒ）関数の位相画像である。
【００８５】
実際の実験において、Ｍ関数は検知モジュールの構造および形状により決められ、コンピュータで模擬する方式により、異なるＭ関数を用いて再構成を行うことができ、Ｍ関数は二次元（平面）実関数であり、通常、その平面での大きさ（画素解像度）は、実際の検知モジュールの物理的（画素）解像度と同じであり、この実施例において、Ｍ関数の画素解像度は５１２＊５１２であり、当該実施例中におけるＭ関数の特徴は、中央半径９１個の画素の円形領域内の値が１であり、白色で示し、他の領域の値が０であり、黒色で示し、対応する検知モジュールに開設される穴が円形であり、回折パターンを取得する集光角が２２ｍｒａｄであり、即ち荷電粒子線の光軸に沿う発散角が２２ｍｒａｄの領域の値が１であり、他の領域が０であり、所用のＭ関数（マトリクス関数）に対応する画像を図８に示す。
【００８６】
当該実施例において、反復を実行する循環数は１０回の循環である。コンピュータにより模擬してシミュレートして得られた最終的な再構成画像の振幅画像を図９に示し、再構成の位相画像を図１０に示す。
【００８７】
実施例１０：
本実施例は実施例９に基づき、反復循環を実行する回数は５０回の循環である以外、他のパラメータ条件はいずれも実施例９に一致し、実験して画像を再構成するために用いられ、再構成画像の結果は図１１、図１２に示すように、図１１は再構成画像の振幅画像であり、図１２は再構成した位相画像である。
【００８８】
実施例１１：
本実施例は実施例９に基づき、反復循環を実行する回数は５０回の循環であり、Ｍ関数の特徴は、中央の内半径２０画素、外半径９０画素の領域の値が１であり、白色で示し、他の領域の値が０であり、黒色で示し、図１３を参照し、対応する検知モジュールに開設された穴は環状孔であり、円環部分は、発散角が５ｍｒａｄ〜２２ｍｒａｄであることに相当する。
【００８９】
他のパラメータ条件はいずれも実施例６に一致し、実験して画像を再構成するために用いられ、再構成画像の結果は図１４、図１５に示すように、図１４は再構成画像の振幅画像であり、図１５は再構成した位相画像である。
【００９０】
実施例１２：
本実施例は実施例９に基づき、反復循環を実行する回数は５０回の循環であり、Ｍ関数の特徴は、中央辺長が１００ｐｉｘｅｌの正方形領域が１であり、白色で示し、他の領域が０であり、黒色で示し、図１６を参照し、対応する検知モジュールに開設された穴が方形孔である。他のパラメータ条件はいずれも実施例６に一致し、実験して画像を再構成するために用いられ、再構成画像の結果は図１７、図１８に示すように、図１７は再構成画像の振幅画像であり、図１８は再構成した位相画像である。
【００９１】
実施例１３：
結像装置および結像方法を含む結像システムであり、結像方法は中空積層結像アルゴリズムを採用し、結像装置は穴が設けられた結像装置を採用し、まず、穴が設けられた結像装置で回折パターンを取得し、さらに中空積層結像方法を利用して当該回折パターンを処理して積層結像を取得する。
【００９２】
以上の記述は、本発明の好適な実施形態に過ぎず、本分野の技術者にとって、本発明の原理を逸脱しない前提で、いくつかの改良および修飾を行うことができ、これらの改良および修飾も本発明の保護範囲と見なすべきであることを注意されたい。
【符号の説明】
【００９３】
１：荷電粒子源
２：収束システム
３：走査制御システム
４：サンプル
５：検知モジュール
６：スペクトル分析モジュール
７：画素化探知器ユニット
８：穴
【図１】
2020537123.tif 000032
【図２】
2020537123.tif 000033
【図３】
2020537123.tif 000034
【図４】
2020537123.tif 000035
【図５】
2020537123.tif 000036
【図６】
2020537123.tif 000037
【図７】
2020537123.tif 000038
【図８】
2020537123.tif 000039
【図９】
2020537123.tif 000040
【図１０】
2020537123.tif 000041
【図１１】
2020537123.tif 000042
【図１２】
2020537123.tif 000043
【図１３】
2020537123.tif 000044
【図１４】
2020537123.tif 000045
【図１５】
2020537123.tif 000046
【図１６】
2020537123.tif 000047
【図１７】
2020537123.tif 000048
【図１８】
2020537123.tif 000049
【国際調査報告】
2020537123.tif 000050
2020537123.tif 000051
2020537123.tif 000052
2020537123.tif 000053