제1원리 시뮬레이션으로 밝힌 Si, Ge, Si/Ge 초격자의 저에너지 방사선 응답 메커니즘

요약

제1원리 분자동역학(AIMD) 시뮬레이션을 통해 벌크 실리콘(Si), 벌크 게르마늄(Ge), Si/Ge 초격자(SL)에서 저에너지 이온 조사에 의한 변위 손상 거동을 규명했다. 연구 결과, SL 계면의 Ge 원자는 벌크 Ge와 유사한 변위 저항을 보인 반면, 계면 근처 Si 원자는 벌크 Si보다 높은 변위 임계값을 보여 방사선 내성이 향상됨을 확인했다. 벌크와 SL 구조 간 결함 생성 메커니즘 차이가 존재하며, SL에서는 더 복잡한 결함 클러스터가 형성된다. 결함 형성 에너지와 이동 장벽 계산 결과, SL에서 점결함 형성이 더 어렵고 공공(vacancy) 이동성이 낮았다. 이러한 결과는 Si/Ge SL이 방사선 환경에서 동작하는 전자·광전자 소재로 유망함을 시사한다.

배경

최근 Si/Ge 초격자는 조절 가능한 밴드 구조와 향상된 캐리어 이동도 덕분에 차세대 전자·광전자 소자 플랫폼으로 부상했다. 응용 분야는 광통신용 고속 포토다이오드, 마이크로일렉트로닉스, 태양전지, 우주용 부품까지 다양하다. 그러나 우주 공간의 고에너지 이온 조사 시 이러한 소재의 광학·전기적 특성이 변질되어 소자 성능이 저하된다. 따라서 극한 조건에서 Si/Ge SL의 방사선 응답 메커니즘을 원자 수준에서 이해하는 것이 필수적이다.

실험 연구들은 Si/Ge SL이 벌크 Si보다 우수한 방사선 저항성을 가짐을 보고해 왔다. 예를 들어, 단일층 Ge 내장 SL에서 전자 조사 시 벌크 Si 대비 발광 효율 저하가 적었고, 2MeV 양성자 조사한 Ge 양자점 SL 구조에서도 뛰어난 내성이 관측되었다[참조]. Ge 양자우물 소자와 열전 측정에서도 이온 플루언스에 따라 포논 산란 증가와 전자 수송 향상으로 성능 지표가 개선되는 경향이 확인되었다. 이러한 통찰에도 불구하고 Si/Ge SL에서 결함 형성과 진화의 포괄적인 원자론적 이해는 여전히 부족한 실정이다.

제1원리 분자동역학(AIMD)은 원자 스케일에서 방사선 손상을 탐구하는 강력한 도구로, 제1원리 전자구조 계산에서 유도된 원자간 포텐셜을 사용한다. 고전적 MD와 달리 AIMD는 전하 이동과 전자 차폐 효과를 포착해 변위 임계 에너지(E_d)와 결함 구성을 정확히 결정할 수 있다. 본 연구에서는 AIMD를 이용해 벌크 Si, 벌크 Ge, Si₂/Ge₂ SL에서 저에너지 변위 이벤트를 비교 분석하고, 이들의 독특한 방사선 거동의 근본 메커니즘을 규명하며, 결함 형성 및 이동 에너지학을 정량화했다.

연구 방법

변위 캐스케이드 시뮬레이션은 SIESTA 코드로 수행했다. 이온-전자 상호작용에는 노름 보존 트롤리에-마틴(Norm-conserving Troullier-Martins) 슈도포텐셜을, 교환-상관 효과에는 세펄리-알더 매개변수화의 국소 밀도 근사(LDA)를 적용했다. 원자가 상태는 단일 제타 플러스 편극(SZP) 기저 세트로 전개했으며, 1×1×1 k-포인트 메쉬와 60 Ry 컷오프를 사용했다. SL 모델 슈퍼셀은 288개 원자로 구성되며, 두 Si 층과 두 Ge 층이 번갈아 적층된 Si₂/Ge₂ 구조다. Fig. 1은 벌크 Si와 Si/Ge SL의 기하 구조를 보여준다.

일차 녹아웃 원자(PKA)를 선택해 운동 에너지를 부여하며 반동 이벤트를 개시했다. PKA가 격자 위치로 복귀하면 에너지를 5eV씩 높여 반복하고, 영구적으로 변위되면 0.5eV 정밀도로 임계값을 세밀 조정했다. 각 원소별로 벌크와 SL에서 4-5개 결정학적 방향을 조사했다. 모든 실행은 NVE 앙상블로 1.2ps 동안 수행해 수치적 불안정성을 방지했다.

결과 및 고찰

벌크 실리콘과 게르마늄에서의 변위 이벤트

AIMD 계산 결과 벌크 Si 격자 상수는 5.50Å(실험값 5.43Å, 이론값 5.48Å와 일치), 벌크 Ge는 5.71Å(실험값 5.77Å, 이론값 5.65Å와 일치)였다. Table 1은 변위 임계 에너지(E_d)와 결과 결함 구조를 요약하며, Fig. 2-3은 Si와 Ge의 손상 구성을 도시한다.

