Comparison of Silvaco and Synopsys TCAD Predictions Including the Perugia Radiation Damage Model in Silicon Pixel Detectors for the HL-LHC

本論文は、高光度LHCで予想される極端な放射線環境下におけるシリコンピクセル検出器の主要性能指標の予測精度を評価するため、ペルージャ放射線損傷モデルを組み込んだSilvacoおよびSynopsys TCADシミュレーションを比較する。

要約

大型ハドロン衝突型加速器（LHC）を、巨大で高速な粒子の競走トラックだと想像してください。科学者たちはこれを「高輝度」バージョン（HL-LHC）にアップグレードしており、これは粒子同士をこれまでよりもはるかに頻繁に衝突させることを意味します。問題は何かというと、この激しい交通量によって大量の「放射線塵」が発生し、衝突の様子を撮影しようとするシリコンセンサー（カメラ）を損傷させてしまうことです。

時間の経過とともに、この放射線塵はシリコンセンサーを「漏れやすく」「硬く」なるデバイスに変えてしまいます。センサーは信号を収集する能力を失い始めます（カメラがピントを失うようなもの）し、動作にははるかに高い電圧が必要となり、それが破損のリスクとなります。

これを起こる前に修正するため、科学者たちはコンピュータシミュレーションを用いて、長年の放射線曝露後にセンサーがどのように振る舞うかを予測します。彼らは以下のことを知る必要があります：必要な電圧はどれくらいか？漏れる電流はどれくらいか？センサーは依然として機能するか？

2 人の「天気予報士」

この論文では、研究者たちはこれらのセンサーの天気予報士のように機能する 2 つの異なるコンピュータプログラム（TCAD ツール）をテストしています：

1. Synopsys
2. Silvaco

両方のプログラムは、「ペルージャ放射線損傷モデル」と呼ばれる特定の規則セットを使用します。このモデルを、放射線塵がシリコンをどのように損傷し、電気の流れを乱す微小なトラップやホールを生成するかをコンピュータに正確に伝える詳細な取扱説明書だと考えてください。

この論文の目的は、同じ取扱説明書を使用した場合、これら 2 人の異なる「天気予報士」が同じ予測を与えるかどうかを確認することです。もし合致すれば、そのマニュアルは信頼性があり、科学者たちはどのソフトウェアを使用しても予測を信頼できることを意味します。

実験：微小なシリコンダイオード

研究者たちは、50 マイクロメートル（人間の髪の毛の幅程度）の厚さを持つ微小なシリコンセンサー（ダイオード）の仮想 2 次元モデルを構築しました。彼らは 2 つのシナリオをシミュレーションしました：

1. 新鮮なセンサー： 放射線が当たる前。
2. 放射線被曝したセンサー： 大量の放射線にさらされた後（HL-LHC の過酷な環境をシミュレート）。

これらのセンサーを 2 つの温度でテストしました：涼しい 248 K（約 -25°C）と暖かい 300 K（室温）。

結果：天気予報士は合意するか？

1. 新鮮なセンサー（非被曝）
センサーが新品のとき、2 つのコンピュータプログラムは約 500 ボルトまで、センサーを流れる電流量と電荷の蓄え方についてほぼ完全に一致しました。

• 不一致： 電圧を非常に高く押し上げると（700 ボルト付近）、センサーが「破損」する（バークダウンする）タイミングについて、プログラム間でわずかな不一致が生じ始めました。著者らは、これはおそらく 2 つのプログラムがセンサーを描画するためにわずかに異なるデジタル「グリッド」（メッシュ）を使用しているためだと示唆しています。これは、2 つの異なる地図アプリが道路をわずかに異なって描画するのと似ています。

2. 放射線被曝したセンサー（真のテスト）
ここで真の魔法が起きました。彼らは放射線にさらされた後のセンサーをシミュレーションしました。

