(57)【要約】【課題】走査型プローブ顕微鏡が表面に追従する能力を高め、表面トポグラフィの描画精度を向上させること。【解決手段】微小電気機械システム（ＭＥＭＳ）アクチュエータ（１１０）によって垂直方向に制御されるプローブ先端の位置を推定する推定器（１８０）を設けた。【選択図】図１ｂ
2008014943.tif 000002
【特許請求の範囲】
【請求項１】
走査型プローブ顕微鏡システムであって、
プローブ先端（１３０）に機械的に結合される可撓性カンチレバーを動かすように動作することができる微小電気機械システムアクチュエータ（１１０）と、
前記可撓性カンチレバーに光学的に結合される光検出器（１６０）と、
前記微小電気機械システムアクチュエータに電気的に接続され、且つ前記光検出器に電気的に接続されるコントローラ（１７０）と、
推定器フィードバック利得を有する推定器（１８０）と
を備え、前記推定器は前記光検出器に電気的に接続され、且つ前記コントローラに電気的に接続され、該推定器は、前記光検出器から誤差信号を受信し、該誤差信号を用いて、前記コントローラに出力するための状態推定値を生成するように動作することができ、
前記コントローラは、前記状態推定値に応答して、前記微小電気機械システムアクチュエータに制御信号を与え、前記可撓性カンチレバーの撓みを概ね一定の値に保持するように動作することができる、走査型プローブ顕微鏡システム。
【請求項２】
前記微小電気機械システムアクチュエータ（２１０）は、該微小電気機械システムアクチュエータのための絶対位置情報信号を与えるように動作することができる位置センサを備え、
前記コントローラは制御入力信号を生成するように動作することができ、
前記推定器（２８０）は、前記絶対位置情報信号及び前記制御入力信号を用いて、前記状態推定値を生成するようにさらに動作することができる、請求項１に記載の走査型プローブ顕微鏡システム。
【請求項３】
前記推定器は予測推定器である、請求項１又は２に記載の走査型プローブ顕微鏡システム。
【請求項４】
前記推定器は電流推定器である、請求項１又は２に記載の走査型プローブ顕微鏡システム。
【請求項５】
前記推定器は離散時間型推定器である、請求項１又は２に記載の走査型プローブ顕微鏡システム。
【請求項６】
前記推定器は広帯域幅推定器である、請求項１又は２に記載の走査型プローブ顕微鏡システム。
【請求項７】
前記推定器フィードバック利得は、極配置の方法によって選択される、請求項１又は２に記載の走査型プローブ顕微鏡システム。
【請求項８】
前記推定器フィードバック利得は、Ｈ∞設計の方法によって選択される、請求項１又は２に記載の走査型プローブ顕微鏡システム。
【請求項９】
前記推定器フィードバック利得は、混在Ｈ２／Ｈ∞設計の方法によって選択される、請求項１又は２に記載の走査型プローブ顕微鏡システム。
【請求項１０】
前記推定器フィードバック利得は、カルマンフィルタを用いることによって選択される、請求項１又は２に記載の走査型プローブ顕微鏡システム。
【請求項１１】
前記カルマンフィルタは定常状態カルマンフィルタである、請求項１０に記載の走査型プローブ顕微鏡システム。
【請求項１２】
前記推定器は表面推定値を生成するように動作することができる、請求項１又は２に記載の走査型プローブ顕微鏡システム。
【請求項１３】
前記制御入力から前記表面推定値への直接的なフィードスルーは存在しない、請求項１２に記載の走査型プローブ顕微鏡システム。
【発明の詳細な説明】
【背景技術】
【０００１】
［背景］
走査型プローブ顕微鏡は、局所的なプローブが存在し、そのプローブが表面と相互作用するということが特徴的である。走査型プローブ顕微鏡の例には、走査型トンネル顕微鏡、磁力顕微鏡、静電力顕微鏡、走査力顕微鏡及び原子間力顕微鏡がある。原子間力顕微鏡（ＡＦＭ）は典型的には、可撓性カンチレバーの端部に取り付けられる細いプローブ先端を、結像されるべき表面に接触させる（apply）ことによって表面を結像する。プローブ先端が表面上を動くのに応じて、典型的にはカンチレバー上のレーザスポットを追跡することによって、プローブ先端の高さの変化が検出される。プローブ先端の高さの変化に起因して、カンチレバーに撓みが生じる結果として、光検出器上の反射レーザスポットが動く。一定のフォースモードで利用する場合、カンチレバーを制御して、カンチレバーの撓みを概ね一定に保持するアクチュエータを動かすことによって、光検出器上の反射レーザスポットの動きを最小限に抑えるために、サーボシステムが用いられる。カンチレバーの撓みはプローブ先端と表面との相互作用によるので、撓みを概ね一定に保持することは、力を概ね一定に保持することと同じである。典型的には、従来技術では、カンチレバーは圧電アクチュエータによって制御される。
【０００２】
光検出器は、カンチレバーの撓みにだけ応答し、先端の絶対位置には応答しないので、サーボシステムのための典型的な制御ループは、出力誤差のみのループである。それゆえ、サーボシステムのコントローラからの出力信号が、表面トポグラフィの推定値として用いられる。これにより典型的には、位置推定値の帯域幅及び精度は、制御ループそのものの帯域幅及び精度に制限される。
【０００３】
【特許文献１】米国特許第５，９８６，３８１号
【発明の開示】
【発明が解決しようとする課題および課題を解決するための手段】
【０００４】
［発明の概要］
【０００５】
本発明によれば、微小電気機械システム（ＭＥＭＳ）アクチュエータによって垂直方向に制御されるプローブ先端の正確な位置推定値を与えるために、走査型プローブ顕微鏡のコントローラ部分において推定器が用いられる。
【０００６】
正確な位置推定値は典型的には、走査型プローブ顕微鏡が表面に追従する能力を高め、典型的には、表面トポグラフィの改善された測定値を与える。
【発明を実施するための最良の形態】
【０００７】
［詳細な説明］
図１aは、本発明による、対応するＡＦＭの物理的な装置の概略図であり、光検出器１６０、ＭＥＭＳアクチュエータ１１０、レーザ源１９０、カンチレバー／プローブ先端１３０及び表面１５０が示される。
【０００８】
図１ｂは、本発明による、垂直制御のためのＭＥＭＳアクチュエータ１１０及び推定器１８０を用いるＡＦＭのブロック図を示す。カンチレバー／プローブ先端１３０が表面１５０上を動くのに応じて、光検出器１６０が、カンチレバー／プローブ先端１３０の撓みを感知する。光検出器１６０への入力は、表面１５０上の表面高（ｒ）と、カンチレバー／プローブ先端１３０の位置（ｐ）との間の相互作用であり、その相互作用の結果として、光検出器１６０上の反射レーザスポットが動く。図示される制御システムは、１入力１出力システム、すなわちＳＩＳＯである。
【０００９】
圧電アクチュエータの代わりにＭＥＭＳアクチュエータ１１０を用いることによって、典型的には、アクチュエータの応答を速くできるようになり、圧電アクチュエータに典型的に関連するヒステリシス効果が避けられる。それぞれのアクチュエータの基本共振周波数を比較すると、圧電アクチュエータが典型的には、約２００Ｈｚ〜約５００Ｈｚの範囲の基本共振周波数を有するのに対して、ＭＥＭＳアクチュエータ１１０の場合には典型的な選択によって、約３ｋＨｚ〜約５０ｋＨｚ又はそれ以上の範囲の基本共振周波数を有することがわかる。さらに、ＭＥＭＳアクチュエータ１１０の場合には典型的な選択によって、基本共振周波数よりもはるかに高い閉ループ周波数において動作することができる。本発明によるＭＥＭＳアクチュエータ１１０の例は、たとえば、参照により本明細書に援用される、米国特許第５，９８６，３８１号において説明される。
【００１０】
図１ｂは、光検出器１６０だけを用いる原子間力顕微鏡のためのＭＥＭＳアクチュエータ１１０を示す。光検出器１６０は出力誤差フィードバックだけを与える。生成される出力誤差は、用いられるシステムモデルに依存するが、表面１５０が表面特徴（高さの変化）を有するときには必ず誤差が生じる。ＭＥＭＳアクチュエータ１１０は典型的には、典型的な圧電アクチュエータよりも正確なシステムモデルを可能にする。この結果として、推定器の性能が高められる。推定器として、予測推定器又は電流推定器を用いることができる。推定器フィードバック利得Ｌを選択する方法は典型的には、極配置、線形二次（ＬＱ）又はカルマンフィルタリング及びＨ∞設計である。ＭＥＭＳアクチュエータ１１０のようなＭＥＭＳアクチュエータを用いることで可能になる、より正確なモデルによれば、精度が悪いシステムモデルしか利用できないときに典型的に必要とされる、より従来的で且つロバストなＨ∞設計の代わりに、カルマンフィルタを用いることができるようになる。
【００１１】
推定器を用いる概念は、最新の制御理論の一部であり、システムの内部エネルギー蓄積又は「状態」を表すためのシステムモデルを有することによって、さらに良好な制御設計が可能になるという概念を前提としている。典型的には、モデルに基づく制御が時間領域において実行され、システムの数学的表現が確立される。
【数１】
2008014943.tif 000003
式（１）及び（２）において、
【数２】
2008014943.tif 000004
はモデル化されるシステム内の全てのエネルギー蓄積要素を表すベクトルであり、ｕはシステムへの入力を表し、ｘはシステムエネルギー蓄積要素の時間導関数であり、
【数３】
2008014943.tif 000005
は、エネルギー蓄積要素をそれらの時間導関数に関連付ける関数であり、
【数４】
2008014943.tif 000006
は、入力ｕをエネルギー蓄積要素の時間導関数に関連付ける関数である。式（２）では、出力ｙは、エネルギー蓄積要素
【数５】
2008014943.tif 000007
と入力ｕとの関数である。
【００１２】
典型的には、モデル化されるシステムは線形化される。線形モデルに類似しているが、いくかの非線形成分を有する中間形が存在することもある。式（１）及び（２）の場合の線形化されたモデルは、以下の式（３）及び（４）によって与えられる。
【数６】
2008014943.tif 000008
図１ｂのＳＩＳＯシステムの場合、ｙ及びｕ及びｄはスカラー量であり、
【数７】
2008014943.tif 000009
及び
【数８】
2008014943.tif 000010
は列ベクトルであり、一方、
【数９】
2008014943.tif 000011
は行ベクトルであり、
【数１０】
2008014943.tif 000012
は正方行列である。推定器が有効であるかは、典型的には、種々の要因によって決まり、その要因のうちの１つは、式が物理系をモデル化することができる精度である。ＭＥＭＳアクチュエータ１１０は典型的には、圧電アクチュエータよりも正確に、且つ一貫してモデル化することができるので、ＭＥＭＳアクチュエータ１１０によれば、より容易に、且つ正確にシステムをモデル化できるようになる。実例となるデータシステムの場合、式（３）及び（４）は以下のようになる。
ｘ（ｋ＋１）＝Φｘ（ｋ）＋Γｕ（ｋ）        （５）
ｙ（ｋ）＝Ｈｘ（ｋ）＋Ｄｕ（ｋ）            （６）
ただし、ｘ（ｋ）は離散時間状態であり、ｋ＝ｋＴｓであり、Ｔｓはサンプル周期であり、ｕ（ｋ）は時間ステップｋにおいて入力である。行列Φは、時間ステップ間の離散時間状態の伝搬を指している。それは典型的には、連続時間導関数行列Ｆ及びサンプル周期Ｔｓを用いて計算される。行列Γは、時間ステップｋにおける入力ｕを、次の時間ステップｘ（ｋ＋１）における状態ｘに関連付ける。状態ｘ及び入力ｕを出力ｙに関連付ける行列Ｈ及びＤは典型的には、連続線形システムから離散線形システムに移行する際に変化しない。ボールド体の取り決めは、制御理論における共通の表記法に従って書かれていることに留意し、それぞれの変数はベクトル及び行列であることを理解されたい。これまでの包括的な説明によって網羅される多数の打切り方法があり、これらは、G. Franklin、D. Powell及びM. Workmanによる「Digital Control of Dynamic Systems」並びにK. Astrom及びB. Wittenmarkによる「Computer Controlled Systems」のような、制御理論に関する標準的な教科書において見ることができる。
【００１３】
以下の説明は、離散時間における推定器を取り扱うが、それらの推定器は、連続時間又はサンプリングされたデータモデルの場合にも用いることができることは理解されたい。いくつかの事例における理論、解析及び設計は連続時間において行われることがあるが、離散時間と等価な形態が得られるほどサンプリングレートが十分に高く、且つ打切りの影響が十分に小さい場合には、実際の実施態様はサンプリングされたデータ方法又は離散時間方法で達成されることに留意されたい。連続形式では、選択されたフィルタを実装するために、アナログ回路が用いられる。サンプリングされたデータ形式の場合、データはサンプリングされるが、アナログ及びデジタルの方法の組み合わせを用いて、選択されたフィルタを実装することができる。離散形式では、典型的には、コンピュータ、デジタルシグナルプロセッサ（ＤＳＰ）又はプログラマブルロジック、たとえばフィールドプログラマブルゲートアレイ（ＦＰＧＡ）のようなデジタルプロセッサが用いられる。
【００１４】
式（５）及び（６）によって記述されるシステムの場合に、離散時間推定器は、システムモデルで開始して、推定器状態を１時間ステップだけ前方に伝搬させることによって形成される。しかしながら、これは、「開ループ」手法であり、式（５）及び（６）によってモデル化されるシステムが開ループ動作において安定している場合にのみ、その誤差が０まで減衰する。誤差が０に接近するのを保証するために、システム出力とモデル出力との間で測定された差を用いて、状態推定値が補正される。一般的な形における推定器の本発明による一実施形態は、以下の式によって与えられる。
【数１１】
2008014943.tif 000013
式（７）は、フィードバック利得行列Ｌを介して補正が追加された状態伝搬を示す。ただし、ｙ（ｋ）はシステム出力である。フィードバック利得行列Ｌの選択は、制御理論の中心的な問題であり、カルマンフィルタリング及びＨ∞設計及び混在Ｈ２／Ｈ∞設計の問題を含む。ｙ及びｕがスカラーである１入力１出力システムの場合、Ｇ、Ｈ及びＤは典型的には、小文字ｇ、ｈ及びｄで置き換えられる。
【００１５】
典型的には、任意の推定器において、重要な１組の数は、測定フィードバック利得行列Ｌの中にある。本発明による制御理論においてＬを求めるために数多くの方法がある。極配置法を用いることは、推定器の閉ループ極を設定することによって、推定器の動態を得ることを含み、閉ループ極は、離散時間定式化におけるΦ−ＬＨの極、又は連続時間定式化におけるＦ−ＬＨの極である。それらの極は、推定器の所望の特性に基づいて設定される。線形二次法を用いることは、システム及びセンサに入る雑音を、既知の平均及び分散を有する白色ガウス雑音としてモデル化することを含む。そのモデルに基づいて、最小二乗推定値が生成される。この方法も、離散時間定式化ではΦ−ＬＨの極、又は連続時間定式化ではＦ−ＬＨの極を配置するが、システムを駆動する雑音の量と、測定における雑音の量との間の最小二乗重み付けを用いる。推定器動態は、最小二乗コスト関数を最小にするように選択される。さらに、時間とともに一定であるのではなく、Ｌは時間とともに徐々に展開する。最後に、Ｐ行列を介して、推定値内の不確定性の測定が含まれる。さらに、最適化判定基準Ｈ∞によって典型的に知られているロバストな推定器を用いることもできるが、たとえば、上述した線形二次法よりも性能が落ちることがある。さらに、混在Ｈ２／Ｈ∞設計は典型的には、カルマンフィルタの性能（Ｈ２）とＨ∞設計のロバストネスとの間のバランスを求める。
【００１６】
本発明による別の実施形態は予測器であり、ｙの最新の測定値が用いられる推定値の形式は以下の式によって与えられる。
【数１２】
2008014943.tif 000014
式（９）は、状態推定値を１時間ステップだけ前方に伝搬させる。式（１０）は、状態推定値に基づいて、出力
【数１３】
2008014943.tif 000015
を予測する。式（１１）は、先行する時間ステップにおけるｙ（ｋ）の測定値に基づいて、式（９）の状態推定値の伝搬を補正する。
【００１７】
対照的に、本発明による別の実施形態は、下記の電流推定器形式を用いる。
【数１４】
2008014943.tif 000016
ただし、式（１２）、（１３）及び（１４）は、式（１４）は式（１１）とは異なるが、式（９）及び（１２）、ならびに式（１０）及び（１３）がそれぞれ同じであるように構成されている。電流推定器の手法では、状態推定値に対する補正は、最新の測定値ｙ（ｋ＋１）を用いて行われる。理想的には、この結果として、フィルタの中の伝搬遅延が１時間ステップだけ小さくなり、それにより、推定器の位相劣化が小さくなる。これは典型的には、推定器がフィードバック制御の場合に用いられるときに重要である。しかしながら、フィルタを計算するためのあらゆる実用的な手段は、デジタルロジックが式（１４）を計算するために０ではない、或る時間を必要とするので、結果として有限のサンプル時間遅延を生成する。予測推定器の場合、フルサンプル遅延が常に存在するが、これは、理論的なモデルにより良好に一致する。したがって、典型的には、予測推定器はモデル化するのが容易である。連続推定器の設計では、全てが連続であるので、システム状態を前方に伝搬させること及び測定補正を含むことの間に別個の時間ステップは存在しない。
【００１８】
式（７）〜式（１４）によれば、システムへの多入力が可能なことに留意されたい。システムモデルの観点から、Ｈ及びＬだけが変化する。入力を増やせることにより、本発明によるさらに別の実施形態が提供される。たとえば、いくつかのＭＥＭＳアクチュエータによれば、モータ位置を直に測定できるようになり、それにより、表面高に対して求められ、相対的な誤差測定値を与える相対的なＡＦＭプローブ先端位置のほかに、別の入力が与えられる。
【００１９】
推定器が、システムへの制御入力ｕ及び出力ｙの両方にアクセスする場合には、推定器が必要とするのは、物理系を推定することだけであり、物理系及びコントローラの閉ループの組み合わせを推定する必要はない。典型的には、ＡＦＭの場合、出力誤差だけが入手可能であり、基準値及び出力は入手することはできない。これは、推定器の形を変更する。
【００２０】
図１ｂに示されるような本発明による一実施形態は、先に言及されたように、推定器フィードバック利得を選択するためにカルマンフィルタを用いることができる。モデル化されるシステムに入るプロセス雑音及び測定雑音がガウス白色雑音であり、無相関である場合には、カルマンフィルタは、推定器フィードバック利得を選択するという問題に対する最小二乗解を生成する。カルマンフィルタ定式化のための出発点は、時間ステップ間の物理系の展開を記述するプロセスモデル方程式で開始する。
ｘ（ｋ＋１）＝Φｘ（ｋ）＋Γｕ（ｋ）＋Γｗ  ｗ（ｋ）        （１５）
ｙ（ｋ）＝Ｈｘ（ｋ）＋Ｄｕ（ｋ）                          （１６）
ｚ（ｋ）＝Ｈｘ（ｋ）＋Ｄｕ（ｋ）＋ｖ（ｋ）                （１７）
上記の式は、式（５）及び（６）を変更したものを表す。式（５）は、ｗ（ｋ）を加えることによって、式（１５）を生成するように変更されており、ｗ（ｋ）は、システムに入るプロセス雑音を表しており、入力行列Γｗを伴っているが、その行列は制御入力行列Γと同じであることも、同じでないこともある。式（１６）のシステム出力ｙ（ｋ）は、状態ｘ（ｋ）を通る場合を除いて、雑音によって汚染されないことに留意されたい。しかしながら、式（１７）の測定されたシステム出力ｚ（ｋ）は、出力雑音、及び測定雑音ｖ（ｋ）の両方を含む。以下の説明を簡単にするために、ｖ（ｋ）及びｗ（ｋ）はいずれも、互いに無相関である、０平均白色ガウス雑音とする。その説明は、０以外の平均を取り扱うように、且つ有色雑音、及びｗ（ｋ）とｖ（ｋ）との間の相関を含むように拡張することができることに留意されたい。プロセス雑音ｗ（ｋ）及び測定雑音ｖ（ｋ）は、共分散Ｒｗ及びＲｖを有するものとする。ただし、その次元は、それぞれ、システム入力及びシステム出力の数によって設定される。
【００２１】
推定器のカルマンフィルタバージョンは２つのステップに分割される。最初に、推定値
【数１５】
2008014943.tif 000017
を時間更新のために時間的に前方に伝搬させる。
【数１６】
2008014943.tif 000018
式（１８）において、
【数１７】
2008014943.tif 000019
はシステム状態
【数１８】
2008014943.tif 000020
及び入力ｕ（ｋ）の最良の推定値を用いる、次の時間ステップにおける状態の最良の予測値である。式（１９）において、Ｍ（ｋ＋１）は、状態推定値の不確定性の予測値である。測定ステップにおいて、推定値は、電流表面測定値を用いて補正される。代替的に、補正は、予測推定器の形で、先行する測定値を用いて行うことができる。
【００２２】
表面測定更新は以下の式によって与えられる。
【数１９】
2008014943.tif 000021
表面測定更新は、測定値からの新たなデータに合わせて調整される。ここで、
【数２０】
2008014943.tif 000022
は時刻ｋにおける表面測定後の状態の最良の推定値であり、式（２０）は、最後の測定値ｚ（ｋ）の結果を含むように状態推定値を調整する。Ｌ（ｋ）は推定器利得行列であり、予測される状態を補正し、状態推定値内の不確定性と測定値内の不確定性とを相殺するために、どの程度の測定値が用いられるかを判定する利得補正を表す。Ｐ（ｋ）は、測定が行われた後の状態推定値内の測定不確定性を表しており、それゆえ、式（２３）は、状態共分散行列Ｐ（ｋ）の推定値を調整する。
【００２３】
出力誤差フィードバックを与えるために光検出器１６０だけを用いる図１ｂに示されるような一実施形態では、システム出力ｙ（ｋ）の実際の測定値ｚ（ｋ）は入手できない。測定誤差だけを入手することができ、これは、表面変化、及びカンチレバー／プローブ先端１３０の任意の基準撓みを含む基準と、システム出力との間の差に相当する。
ｅ（ｋ）＝ｒ（ｋ）−ｙ（ｋ）＋ｖ（ｋ）          （２４）
それゆえ、出力が摂動するとき、その摂動が誤差推定器動態の影響を受けるまで、その摂動は推定器内に現われない。これは、典型的には、Ｈ∞のようなロバストな設計で典型的に利用できる推定器よりも広い帯域幅を有する推定器を有することが好都合であることを指示する。ロバストな推定器の場合、帯域幅は典型的には比較的狭いので、誤差は典型的には、安定するのに時間がかかるのに対して、より高い性能の推定器の場合、帯域幅は典型的には相対的に広いので、誤差は典型的には短い時間で抑圧される。図１ｂに示される実施形態の場合、式（２０）は以下のようになる。
【数２１】
2008014943.tif 000023
そして、測定更新は、式（２５）、（２１）、（２２）及び（２３）から成る。初期行列Ｒｗ、Ｒｖ、Ｍ（０）、Ｐ（０）及びＬ（０）で開始するとき、表面測定更新の式（２５）、（２１）、（２２）及び（２３）とともに、時間更新の式（１８）及び（１９）を用いて、時間ステップ毎にシステム状態推定値が与えられる。
【００２４】
システム状態推定値は、多様に用いることができる。図１ｂに示されるような本発明による一実施形態の場合、ＡＦＭのカンチレバー／プローブ先端１３０の位置の最小二乗推定値は以下の式によって与えられる。
【数２２】
2008014943.tif 000024
上記の式から、表面１５０の最小二乗推定値を求めることができる。制御入力ｕ（ｋ）から出力
【数２３】
2008014943.tif 000025
への直接的なフィードスルーが存在しない場合には、Ｄ＝０であり、式（２６）は以下のようになる。
【数２４】
2008014943.tif 000026
これは典型的には、入力が出力に直に作用する可能性がない多種多様な物理系に相当する。代わりに、入力ｕはシステム状態ｘに作用し、システム状態がさらに出力ｙに作用する。圧電アクチュエータ、ボイスコイルアクチュエータ及びＭＥＭＳアクチュエータのような典型的なアクチュエータは、この特徴を有する。
【００２５】
さらに、図１ｂに示されるような本発明による一実施形態の場合、ＡＦＭのカンチレバー／プローブ先端１３０の位置の最小二乗推定値は以下の式によっても与えられる。
【数２５】
2008014943.tif 000027
上記の式から、表面１５０の最小二乗推定値を求めることができる。制御入力ｕ（ｋ）から出力
【数２６】
2008014943.tif 000028
への直接的なフィードスルーが存在しない場合には、Ｄ＝０であり、式（２８）は以下のようになる。
【数２７】
2008014943.tif 000029
ここで、
【数２８】
2008014943.tif 000030
を生成する際に
【数２９】
2008014943.tif 000031
を用いることが因果的に関係しない限り、
【数３０】
2008014943.tif 000032
を表面推定値として用いることができる。この場合、
【数３１】
2008014943.tif 000033
は最新の測定値を利用するので、
【数３２】
2008014943.tif 000034
は典型的には、
【数３３】
2008014943.tif 000035
よりも改善された表面１５０の推定値である。
【００２６】
Ｐ（ｋ）は、
【数３４】
2008014943.tif 000036
における不確定性の推定値であり、一方、Ｍ（ｋ）は、
【数３５】
2008014943.tif 000037
における不確定性の推定値である。Ｐ（ｋ）又はＭ（ｋ）のいずれか、又はＰ（ｋ）及びＭ（ｋ）の両方を用いて、表面測定値における不確定性の推定値を与えることができ、
【数３６】
2008014943.tif 000038
又は
【数３７】
2008014943.tif 000039
のいずれかをフィードバックコントローラ１７０において用いることができる。推定値が、次の測定の前に１時間ステップだけ前方に伝搬されるときに、量
【数３８】
2008014943.tif 000040
及びＭ（ｋ）は時間更新に由来する。Ｍ（ｋ）は
【数３９】
2008014943.tif 000041
の共分散のための推定値である。量
【数４０】
2008014943.tif 000042
及びＰ（ｋ）は測定後に得られる。典型的には、全ての測定値が或る情報を含むものと仮定し、且つ典型的には、Ｐ（ｋ）、すなわち
【数４１】
2008014943.tif 000043
の共分散のための推定値が、Ｍ（ｋ）、すなわち
【数４２】
2008014943.tif 000044
の共分散のための推定値よりも小さくなるように、モデルによって伝搬される量に対して、測定値に適当な重みが割り当てられているものと仮定される。
【００２７】
本発明による一実施形態は、定常状態カルマンフィルタを用いることができる。定常状態カルマンフィルタを用いることは、初期条件からの不確定性が減衰する場合にＭ（ｋ）が定常状態値に収束するときに適しており、これは、Ｐ（ｋ）及びＬ（ｋ）がいずれも収束することを意味する。時間更新は以下のようになる。
【数４３】
2008014943.tif 000045
一方、表面測定値更新は以下のようになる。
【数４４】
2008014943.tif 000046
ただし、ＰＳＳはＰの定常状態値を示す。式（３３）によって与えられる推定値行列Ｌは、予め計算し、格納することができることに留意されたい。
【００２８】
図２は、ＡＦＭに基づくＭＥＭＳにおいて多数のセンサを使用することを示す。本発明の一実施形態によれば、ＭＥＭＳアクチュエータ２１０は、ＭＥＭＳアクチュエータ２１０のための絶対位置情報を与える位置センサを含む。これは、表面１５０の推定を補うために用いられる。アクチュエータ位置センサは、カンチレバー動態の作用を通じてプローブ先端位置に関連付けられるアクチュエータの絶対位置情報を与えるので、その絶対位置情報を用いて、光検出器１６０からの入力だけを用いて表面トポグラフィを推定するＡＦＭにおいて典型的に生じるドリフト問題を除去することができる。相対的な位置誤差しか与えられない場合には、光検出器１６０は、ＭＥＭＳアクチュエータ２１０の絶対的な空間位置を知ることができないので、ドリフト問題が生じる。ＭＥＭＳアクチュエータ２１０の絶対位置は、物理的及び動的にカンチレバー／プローブ先端１３０の極めて近くで測定される。それゆえ、表面１５０の測定値であるカンチレバー／プローブ先端１３０の位置と、ＭＥＭＳアクチュエータ２１０の測定された位置との間に存在する動態は制限される。さらに、さらに多くの情報を推定器２８０に追加することによって、典型的には、信号対雑音比が改善される。
【００２９】
推定器２８０のための式は、式（２５）、（２１）、（２２）及び（２３）から次元を変更しているが、一方、時間更新は変更されない。
【数４５】
2008014943.tif 000047
又は
【数４６】
2008014943.tif 000048
によって推定されるシステム状態は前と同じ次元を有するが、出力行列Ｈ及び推定器利得行列Ｌ（ｋ）は次元を変更する。システム状態が次元数ｎを有するものと仮定すると、Ｒｖは２×２行列であり、Ｈは２×ｎ行列であり、Ｌはｎ×２行列である。式（２５）、（２１）、（２２）及び（２３）は、それぞれ以下のようになる。
【数４７】
2008014943.tif 000049
ただし、Ｈ（２，：）はＨ行列の第２の行を表し、ｅ０（ｋ）は光センサ１６０から導出される誤差である。行列Ｌは、２つの入手可能な測定値を用いて、システム状態推定値
【数４８】
2008014943.tif 000050
を調整する。システム状態推定値
【数４９】
2008014943.tif 000051
又は
【数５０】
2008014943.tif 000052
が依然としてｎ×１列ベクトルであるとき、不確定性を表すそれらの共分散は典型的には、２つの測定値を入手可能であることに起因して低減される。ｅ０（ｋ）が光センサ１６０から導出される誤差項であるとき、第２の誤差項が、ＭＥＭＳアクチュエータ２１０の位置ｚｍをとり、その値から、推定されたモータ位置
【数５１】
2008014943.tif 000053
を減算することによって形成される。第２の誤差項
【数５２】
2008014943.tif 000054
はＭＥＭＳアクチュエータ２１０に由来する。それは、制御入力ｕ（ｋ）よりも、カンチレバー／プローブ先端１３０の位置Ｐに物理的且つ動的にはるかに近いので、第２の誤差項は、推定器に入る、さらに多くの位置情報を与え、結果として、システム状態及び表面１５０の推定値がさらに良好になる。初期行列Ｒｗ、Ｒｖ、Ｍ（０）、Ｐ（０）及びＬ（０）で開始するとき、表面測定値更新の式（３４）、（３５）、（３６）及び（３７）とともに、時間更新の式（１８）及び（１９）を用いて、時間ステップ毎にシステム状態推定値が与えられる。
【００３０】
状態推定値は多様に用いることができる。図２に示されるような本発明による一実施形態の場合、ＡＦＭのカンチレバー／プローブ先端１３０の位置の最小二乗推定値は、以下の式によって与えられる。
【数５３】
2008014943.tif 000055
その式から、表面１５０の最小二乗推定値を求めることができる。制御入力ｕ（ｋ）から出力
【数５４】
2008014943.tif 000056
への直接的なフィードスルーが存在しない場合には、Ｄ＝０であり、式（３９）は以下のようになる。
【数５５】
2008014943.tif 000057
【００３１】
さらに、図２に示されるような本発明による一実施形態の場合に、ＡＦＭのカンチレバー／プローブ先端１３０の位置の最小二乗推定値は、以下の式によっても与えられる。
【数５６】
2008014943.tif 000058
その式から、表面１５０の最小二乗推定値を求めることができる。制御入力ｕ（ｋ）から出力
【数５７】
2008014943.tif 000059
への直接的なフィードスルーが存在しない場合には、Ｄ＝０であり、式（４１）は以下のようになる。
【数５８】
2008014943.tif 000060
ここで、
【数５９】
2008014943.tif 000061
を生成する際に
【数６０】
2008014943.tif 000062
を用いることが因果的に関係しない限り、
【数６１】
2008014943.tif 000063
を表面推定値として用いることができる。この場合、
【数６２】
2008014943.tif 000064
は最新の測定値を利用するので、
【数６３】
2008014943.tif 000065
は典型的には、
【数６４】
2008014943.tif 000066
よりも改善された表面１５０の推定値である。
【００３２】
再び、本発明による一実施形態は、定常状態カルマンフィルタを用いることができる。定常状態カルマンフィルタを用いることは、初期条件からの不確定性が減衰する場合にＭ（ｋ）が定常状態値に収束するときに適しており、これは、Ｐ（ｋ）及びＬ（ｋ）がいずれも収束することを意味する。時間更新は以下のようになる。
【数６５】
2008014943.tif 000067
一方、表面測定値更新は以下のようになる。
【数６６】
2008014943.tif 000068
【００３３】
先の定常状態カルマンフィルタの実施形態の場合のように、Ｐ行列及びＬ行列は一定であり、Ｌは予め計算して、格納することができる。
【００３４】
本発明は、具体的な実施形態に関連して説明されてきたが、これまでの説明に鑑みて、数多くの改変、変更及び変形が明らかであろうことが当業者には自明である。したがって、本発明は、添付の特許請求の範囲の精神及び範囲内に入る全ての他の改変、変更及び変形を含むことを意図している。
【図面の簡単な説明】
【００３５】
【図１ａ】本発明による、典型的な原子間力顕微鏡の物理的な装置の概略図である。
【図１ｂ】本発明による一実施形態のブロック図である。
【図２】本発明による一実施形態のブロック図である。
【図１ａ】
2008014943.tif 000069
【図１ｂ】
2008014943.tif 000070
【図２】
2008014943.tif 000071
【外国語明細書】
2008014943000001.pdf