벌크 Si의 경우 계산된 E_d 값이 실험 측정치와 잘 일치했다: [001] 21eV, [110] ~47.6eV, [111] ~12.9eV. 모든 방향에서 주된 결함 유형은 프렌켈 쌍(공공+간격 원자)이었다. [110] 방향의 E_d(47eV)는 고전적 MD 값(24eV)을 크게 상회했는데, 이는 AIMD가 전하 이동 효과를 포함하기 때문이다. 변위 메커니즘은 방향에 따라 달랐다: [111]에서는 PKA가 간격 위치로 산란되고 밀려난 원자가 격자 위치로 복귀; [001]에서는 복잡한 교환 시퀀스로 2개 공공과 2개 간격 원자 생성; [110]에서는 다중 충돌로 프렌켈 쌍 클러스터 형성. 이는 벌크 Si에서 [110] 방향이 변위에 가장 강함 확인한다.

벌크 Ge는 실험·이론과 잘 부합하는 E_d 값을 보였다: [001] ~18eV, [111]과 [1 1 1] ~9.5eV. [111] 방향에서는 단순 간격 원자 형성, [110]에서는 PKA가 인접 Ge 원자와 자리 바꿈해 1개 공공과 1개 간격 원자 생성. Si와 마찬가지로 [110] 방향이 가장 변위하기 어려웠다. 벌크 Si와 Ge의 결함 지형은 프렌켈 쌍이 지배적이었다.

Si/Ge 초격자에서의 변위 이벤트

Si/Ge 계면 인접 원자를 PKA로 선정해 Si₂/Ge₂ SL을 조사했다. Table 2에 E_d 값과 관련 결함이 정리되어 있으며, Fig. 4-5에 결과 구성이 도시되어 있다.

SL에서 Si 원자는 [1 1 1] 방향에 대해 높은 저항성(E_d=10eV)을 보였으며, Si[001]과 Si[00 1̅]은 각각 46.5eV, 42.5eV로 복잡한 반자리(antisite) 클러스터를 생성했다. SL 내 Ge 원자는 [111]과 [1 1 1̅] 방향에서 벌크 Ge와 유사한 낮은 E_d(≈16eV)을 유지했으나, [001](≈17eV)과 [00 1̅](≈18eV)에서는 임계값이 상승했다. 주목할 점은 Ge[110]의 E_d가 8.5eV로 감소하며 혼합 결함 클러스터(공공, 반자리, 간격 원자)를 생성한다는 것이다. 이러한 결과는 SL이 Si의 방사선 저항성을 향상시키고 Ge에도 어느 정도 기여함을 확인시켜 주며, 기존 실험 관측과 일치한다.

벌크 Si, Ge, Si/Ge 초격자에서의 결함 형성 에너지와 이동 장벽

64원자 슈퍼셀, 6×6×6 k-포인트 샘플링, 500eV 컷오프 조건으로 벌크 Si, 벌크 Ge, SL에서 공공, 간격 원자, 반자리의 결함 형성 에너지를 계산했다(Table 3). 결과적으로 벌크 Si에서 V_Si(3.60eV), 벌크 Ge에서 V_Ge(2.23eV)가 가장 안정한 결함이었으며, SL에서는 대부분의 결함 형성 에너지가 약간 높게 나타났으나 V_Ge(2.73eV)와 V_Si(2.85eV)가 여전히 최저였다.

넛지드 일래스틱 밴드(NEB) 법으로 SL 계면 인접 V_Ge와 V_Si의 이동 장벽을 평가했다(Table 4). 공공은 벌크와 SL 모두에서 [111] 방향으로 가장 쉽게 이동했으나, 장벽은 SL에서 일관되게 더 높았다(예: V_Si [111] 방향: 벌크 0.11eV vs SL 0.12eV). 반면 [100]과 [110] 방향에서는 차이가 미미했다. 이는 SL에서 공공 이동성이 억제되어 조사 하에서 결함 클러스터링과 부풀림을 감소시킬 수 있음을 시사한다.

결론

본 AIMD 연구는 벌크 Si와 Ge의 변위 임계 에너지가 강하게 이방적이며 [110] 방향이 가장 저항적임을 보였다. Si/Ge SL에서는 Si 원자가 [111]과 [001] 방향에서 저항성이 향상되고, Ge 원자는 [110] 방향의 감소분을 제외하면 벌크와 유사한 임계값을 유지한다. 결함 형성 에너지는 SL에서 점결함 생성이 더 어려움을, 이동 장벽은 공공 이동성이 감소함을 입증한다. 특히 가장 유리한 [111] 방향에서 이동 장벽 증가가 두드러진다. 이러한 복합적 효과가 Si/Ge 초격자의 우수한 방사선 내성의 근간이며, 가혹한 방사선 환경에서 동작해야 하는 소자에 매력적인 후보가 된다.