• リーク電流： 2 つのプログラムは「リーク」（不要な電流）をほぼ同一に予測しました。
• 空乏電圧： センサーを再び動作させるために必要な電圧について、両者は完全に合意しました。
• 電界： 彼らはシリコン内部の目に見えない電気的な力をマッピングしました。センサーの中央部（バルク）では、2 つのプログラムはほぼ完全に一致しました（互いの 1% 以内）。
• 「トラップ」： また、放射線によって生成され電子を捕獲する微小な「トラップ」も調べました。2 つのプログラムは、これらのトラップの挙動について非常に妥当な範囲（約 20%）内で合意しました。

温度のひねり：
室温（300 K）では、プログラムは最高レベルの放射線被曝において少し不一致が見られました。しかし、著者らはこれは大きな懸念事項ではないと指摘しています。なぜなら、現実世界では、これらの損傷したセンサーはほぼ決して室温で動作させられず、生存のために非常に冷たく保たれているからです。したがって、重要なのは低温（248 K）での合意であり、そこでは 2 つのプログラムは完璧に同期していました。

結論

この論文は、Synopsys と Silvaco は、全く同じレシピ（ペルージャモデル）に従う 2 人の異なるシェフのように、同じ美味しい料理を作り出すと結論付けています。

ソフトウェアツールが異なっていても、ペルージャ放射線損傷モデルを使用すれば、将来の HL-LHC の過酷な放射線環境においてシリコンセンサーがどのように生存するかについての予測はほぼ同一になります。これにより、科学者たちは彼らのモデルが確固たるものであると確信し、次世代の粒子検出器を設計するためにどちらのツールでも使用できるという信頼を得ています。

注：著者らは将来「収集された電荷」を検討する計画を述べていますが、この論文は電圧、電流、および電界に厳密に焦点を当てました。

技術概要

論文「HL-LHC 用シリコンピクセル検出器におけるペルージャ放射線損傷モデルを含む Silvaco および Synopsys TCAD 予測の比較」の詳細な技術的要約を以下に示す。

1. 問題提起

高光度大型ハドロン衝突型加速器（HL-LHC）は、現在のレベルの 5 倍から 10 倍の放射線フラックスにシリコンピクセル検出器をさらすことになる。この激しい放射線は、シリコンセンサーに体積および表面損傷を引き起こし、以下の結果をもたらす。

• リーク電流の増加。
• 動作電圧（空乏電圧）のシフト。
• 信号損失および電荷収集の劣化。
• フロントエンド電子回路の動作変化。

これらの影響を正確に予測することは、動作電圧の見積もりや再構成アルゴリズムの訓練に不可欠である。特定の体積および表面欠陥物理学を組み込んだペルージャ大学放射線損傷モデルは、すでにSynopsys Sentaurus TCADで成功裏に実装されているが、そのSilvaco TCADへの移植性と予測能力を検証する必要があった。これらの 2 つの主要な業界標準シミュレーションツール間での相互検証が行われなければ、放射線損傷モデリングの堅牢性は、より広範な採用に対して未確認のままとされる。

2. 手法

著者らは、SynopsysおよびSilvacoの両 TCAD ツールを用いたデバイスシミュレーションの比較研究を実施した。

• デバイス構造：n-on-p ダイオードの 2 次元シミュレーションを使用。
  • 厚さ：50 µm。
  • 体積濃度：$5 \times 10^{12} \text{ cm}^{-3}$（ホウ素）。
  • ピッチ：20 µm（小さな収集電極付き）。
  • 表面：1 µm の酸化膜層。
• 物理モデル：公平な比較を確保するため、両ツールは同一の物理モデルを使用。
  • 統計：フェルミ - ディラック統計。
  • 再結合：Scharfetter - Shockley Read Hall。
  • バンド間トンネル効果：Schenk。
  • バンドギャップ縮小：Slotboom。
  • 衝撃電離：van Overstraeten および de Man。
• 放射線損傷の実装：ペルージャモデルを Silvaco へ移植。これには以下が含まれる。
  • 体積モデル：フラックス（$\Phi$）に比例する密度を持つ 2 つのアクセプター欠陥と 1 つのドナー欠陥。高フラックス（$> 3 \times 10^{15} \text{ neq/cm}^2$）を処理するための、アクセプター生成（$N_A$）のフラックス依存関数を利用。
  • 表面モデル：欠陥の連続的なエネルギー分布（上部バンドギャップではアクセプター様、価数帯近傍ではドナー様）および線形に増加する酸化膜電荷（線量に比例）。
