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

본 논문은 고광도 LHC 에서 예상되는 극한의 방사선 조건 하에서 실리콘 픽셀 검출기의 주요 성능 지표에 대한 예측 정확도를 평가하기 위해 페루기아 방사선 손상 모델을 통합한 실바코와 신옵시스 TCAD 시뮬레이션을 비교한다.

핵심

대형 강입자 충돌기 (LHC) 를 거대하고 초고속인 입자 레이스 트랙으로 상상해 보세요. 과학자들은 이를 "고광도 (High-Luminosity)" 버전 (HL-LHC) 으로 업그레이드하고 있는데, 이는 입자들을 훨씬 더 빈번하게 충돌시킨다는 것을 의미합니다. 문제는 무엇일까요? 이 치열한 교통 흐름은 실리콘 센서 (즉, 충돌 현상을 촬영하는 카메라) 를 손상시키는 많은 양의 '방사선 먼지'를 생성합니다.

시간이 지남에 따라 이 방사선 먼지는 실리콘 센서를 '누수'가 발생하고 '뻣뻣한' 장치로 변모시킵니다. 이들은 신호를 수집하는 능력을 잃기 시작합니다 (카메라가 초점을 잃는 것과 유사하게). 또한 작동하려면 훨씬 더 높은 전압이 필요해지는데, 이는 장치를 파괴할 위험을 초래합니다.

이 문제가 발생하기 전에 이를 해결하기 위해 과학자들은 컴퓨터 시뮬레이션을 사용하여 수년 간의 방사선 노출 후 센서가 어떻게 행동할지 예측합니다. 그들은 다음과 같은 사항을 알아야 합니다: 우리는 얼마나 많은 전압이 필요한가? 누출되는 전류는 얼마나 될까? 센서가 여전히 작동할 것인가?

두 명의 "기상 예보관"

이 논문에서 연구자들은 이러한 센서들을 위한 기상 예보관 역할을 하는 두 가지 다른 컴퓨터 프로그램 (TCAD 도구) 을 테스트하고 있습니다:

1. Synopsys
2. Silvaco

두 프로그램 모두 **"페루기아 방사선 손상 모델 (Perugia Radiation Damage Model)"**이라는 특정 규칙 세트를 사용합니다. 이 모델을 방사선 먼지가 실리콘을 손상시키고 전류 흐름을 방해하는 미세한 함정과 구멍을 생성하는 방식을 컴퓨터에 정확히 알려주는 상세한 사용 설명서로 생각할 수 있습니다.

이 논문의 목표는 동일한 사용 설명서를 사용할 때 이 두 가지 다른 "예보관"이 동일한 예측을 제공하는지 확인하는 것입니다. 만약 그들이 동의한다면, 그 설명서는 신뢰할 수 있으며 과학자들은 어떤 소프트웨어를 사용하든 예측을 신뢰할 수 있음을 의미합니다.

실험: 미세한 실리콘 다이오드

연구자들은 50 마이크로미터 두께 (인간 머리카락 너비 정도) 의 미세한 실리콘 센서 (다이오드) 의 가상 2 차원 모델을 구축했습니다. 그들은 두 가지 시나리오를 시뮬레이션했습니다:

1. 신선한 센서: 어떤 방사선도 충돌하기 전.
2. 방사선 조사된 센서: HL-LHC 의 가혹한 환경을 시뮬레이션하는 막대한 양의 방사선에 충돌한 후.

그들은 이 센서들을 두 가지 온도에서 테스트했습니다: 차가운 248 K(약 -25°C) 와 따뜻한 300 K(실내 온도).

결과: 예보관들이 동의할까요?

1. 신선한 센서 (비조사)
센서가 갓 제작되었을 때, 두 컴퓨터 프로그램은 약 500 볼트까지 센서를 통해 흐르는 전류량과 전하 저장 방식에 대해 거의 완벽하게 동의했습니다.

