Nonequilibrium Green's Function Formalism Applicable to Discrete Impurities in Semiconductor Nanostructures

本論文は、半導体ナノ構造における不純物の離散的性質を明示的に考慮するため、短距離散乱と長距離ハートリーポテンシャルを分離する新しい非平衡グリーン関数枠組みを提示し、不純物散乱の内在的な非局所性を明らかにするとともに、準一次元系における輸送特性へのその顕著な影響を実証する。

要約

半導体ナノ構造（将来のコンピュータチップに使用される微小なワイヤなど）を、長く細い廊下だと想像してください。この廊下の中を、電子が電気電流を運ぶために一端から他端へと走ろうとしています。しかし、この廊下は空っぽではなく、製造過程で偶然残された「不純物」、つまり汚れや破片で満たされています。これらの不純物は、実際には個々の原子（ドナー）であり、障害物として機能します。

長年にわたり、科学者たちはこれらの障害物を、廊下全体に均等に広がった滑らかで目に見えない霧だと仮定してモデル化してきました。障害物が非常に多いため、電子は平均的な「抵抗の雲」しか見えないと想定されたのです。これは、広くて大きな廊下では機能しました。

しかし、現代技術における微小で極薄のワイヤでは、この「霧」という考え方は破綻します。廊下があまりにも狭いため、各々の汚れの粒子の具体的な「位置」が重要になります。もし粒子が経路の真ん中にあれば、電子をブロックします。もし横にあれば、電子はすり抜けるかもしれません。古い「霧」モデルはこの決定的な詳細を見落としています。

新しい枠組み：二種類のトラブル

この論文は、野中伸幸氏によって、不純物を霧ではなく、個別の点として扱うことで、これらの微小ワイヤを電子がどのように移動するかを計算する新しい方法を提案しています。著者は、巧妙なアナロジーを用いて、不純物の問題を二つの部分に分けます。

1. 「長距離」部分（近隣効果）： 不純物を廊下に立っている人だと想像してください。たとえ触れていなくても、その存在はわずかに雰囲気を 변화させます。遠くから人を押し返したり、引き寄せたりするかもしれません。物理学では、これが「長距離」の電場です。論文では、これを全員に影響を与える滑らかな自己無撞着な背景ポテンシャル（廊下内の緩やかな傾斜のようなもの）として扱います。
2. 「短距離」部分（つまずきの危険）： これは、不純物の上に直接足を踏み入れたときに転倒する、即座に鋭い段差です。これが「短距離」散乱です。論文では、これを電子が特定の不純物原子に非常に接近したときにのみ起こる、特定の局所的な衝突として扱います。

「ゴースト」座標系

この論文における最も驚くべき発見は、これらの衝突が「どこ」で起こるかに関するものです。

従来の物理学では、衝突は地図上の特定の場所（実空間）で起こると考えられています。不純物が位置 X にあれば、衝突も X で起こります。

しかし、この論文は、これらの微小ワイヤの量子世界では、衝突の「位置」は実際には電子が「どこにいたか」と「どこに向かっているか」の混合であると示しています。著者は、電子の経路の「重心」を記述するために、ウィグナー座標と呼ばれる数学的ツールを使用しています。

アナロジー：
運動のぼやけと考えてください。高速で移動する車の写真を撮ると、それは一つの正確な点としてではなく、滲みとして見えます。論文は、「散乱率」（電子が不純物に跳ね返る確率）は地図上の単一の点に縛られているのではなく、電子の旅の「平均的な位置」に結びついていると主張します。

これは、散乱が非局所的であることを意味します。電子は、不純物に触れたときだけでなく、過去と未来の位置の間のより広範で曖昧な関係に基づいて不純物を感じ取ります。まるで電子が、物理的な接触点を超えて延びる障害物に対する「幽霊のような」感覚を持っているかのようです。

新しいモデルを使用するとどうなるか？

著者はこの新しい数学を適用して円筒形のワイヤ（ナノワイヤ）をシミュレーションし、古い「霧」モデルと比較しました。

• 古いモデル（局所的/対角）： 散乱が単一の点で起こり、単純な壁のように作用すると仮定しています。このモデルは、電子の移動度（移動速度）を過大評価する傾向があります。電子は実際よりも障害物に「混乱」していないと考えるからです。
• 新しいモデル（非局所的/非対角）： 衝突の「曖昧」な性質を考慮しているため、電子が「位相コヒーレンス」（同期したリズム）を失う速度がはるかに速いことを示しています。電子はより簡単に混乱し、散乱されます。
• 結果： 新しいモデルは、不純物の数が中程度（少なすぎず、多すぎず）の場合、特に電気電流と移動度が古いモデルが示唆したものよりも実際には低いと予測しています。