• シミュレーションシナリオ：
  • フラックス：非照射、$\Phi = 1.5 \times 10^{15} \text{ neq/cm}^2$、および$\Phi = 1.0 \times 10^{16} \text{ neq/cm}^2$。
  • 温度：$T = 248 \text{ K}$および$T = 300 \text{ K}$。
• 観測量：電流 - 電圧（IV）、静電容量 - 電圧（CV）、電界プロファイル、およびトラップ電離確率。

3. 主要な貢献

• モデルの移植性：複雑なペルージャ体積および表面放射線損傷モデルを、Synopsys から Silvaco TCAD へ成功裏に移植。
• 相互検証：HL-LHC に関連する高放射線環境をシミュレートする際、これら 2 つの異なる TCAD エンジンの直接的な比較を初めて提供。
• 欠陥統計の分析：巨視的な電気的特性だけでなく、微視的なトラップ電離確率も分析し、欠陥物理学の実装の一貫性に関する深い洞察を提供。

4. 結果

非照射デバイス

• IV 特性：約 500 V までツール間で優れた一致。高電圧では不一致が見られた（Synopsys は約 700 V でブレークダウンを予測、Silvaco は約 760 V で予測）、これはおそらくメッシュ生成の違いによるもの。
• CV 特性：両ツールとも約 10 V の空乏電圧を正確に予測し、約 55 V における「ステップ」特徴（酸化膜下の空乏に起因）と一致。

照射デバイス

• リーク電流（IV）：
  • 248 K：すべてのフラックスで驚くほど良い一致。わずかな偏差は、非常に高い電圧および最高フラックスでのみ発生し、メッシュの違いに起因。
  • 300 K：最高フラックスでは定性的な一致のみ。ただし、著者らは、このような高照射後には 0°C 以上でデバイスを動作させることは稀であるため、これは重要度が低いと指摘。
• 空乏電圧（CV）：248 K において、両フラックスとも 2 ツール間の一致は「実質的に完璧」であった。
• 電界プロファイル：
  • 体積領域において、電界プロファイルは両フラックスおよびすべてのバイアス電圧（100 V、300 V、500 V）で**1%**以内で一致。
  • 前面および背面のイオン注入近傍でわずかな不一致が観測されたが、これもメッシュの不一致およびイオン注入モデリングの違いに起因。
• トラップ電離確率：
  • 体積欠陥（アクセプタートラップ 1 および 2）について、電離確率は 1 桁以内で一致し、体積領域では**20%**以内で一致。
  • この一致レベルは、電界におけるパーセントレベルの一致を考慮すると、満足すべきものとみなされる。

5. 意義

• ペルージャモデルの検証：本研究は、ペルージャ大学放射線損傷モデルが堅牢であり、プラットフォームに依存せず、異なる TCAD ソフトウェア環境で一貫した結果をもたらすことを確認。
• HL-LHC への信頼性：リーク電流、空乏電圧、電界といった重要な観測量における高い一致度は、HL-LHC 条件下でのセンサー性能予測にこれらのシミュレーションを使用する信頼性を提供。
• 運用ガイダンス：これらのモデルを用いて安全な動作電圧を推定し、信号損失メカニズムを理解することを支持する結果は、アップグレードされた ATLAS および CMS 検出器の設計と再構成アルゴリズムに不可欠。
• 将来の作業：著者らは、この検証を収集電荷のシミュレーションへ拡張する計画。収集電荷は、追跡性能の究極の指標である。

結論：本論文は、ペルージャ放射線損傷モデルを使用する場合、Silvaco および Synopsys TCAD ツールが照射されたシリコンピクセル検出器に対して相互に一貫した予測を生み出すことを実証し、将来の高光度コライダー応用におけるモデルの妥当性を検証した。