• 불일치: 전압을 매우 높게 (약 700 볼트 근처) 올렸을 때, 프로그램들은 센서가 언제 "파손 (breakdown)"될 것인지에 대해 약간씩 다르게 예측하기 시작했습니다. 저자들은 이것이 두 프로그램이 센서를 그리기 위해 약간 다른 디지털 "격자 (meshes)"를 사용하기 때문일 가능성이 높다고 제안합니다. 이는 두 가지 다른 지도 앱이 도로를 약간 다르게 그리는 것과 유사합니다.

2. 방사선 조사된 센서 (실제 테스트)
이곳에서 진정한 마법이 일어났습니다. 그들은 센서가 방사선에 폭격당한 후를 시뮬레이션했습니다.

• 누설 전류: 두 프로그램 모두 "누설" (원치 않는 전류) 을 거의 동일하게 예측했습니다.
• 공핍 전압: 두 프로그램은 센서가 다시 작동하도록 하는 데 필요한 전압량에 대해 완벽하게 동의했습니다.
• 전기장: 그들은 실리콘 내부의 보이지 않는 전기력을 매핑했습니다. 센서의 중앙 ("벌크") 에서 두 프로그램은 거의 완벽하게 일치했습니다 (서로 1% 이내).
• "함정": 그들은 전자를 잡는 방사선에 의해 생성된 미세한 "함정"들도 살펴보았습니다. 두 프로그램은 이러한 함정의 행동에 대해 매우 합리적인 오차 범위 (약 20%) 내에서 동의했습니다.

온도의 반전:
실내 온도 (300 K) 에서 프로그램들은 가장 높은 방사선 수준에서 약간 더 불일치했습니다. 그러나 저자들은 이것이 큰 우려 사항은 아니라고 지적합니다. 왜냐하면 실제 세계에서는 이러한 손상된 센서가 거의 실내 온도에서 작동하지 않기 때문입니다. 대신 생존을 위해 매우 차갑게 유지됩니다. 따라서 차가운 온도 (248 K) 에서의 일치가 실제로 중요하며, 그곳에서 두 프로그램은 완벽하게 동기화되었습니다.

결론

이 논문은 Synopsys 와 Silvaco 가 정확히 동일한 레시피 (페루기아 모델) 를 따르는 두 명의 다른 요리사처럼, 결국 동일한 맛있는 요리를 만들어낸다고 결론 내립니다.

소프트웨어 도구가 다르더라도 페루기아 방사선 손상 모델을 사용할 때, 그들은 미래 HL-LHC 의 가혹한 방사선 환경에서 실리콘 센서가 어떻게 생존할 것인지에 대한 거의 동일한 예측을 생성합니다. 이는 과학자들이 그들의 모델이 견고하다는 확신을 갖게 하며, 차세대 입자 검출기를 설계하기 위해 두 도구 중 어느 것을 사용해도 된다는 신뢰를 줍니다.

참고: 저자들은 향후 "수집된 전하"를 조사할 계획이라고 언급했지만, 이 논문은 전압, 전류, 전기장에만 엄격하게 초점을 맞추었습니다.

기술 요약

"HL-LHC용 실리콘 픽셀 검출기에서 페루기아 방사선 손상 모델을 포함한 Silvaco 와 Synopsys TCAD 예측 비교"에 대한 상세 기술 요약은 다음과 같습니다.

1. 문제 제기

고광도 대형 강입자 충돌기 (HL-LHC) 는 기존 수준보다 5 배에서 10 배 높은 방사선 선속 (fluence) 을 실리콘 픽셀 검출기에 가하게 됩니다. 이러한 강력한 방사선은 실리콘 센서 내부와 표면에 손상을 일으켜 다음과 같은 결과를 초래합니다.

• 누설 전류 증가.
• 작동 (공핍) 전압의 변화.
• 신호 손실 및 전하 수집 열화.
• 프론트엔드 전자 장치 동작의 변경.