「自己平均化」の驚き

論文はまた、平均に関する興味深い発見もしています。不純物のさまざまなランダムな配置を多数取り、それらを平均化する場合（遠くから群衆を見るようなもの）、新しい「非局所的」モデルは、驚くほど古い「霧」モデルの結果と一致します。

しかし、不純物の特定の配置を持つ「単一の」特定のワイヤを見ると、古いモデルは完全に失敗します。汚れの粒子がわずかに異なる場所に落ちただけで、一つ一つの微小なワイヤで発生する性能の激しい変動を見逃してしまうのです。

まとめ

この論文は、微小ワイヤの量子世界をナビゲートするためのより正確な「地図」を提供します。不純物を滑らかな霧や単純な点の段差として扱うだけではならないことを教えてくれます。量子領域では、衝突の「位置」は少し曖昧であり、電子の経路全体に依存することを認める必要があります。これを行うことで、次世代のコンピュータチップにおける電気が実際にどれほど速く流れることができるかについての真実の姿が得られ、それらが以前考えられていたよりもわずかに遅く（そして変動が大きい）であることが明らかになります。

技術概要

Nobuyuki Sano による論文「半導体ナノ構造における離散不純物に適用可能な非平衡グリーン関数形式」の詳細な技術的要約を以下に示す。

1. 問題提起

ドープされた半導体ナノ構造において、電子輸送特性は室温であってもデバイス間で顕著なばらつきを示す。このばらつきは、従来のモデリング手法では正確に捉えられない不純物（ドナー）の離散的な性質に起因する。

• 従来の NEGF の限界: 標準的な非平衡グリーン関数（NEGF）スキームは、通常、不純物を連続的な「ジェリウム」背景として扱う「自己平均化」を仮定している。これは、デバイスが位相コヒーレンス長よりもはるかに大きいと仮定し、微視的な詳細を平均化するものである。
• 代替手法の限界: ハミルトニアン（未摂動）において不純物を明示的な点電荷として扱うと、特異なポテンシャルが生じ、有限差分法において数値的不安定性を引き起こす。
• ギャップ: 数値的特異性を伴わずに、長距離ポテンシャル変動と短距離散乱を正しく区別しつつ、不純物の離散的空間分布を NEGF 自己エネルギーに組み込むための厳密な理論的枠組みが必要である。

2. 手法

著者は、離散不純物を処理するように特別に設計された、NEGF スキームとポアソン方程式を自己的一貫的に結合する新しい理論的枠組みを提案する。

A. 不純物ポテンシャルの分離

離散不純物の静電ポテンシャルは、2 つの成分に分割される：

1. 長距離成分（$V_{long}$）: 自己一貫的なハートリーポテンシャルとして扱われる。キャリア密度と結合したポアソン方程式を解くことで、不純物の集団的遮蔽によって引き起こされるポテンシャル変動を考慮する。
2. 短距離成分（$V_{short}$）: 自己エネルギー（$\Sigma_{imp}$）に含まれる散乱ポテンシャルとして扱われる。これは局所的な散乱事象を表す。

B. 新しい自己エネルギー定式化

核心的な革新は、不純物散乱の自己エネルギー（$\Sigma_{imp}$）の導出にある：

• ウィグナー座標: 自己エネルギーは、標準的な実空間座標ではなく、ウィグナー座標（$X = (r+r')/2$, $\xi = r-r'$）を用いて導出される。
• 非局所性: 得られる自己エネルギー行列は、実空間において本質的に非対角である。散乱率の位置依存性は、局所的な実空間座標（$r$）ではなく、ウィグナー表現における「重心」座標（$X$）によって支配される。
• 数学的形態: 自己エネルギーは、遮蔽されたポテンシャルの二乗のフーリエ変換である相関関数 $C_X(\xi)$ を含む、離散不純物位置（$R_l$）の和として表される。
  $$\Sigma^\alpha_{imp}(r, r'; E) \approx \sum_l \delta(X - R_l) C_X(\xi) G^\alpha(X, \xi; E)$$

C. 散乱率への還元

本論文は、この新しい形式に対して不純物分布上で「自己平均化」（アンサンブル平均）を行うと、フェルミ黄金律から導出された標準的な散乱率が回復することを示している。しかし、散乱率における位置依存性は、局所座標（$r$）ではなく、本質的にウィグナー座標（$X$）に結びついていることを明確にしている。

D. 適用モデル

この枠組みは、準 1 次元近似の下で円筒形細線（半径 $\rho_s = 4$ nm、長さ $L_c = 30$ nm）に適用される：

• 3 次元ポアソン方程式は、ワイヤー軸に沿った「凍結」された長距離ポテンシャルを提供する。
• NEGF 計算は、最低サブバンド近似を用いて 1 次元で行われる。
• 散乱率行列は、離散不純物位置に基づいて明示的に構築される。

3. 主要な貢献

1. 位置依存散乱の厳密な導出: 不均一な不純物プロファイル下における不純物散乱率の位置依存性を体系的に導出し、散乱率が局所的な実空間座標ではなく、重心座標（ウィグナー）に依存することを証明した。
2. 本質的な非局所性: 不純物散乱が空間的に本質的に非局所的であることを確立した。散乱率行列は非対角であり、量子位相干渉を反映しており、従来の「局所」近似で用いられる対角行列とは対照的である。
3. 統合された枠組み: 離散不純物モデリングと従来のアンサンブル平均モデルの間のギャップを成功裡に埋めた。新しい形式は自己平均化により標準的な「ジェリウム」モデルに帰着するが、特定の離散構成をシミュレートする能力を保持している。
4. 特異性の解決: 短距離部分をハミルトニアン内の特異ポテンシャルではなく、自己エネルギー内の散乱ポテンシャルとして扱うことで、点電荷シミュレーションに伴う数値的不安定性を回避した。

4. 結果と考察

この枠組みは、異なるドナー密度（$\bar{N}_d^+$ が $10^{18}$ から $10^{20}$ cm$^{-3}$）を持つシリコンナノワイヤーに適用された。

• ポテンシャル変動: 離散不純物は、数十 mV のオーダー（高ドープでは最大約 100 mV）のポテンシャル変動を生み出す。これらの変動は、室温であっても局所的な電子の局在化を引き起こすのに十分な大きさである。
• 散乱率行列:
  • 離散の場合: 行列は非対角、対称であり、不純物位置の周に局在化している。非対角要素は負になり得て、強い位相干渉（ウィグナー関数に類似）を示している。
  • 自己平均化（非局所）の場合: 平均化後であっても、ウィグナー座標依存性により行列は非対角のまま残る。
  • 局所近似（従来）の場合: 対角行列を仮定しており、位相の乱れ効果を捉え損なう。
• 輸送特性（電流と移動度）:
  • 従来モデルでの過大評価: 従来の「局所」自己平均化モデルは、弾性散乱による位相の乱れを過小評価するため、電子移動度と電流を著しく過大評価する。
  • 離散対アンサンブル: 中程度のドープ領域（$10^{18} \leq \bar{N}_d^+ \leq 10^{19}$ cm$^{-3}$）では、特定のデバイスの輸送特性は、正確な不純物構成に依存して激しく変動する。「非局所自己平均化」モデルはアンサンブル平均の良い近似を提供するのに対し、「局所」モデルはそうではない。
  • 定量的データ: $\bar{N}_d^+ = 10^{19}$ cm$^{-3}$ の場合、局所モデルは移動度を215 cm$^2$/V/sと予測したが、より正確な非局所自己平均化モデルは100 cm$^2$/V/s、離散シミュレーションは95 cm$^2$/V/sを予測した。

5. 意義

• デバイスばらつき: この研究は、ドナー数が少なく位置がランダムである次世代半導体デバイス（FinFET、ナノワイヤー FET など）における統計的ばらつきを理解し予測するための理論的基盤を提供する。
• モデリングの精度: 散乱率に対する「局所」近似がナノ構造には不十分であることを浮き彫りにした。散乱の非局所性（ウィグナー座標依存性）を無視することは、移動度と電流の誤った予測につながる。
• 将来の方向性: この枠組みは、他の散乱機構（フォノン散乱、表面粗さなど）にも一般化可能であり、不規則系における量子輸送の本質的な特徴として非局所性を捉え、正確なデバイスシミュレーションで考慮する必要があることを示唆している。

結論として、Sano は不純物の離散的性質を正しく扱う堅牢な自己一貫的 NEGF 形式を提示し、不純物散乱が本質的に非局所的であり、従来の局所近似がナノスケールデバイスにおける輸送性能を著しく過大評価していることを明らかにした。