이러한 효과의 정확한 예측은 작동 전압을 추정하고 재구성 알고리즘을 훈련시키는 데 필수적입니다. 페루기아 대학교 방사선 손상 모델(특정 벌크 및 표면 결함 물리를 포함) 이 Synopsys Sentaurus TCAD에 성공적으로 구현되었음에도 불구하고, Silvaco TCAD에서의 이식성과 예측 능력을 검증할 필요가 있었습니다. 이러한 두 가지 주요 산업 표준 시뮬레이션 도구 간의 교차 검증이 없다면, 방사선 손상 모델링의 견고성은 광범위한 채택을 위해 검증되지 않은 채 남아있게 됩니다.

2. 방법론

저자들은 Synopsys와 Silvaco TCAD 도구를 모두 사용하여 장치 시뮬레이션에 대한 비교 연구를 수행했습니다.

• 장치 구조: n-on-p 다이오드의 2 차원 시뮬레이션을 사용했습니다.
  • 두께: 50 µm.
  • 벌크 농도: $5 \times 10^{12} \text{ cm}^{-3}$ (보론).
  • 피치: 작은 수집 전극을 가진 20 µm.
  • 표면: 1 µm 산화막 층.
• 물리 모델: 공정한 비교를 위해 두 도구 모두 동일한 물리 모델을 사용했습니다.
  • 통계: 페르미 - 디랙.
  • 재결합: 샤르페터 - 쇼클리 - 리드 - 홀.
  • 밴드 - 밴드 터널링: 쉔크.
  • 밴드갭 축소: 슬롯붐.
  • 충돌 이온화: 반 오버스트라텐과 드 만.
• 방사선 손상 구현: 페루기아 모델을 Silvaco 로 이식했습니다. 여기에는 다음이 포함됩니다.
  • 벌크 모델: 선속 ($\Phi$) 에 비례하는 밀도를 가진 두 개의 수용체 결함과 하나의 공여체 결함. 고선속 ($> 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 와 관련된 고방사선 환경을 시뮬레이션할 때 두 가지 다른 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 에서 두 도구 간의 일치는 두 선속 모두에서 "실질적으로 완벽"했습니다.
• 전기장 프로파일:
  • 벌크 영역에서 전기장 프로파일은 두 선속과 모든 바이어스 전압 (100 V, 300 V, 500 V) 에서 1% 이내로 일치했습니다.
  • 전면 및 후면 임플란트 근처에서는 약간의 불일치가 관찰되었으며, 이는 다시 메쉬 불일치 및 임플란트 모델링 차이로 귀결되었습니다.
• 트랩 이온화 확률:
  • 벌크 결함 (수용체 트랩 1 및 2) 의 경우, 이온화 확률은 한 자릿수 범위 내에서 일치했으며, 벌크 영역에서는 20% 이내로 일치했습니다.
  • 전기장에서 퍼센트 수준의 일치를 고려할 때, 이러한 수준의 일치는 만족스러운 것으로 간주됩니다.

5. 중요성

• 페루기아 모델 검증: 이 연구는 페루기아 대학교 방사선 손상 모델이 견고하고 플랫폼에 독립적이며 서로 다른 TCAD 소프트웨어 환경에서 일관된 결과를 산출함을 확인했습니다.
• HL-LHC 를 위한 신뢰성: 누설 전류, 공핍 전압, 전기장과 같은 중요한 관측 가능량에서의 높은 일치도는 HL-LHC 조건 하에서 센서 성능을 예측하는 데 이러한 시뮬레이션을 사용하는 것에 대한 신뢰를 제공합니다.
• 작동 지침: 이 결과는 안전한 작동 전압을 추정하고 신호 손실 메커니즘을 이해하는 데 이러한 모델의 사용을 지지하며, 이는 업그레이드된 ATLAS 및 CMS 추적기의 설계 및 재구성 알고리즘에 필수적입니다.
• 향후 작업: 저자들은 이 검증을 추적 성능의 궁극적인 척도인 수집된 전하 시뮬레이션으로 확장할 계획을 가지고 있습니다.

결론: 이 논문은 페루기아 방사선 손상 모델을 사용할 때 Silvaco 와 Synopsys TCAD 도구가 조사된 실리콘 픽셀 검출기에 대해 상호 일관된 예측을 산출할 수 있음을 보여주며, 향후 고광도 충돌기 응용을 위한 모델의 건전성을 검증합니다.